Bayesian indirect inference using a parametric auxiliary model

Indirect inference (II) is a methodology for estimating the parameters of an intractable (generative) model on the basis of an alternative parametric (auxiliary) model that is both analytically and computationally easier to deal with. Such an approach has been well explored in the classical literature but has received substantially less attention in the Bayesian paradigm. The purpose of this paper is to compare and contrast a collection of what we call parametric Bayesian indirect inference (pBII) methods. One class of pBII methods uses approximate Bayesian computation (referred to here as ABC II) where the summary statistic is formed on the basis of the auxiliary model, using ideas from II. Another approach proposed in the literature, referred to here as parametric Bayesian indirect likelihood (pBIL), we show to be a fundamentally different approach to ABC II. We devise new theoretical results for pBIL to give extra insights into its behaviour and also its differences with ABC II. Furthermore, we examine in more detail the assumptions required to use each pBII method. The results, insights and comparisons developed in this paper are illustrated on simple examples and two other substantive applications. The first of the substantive examples involves performing inference for complex quantile distributions based on simulated data while the second is for estimating the parameters of a trivariate stochastic process describing the evolution of macroparasites within a host based on real data. We create a novel framework called Bayesian indirect likelihood (BIL) which encompasses pBII as well as general ABC methods so that the connections between the methods can be established.

Fate and transport of bensulfuron-methyl and imazosulfuron in paddy fields: Experiments and model simulation

Experiments were conducted to determine the fate of bensulfuron-methyl (BSM) and imazosulfuron (IMS) under paddy conditions. Initially, laboratory experiments were conducted and the photolysis half-lives of the two herbicides were found to be much shorter than their hydrolysis half-lives in aqueous solutions. In the aerobic water–soil system, dissipation followed first-order kinetics with water half-lives of 9.1 and 11.0 days and soil half-lives of 12.4 and 18.5 days (first phase) and 35.0 and 44.1 days (second phase) for bensulfuron-methyl and imazosulfuron, respectively. However, the anaerobic soil half-lives were only 12.7 and 9.8 days for BSM and IMS, respectively. The values of K d were determined to be 16.0 and 13.8 for BSM and IMS, respectively. Subsequent field measurements for the two herbicides revealed that dissipation of both herbicides in paddy water involved biphasic first-order kinetics, with the dissipation rates in the first phase being much faster than those in the second phase. The dissipation of bensulfuron-methyl and imazosulfuron in the paddy surface soil were also followed biphasic first-order kinetics. These results were then used as input parameters for the PCPF-1 model to simulate the fate and transport of BSM and IMS in the paddy environment (water and 1-cm surface soil layer). The measured and simulated values agreed well and the mass balance error during the simulation period was −1.2 and 2.8% of applied pesticide, respectively, for BSM and IMS.

An approximate method for the solution of an influence function foil rolling model

Fleck and Johnson (Int. J. Mech. Sci. 29 (1987) 507) and Fleck et al. (Proc. Inst. Mech. Eng. 206 (1992) 119) have developed foil rolling models which allow for large deformations in the roll profile, including the possibility that the rolls flatten completely. However, these models require computationally expensive iterative solution techniques. A new approach to the approximate solution of the Fleck et al. (1992) Influence Function Model has been developed using both analytic and approximation techniques. The numerical difficulties arising from solving an integral equation in the flattened region have been reduced by applying an Inverse Hilbert Transform to get an analytic expression for the pressure. The method described in this paper is applicable to cases where there is or there is not a flat region.