We develop a Markov chain Monte Carlo method for a linear regression model with an ARMA(p, q)-GARCH(r, s) error. To generate a Monte Carlo sample from the joint posterior distribution, we employ a Markov chain samplin...
详细信息
We develop a Markov chain Monte Carlo method for a linear regression model with an ARMA(p, q)-GARCH(r, s) error. To generate a Monte Carlo sample from the joint posterior distribution, we employ a Markov chain sampling with the metropolis-hastings algorithm. As illustration, we estimate an ARMA-GARCH model of simulated time series data. (C) 2000 Elsevier Science S.A. All righter reserved. JEL classification: C11;C22.
A Bayesian statistical framework is presented for Zimmerman and Weissenburger flutter margin method which considers the uncertainties in aeroelastic modal parameters. The proposed methodology overcomes the limitations...
详细信息
A Bayesian statistical framework is presented for Zimmerman and Weissenburger flutter margin method which considers the uncertainties in aeroelastic modal parameters. The proposed methodology overcomes the limitations of the previously developed least square based estimation technique which relies on the Gaussian approximation of the flutter margin probability density function (pdf). Using the measured free-decay responses at subcritical (preflutter) airspeeds, the joint non-Gaussain posterior pdf of the modal parameters is sampled using the metropolis-hastings (MH)Markov chain Monte Carlo (MCMC) algorithm. The posterior MCMC samples of the modal parameters are then used to obtain the flutter margin pdfs and finally the flutter speed pdf. The usefulness of the Bayesian flutter margin method is demonstrated using synthetic data generated from a two-degree-of-freedom pitch-plunge aeroelastic model. The robustness of the statistical framework is demonstrated using different sets of measurement data. It will be shown that the probabilistic (Bayesian) approach reduces the number of test points required in providing a flutter speed estimate for a given accuracy and precision. (C) 2016 Elsevier Ltd. All rights reserved.
In this paper we address the problem of the prohibitively large computational cost of existing Markov chain Monte Carlo methods for large-scale applications with high-dimensional parameter spaces, e.g., in uncertainty...
详细信息
In this paper we address the problem of the prohibitively large computational cost of existing Markov chain Monte Carlo methods for large-scale applications with high-dimensional parameter spaces, e.g., in uncertainty quantification in porous media flow. We propose a new multilevel metropolis-hastings algorithm and give an abstract, problem-dependent theorem on the cost of the new multilevel estimator based on a set of simple, verifiable assumptions. For a typical model problem in subsurface flow, we then provide a detailed analysis of these assumptions and show significant gains over the standard metropolis-hastings estimator. Numerical experiments confirm the analysis and demonstrate the effectiveness of the method with consistent reductions of more than an order of magnitude in the cost of the multilevel estimator over the standard metropolis-hastings algorithm for tolerances epsilon < 10(-2).
We examine some Markov chain Monte Carlo (MCMC) methods for a generalized non-linear regression model, the Logit model. It is first shown that MCMC algorithms may be used since the posterior is proper under the choice...
详细信息
We examine some Markov chain Monte Carlo (MCMC) methods for a generalized non-linear regression model, the Logit model. It is first shown that MCMC algorithms may be used since the posterior is proper under the choice of non-informative priors. Then two non-standard MCMC methods are compared: a metropolis-hastings algorithm with a bivariate normal proposal resulting from an approximation, and a metropolis-hastings algorithm with an adaptive proposal. The results presented here are illustrated by simulations, and show the good behavior of both methods, and superior performances of the method with an adaptive proposal in terms of convergence to the stationary distribution and exploration of the posterior distribution surface. (C) 2002 Elsevier Science B.V. All rights reserved.
The design of a blind receiver for coded orthogonal frequency-division multiplexing communication systems in the presence of frequency offset and frequency-selective fading is investigated. The proposed blind receiver...
详细信息
The design of a blind receiver for coded orthogonal frequency-division multiplexing communication systems in the presence of frequency offset and frequency-selective fading is investigated. The proposed blind receiver iterates between a Bayesian demodulation stage and a maximum a posteriori channel decoding stage. The extrinsic a posteriori probabilities of data symbols are iteratively exchanged between these two stages to achieve successively improved performance. The Bayesian demodulator computes the a posteriori data symbol probabilities, based on the received signals (without knowing or explicitly estimating the frequency offset and the fading channel states), by using Markov chain Monte Carlo (MCMC) techniques. In particular, two MCMC methods-the metropolis-hastings algorithm and the Gibbs sampler-are studied for this purpose. Computer simulation results show that the proposed Bayesian blind turbo receiver can achieve good performance and is robust against modeling mismatch.
This work investigates joint estimation of symbol timing synchronization and channel response in two-way relay networks (TWRN) that utilize amplify-and-forward (AF) relay strategy. With unknown relay channel gains and...
详细信息
This work investigates joint estimation of symbol timing synchronization and channel response in two-way relay networks (TWRN) that utilize amplify-and-forward (AF) relay strategy. With unknown relay channel gains and unknown timing offset, the optimum maximum likelihood (ML) algorithm for joint timing recovery and channel estimation can be overly complex. We develop a new Bayesian based Markov chain Monte Carlo (MCMC) algorithm in order to facilitate joint symbol timing recovery and effective channel estimation. In particular, we present a basic metropolis-hastings algorithm (BMH) and a metropolis-hastings-ML (MH-ML) algorithm for this purpose. We also derive the Cramer-Rao lower bound (CRLB) to establish a performance benchmark. Our test results of ML, BMH, and MH-ML estimation illustrate near-optimum performance in terms of mean-square errors (MSE) and estimation bias. We further present bit error rate (BER) performance results.
In this paper we develop new Markov chain Monte Carlo schemes for the estimation of Bayesian models. One key feature of our method, which we call the tailored randomized block metropolis-hastings (TaRB-MH) method, is ...
详细信息
In this paper we develop new Markov chain Monte Carlo schemes for the estimation of Bayesian models. One key feature of our method, which we call the tailored randomized block metropolis-hastings (TaRB-MH) method, is the random clustering of the parameters at every iteration into all arbitrary number of blocks. Then each block is sequentially updated through an M-H step. Another feature is that the proposal density for each block is tailored to the location and curvature of the target density based on the output of simulated annealing, following Chib and Greenberg ( 1994, 1995) and Chib and Ergashev (in press). We also provide all extended version of our method for sampling multi-modal distributions in which at a pre-specified mode jumping iteration, a single-block proposal is generated from one of the modal regions using a Mixture proposal density, and this proposal is then accepted according to all M-H probability of move. At the non-mode jumping iterations, the draws are obtained by applying the TaRB-MH algorithm. We also discuss flow the approaches of Chib (1995) and Chib and Jeliazkov (2001) call be adapted to these sampling schemes for estimating the model marginal likelihood. The methods are illustrated in several problems. In the DSGE model of Smets and Wouters (2007), for example, which involves a 36-dimensional posterior distribution, we show that the autocorrelations of the sampled draws from the TaRB-MH algorithm decay to zero within 30-40 lags for most parameters. In contrast, the sampled draws from the random-walk M-H method, the algorithm that has been used to date in the context of DSGE models, exhibit significant autocorrelations even at lags 2500 and beyond. Additionally, the RW-MH does not explore the same high density regions of the posterior distribution as the TaRB-MH algorithm. Another example concerns the model of An and Schorfheide (2007) where the posterior distribution is Multimodal. While the RW-MH algorithm is Unable to jump from the low
We introduce a new adaptive MCMC algorithm, based on the traditional single component metropolis-hastings algorithm and on our earlier adaptive metropolisalgorithm (AM). In the new algorithm the adaption is performed...
详细信息
We introduce a new adaptive MCMC algorithm, based on the traditional single component metropolis-hastings algorithm and on our earlier adaptive metropolisalgorithm (AM). In the new algorithm the adaption is performed component by component. The chain is no more Markovian, but it remains ergodic. The algorithm is demonstrated to work well in varying test cases up to 1000 dimensions.
In recent years much effort has been devoted to maximum likelihood estimation of generalized linear mixed models. Most of the existing methods use the EM algorithm, with various techniques in handling the intractable ...
详细信息
In recent years much effort has been devoted to maximum likelihood estimation of generalized linear mixed models. Most of the existing methods use the EM algorithm, with various techniques in handling the intractable E-step. In this paper, a new implementation of a stochastic approximation algorithm with Markov chain Monte Carlo method is investigated. The proposed algorithm is computationally straightforward and its convergence is guaranteed. A simulation and three real data sets, including the challenging salamander data, are used to illustrate the procedure and to compare it with some existing methods. The results indicate that the proposed algorithm is an attractive alternative for problems with a large number of random effects or with high dimensional intractable integrals in the likelihood function.
We propose a fully Bayesian inference for semiparametric joint mean and variance models on the basis of B-spline approximations of nonparametric components. An efficient MCMC method which combines Gibbs sampler and Me...
详细信息
We propose a fully Bayesian inference for semiparametric joint mean and variance models on the basis of B-spline approximations of nonparametric components. An efficient MCMC method which combines Gibbs sampler and metropolis-hastings algorithm is suggested for the inference, and the methodology is illustrated through a simulation study and a real example. (c) 2013 Elsevier B.V. All rights reserved.
暂无评论