This study proposes a stochastic model updating approach using an improved multi-population migrant genetic algorithm (MMGA) with the metropolis-hastings (MH) algorithm based on the Bayesian inference approach. To enh...
详细信息
This study proposes a stochastic model updating approach using an improved multi-population migrant genetic algorithm (MMGA) with the metropolis-hastings (MH) algorithm based on the Bayesian inference approach. To enhance the search capacity of the genetic algorithm (GA), the MMGA introduces heuristic: the promotion and elimination mechanisms of urban population migration. The stochastic model updating approach includes two stages. The optimal values identified by the MMGA without considering the parameter uncertainties in the first stage are used as the initial values for the MH algorithm to accelerate Bayesian inference for uncertain parameters in the second stage. A practical application, the Yellow River Bridge with a heavy-haul railway line, demonstrates that the proposed approach offers an increased probability of gaining the global optimum, and predictions of displacements and frequencies provided by the updated finite element model are comparable to those obtained from the field tests.
The metropolis-hastings algorithm has been important in the recent development of Bayes methods. This algorithm generates random draws from a target distribution utilizing a sampling (or proposal) distribution. This a...
详细信息
The metropolis-hastings algorithm has been important in the recent development of Bayes methods. This algorithm generates random draws from a target distribution utilizing a sampling (or proposal) distribution. This article compares the properties of three sampling distributions-the independence chain, the random walk chain, and the Taylored chain suggested by Geweke and Tanizaki (Geweke, J., Tanizaki, H. (1999). On Markov Chain Monte-Carlo methods for nonlinear and non-Gaussian state-space models. Communications in Statistics, Simulation and. Computation 28(4):867-894, Geweke, J., Tanizaki, H. (2001). Bayesian estimation of state-space model using the metropolis-hastings algorithm within Gibbs sampling. Computational Statistics and Data Analysis 37(2):151-170).
In this paper, an attempt is made to show a general solution to nonlinear and/or non-Gaussian state-space modeling in a Bayesian framework, which corresponds to an extension of Carlin et al. (J. Amer. Statist. Assoc. ...
详细信息
In this paper, an attempt is made to show a general solution to nonlinear and/or non-Gaussian state-space modeling in a Bayesian framework, which corresponds to an extension of Carlin et al. (J. Amer. Statist. Assoc. 87(418) (1992) 493-500) and Carter and Kohn (Biometrika 81(3) (1994) 541-553;Biometrika 83(3) (1996) 589-601). Using the Gibbs sampler and the metropolis-hastings algorithm, an asymptotically exact estimate of the smoothing mean is obtained from any nonlinear and/or non-Gaussian model. Moreover, taking several candidates of the proposal density function, we examine precision of the proposed Bayes estimator. (C) 2001 Elsevier Science B.V. All rights reserved.
This paper proposes a new method to identify local damages in frame structures based on approximate metropolis-hastings (AMH) algorithm and statistical moment. By analyzing the sensitivity of different statistical mom...
详细信息
This paper proposes a new method to identify local damages in frame structures based on approximate metropolis-hastings (AMH) algorithm and statistical moment. By analyzing the sensitivity of different statistical moment-based damage indices, the fusion index of fourth-order displacement moment and eighth-order acceleration moment is selected. Then the local damages in frame structures are primarily evaluated by AMH algorithm, where Gibbs sampling is adopted. Finally, the uncertainty of identified local damages is analyzed by using probability density evolution method (PDEM). Numerical simulations have been conducted to compare the proposed method with other similar damage detection methods, showing that the proposed method is more time-saving due to the involvement of Gibbs sampling and more accurate in assessing the damage severity. Experimental study of a 12-story benchmark frame testing has also been carried out, further validating the effectiveness of the proposed method.
The development of an efficient MCMC strategy for sampling from complex distributions is a difficult task that needs to be solved for calculating the small failure probabilities encountered in the high-dimensional rel...
详细信息
The development of an efficient MCMC strategy for sampling from complex distributions is a difficult task that needs to be solved for calculating the small failure probabilities encountered in the high-dimensional reliability analysis of engineering systems. Usually different variations of the metropolis-hastings algorithm (MH) are used. However, the standard MH algorithm does not generally work in high dimensions, since it leads to very frequent repeated samples. In order to overcome this deficiency one can use the Modified metropolis-hastings algorithm (MMH) proposed in Au and Beck (2001)[1]. Another variation of the MH algorithm, called the metropolis-hastings algorithm with delayed rejection (MHDR) has been proposed by Tierney and Mira (1999)[7]. The key idea behind the MHDR algorithm is to reduce the correlation between states of the Markov chain. In this paper we combine the ideas of MMH and MHDR and propose a novel modification of the MH algorithm, called the Modified metropolis-hastings algorithm with delayed rejection (MMHDR). The efficiency of the new algorithm is demonstrated with a numerical example where MMHDR is used together with Subset simulation for computing small failure probabilities in high dimensions. (C) 2010 Elsevier Ltd. All rights reserved.
Kernel density estimation is an important tool in visualizing posterior densities from Markov chain Monte Carlo output. It is well known that when smooth transition densities exist, the asymptotic properties of the es...
详细信息
Kernel density estimation is an important tool in visualizing posterior densities from Markov chain Monte Carlo output. It is well known that when smooth transition densities exist, the asymptotic properties of the estimator agree with those for independent data. In this paper, we show that because of the rejection step of the metropolis-hastings algorithm, this is no longer true and the asymptotic variance will depend on the probability of accepting a proposed move. We find an expression for this variance and apply the result to algorithms for automatic bandwidth selection.
Simulation of Markov chain samples using the metropolis-hastings algorithm is useful for reliability estimation. Subset simulation is an example of the reliability estimation method utilizing this algorithm. The effic...
详细信息
Simulation of Markov chain samples using the metropolis-hastings algorithm is useful for reliability estimation. Subset simulation is an example of the reliability estimation method utilizing this algorithm. The efficiency of the simulation is governed by the correlation between the simulated Markov chain samples. The objective of this study is to propose a modified metropolis-hastings algorithm with reduced chain correlation. The modified algorithm differs from the original in terms of the transition probability. It has been verified that the modified algorithm satisfies the reversibility condition and therefore the simulated samples follow the target distribution for the correct theoretical reasons. When applied to subset simulation, the modified algorithm produces a more accurate estimate of failure probability as indicated by a lower coefficient of variation and a lower mean square error. The advantage is more significant for small failure probability. Examples of soil slope with spatially variable properties were presented to demonstrate the applicability of the proposed modification to reliability estimation of engineering problems. It was found that the modified algorithm produces a more accurate estimator over the range of random dimensions studied. (C) 2010 Elsevier Ltd. All rights reserved.
Nonlinear latent variable models are specified that include quadratic forms and interactions of latent regressor variables as special cases. To estimate the parameters, the models are put in a Bayesian framework with ...
详细信息
Nonlinear latent variable models are specified that include quadratic forms and interactions of latent regressor variables as special cases. To estimate the parameters, the models are put in a Bayesian framework with conjugate priors for the parameters. The posterior distributions of the parameters and the latent variables are estimated using Markov chain Monte Carlo methods such as the Gibbs sampler and the metropolis-hastings algorithm. The proposed estimation methods are illustrated by two simulation studies and by the estimation of a non-linear model for the dependence of performance on task complexity and goal specificity using empirical data.
The waste-recycling Monte Carlo (WRMC) algorithm introduced by physicists is a modification of the (multi-proposal) metropolis-hastings algorithm, which makes use of all the proposals in the empirical mean, whereas th...
详细信息
The waste-recycling Monte Carlo (WRMC) algorithm introduced by physicists is a modification of the (multi-proposal) metropolis-hastings algorithm, which makes use of all the proposals in the empirical mean, whereas the standard (multi-proposal) metropolis-hastings algorithm uses only the accepted proposals. In this paper we extend the WRMC algorithm to a general control variate technique and exhibit the optimal choice of the control variate in terms of the asymptotic variance. We also give an example which shows that, in contradiction to the intuition of physicists, the WRMC algorithm can have an asymptotic variance larger than that of the metropolis-hastings algorithm. However, in the particular case of the metropolis-hastings algorithm called the Boltzmann algorithm, we prove that the WRMC algorithm is asymptotically better than the metropolis-hastings algorithm. This last property is also true for the multiproposal metropolis-hastings algorithm. In this last framework we consider a linear parametric generalization of WRMC, and we propose an estimator of the explicit optimal parameter using the proposals.
Judicious choice of candidate generating distributions improves efficiency of the metropolis-hastings algorithm. In Bayesian applications, it is sometimes possible to identify an approximation to the target posterior ...
详细信息
Judicious choice of candidate generating distributions improves efficiency of the metropolis-hastings algorithm. In Bayesian applications, it is sometimes possible to identify an approximation to the target posterior distribution;this approximate posterior distribution is a good choice for candidate generation. These observations are applied to analysis of the Cormack-Jolly-Seber model and its extensions.
暂无评论