1 MCMC Methodology. ISyE8843A, Brani Vidakovic Handout Theoretical Background and Notation
|
|
- Marcus Chambers
- 6 years ago
- Views:
Transcription
1 ISyE8843A, Brani Vidakovic Handout MCMC Methodology. Independence of X,..., X n is not critical for an approximation of the form E θ x h(x) = n n i= h(x i), X i π(θ x). In fact, when X s are dependent, the ergodic theorems describe the approximation. An easy and convenient form of dependence is Markov chain dependence. The Markov dependence is perfect for computer simulations since for producing a future realization of the chain, only the current state is needed.. Theoretical Background and Notation Random variables X, X 2,..., X n,... constitute a Markov Chain on continuous state space if they possess a Markov property, P (X n+ A X,..., X n ) = P (X n+ A X,..., X n ) = Q(X n, A) = Q(A X n ), for some probability distribution Q. Typically, Q is assumed a time-homogeneous, i.e., independent on n ( time ). The transition (from the state n to the state n + ) kernel defines a probability measure on the state space and we will assume that the density q exists, i.e., Q(A X n = x) = q(x, y)dy = q(y x)dy. Distribution Π is invariant, if for all measurable sets A Π(A) = Q(A x)π(dx). A If the transition density π exists, it is stationary if q(x y)π(y) = q(y x)π(x). Here and in the sequel we assume that the density for Π exists, Π(A) = A π(x)dx. A distribution Π is an equilibrium distribution if for Q n (A x) = P (X n A X = x), lim n Qn (A x) = Π(A). In plain terms, the Markov chain will forget the initial distribution and will converge to the stationary distribution. The Markov Chain is irreducible if for each A for which Π(A) >, and for each x, one can find n, so that Q n (A x) >. The Markov Chain X,..., X n,... is recurrent if for each B such that Π(B) >, P (X n B i.o. X = x) =, a.s.(in distribution of X ) A It is Harris recurrent if P (X n B i.o. X = x) =, ( x). The acronym i.o. stands for infinitely often.
2 Figure : Nicholas Constantine Metropolis, Metropolis Algorithm Metropolis algorithm is the fundamental to MCMC development. Assume that the target distribution is known up to a normalizing constant. We would like to construct a chain with π as its stationary distribution. As in ARM, we take a proposal distribution q(x, y) = q(y x), where the proposal for a new value of a chain is y, given that the chain is at value x. Thus q defines transition kernel Q(A, x) = A q(y x)dx which is the probability of transition to some y A. Detailed Balance Equation. A Markov Chain with transition density q(x, y) = q(y x) satisfies detailed balance equation if there exists a distribution f such that q(y x)f(x) = q(x y)f(y). () The distribution f is stationary (invariant) and the chain is reversible. Indeed, if () holds, q(x y)f(y)dy = q(y x)f(x)dy = f(x) q(y x)dy = f(x), which is the definition of invariant distribution. For a given target distribution π, the proposal q is admissible if supp π(x) x supp q( x). Metropolis-Hastings Algorithm is universal. One can select an arbitrary proposal distribution that is admissible. Of course such arbitrary distribution/kernel cannot be expected to satisfy the detailed balance equation () for the target distribution π, i.e, Suppose (wlog) q(y x)π(x) q(x y)π(y). q(y x)π(x) > q(x y)π(y). Then there is a factor ρ(x, y) such that the above inequality is balanced, q(y x) ρ(x, y) π(x) = q(x y)π(y). 2
3 By solving with respect to ρ(x, y) one obtains, ρ(x, y) = q(x y)π(y) q(y x)π(x), where a b denotes min{a, b}. What is the transition kernel corresponding to modified equation? q M (y x) = q(y x)ρ(x, y) + (y = x)( q(y x)ρ(x, y)dy). Metropolis-Hastings Algorithm. Assume that target distribution π is known up to the normalizing constant. This may be the case of posteriors which are always known up to the proportionality constant as products of the likelihood and a prior. STEP Start with arbitrary x from the support of target distribution. STEP 2 At stage n, generate proposal y from q(y x n ). Take x n+ = y with probability ρ(x n, y) = q(xn y)π(y) q(y x n )π(x n ). Otherwise, take x n+ = x n. This random acceptance is done by gen- STEP 3 erating a uniform on (,) random variable U and accepting the proposal y if U ρ(x n, y). STEP 4 Increase n and return to STEP 2. Some Common Choices for q. If q(x y) = q(y x), i.e. if the kernel is symmetric, the acceptance ratio ρ(x, y) simplifies to π(y) π(x), since the proposal kernels from the numerator and denominator cancel. If in addition q depends on (x, y) via y x, i.e., q(x, y) = q ( y x ), for some distribution q, the algorithm is called the Metropolis random walk. A symmetric kernel is the original proposal from Metropolis et al. (953). If the proposal q(x, y) does not depend on x, i.e., q(y x) = q(y), the algorithm is called theindependence Metropolis. It is similar to the aceptance/rejection method (ARM) but unlike the ARM, every step produces a realization from the target distribution. That realization may be repeated many times which is the case when proposal is not accepted and current state is repeatedly taken to be the new state..2. Examples [From Johnson and Albert (999)] A small company improved a product and wants to infer about the proportion of potential customers who will buy the product if the new product is preferred to the old one. The company is certain that this proportion will exceed.5, i.e. and uses the uniform prior on [.5, ]. Out of 2 customers surveyed, 2 prefer the new product. Find the posterior for p. 3
4 (a) x 4 (b) Figure 2: s =.5 Since the support of p is [.5,], we transform the data by θ = log p.5 p, so that θ (, ). For such θ it is easier to specify the proposal, although one can construct Metropolis chain for the untransformed parameter. The inverse transformation is p = /2 + exp{θ} + exp{θ}, with Jacobian /2 exp{θ} (+exp{θ}) 2, and the density for θ is proportional to (/2 + exp{θ}) 2 exp{θ} ( + exp{θ}) 22. The proposal distribution is normal N (θ n, s 2 ), where θ n is current state of the chain and s 2 is to be specified. Here is matlab program illustrating the sampling (albertmc.m at the course web page). % nn = 4; % nn=number of metropolis iterations s=; % s = std of normal proposal density burn=2; % burn = burnin amount % ps=[]; thetas=[]; %transformed p s old = ; % start, theta_ for i = :nn prop = old + s*randn(,); %proposal from N(theta_old, sˆ2) u = rand(,); ep=exp(prop); eo=exp(old); post_p=((/2 + ep)ˆ2 * ep)/((+ep)ˆ22); post_o=((/2 + eo)ˆ2 * eo)/((+eo)ˆ22); new = old; if u <= min(post_p/post_o, ) 4
5 new = prop; %accept proposal as new old = new; % and set old to be the new % for the next iteration; end thetas = [thetas, new]; %collect all theta s ps=[ps, (/2+exp(new))/(+exp(new))]; %back-transformation to p s. end (a) x 4 (b) Figure 3: s = Figures 2 and 3 illustrate the simulation for s =.5 and s = in the proposal distribution. Panels (a) depict the histograms for θ and p, while the panel (b) depicts the last 5 simulations of the chain. Notice that the chain in Figure 2(b) mixes well and indeed reminds a random walk, while its counterpart in Figure 3(b) shows poor mixing. Weibull Example. The Weibull distribution is used extensively in reliability, queueing theory, and many other engineering applications, partly for its ability to describe different hazard rate behavior and partly for historic reasons. The Weibull distribution parameterized by α - the shape or slope, and η /α - the scale, f(x α, η) = αηx α e xαη, is not a member of the exponential family of distributions and explicit posteriors for α and η are impossible. Consider the prior π(α, η) e α ηβ e βη, and observations data = [.2..25];. Imagine these data are extremely expensive obtained by performing a destructive inspection of the pricey products. Construct MCMC based on the Metropolis-Hastings algorithm and approximate posteriors for α and η. Assume the hyperparameter beta = 2; and proposal distribution q(α, η α, η) = } { αη exp α α η η 5
6 (product of two exponentials with means α and η). Note that q(α, η α, η) q(α, η α, η ) and q does not cancel in the expression for ρ. Some hints that are checked and should work well: (i) Start with arbitrary initial values, say: alpha = 2; eta = 2; (ii) Set alphas = []; etas = []; and do: for i = : alpha_prop = - alpha * log(rand); eta_prop = - eta * log(rand); % prod = prod(data); prod2 = prod( exp( eta * data.ˆalpha - eta_prop * data.ˆalpha_prop)); % rr = (eta_prop/eta)ˆ(beta-) * exp(alpha - alpha_prop - beta *... (eta_prop - eta)) * exp(- alpha/alpha_prop - eta/eta_prop +... alpha_prop/alpha + eta_prop/eta)*prod.ˆ(alpha_prop - alpha) *... prod2 *((alpha_prop * eta_prop)/(alpha * eta))ˆ(n-); % r = min( rr,); if (rand < r) alpha = alpha_prop; eta = eta_prop; end alphas = [alphas alpha]; etas = [etas eta]; end Values alpha prop and eta prop are proposals from independent exponential distributions with means alpha and eta. Burn in 5 out of simulations (usually -5 is enough) to make sure that there is no influence of the initial values for α and η, and plot the histograms of their posterior distributions. figure() subplot(,2,) hist(alphas(5:end),) subplot(,2,2) hist(etas(5:end),) Finally, report the mean and variance of alphas and etas. These are desired Bayes estimators with their posterior precisions..3 Gibbs Sampler The Gibbs sampler, introduced in [], is a special case of A Single Component Metropolis algorithm. Define X = (X, X 2,..., X p ). Each step of the algorithm will consist of p coordinatewise updates. Define X i = (X,..., X i, X i+,..., X p ), and X n i = (X n+,..., Xi n+, Xn i+,..., Xn p ), i =,..., p. X n i is the update at step n, where first i coordinates are updated to their values at step n + and the coordinates at positions i +, i + 2,..., p are not updated and still are at step n. Let q i (Y i Xi n, Xn i) be the proposal distribution that generates proposals for ith coordinate only. Then the single component Metropolis proceeds as follows: Generate a candidate Y i q i (Y i Xi n, Xn i). 6
7 Accept the candidate Y i with probability as Xi n+ = Y i. Otherwise set Xi n+ = Xi n. Go to the next coordinate at step n. Let π(x i X i ) = ρ(x n i, X n i, Y i ) = π(y i X n i)q i (X n i Y i, X n i) π(x i X n i)q i (Y i X i, X n i). π(x) R π(x)dx i be the full conditional. The Gibbs sampler is a single step Metropolis algorithm with q i (Y i X i, X i ) = π(y i X i ). Obviously, at each step n, ρ(x n i, Xn i, Y i ) =, and each update in Gibbs algorithm is accepted. In more familiar notation, suppose that θ = (θ,..., θ p ) is multidimensional parameter of interest. Each component can be univariate of multivariate. Suppose that we can simulate from the conditional densities π(θ i θ i ), where as the above, θ i denotes the parameter vector θ without the component i. If the current state of θ is θ n = (θ n, θn 2,..., θn p ), the Gibbs sampler produces θ n+ in the following way: Draw θ n+ from π(θ θ2 n, θn 3,..., θn p ) Draw θ2 n+ from π(θ 2 θ n+, θ3 n,..., θn p ) Draw θ3 n+ from π(θ 3 θ n+, θ2 n+, θ4 n,..., θn p )... Draw θp n+ from π(θ p θ n+, θ2 n+,..., θp 2 n+, θn p ) Draw θp n+ from π(θ p θ n+, θ2 n+,..., θp n+ ) Above, we have assumed a fixed updating order. This may not always be the case, since it is possible to generalize the Gibbs Sampler in a number of ways. It is possible to assume the random updating order, i.e. pick the block to update randomly. It is also possible to update only one block per iteration and to choose the block to update with some preassigned probability, see Gilks et al. [6] for discussion..3. Finding the Full Conditionals The full conditionals, needed for implementation of the Gibbs sampler are conceptually easy to find. From the joint distribution of all variables, only expressions that contain the particular variable are entering to the conditional distribution. The difficulty is (as always) in finding normalizing constants. Suppose θ = (θ s, θ s ) and we are interested in the full conditional for θ s. The full conditional is π(θ s θ s ) = π(θ s, θ s ) π(θs, θ s )dθ s π(θ s, θ s ). Example. A popular simple model to illustrate Gibbs sampler and finding the full conditionals the following [Gilks (from [2]), Chapter 5, page 76.] Y, Y 2,..., Y n N (µ, /τ) µ N (, ) τ Gamma(2, ). 7
8 The joint distribution is { n } f(y, µ, τ) = f(y i µ, τ) π(µ)π(τ) i= { = (2π) (n+)/2 τ n/2 exp τ/2 } n (y i µ) 2 i= exp{ /2µ 2 } τ exp{ τ}. To find the full conditional for µ we select the terms from f(y, µ, τ) that contain µ and normalize. Indeed, π(µ τ, y) = = π(µ, τ y) π(τ y) π(µ, τ, y) π(τ, y) π(µ, τ, y). Thus, { } π(µ τ, y) exp τ n (y i µ) 2 exp{ /2µ 2 } 2 i= exp { 2 ( ( + nτ) µ τ ) } y 2 i, + nτ which is normal N ( τ P y i +nτ, +nτ ) distribution. Similarly, { } n π(τ µ, y) τ n/2 exp τ/2 (y i µ) 2 τ exp{ τ} = τ n/2+ exp { τ i= [ + 2 ]} n (y i µ) 2, which is unnurmalized gamma Gamma(2 + n/2, + n 2 i= (y i µ) 2 ). The matlab code mcmc.m implements the sampler. Function rand gamma.m generates random gamma variates and is a part of BayesLab. We simulated n = 2 observations from N (, 4 2 ) distribution and started a gibs with µ = and τ = 2. n=3; % sample size randn( state, ); y = 4 * randn(,n) + ; % NN = ; mus = []; taus = []; suma = sum(y); mu = ; % set the parameters as prior means tau = 2; % for i = : NN new_mu = sqrt(/(+n*tau)) * randn + (tau * suma)/(+n*tau); 8 i=
9 par = +/2 * sum ( (y - mu).ˆ2 ); new_tau = rand_gamma(2 + n/2, par,,); mus = [mus new_mu]; taus = [taus new_tau]; mu=new_mu; tau=new_tau; end (a) (b) (c) Figure 4: (a) Last 5 simulations (out of ) for µ (top) and τ (below); (b) histograms of µ and τ; (c) joint simulation for µ and τ. Figure 4(a) depicts last 5 simulations (out of ) for µ (top) and τ (below). Panel (b) gives histograms of µ and τ, and panel (c) presents a joint simulation for µ and τ. The burn-in period was, so only 9 variates have been used to approximate the posteriors. The MCMC estimators are ˆµ =.938 and ˆτ =.652. The performance is quite good since the theoretical parameters are and /6=.625, respectively. Exercises. Beetle Mortality. The data comes from Bliss [] (cited in Dobson [2]) and is shown in Table. The data involves counting the number of beetles killed after five hours of exposure to various concentrations of gaseous carbon disulphide (CS 2 ). The analysis concerns estimating the proportion r i /n i of beetles that are killed by the gas. Consider the model ( ) exp{xi } m P (death w i ) = h(w i ) =, + exp{x i } where m >, and w i is a covariate (dose), and x i = w i µ σ, µ R, σ2 >. 9
10 Dosage (log CS 2 mg/litre) Beetles Killed Table : Data on Beetle Mortality from Bliss (935). Batches of adult beetles were exposed to gaseous carbon disulphide for five hours. The priors are m Gamma(a, b ), µ N (c, d ), σ 2 IG(e, f ) [τ = σ 2 Gamma(e, f ). The joint posterior π(µ, σ 2, m y) is proportional to f(y µ, σ 2, m )π(µ, σ 2, m ) ( k ) [h(w i )] y i [ h(w i )] n i y i i= { ma (σ 2 ) e + exp 2 ( µ c d ) } 2 m b f σ 2. The transformation θ = (θ, θ 2, θ 3 ) = (µ, 2 log(σ2 ), log m ) is supported by R 2 and multivariate normal proposal for θ is possible. In the new variables, ( k ) π(θ y) [h(w i )] y i [ h(w i )] n i y i exp{a θ 3 2e θ 2 ) i= exp { 2 ( θ c d ) } 2 exp{θ 3} exp{ 2θ 2}. b f The acceptance probability is more calculationally stable if it is represented as ρ = exp{log π (θ y) log π (θ n y)}, for unnormalized posterior π, θ proposal, and θ n current state. The choice of the hyperparameters is a =.25, b = 4 (m has prior mean, as in the standard logit model), c = 2, d =, e = 2, f =. Proposal density is MVN 3 (θ n, Σ) with Σ = diag(.2,.33,.). 2. Contingency Table. See the matlab file albertmc2.m on the Bayes page. Figure 5 depicts the output.
11 α η p C p L Histogram of differences: p L p C (a) (b) (c) References Figure 5: Graphical output from albertmc2.m. [] Besag, J. (974). Spatial interaction and the statistical analysis of lattice systems. J. Roy. Statist. Soc. Ser. B, 36, [2] Bliss, C. I. (935). The calculation of the dosage-mortality curve. The Annals of Applied Biology, 22, [3] Brooks, S.P. (998). Monte Carlo Methods and its application. The Statistician, 47, 69. [4] Brooks, S.P. (998). Quantitative convergence assessment for Markov chain Monte Carlo via cusums. Statistics and Computing, 8, [5] Casella, G., George, E.I. (992). Explaining the Gibbs Sampler. The American Statistician, 46, (Celebrated Gibbs for Kids ). [6] Chib, S., Greenberg, E. (995). Understanding the Metropolis-Hastings Algorithm. The American Statistician, 49, [7] Gelfand, A.E., Smith, A.F.M. (99). Sampling based approaches to calculating marginal densities. Journal of the American Statistical Association, 85, [8] Gelfand, A.E. (2). Gibbs Sampling. Journal of the American Statistical Association, 95, [9] Robert, C. (2). Bayesian Choice, Second Edition, Springer Verlag. [] Dobson. A. J. (983). An Introduction to Statistical Modelling. Chapman and Hall. [] Geman, S. and Geman, D. (984). Stochastic Relaxation, Gibbs Distributions, and the Bayesian Restoration of Images. IEEE Transactions on Pattern Analysis and Machine Intelligence, Vol. Pami-6, No. 6, [2] Gilks, W.R., Richardson, S., Spiegelhalter, D.J. (996). Markov Chain Monte Carlo in Practice. Chapman and Hall.
12 [3] Hastings, W.K. (97). Monte Carlo sampling methods using Markov chains and their applications. Biometrika, 57, [4] Johnson, V. and Albert, J. (999). Ordinal Data Modeling. Springer Verlag, NY. [5] Metropolis, N., Rosenbluth, M.N., Teller, A.H., Teller, E. (953). Equations of State Calculations by Fast Computing Machines. The Journal of Chemical Physics, 2, [6] Robert, C.P., Casella, G. (999). Monte Carlo Statistical Methods. Springer-Verlag, New York. [7] Roberts, G.O., Gelman, A., Gilks, W.R. (997). Weak Convergence and Optimal Scaling of Random Walk Metropolis Algorithms. Annals of Applied Probability, 7, 2. [8] Tanner, M.A. (996). Tools for Statistical Inference. Springer-Verlag, New York. [9] Tierney, L. (994). Markov Chains for exploring posterior distributions( with discussion). Annals of Statistics, 22,
eqr094: Hierarchical MCMC for Bayesian System Reliability
eqr094: Hierarchical MCMC for Bayesian System Reliability Alyson G. Wilson Statistical Sciences Group, Los Alamos National Laboratory P.O. Box 1663, MS F600 Los Alamos, NM 87545 USA Phone: 505-667-9167
More informationSession 3A: Markov chain Monte Carlo (MCMC)
Session 3A: Markov chain Monte Carlo (MCMC) John Geweke Bayesian Econometrics and its Applications August 15, 2012 ohn Geweke Bayesian Econometrics and its Session Applications 3A: Markov () chain Monte
More informationA quick introduction to Markov chains and Markov chain Monte Carlo (revised version)
A quick introduction to Markov chains and Markov chain Monte Carlo (revised version) Rasmus Waagepetersen Institute of Mathematical Sciences Aalborg University 1 Introduction These notes are intended to
More informationMarkov Chain Monte Carlo methods
Markov Chain Monte Carlo methods Tomas McKelvey and Lennart Svensson Signal Processing Group Department of Signals and Systems Chalmers University of Technology, Sweden November 26, 2012 Today s learning
More informationHastings-within-Gibbs Algorithm: Introduction and Application on Hierarchical Model
UNIVERSITY OF TEXAS AT SAN ANTONIO Hastings-within-Gibbs Algorithm: Introduction and Application on Hierarchical Model Liang Jing April 2010 1 1 ABSTRACT In this paper, common MCMC algorithms are introduced
More informationMarkov Chain Monte Carlo (MCMC)
Markov Chain Monte Carlo (MCMC Dependent Sampling Suppose we wish to sample from a density π, and we can evaluate π as a function but have no means to directly generate a sample. Rejection sampling can
More informationBayesian Inference in GLMs. Frequentists typically base inferences on MLEs, asymptotic confidence
Bayesian Inference in GLMs Frequentists typically base inferences on MLEs, asymptotic confidence limits, and log-likelihood ratio tests Bayesians base inferences on the posterior distribution of the unknowns
More informationWhat is the most likely year in which the change occurred? Did the rate of disasters increase or decrease after the change-point?
Chapter 11 Markov Chain Monte Carlo Methods 11.1 Introduction In many applications of statistical modeling, the data analyst would like to use a more complex model for a data set, but is forced to resort
More informationComputational statistics
Computational statistics Markov Chain Monte Carlo methods Thierry Denœux March 2017 Thierry Denœux Computational statistics March 2017 1 / 71 Contents of this chapter When a target density f can be evaluated
More informationMONTE CARLO METHODS. Hedibert Freitas Lopes
MONTE CARLO METHODS Hedibert Freitas Lopes The University of Chicago Booth School of Business 5807 South Woodlawn Avenue, Chicago, IL 60637 http://faculty.chicagobooth.edu/hedibert.lopes hlopes@chicagobooth.edu
More informationBayesian Methods for Machine Learning
Bayesian Methods for Machine Learning CS 584: Big Data Analytics Material adapted from Radford Neal s tutorial (http://ftp.cs.utoronto.ca/pub/radford/bayes-tut.pdf), Zoubin Ghahramni (http://hunch.net/~coms-4771/zoubin_ghahramani_bayesian_learning.pdf),
More informationMarkov Chain Monte Carlo
Markov Chain Monte Carlo Recall: To compute the expectation E ( h(y ) ) we use the approximation E(h(Y )) 1 n n h(y ) t=1 with Y (1),..., Y (n) h(y). Thus our aim is to sample Y (1),..., Y (n) from f(y).
More informationSimulation of truncated normal variables. Christian P. Robert LSTA, Université Pierre et Marie Curie, Paris
Simulation of truncated normal variables Christian P. Robert LSTA, Université Pierre et Marie Curie, Paris Abstract arxiv:0907.4010v1 [stat.co] 23 Jul 2009 We provide in this paper simulation algorithms
More informationMetropolis-Hastings Algorithm
Strength of the Gibbs sampler Metropolis-Hastings Algorithm Easy algorithm to think about. Exploits the factorization properties of the joint probability distribution. No difficult choices to be made to
More informationA note on Reversible Jump Markov Chain Monte Carlo
A note on Reversible Jump Markov Chain Monte Carlo Hedibert Freitas Lopes Graduate School of Business The University of Chicago 5807 South Woodlawn Avenue Chicago, Illinois 60637 February, 1st 2006 1 Introduction
More informationIntroduction to Markov Chain Monte Carlo & Gibbs Sampling
Introduction to Markov Chain Monte Carlo & Gibbs Sampling Prof. Nicholas Zabaras Sibley School of Mechanical and Aerospace Engineering 101 Frank H. T. Rhodes Hall Ithaca, NY 14853-3801 Email: zabaras@cornell.edu
More informationMarkov Chain Monte Carlo methods
Markov Chain Monte Carlo methods By Oleg Makhnin 1 Introduction a b c M = d e f g h i 0 f(x)dx 1.1 Motivation 1.1.1 Just here Supresses numbering 1.1.2 After this 1.2 Literature 2 Method 2.1 New math As
More informationSampling Methods (11/30/04)
CS281A/Stat241A: Statistical Learning Theory Sampling Methods (11/30/04) Lecturer: Michael I. Jordan Scribe: Jaspal S. Sandhu 1 Gibbs Sampling Figure 1: Undirected and directed graphs, respectively, with
More informationMarkov chain Monte Carlo
Markov chain Monte Carlo Karl Oskar Ekvall Galin L. Jones University of Minnesota March 12, 2019 Abstract Practically relevant statistical models often give rise to probability distributions that are analytically
More informationMinicourse on: Markov Chain Monte Carlo: Simulation Techniques in Statistics
Minicourse on: Markov Chain Monte Carlo: Simulation Techniques in Statistics Eric Slud, Statistics Program Lecture 1: Metropolis-Hastings Algorithm, plus background in Simulation and Markov Chains. Lecture
More informationThe simple slice sampler is a specialised type of MCMC auxiliary variable method (Swendsen and Wang, 1987; Edwards and Sokal, 1988; Besag and Green, 1
Recent progress on computable bounds and the simple slice sampler by Gareth O. Roberts* and Jerey S. Rosenthal** (May, 1999.) This paper discusses general quantitative bounds on the convergence rates of
More informationA short diversion into the theory of Markov chains, with a view to Markov chain Monte Carlo methods
A short diversion into the theory of Markov chains, with a view to Markov chain Monte Carlo methods by Kasper K. Berthelsen and Jesper Møller June 2004 2004-01 DEPARTMENT OF MATHEMATICAL SCIENCES AALBORG
More information16 : Approximate Inference: Markov Chain Monte Carlo
10-708: Probabilistic Graphical Models 10-708, Spring 2017 16 : Approximate Inference: Markov Chain Monte Carlo Lecturer: Eric P. Xing Scribes: Yuan Yang, Chao-Ming Yen 1 Introduction As the target distribution
More informationKobe University Repository : Kernel
Kobe University Repository : Kernel タイトル Title 著者 Author(s) 掲載誌 巻号 ページ Citation 刊行日 Issue date 資源タイプ Resource Type 版区分 Resource Version 権利 Rights DOI URL Note on the Sampling Distribution for the Metropolis-
More informationMarkov Chain Monte Carlo
Chapter 5 Markov Chain Monte Carlo MCMC is a kind of improvement of the Monte Carlo method By sampling from a Markov chain whose stationary distribution is the desired sampling distributuion, it is possible
More informationLecture 8: The Metropolis-Hastings Algorithm
30.10.2008 What we have seen last time: Gibbs sampler Key idea: Generate a Markov chain by updating the component of (X 1,..., X p ) in turn by drawing from the full conditionals: X (t) j Two drawbacks:
More informationMarkov Chain Monte Carlo Using the Ratio-of-Uniforms Transformation. Luke Tierney Department of Statistics & Actuarial Science University of Iowa
Markov Chain Monte Carlo Using the Ratio-of-Uniforms Transformation Luke Tierney Department of Statistics & Actuarial Science University of Iowa Basic Ratio of Uniforms Method Introduced by Kinderman and
More informationComputer intensive statistical methods
Lecture 13 MCMC, Hybrid chains October 13, 2015 Jonas Wallin jonwal@chalmers.se Chalmers, Gothenburg university MH algorithm, Chap:6.3 The metropolis hastings requires three objects, the distribution of
More informationTheory of Stochastic Processes 8. Markov chain Monte Carlo
Theory of Stochastic Processes 8. Markov chain Monte Carlo Tomonari Sei sei@mist.i.u-tokyo.ac.jp Department of Mathematical Informatics, University of Tokyo June 8, 2017 http://www.stat.t.u-tokyo.ac.jp/~sei/lec.html
More informationBridge estimation of the probability density at a point. July 2000, revised September 2003
Bridge estimation of the probability density at a point Antonietta Mira Department of Economics University of Insubria Via Ravasi 2 21100 Varese, Italy antonietta.mira@uninsubria.it Geoff Nicholls Department
More informationMarkov Chain Monte Carlo and Applied Bayesian Statistics
Markov Chain Monte Carlo and Applied Bayesian Statistics Trinity Term 2005 Prof. Gesine Reinert Markov chain Monte Carlo is a stochastic simulation technique that is very useful for computing inferential
More informationPrecision Engineering
Precision Engineering 38 (2014) 18 27 Contents lists available at ScienceDirect Precision Engineering j o ur nal homep age : www.elsevier.com/locate/precision Tool life prediction using Bayesian updating.
More informationMarkov Chain Monte Carlo A Contribution to the Encyclopedia of Environmetrics
Markov Chain Monte Carlo A Contribution to the Encyclopedia of Environmetrics Galin L. Jones and James P. Hobert Department of Statistics University of Florida May 2000 1 Introduction Realistic statistical
More informationMonte Carlo methods for sampling-based Stochastic Optimization
Monte Carlo methods for sampling-based Stochastic Optimization Gersende FORT LTCI CNRS & Telecom ParisTech Paris, France Joint works with B. Jourdain, T. Lelièvre, G. Stoltz from ENPC and E. Kuhn from
More informationMarkov Chain Monte Carlo in Practice
Markov Chain Monte Carlo in Practice Edited by W.R. Gilks Medical Research Council Biostatistics Unit Cambridge UK S. Richardson French National Institute for Health and Medical Research Vilejuif France
More informationMarkov Chain Monte Carlo
Markov Chain Monte Carlo Michael Johannes Columbia University Nicholas Polson University of Chicago August 28, 2007 1 Introduction The Bayesian solution to any inference problem is a simple rule: compute
More informationMarkov Chain Monte Carlo
Markov Chain Monte Carlo (and Bayesian Mixture Models) David M. Blei Columbia University October 14, 2014 We have discussed probabilistic modeling, and have seen how the posterior distribution is the critical
More informationStat 451 Lecture Notes Markov Chain Monte Carlo. Ryan Martin UIC
Stat 451 Lecture Notes 07 12 Markov Chain Monte Carlo Ryan Martin UIC www.math.uic.edu/~rgmartin 1 Based on Chapters 8 9 in Givens & Hoeting, Chapters 25 27 in Lange 2 Updated: April 4, 2016 1 / 42 Outline
More informationSTA 4273H: Statistical Machine Learning
STA 4273H: Statistical Machine Learning Russ Salakhutdinov Department of Computer Science! Department of Statistical Sciences! rsalakhu@cs.toronto.edu! h0p://www.cs.utoronto.ca/~rsalakhu/ Lecture 7 Approximate
More informationControl Variates for Markov Chain Monte Carlo
Control Variates for Markov Chain Monte Carlo Dellaportas, P., Kontoyiannis, I., and Tsourti, Z. Dept of Statistics, AUEB Dept of Informatics, AUEB 1st Greek Stochastics Meeting Monte Carlo: Probability
More informationIntroduction to Machine Learning CMU-10701
Introduction to Machine Learning CMU-10701 Markov Chain Monte Carlo Methods Barnabás Póczos & Aarti Singh Contents Markov Chain Monte Carlo Methods Goal & Motivation Sampling Rejection Importance Markov
More informationCPSC 540: Machine Learning
CPSC 540: Machine Learning MCMC and Non-Parametric Bayes Mark Schmidt University of British Columbia Winter 2016 Admin I went through project proposals: Some of you got a message on Piazza. No news is
More informationBayesian Nonparametric Regression for Diabetes Deaths
Bayesian Nonparametric Regression for Diabetes Deaths Brian M. Hartman PhD Student, 2010 Texas A&M University College Station, TX, USA David B. Dahl Assistant Professor Texas A&M University College Station,
More informationAdaptive Monte Carlo methods
Adaptive Monte Carlo methods Jean-Michel Marin Projet Select, INRIA Futurs, Université Paris-Sud joint with Randal Douc (École Polytechnique), Arnaud Guillin (Université de Marseille) and Christian Robert
More informationComputer Practical: Metropolis-Hastings-based MCMC
Computer Practical: Metropolis-Hastings-based MCMC Andrea Arnold and Franz Hamilton North Carolina State University July 30, 2016 A. Arnold / F. Hamilton (NCSU) MH-based MCMC July 30, 2016 1 / 19 Markov
More informationBRIDGE ESTIMATION OF THE PROBABILITY DENSITY AT A POINT
Statistica Sinica 14(2004), 603-612 BRIDGE ESTIMATION OF THE PROBABILITY DENSITY AT A POINT Antonietta Mira and Geoff Nicholls University of Insubria and Auckland University Abstract: Bridge estimation,
More informationPrinciples of Bayesian Inference
Principles of Bayesian Inference Sudipto Banerjee University of Minnesota July 20th, 2008 1 Bayesian Principles Classical statistics: model parameters are fixed and unknown. A Bayesian thinks of parameters
More informationThe Mixture Approach for Simulating New Families of Bivariate Distributions with Specified Correlations
The Mixture Approach for Simulating New Families of Bivariate Distributions with Specified Correlations John R. Michael, Significance, Inc. and William R. Schucany, Southern Methodist University The mixture
More informationMCMC algorithms for fitting Bayesian models
MCMC algorithms for fitting Bayesian models p. 1/1 MCMC algorithms for fitting Bayesian models Sudipto Banerjee sudiptob@biostat.umn.edu University of Minnesota MCMC algorithms for fitting Bayesian models
More informationMarkov chain Monte Carlo
Markov chain Monte Carlo Peter Beerli October 10, 2005 [this chapter is highly influenced by chapter 1 in Markov chain Monte Carlo in Practice, eds Gilks W. R. et al. Chapman and Hall/CRC, 1996] 1 Short
More informationOn the Optimal Scaling of the Modified Metropolis-Hastings algorithm
On the Optimal Scaling of the Modified Metropolis-Hastings algorithm K. M. Zuev & J. L. Beck Division of Engineering and Applied Science California Institute of Technology, MC 4-44, Pasadena, CA 925, USA
More informationOn Reparametrization and the Gibbs Sampler
On Reparametrization and the Gibbs Sampler Jorge Carlos Román Department of Mathematics Vanderbilt University James P. Hobert Department of Statistics University of Florida March 2014 Brett Presnell Department
More informationLECTURE 15 Markov chain Monte Carlo
LECTURE 15 Markov chain Monte Carlo There are many settings when posterior computation is a challenge in that one does not have a closed form expression for the posterior distribution. Markov chain Monte
More informationStatistical Inference for Stochastic Epidemic Models
Statistical Inference for Stochastic Epidemic Models George Streftaris 1 and Gavin J. Gibson 1 1 Department of Actuarial Mathematics & Statistics, Heriot-Watt University, Riccarton, Edinburgh EH14 4AS,
More informationOverlapping block proposals for latent Gaussian Markov random fields
NORGES TEKNISK-NATURVITENSKAPELIGE UNIVERSITET Overlapping block proposals for latent Gaussian Markov random fields by Ingelin Steinsland and Håvard Rue PREPRINT STATISTICS NO. 8/3 NORWEGIAN UNIVERSITY
More information6 Markov Chain Monte Carlo (MCMC)
6 Markov Chain Monte Carlo (MCMC) The underlying idea in MCMC is to replace the iid samples of basic MC methods, with dependent samples from an ergodic Markov chain, whose limiting (stationary) distribution
More informationMARKOV CHAIN MONTE CARLO
MARKOV CHAIN MONTE CARLO RYAN WANG Abstract. This paper gives a brief introduction to Markov Chain Monte Carlo methods, which offer a general framework for calculating difficult integrals. We start with
More informationInverse Problems in the Bayesian Framework
Inverse Problems in the Bayesian Framework Daniela Calvetti Case Western Reserve University Cleveland, Ohio Raleigh, NC, July 2016 Bayes Formula Stochastic model: Two random variables X R n, B R m, where
More informationBayesian Computational Methods: Monte Carlo Methods
Bayesian Computational Methods: Monte Carlo Methods c These notes are copyrighted by the authors. Unauthorized use is not permitted. Bayesian Computational Methods p.1/?? Noniterative Monte Carlo Methods
More informationMultivariate Slice Sampling. A Thesis. Submitted to the Faculty. Drexel University. Jingjing Lu. in partial fulfillment of the
Multivariate Slice Sampling A Thesis Submitted to the Faculty of Drexel University by Jingjing Lu in partial fulfillment of the requirements for the degree of Doctor of Philosophy June 2008 c Copyright
More informationMarkov Chain Monte Carlo Methods
Markov Chain Monte Carlo Methods John Geweke University of Iowa, USA 2005 Institute on Computational Economics University of Chicago - Argonne National Laboaratories July 22, 2005 The problem p (θ, ω I)
More informationNonlinear Inequality Constrained Ridge Regression Estimator
The International Conference on Trends and Perspectives in Linear Statistical Inference (LinStat2014) 24 28 August 2014 Linköping, Sweden Nonlinear Inequality Constrained Ridge Regression Estimator Dr.
More informationA Search and Jump Algorithm for Markov Chain Monte Carlo Sampling. Christopher Jennison. Adriana Ibrahim. Seminar at University of Kuwait
A Search and Jump Algorithm for Markov Chain Monte Carlo Sampling Christopher Jennison Department of Mathematical Sciences, University of Bath, UK http://people.bath.ac.uk/mascj Adriana Ibrahim Institute
More informationDavid Giles Bayesian Econometrics
David Giles Bayesian Econometrics 5. Bayesian Computation Historically, the computational "cost" of Bayesian methods greatly limited their application. For instance, by Bayes' Theorem: p(θ y) = p(θ)p(y
More informationIntroduction to Bayesian methods in inverse problems
Introduction to Bayesian methods in inverse problems Ville Kolehmainen 1 1 Department of Applied Physics, University of Eastern Finland, Kuopio, Finland March 4 2013 Manchester, UK. Contents Introduction
More informationCSC 2541: Bayesian Methods for Machine Learning
CSC 2541: Bayesian Methods for Machine Learning Radford M. Neal, University of Toronto, 2011 Lecture 3 More Markov Chain Monte Carlo Methods The Metropolis algorithm isn t the only way to do MCMC. We ll
More informationLikelihood-free MCMC
Bayesian inference for stable distributions with applications in finance Department of Mathematics University of Leicester September 2, 2011 MSc project final presentation Outline 1 2 3 4 Classical Monte
More informationThe Metropolis-Hastings Algorithm. June 8, 2012
The Metropolis-Hastings Algorithm June 8, 22 The Plan. Understand what a simulated distribution is 2. Understand why the Metropolis-Hastings algorithm works 3. Learn how to apply the Metropolis-Hastings
More informationBrief introduction to Markov Chain Monte Carlo
Brief introduction to Department of Probability and Mathematical Statistics seminar Stochastic modeling in economics and finance November 7, 2011 Brief introduction to Content 1 and motivation Classical
More informationThe Recycling Gibbs Sampler for Efficient Learning
The Recycling Gibbs Sampler for Efficient Learning L. Martino, V. Elvira, G. Camps-Valls Universidade de São Paulo, São Carlos (Brazil). Télécom ParisTech, Université Paris-Saclay. (France), Universidad
More informationBayesian Estimation of Expected Cell Counts by Using R
Bayesian Estimation of Expected Cell Counts by Using R Haydar Demirhan 1 and Canan Hamurkaroglu 2 Department of Statistics, Hacettepe University, Beytepe, 06800, Ankara, Turkey Abstract In this article,
More informationMarkov chain Monte Carlo methods for visual tracking
Markov chain Monte Carlo methods for visual tracking Ray Luo rluo@cory.eecs.berkeley.edu Department of Electrical Engineering and Computer Sciences University of California, Berkeley Berkeley, CA 94720
More informationarxiv: v1 [stat.co] 18 Feb 2012
A LEVEL-SET HIT-AND-RUN SAMPLER FOR QUASI-CONCAVE DISTRIBUTIONS Dean Foster and Shane T. Jensen arxiv:1202.4094v1 [stat.co] 18 Feb 2012 Department of Statistics The Wharton School University of Pennsylvania
More informationBayesian Estimation of DSGE Models 1 Chapter 3: A Crash Course in Bayesian Inference
1 The views expressed in this paper are those of the authors and do not necessarily reflect the views of the Federal Reserve Board of Governors or the Federal Reserve System. Bayesian Estimation of DSGE
More informationMCMC: Markov Chain Monte Carlo
I529: Machine Learning in Bioinformatics (Spring 2013) MCMC: Markov Chain Monte Carlo Yuzhen Ye School of Informatics and Computing Indiana University, Bloomington Spring 2013 Contents Review of Markov
More information17 : Markov Chain Monte Carlo
10-708: Probabilistic Graphical Models, Spring 2015 17 : Markov Chain Monte Carlo Lecturer: Eric P. Xing Scribes: Heran Lin, Bin Deng, Yun Huang 1 Review of Monte Carlo Methods 1.1 Overview Monte Carlo
More informationBayesian Methods with Monte Carlo Markov Chains II
Bayesian Methods with Monte Carlo Markov Chains II Henry Horng-Shing Lu Institute of Statistics National Chiao Tung University hslu@stat.nctu.edu.tw http://tigpbp.iis.sinica.edu.tw/courses.htm 1 Part 3
More informationBayesian modelling. Hans-Peter Helfrich. University of Bonn. Theodor-Brinkmann-Graduate School
Bayesian modelling Hans-Peter Helfrich University of Bonn Theodor-Brinkmann-Graduate School H.-P. Helfrich (University of Bonn) Bayesian modelling Brinkmann School 1 / 22 Overview 1 Bayesian modelling
More information(5) Multi-parameter models - Gibbs sampling. ST440/540: Applied Bayesian Analysis
Summarizing a posterior Given the data and prior the posterior is determined Summarizing the posterior gives parameter estimates, intervals, and hypothesis tests Most of these computations are integrals
More informationMarkov chain Monte Carlo
1 / 26 Markov chain Monte Carlo Timothy Hanson 1 and Alejandro Jara 2 1 Division of Biostatistics, University of Minnesota, USA 2 Department of Statistics, Universidad de Concepción, Chile IAP-Workshop
More informationBayesian Inference. Chapter 1. Introduction and basic concepts
Bayesian Inference Chapter 1. Introduction and basic concepts M. Concepción Ausín Department of Statistics Universidad Carlos III de Madrid Master in Business Administration and Quantitative Methods Master
More informationINTRODUCTION TO BAYESIAN STATISTICS
INTRODUCTION TO BAYESIAN STATISTICS Sarat C. Dass Department of Statistics & Probability Department of Computer Science & Engineering Michigan State University TOPICS The Bayesian Framework Different Types
More informationBayesian Prediction of Code Output. ASA Albuquerque Chapter Short Course October 2014
Bayesian Prediction of Code Output ASA Albuquerque Chapter Short Course October 2014 Abstract This presentation summarizes Bayesian prediction methodology for the Gaussian process (GP) surrogate representation
More informationBayesian data analysis in practice: Three simple examples
Bayesian data analysis in practice: Three simple examples Martin P. Tingley Introduction These notes cover three examples I presented at Climatea on 5 October 0. Matlab code is available by request to
More informationApril 20th, Advanced Topics in Machine Learning California Institute of Technology. Markov Chain Monte Carlo for Machine Learning
for for Advanced Topics in California Institute of Technology April 20th, 2017 1 / 50 Table of Contents for 1 2 3 4 2 / 50 History of methods for Enrico Fermi used to calculate incredibly accurate predictions
More informationBagging During Markov Chain Monte Carlo for Smoother Predictions
Bagging During Markov Chain Monte Carlo for Smoother Predictions Herbert K. H. Lee University of California, Santa Cruz Abstract: Making good predictions from noisy data is a challenging problem. Methods
More informationMarkov chain Monte Carlo
Markov chain Monte Carlo Markov chain Monte Carlo (MCMC) Gibbs and Metropolis Hastings Slice sampling Practical details Iain Murray http://iainmurray.net/ Reminder Need to sample large, non-standard distributions:
More informationBootstrap Approximation of Gibbs Measure for Finite-Range Potential in Image Analysis
Bootstrap Approximation of Gibbs Measure for Finite-Range Potential in Image Analysis Abdeslam EL MOUDDEN Business and Management School Ibn Tofaïl University Kenitra, Morocco Abstract This paper presents
More information1 Priors, Contd. ISyE8843A, Brani Vidakovic Handout Nuisance Parameters
ISyE8843A, Brani Vidakovic Handout 6 Priors, Contd.. Nuisance Parameters Assume that unknown parameter is θ, θ ) but we are interested only in θ. The parameter θ is called nuisance parameter. How do we
More informationMCMC Review. MCMC Review. Gibbs Sampling. MCMC Review
MCMC Review http://jackman.stanford.edu/mcmc/icpsr99.pdf http://students.washington.edu/fkrogsta/bayes/stat538.pdf http://www.stat.berkeley.edu/users/terry/classes/s260.1998 /Week9a/week9a/week9a.html
More informationBayesian Inference for DSGE Models. Lawrence J. Christiano
Bayesian Inference for DSGE Models Lawrence J. Christiano Outline State space-observer form. convenient for model estimation and many other things. Bayesian inference Bayes rule. Monte Carlo integation.
More informationMonte Carlo Methods in Bayesian Inference: Theory, Methods and Applications
University of Arkansas, Fayetteville ScholarWorks@UARK Theses and Dissertations 1-016 Monte Carlo Methods in Bayesian Inference: Theory, Methods and Applications Huarui Zhang University of Arkansas, Fayetteville
More informationMCMC Methods: Gibbs and Metropolis
MCMC Methods: Gibbs and Metropolis Patrick Breheny February 28 Patrick Breheny BST 701: Bayesian Modeling in Biostatistics 1/30 Introduction As we have seen, the ability to sample from the posterior distribution
More informationThe Pennsylvania State University The Graduate School RATIO-OF-UNIFORMS MARKOV CHAIN MONTE CARLO FOR GAUSSIAN PROCESS MODELS
The Pennsylvania State University The Graduate School RATIO-OF-UNIFORMS MARKOV CHAIN MONTE CARLO FOR GAUSSIAN PROCESS MODELS A Thesis in Statistics by Chris Groendyke c 2008 Chris Groendyke Submitted in
More informationSampling Algorithms for Probabilistic Graphical models
Sampling Algorithms for Probabilistic Graphical models Vibhav Gogate University of Washington References: Chapter 12 of Probabilistic Graphical models: Principles and Techniques by Daphne Koller and Nir
More informationHmms with variable dimension structures and extensions
Hmm days/enst/january 21, 2002 1 Hmms with variable dimension structures and extensions Christian P. Robert Université Paris Dauphine www.ceremade.dauphine.fr/ xian Hmm days/enst/january 21, 2002 2 1 Estimating
More informationIntroduction to Computational Biology Lecture # 14: MCMC - Markov Chain Monte Carlo
Introduction to Computational Biology Lecture # 14: MCMC - Markov Chain Monte Carlo Assaf Weiner Tuesday, March 13, 2007 1 Introduction Today we will return to the motif finding problem, in lecture 10
More informationSurveying the Characteristics of Population Monte Carlo
International Research Journal of Applied and Basic Sciences 2013 Available online at www.irjabs.com ISSN 2251-838X / Vol, 7 (9): 522-527 Science Explorer Publications Surveying the Characteristics of
More informationMarkov Chain Monte Carlo, Numerical Integration
Markov Chain Monte Carlo, Numerical Integration (See Statistics) Trevor Gallen Fall 2015 1 / 1 Agenda Numerical Integration: MCMC methods Estimating Markov Chains Estimating latent variables 2 / 1 Numerical
More informationAppendix 3. Markov Chain Monte Carlo and Gibbs Sampling
Appendix 3 Markov Chain Monte Carlo and Gibbs Sampling A constant theme in the development of statistics has been the search for justifications for what statisticians do Blasco (2001) Draft version 21
More information