This article describes an inhomogeneous Poisson point process in the plane with an intensity function based on a Dirichlet tessellation process and a method for using observations on the point process to make fully Ba...
详细信息
This article describes an inhomogeneous Poisson point process in the plane with an intensity function based on a Dirichlet tessellation process and a method for using observations on the point process to make fully Bayesian inferences about the underlying tessellation. The method is implemented using a Markov chain Monte Carlo approach. An application to modeling the territories of clans of badgers, Meles meles, is described.
We study the appraisal problem for the joint inversion of seismic and controlled source electro-magnetic (CSEM) data and utilize rock-physics models to integrate these two disparate data sets. The appraisal problem is...
详细信息
We study the appraisal problem for the joint inversion of seismic and controlled source electro-magnetic (CSEM) data and utilize rock-physics models to integrate these two disparate data sets. The appraisal problem is solved by adopting a Bayesian model and we incorporate four representative sources of uncertainty. These are uncertainties in 1) seismic wave velocity, 2) electric conductivity, 3) seismic data and 4) CSEM data. The uncertainties in porosity and water saturation are quantified by a posterior random sampling in the model space of porosity and water saturation in a marine one-dimensional structure. We study the relative contributions from the four individual sources of uncertainty by performing several statistical experiments. The uncertainties in the seismic wave velocity and electric conductivity play a more significant role on the variation of posterior uncertainty than do the seismic and CSEM data noise. The numerical simulations also show that the uncertainty in porosity is most affected by the uncertainty in the seismic wave velocity and that the uncertainty in water saturation is most influenced by the uncertainty in electric conductivity. The framework of the uncertainty analysis presented in this study can be utilized to effectively reduce the uncertainty of the porosity and water saturation derived from the integration of seismic and CSEM data.
In this paper we propose a new Monte Carlo EM algorithm to compute maximum likelihood estimates in the context of random effects models. The algorithm involves the construction of efficient sampling distributions for ...
详细信息
In this paper we propose a new Monte Carlo EM algorithm to compute maximum likelihood estimates in the context of random effects models. The algorithm involves the construction of efficient sampling distributions for the Monte Carlo implementation of the E-step, together with a reweighting procedure that allows repeatedly using a same sample of random effects. In addition, we explore the use of stochastic approximations to speed up convergence once stability has been reached. Our algorithm is compared with that of McCulloch (1997). Extensions to more general problems are discussed. (C) 1999 Elsevier Science B.V. All rights reserved.
Stochastic search algorithms inspired by physical and biological systems are applied to the problem of learning directed graphical probability models in the presence of missing observations and hidden variables. For t...
详细信息
Stochastic search algorithms inspired by physical and biological systems are applied to the problem of learning directed graphical probability models in the presence of missing observations and hidden variables. For this class of problems, deterministic search algorithms tend to halt at local optima, requiring random restarts to obtain solutions of acceptable quality. We compare three stochastic search algorithms: a metropolis-hastings Sampler (MHS), an Evolutionary algorithm (EA), and a new hybrid algorithm called Population Markov Chain Monte Carlo, or popMCMC. PopMCMC uses statistical information from a population of MHSs to inform the proposal distributions for individual samplers in the population. Experimental results show that popMCMC and EAs learn more efficiently than the MHS with no information exchange. Populations of MCMC samplers exhibit more diversity than populations evolving according to EAs not satisfying physics-inspired local reversibility conditions.
We present explicit methods for simulating diffusions whose generator is self-adjoint with respect to a known (but possibly not normalizable) density. These methods exploit this property and combine an optimized Runge...
详细信息
We present explicit methods for simulating diffusions whose generator is self-adjoint with respect to a known (but possibly not normalizable) density. These methods exploit this property and combine an optimized Runge Kutta algorithm with a metropolishastings Monte Carlo scheme. The resulting numerical integration scheme is shown to be weakly accurate at finite noise and to gain higher order accuracy in the small noise limit. It also permits the user to avoid computing explicitly certain terms in the equation, such as the divergence of the mobility tensor, which can be tedious to calculate. Finally, the scheme is shown to be ergodic with respect to the exact equilibrium probability distribution of the diffusion when it exists. These results are illustrated in several examples, including a Brownian dynamics simulation of DNA in a solvent. In this example, the proposed scheme is able to accurately compute dynamics at time step sizes that are an order of magnitude (or more) larger than those permitted with commonly used explicit predictor-corrector schemes.
In this paper, we propose a framework for exploring the statistical advantages of nondestructive evaluation (NDE) over destructive testing (DT). Two cases are considered: (1) 0-1 or pass/fail data, and (2) continuous ...
详细信息
In this paper, we propose a framework for exploring the statistical advantages of nondestructive evaluation (NDE) over destructive testing (DT). Two cases are considered: (1) 0-1 or pass/fail data, and (2) continuous measurement data. While NDE data are less expensive to collect, they are less precise than DT data. However, taking more NDE data provides more precision. The proposed framework provides a way to calculate equivalent sample sizes of NDE data for given sample sizes of DT data. We illustrate the proposed framework with a radiographic example.
Among different numerical methods for modeling turbulent flow, Reynolds-averaged Navier-Stokes (RANS) is the most commonly used and computationally reasonable. However, the accuracy of RANS is lower than that of other...
详细信息
Among different numerical methods for modeling turbulent flow, Reynolds-averaged Navier-Stokes (RANS) is the most commonly used and computationally reasonable. However, the accuracy of RANS is lower than that of other high-fidelity numerical methods. In this work, the uncertainties associated with the coefficients of the standard k - epsilon RANS turbulence model are estimated and calibrated to improve the accuracy. The calibration is performed by considering the coefficients individually as well as collectively. The first three coefficients of the standard k - epsilon turbulence model are calibrated among the five coefficients (C-mu, C-epsilon 1, C-epsilon 2, sigma(epsilon). and sigma(k)). The Bayesian inference technique using the metropolis-hastings algorithm is applied to quantify uncertainties and calibration. Flow over a periodic hill is selected as a test case. The separation height of the bubble at x/h = 2 and x/h = 4, along with the streamwise velocity at various locations, has been chosen as the quantities of interest for comparing the results with DNS. The calibration is performed using known high-fidelity data (direct numerical simulation) from the available data set. The velocity field is re-calculated from the calibrated closure coefficients and compared with the same calculated with the standard coefficients of k - epsilon turbulence model (baseline). The deviation of calibrated C-mu is almost 50%-60% from baseline and for C-epsilon 1 and C-epsilon 2 it is 3%-12% and 6%-9% respectively. The algorithm is tested for different Reynold numbers and data points. A sensitivity analysis is also performed.
作者:
Lang, JBUniv Iowa
Dept Stat & Actuarial Sci Iowa City IA 52242 USA
An ordinal and binary regression model with parametric link is introduced. The link is a member of a one-parameter family of "mixture links", a family that comprises smooth mixtures of the extreme minimum-va...
详细信息
An ordinal and binary regression model with parametric link is introduced. The link is a member of a one-parameter family of "mixture links", a family that comprises smooth mixtures of the extreme minimum-value, extreme maximum-value, and logistic distributions. A Bayesian version of this flexible model serves as a vehicle for introducing a priori information regarding the choice of link. Owing to non-conjugacy, posterior and predictive distributions are approximated using Markov chain Monte Carlo simulation methods. Link-independent, Bayesian interpretations of covariate effects are described. The method is illustrated through the analyses of several data sets. (C) 1999 Elsevier Science B.V. All rights reserved.
作者:
Wang, ZhongleiXiamen Univ
Wang Yanan Inst Studies Econ MOE Key Lab Econometr Xiamen 361005 Fujian Peoples R China Xiamen Univ
Sch Econ Xiamen 361005 Fujian Peoples R China
Markov chain Monte Carlo methods are widely used to draw a sample from a target distribution which is hard to characterize analytically, and reservoir sampling is developed to obtain a sample from a data stream sequen...
详细信息
Markov chain Monte Carlo methods are widely used to draw a sample from a target distribution which is hard to characterize analytically, and reservoir sampling is developed to obtain a sample from a data stream sequentially in a single pass. A stochastic thinning algorithm using reservoir sampling is proposed, and it can be embedded in most Markov chain Monte Carlo methods to reduce the autocorrelation among the generated sample. The distribution of the sample generated by the proposed sampling algorithm converges in total variation to the target distribution in probability under mild conditions. A practical method is introduced to detect the convergence of the proposed sampling algorithm. Two simulation studies are conducted to compare the proposed sampling algorithm and the corresponding Markov chain Monte Carlo methods without thinning, and results show that the estimation bias of the proposed sampling algorithm is approximately the same as the corresponding Markov chain Monte Carlo method, but the proposed sampling algorithm has a smaller Monte Carlo variance. The proposed sampling algorithm saves computer memory in the sense that the storage of a small portion of the Markov chain is required in each iteration. (C) 2019 Elsevier B.V. All rights reserved.
In the current study we develop the robust Bayesian inference for the generalized inverted family of distributions (GIFD) under an epsilon-contamination class of prior distributions for the shape parameter alpha, with...
详细信息
In the current study we develop the robust Bayesian inference for the generalized inverted family of distributions (GIFD) under an epsilon-contamination class of prior distributions for the shape parameter alpha, with different possibilities of known and unknown scale parameter. We used Type II censoring and Bartholomew sampling scheme (1963) for the following derivations under the squared-error loss function (SELF) and linear exponential (LINEX) loss function : ML-II Bayes estimators of the i) parameters;ii) Reliability function and;iii) Hazard function. We also present simulation study and analysis of a real data set.
暂无评论