Likelihood-free MCMC
|
|
- Brook Blair
- 5 years ago
- Views:
Transcription
1 Bayesian inference for stable distributions with applications in finance Department of Mathematics University of Leicester September 2, 2011 MSc project final presentation
2 Outline
3 Classical Monte Carlo Integration Bayes formula f (θ D) π (θ) P (D θ) f (θ D): posterior π (θ): prior P (D θ):likelihood Evaluating integrals Normalisation Z = π (θ) P (D θ) dθ Marginalisation f (θ D) = f (θ, x D) dx Expectation E f [h (θ)] = h (θ) f (θ D) dθ Suppose we can draw samples θ (j) f (θ D), j = 1,..., m E f [h (θ)] 1 m m j=1 h ( θ (j))
4 Rejection sampling 1 Sample x (i) g (x) 2 Accept x (i) with probability f (x(i) )/Mg(x (i) ), then go to 1
5 Importance sampling A different way to view E f [h (θ)] E f [h (θ)] = ˆ 1 m h (θ) f (θ) [ ] f (θ) g (θ) g (θ) dθ = E g g (θ) h (θ) m f ( θ (j)) ( g ( θ (j))h θ (j)) (1) j=1 for θ (j) drawn from g (θ) Importance sampling does not throw away samples, it gives different weights(importance) f ( θ (j) D ) /g ( θ (j)).
6 Importance sampling In Bayesian context with normalising constant not known: [ ] f (θ D) E f [h (θ)] = E g g (θ) h (θ) [ ] 1 π (θ) P (D θ) = P (D) E g h (θ) g (θ) m j=1 π ( θ (j)) P ( D θ (j)) h ( θ (j)) /g ( θ (j)) = m j=1 P ( D θ (j)) π ( θ (j)) /g ( θ (j)) where θ (j) g (θ) This can also be used in general setting, i.e. use m j=1 h(θ(j) )f (θ (j) )/g(θ (j) ) m as an alternative to (1), with an j=1 f (θ(j) )/g(θ (j) ) improvement in variance.
7 What if likelihoods are unavailable? Approximating posterior, avoid likelihood evaluation known as approximate Bayesian computation Some early literatures LF-RS Tavaré et al., 1997 Inferring Coalescence Times From DNA Sequence Data replacing the full dataset with summary statistics. Fu and Li, 1997 Estimating the age of the common ancestor of a sample of DNA sequences simulating a new dataset, comparing with the observed one. LF-MCMC Marjoram et al., 2003 Markov Chain Monte Carlo without likelihoods MCMC approach generalized from LF-RS
8 Likelihood-free rejection sampling The idea can be seen in the following algorithm. LF-RS 1 Simulate from the prior θ π 2 Generate D under the model with parameter θ 3 Accept θ if D = D ;go to 1 D : observed dataset D : simulated dataset In practice one replace D and D with corresponding summary statistics S and S.The condition can be rewritten as ρ (S, S ) ε for some distance measure ρ (e.g. Euclidean). This will result in an approximate posterior f (θ ρ (S, S ) ε).
9 LF-RS example Example Suppose y 1, y 2,... y n are observations from Exp (θ) with density f (y θ) = θe θy, y > 0. The prior for θ is conjugate gamma distribution θ Gamma (α, β), then the posterior is gamma with altered parameters θ D Gamma (n + α, β/ (β y i + 1)). Let α = 3, β = 1 and n = 5 observations from Exp (2), and choose the sample mean y as a sufficient statistic. We simulate the posterior distribution using LF-RS algorithm with ε = 1, and 0.1
10 LF-RS example Simulation results Results for ε = 1(left) and ε = 0.1(right)
11 Markov Chain Monte Carlo(MCMC) About MCMC Algorithms that realize Markov chain We want the invariant distribution of the chain to be our target distribution Samples can be taken as drawn from the target distribution after running the chain for a long time
12 Markov Chain Monte Carlo(MCMC) About MCMC Algorithms that realize Markov chain We want the invariant distribution of the chain to be our target distribution Samples can be taken as drawn from the target distribution after running the chain for a long time MCMC History Metropolis, et.al.(1953). Equations of state calculations by fast computing machines. J. Chem. Phys Hastings, W. (1970). Monte Carlo sampling methods using Markov chains and their application. Biometrika Gelfand, A. E. and Smith, A. F. M. (1990). Sampling-based approaches to calculating marginal densities. J. Amer. Statist. Assoc
13 Constructing MCMC algorithms Ergodic Theorem gaurantees convergence. From Markov chain theory, general balance implies(in discrete setting) f P = f (2) f : invariant distribution P: transition matrix with elements P ij = P (x t+1 = j x t = i) := P (i j). Sum over each row is one. Detailed balance P ( x x ) f (x) = P ( x x ) f ( x ) Summing both sides over x, we get (2).
14 Metropolis-Hastings algorithm Metropolis-Hastings 1 If now at θ, propose a move to θ according to a proposal distribution q (θ θ ) 2 Accept θ with probability A (θ, θ ) = min{1, f (θ )q(θ θ) f (θ)q(θ θ ) } 3 Go to 1 until desired number of iterations
15 Theorem The invariant distribution of the chain is f (θ). Proof. We show that detailed balance is satisfied. The M-H transition probability is Choose(w.l.o.g) Then P ( θ θ ) = q ( θ θ ) A ( θ, θ ) f (θ ) q (θ θ) f (θ) q (θ θ ) 1 P ( θ θ ) f (θ) = q ( θ θ ) f (θ ) q (θ θ) f (θ) q (θ θ ) f (θ) = f ( θ ) q ( θ θ ) A ( θ, θ ) = P ( θ θ ) f ( θ )
16 (Marjoram et al.,2003) proposed a MCMC method without likelihood evaluation LF-MCMC 1.If now at θ, propose a move to θ according to a proposal distribution q (θ θ ) 2.Generate D under model with θ 3.If D = D, go to 4; otherwise return to 1 4.Accept θ with probability A (θ, θ ) = min{1, π(θ )q(θ θ) π(θ)q(θ θ ) }, then go to 1 One can proof the invariant distribution is f (θ D). Approximate posterior: replacing D = D with ρ (S, S ) ε
17 Stable distributions No closed form densities in general.
18 Stable distributions No closed form densities in general. 4 parameters: α (0, 2] determins tail behavior, β [ 1, 1] the skewness, γ > 0 the scale and δ R the location.
19 Stable distributions No closed form densities in general. 4 parameters: α (0, 2] determins tail behavior, β [ 1, 1] the skewness, γ > 0 the scale and δ R the location. Special cases: Cauchy (α = 1, β = 0) Normal (α = 2, β = 0) Levy (α = 1/2, β = 1)
20 Stable distributions No closed form densities in general. 4 parameters: α (0, 2] determins tail behavior, β [ 1, 1] the skewness, γ > 0 the scale and δ R the location. Special cases: Cauchy (α = 1, β = 0) Normal (α = 2, β = 0) Levy (α = 1/2, β = 1) Infinite variance(except α = 2), mean is existed only if 1 < α 2.
21 Stable distributions No closed form densities in general. 4 parameters: α (0, 2] determins tail behavior, β [ 1, 1] the skewness, γ > 0 the scale and δ R the location. Special cases: Cauchy (α = 1, β = 0) Normal (α = 2, β = 0) Levy (α = 1/2, β = 1) Infinite variance(except α = 2), mean is existed only if 1 < α 2. Generalized CLT
22 Some literatures Bayesian inference for stable models Buckle, D.J., Bayesian inference for stable distributions. Journal of the American Statistical Association 90, Auxiliary variable Gibbs sampler Lombardi, M.J., Bayesian inference for alpha stable distributions: a random walk MCMC approach. Computational Statistics & Data Analysis 51, Evaluating likelihood via inverse Fourier transform combined with a series expansion Peters, G.W., Sisson, S.A., Fan, Y., Likelihood-free Bayesian inference for α-stable models. Computational Statistics and Data Analysis. doi: /j.csda Likelihood-free sequential Monte Carlo sampler
23 Before implementing LF-MCMC Assumptions of simulation: Estimate one parameter with the other three parameters fixed. * Use flat prior for the parameter to be estimated. Use a Gaussian transition kernel centered at current state. If the parameter is within some interval, simply truncate those values that are outside the interval. Use quantiles and Kolmogorov-Smirnov statistic as summary statistics. Use a fixed ε value during computation. * *: These assumptions will be dropped later.
24 Simulation results Fix parameters β, γ, δ Simulation results for α based on 200 observations from Stable (1.5, 0.5, 10, 10) using a fixed ε = 25. (Left) Sample path of α, true value is 1.5. (Right) Trace of sample average.
25 Simulation results Fix parameters β, γ, δ Sample path and ergodic average plot for α. Top: ε = 15, acceptance rate: 1.3% Bottom: ε = 50 acceptance rate:34.2%
26 Modified LF-MCMC Motivation Modification dynamically define ε t as a monotonically decreasing sequence: { max{ε min, min{ε, ε t 1 }} if accept θ ε t = otherwise ε t 1 ε 0 = ρ (S, S 0 ), ε = ρ (S, S ), where S 0 : summary statistics for the dataset generated by the intital value and ε min : target ε value. Before: compare with the target ε value(global comparison) Now: compare with the previous ε value(local comparison) adaptively change the variance of the proposal distribution accelerate/control chain mixing
27 Modified LF-MCMC Motivation Modification dynamically define ε t as a monotonically decreasing sequence: { max{ε min, min{ε, ε t 1 }} if accept θ ε t = otherwise ε t 1 ε 0 = ρ (S, S 0 ), ε = ρ (S, S ), where S 0 : summary statistics for the dataset generated by the intital value and ε min : target ε value. Before: compare with the target ε value(global comparison) Now: compare with the previous ε value(local comparison) adaptively change the variance of the proposal distribution accelerate/control chain mixing
28 Modified LF-MCMC Motivation Modification dynamically define ε t as a monotonically decreasing sequence: { max{ε min, min{ε, ε t 1 }} if accept θ ε t = otherwise ε t 1 ε 0 = ρ (S, S 0 ), ε = ρ (S, S ), where S 0 : summary statistics for the dataset generated by the intital value and ε min : target ε value. Before: compare with the target ε value(global comparison) Now: compare with the previous ε value(local comparison) adaptively change the variance of the proposal distribution accelerate/control chain mixing
29 Simulation results All four parameters unknown Simulation results for Stable (α, β, γ, δ) based on 500 observations from Stable (1.5, 0.5, 10, 10), using iterations and ε min = 15.
30 LF-MCMC for variance gamma (VG) distribution For comparison, we apply the method to VG distribution.
31 LF-MCMC for variance gamma (VG) distribution For comparison, we apply the method to VG distribution. VG process (Madan and Seneta, 1990) (Madan, Carr and Chang, 1998) X (VG) t = θg t + σw Gt G t is a gamma process with mean rate unity and variance rate ν, W t is the standard Brownian motion.
32 LF-MCMC for variance gamma (VG) distribution For comparison, we apply the method to VG distribution. VG process (Madan and Seneta, 1990) (Madan, Carr and Chang, 1998) X (VG) t = θg t + σw Gt G t is a gamma process with mean rate unity and variance rate ν, W t is the standard Brownian motion. Unit period distribution VG (σ, ν, θ) pdf can be written in terms of modified Bessel function of the second kind
33 LF-MCMC for variance gamma (VG) distribution For comparison, we apply the method to VG distribution. VG process (Madan and Seneta, 1990) (Madan, Carr and Chang, 1998) X (VG) t = θg t + σw Gt G t is a gamma process with mean rate unity and variance rate ν, W t is the standard Brownian motion. Unit period distribution VG (σ, ν, θ) pdf can be written in terms of modified Bessel function of the second kind VG distribution has finite moments of all order.
34 LF-MCMC for variance gamma (VG) distribution Simulation results Simulation results for VG (σ, ν, θ, µ) based on 500 observations from VG (0.8, 1, 0.5, 10), using iterations and ε min = 1. Added summary statistics: mean and variance.
35 Application to financial data Fit stable distribution to real financial data. The data is the S&P 500 index from the period of January 2009 to July 2011, with 629 daily log returns and the prices are adjusted close price. Implement iterations of LF-MCMC, discard first 2000 iterations, averaging over the samples gave the values of posterior estimates: α: β: γ: δ:
36 blue=stable fit, green=smoothed data The figure is produced using J.P. Nolan s STABLE program, available at jpnolan
37 Concluding remarks Our results: apply LF-MCMC to the inference for stable models
38 Concluding remarks Our results: apply LF-MCMC to the inference for stable models make the method applicable to general cases
39 Concluding remarks Our results: apply LF-MCMC to the inference for stable models make the method applicable to general cases relatively low computational cost
40 Concluding remarks Our results: Pitfalls: apply LF-MCMC to the inference for stable models make the method applicable to general cases relatively low computational cost need to specify a proper target ε value
41 Concluding remarks Our results: Pitfalls: apply LF-MCMC to the inference for stable models make the method applicable to general cases relatively low computational cost need to specify a proper target ε value don t know when convergence will happen, need more iterations
42 Concluding remarks Our results: Pitfalls: apply LF-MCMC to the inference for stable models make the method applicable to general cases relatively low computational cost need to specify a proper target ε value don t know when convergence will happen, need more iterations choice of summary statistics can crucially affect sampler performance
43 Acknowledgement Thanks to My supervisor Dr. Ray Kawai
F denotes cumulative density. denotes probability density function; (.)
BAYESIAN ANALYSIS: FOREWORDS Notation. System means the real thing and a model is an assumed mathematical form for the system.. he probability model class M contains the set of the all admissible models
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 informationMonte Carlo in Bayesian Statistics
Monte Carlo in Bayesian Statistics Matthew Thomas SAMBa - University of Bath m.l.thomas@bath.ac.uk December 4, 2014 Matthew Thomas (SAMBa) Monte Carlo in Bayesian Statistics December 4, 2014 1 / 16 Overview
More informationMarkov Chain Monte Carlo
1 Motivation 1.1 Bayesian Learning Markov Chain Monte Carlo Yale Chang In Bayesian learning, given data X, we make assumptions on the generative process of X by introducing hidden variables Z: p(z): prior
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 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 informationCalibration of Stochastic Volatility Models using Particle Markov Chain Monte Carlo Methods
Calibration of Stochastic Volatility Models using Particle Markov Chain Monte Carlo Methods Jonas Hallgren 1 1 Department of Mathematics KTH Royal Institute of Technology Stockholm, Sweden BFS 2012 June
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 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 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 informationEco517 Fall 2014 C. Sims MIDTERM EXAM
Eco57 Fall 204 C. Sims MIDTERM EXAM You have 90 minutes for this exam and there are a total of 90 points. The points for each question are listed at the beginning of the question. Answer all questions.
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 informationBayesian Inference and MCMC
Bayesian Inference and MCMC Aryan Arbabi Partly based on MCMC slides from CSC412 Fall 2018 1 / 18 Bayesian Inference - Motivation Consider we have a data set D = {x 1,..., x n }. E.g each x i can be the
More informationApproximate Bayesian Computation: a simulation based approach to inference
Approximate Bayesian Computation: a simulation based approach to inference Richard Wilkinson Simon Tavaré 2 Department of Probability and Statistics University of Sheffield 2 Department of Applied Mathematics
More informationBayesian Estimation with Sparse Grids
Bayesian Estimation with Sparse Grids Kenneth L. Judd and Thomas M. Mertens Institute on Computational Economics August 7, 27 / 48 Outline Introduction 2 Sparse grids Construction Integration with sparse
More informationPseudo-marginal MCMC methods for inference in latent variable models
Pseudo-marginal MCMC methods for inference in latent variable models Arnaud Doucet Department of Statistics, Oxford University Joint work with George Deligiannidis (Oxford) & Mike Pitt (Kings) MCQMC, 19/08/2016
More informationDown by the Bayes, where the Watermelons Grow
Down by the Bayes, where the Watermelons Grow A Bayesian example using SAS SUAVe: Victoria SAS User Group Meeting November 21, 2017 Peter K. Ott, M.Sc., P.Stat. Strategic Analysis 1 Outline 1. Motivating
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 informationComputationalToolsforComparing AsymmetricGARCHModelsviaBayes Factors. RicardoS.Ehlers
ComputationalToolsforComparing AsymmetricGARCHModelsviaBayes Factors RicardoS.Ehlers Laboratório de Estatística e Geoinformação- UFPR http://leg.ufpr.br/ ehlers ehlers@leg.ufpr.br II Workshop on Statistical
More informationTutorial on ABC Algorithms
Tutorial on ABC Algorithms Dr Chris Drovandi Queensland University of Technology, Australia c.drovandi@qut.edu.au July 3, 2014 Notation Model parameter θ with prior π(θ) Likelihood is f(ý θ) with observed
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 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 informationStat 535 C - Statistical Computing & Monte Carlo Methods. Lecture February Arnaud Doucet
Stat 535 C - Statistical Computing & Monte Carlo Methods Lecture 13-28 February 2006 Arnaud Doucet Email: arnaud@cs.ubc.ca 1 1.1 Outline Limitations of Gibbs sampling. Metropolis-Hastings algorithm. Proof
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 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 informationLecture 7 and 8: Markov Chain Monte Carlo
Lecture 7 and 8: Markov Chain Monte Carlo 4F13: Machine Learning Zoubin Ghahramani and Carl Edward Rasmussen Department of Engineering University of Cambridge http://mlg.eng.cam.ac.uk/teaching/4f13/ Ghahramani
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 informationLearning the hyper-parameters. Luca Martino
Learning the hyper-parameters Luca Martino 2017 2017 1 / 28 Parameters and hyper-parameters 1. All the described methods depend on some choice of hyper-parameters... 2. For instance, do you recall λ (bandwidth
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 Phylogenetics:
Bayesian Phylogenetics: an introduction Marc A. Suchard msuchard@ucla.edu UCLA Who is this man? How sure are you? The one true tree? Methods we ve learned so far try to find a single tree that best describes
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 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 informationSAMPLING ALGORITHMS. In general. Inference in Bayesian models
SAMPLING ALGORITHMS SAMPLING ALGORITHMS In general A sampling algorithm is an algorithm that outputs samples x 1, x 2,... from a given distribution P or density p. Sampling algorithms can for example be
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 informationBayesian Estimation of Input Output Tables for Russia
Bayesian Estimation of Input Output Tables for Russia Oleg Lugovoy (EDF, RANE) Andrey Polbin (RANE) Vladimir Potashnikov (RANE) WIOD Conference April 24, 2012 Groningen Outline Motivation Objectives Bayesian
More informationNonparametric Bayesian Methods - Lecture I
Nonparametric Bayesian Methods - Lecture I Harry van Zanten Korteweg-de Vries Institute for Mathematics CRiSM Masterclass, April 4-6, 2016 Overview of the lectures I Intro to nonparametric Bayesian statistics
More informationA Review of Pseudo-Marginal Markov Chain Monte Carlo
A Review of Pseudo-Marginal Markov Chain Monte Carlo Discussed by: Yizhe Zhang October 21, 2016 Outline 1 Overview 2 Paper review 3 experiment 4 conclusion Motivation & overview Notation: θ denotes the
More informationDefault Priors and Effcient Posterior Computation in Bayesian
Default Priors and Effcient Posterior Computation in Bayesian Factor Analysis January 16, 2010 Presented by Eric Wang, Duke University Background and Motivation A Brief Review of Parameter Expansion Literature
More informationBAYESIAN METHODS FOR VARIABLE SELECTION WITH APPLICATIONS TO HIGH-DIMENSIONAL DATA
BAYESIAN METHODS FOR VARIABLE SELECTION WITH APPLICATIONS TO HIGH-DIMENSIONAL DATA Intro: Course Outline and Brief Intro to Marina Vannucci Rice University, USA PASI-CIMAT 04/28-30/2010 Marina Vannucci
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 informationEco517 Fall 2013 C. Sims MCMC. October 8, 2013
Eco517 Fall 2013 C. Sims MCMC October 8, 2013 c 2013 by Christopher A. Sims. This document may be reproduced for educational and research purposes, so long as the copies contain this notice and are retained
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 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 informationCS242: Probabilistic Graphical Models Lecture 7B: Markov Chain Monte Carlo & Gibbs Sampling
CS242: Probabilistic Graphical Models Lecture 7B: Markov Chain Monte Carlo & Gibbs Sampling Professor Erik Sudderth Brown University Computer Science October 27, 2016 Some figures and materials courtesy
More informationPattern Recognition and Machine Learning. Bishop Chapter 11: Sampling Methods
Pattern Recognition and Machine Learning Chapter 11: Sampling Methods Elise Arnaud Jakob Verbeek May 22, 2008 Outline of the chapter 11.1 Basic Sampling Algorithms 11.2 Markov Chain Monte Carlo 11.3 Gibbs
More informationSC7/SM6 Bayes Methods HT18 Lecturer: Geoff Nicholls Lecture 2: Monte Carlo Methods Notes and Problem sheets are available at http://www.stats.ox.ac.uk/~nicholls/bayesmethods/ and via the MSc weblearn pages.
More informationAdvanced Statistical Modelling
Markov chain Monte Carlo (MCMC) Methods and Their Applications in Bayesian Statistics School of Technology and Business Studies/Statistics Dalarna University Borlänge, Sweden. Feb. 05, 2014. Outlines 1
More informationHierarchical models. Dr. Jarad Niemi. August 31, Iowa State University. Jarad Niemi (Iowa State) Hierarchical models August 31, / 31
Hierarchical models Dr. Jarad Niemi Iowa State University August 31, 2017 Jarad Niemi (Iowa State) Hierarchical models August 31, 2017 1 / 31 Normal hierarchical model Let Y ig N(θ g, σ 2 ) for i = 1,...,
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 informationRiemann Manifold Methods in Bayesian Statistics
Ricardo Ehlers ehlers@icmc.usp.br Applied Maths and Stats University of São Paulo, Brazil Working Group in Statistical Learning University College Dublin September 2015 Bayesian inference is based on Bayes
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 informationDeblurring Jupiter (sampling in GLIP faster than regularized inversion) Colin Fox Richard A. Norton, J.
Deblurring Jupiter (sampling in GLIP faster than regularized inversion) Colin Fox fox@physics.otago.ac.nz Richard A. Norton, J. Andrés Christen Topics... Backstory (?) Sampling in linear-gaussian hierarchical
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 Semiparametric GARCH Models
Bayesian Semiparametric GARCH Models Xibin (Bill) Zhang and Maxwell L. King Department of Econometrics and Business Statistics Faculty of Business and Economics xibin.zhang@monash.edu Quantitative Methods
More informationBayesian model selection in graphs by using BDgraph package
Bayesian model selection in graphs by using BDgraph package A. Mohammadi and E. Wit March 26, 2013 MOTIVATION Flow cytometry data with 11 proteins from Sachs et al. (2005) RESULT FOR CELL SIGNALING DATA
More informationBayesian Regression Linear and Logistic Regression
When we want more than point estimates Bayesian Regression Linear and Logistic Regression Nicole Beckage Ordinary Least Squares Regression and Lasso Regression return only point estimates But what if we
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 informationCSC 2541: Bayesian Methods for Machine Learning
CSC 2541: Bayesian Methods for Machine Learning Radford M. Neal, University of Toronto, 2011 Lecture 10 Alternatives to Monte Carlo Computation Since about 1990, Markov chain Monte Carlo has been the dominant
More informationBayesian Semiparametric GARCH Models
Bayesian Semiparametric GARCH Models Xibin (Bill) Zhang and Maxwell L. King Department of Econometrics and Business Statistics Faculty of Business and Economics xibin.zhang@monash.edu Quantitative Methods
More informationSTAT 425: Introduction to Bayesian Analysis
STAT 425: Introduction to Bayesian Analysis Marina Vannucci Rice University, USA Fall 2017 Marina Vannucci (Rice University, USA) Bayesian Analysis (Part 2) Fall 2017 1 / 19 Part 2: Markov chain Monte
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 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 informationBayesian Inference for Discretely Sampled Diffusion Processes: A New MCMC Based Approach to Inference
Bayesian Inference for Discretely Sampled Diffusion Processes: A New MCMC Based Approach to Inference Osnat Stramer 1 and Matthew Bognar 1 Department of Statistics and Actuarial Science, University of
More informationKernel Sequential Monte Carlo
Kernel Sequential Monte Carlo Ingmar Schuster (Paris Dauphine) Heiko Strathmann (University College London) Brooks Paige (Oxford) Dino Sejdinovic (Oxford) * equal contribution April 25, 2016 1 / 37 Section
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 informationOnline appendix to On the stability of the excess sensitivity of aggregate consumption growth in the US
Online appendix to On the stability of the excess sensitivity of aggregate consumption growth in the US Gerdie Everaert 1, Lorenzo Pozzi 2, and Ruben Schoonackers 3 1 Ghent University & SHERPPA 2 Erasmus
More informationKernel adaptive Sequential Monte Carlo
Kernel adaptive Sequential Monte Carlo Ingmar Schuster (Paris Dauphine) Heiko Strathmann (University College London) Brooks Paige (Oxford) Dino Sejdinovic (Oxford) December 7, 2015 1 / 36 Section 1 Outline
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 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 informationInference in state-space models with multiple paths from conditional SMC
Inference in state-space models with multiple paths from conditional SMC Sinan Yıldırım (Sabancı) joint work with Christophe Andrieu (Bristol), Arnaud Doucet (Oxford) and Nicolas Chopin (ENSAE) September
More information(I AL BL 2 )z t = (I CL)ζ t, where
ECO 513 Fall 2011 MIDTERM EXAM The exam lasts 90 minutes. Answer all three questions. (1 Consider this model: x t = 1.2x t 1.62x t 2 +.2y t 1.23y t 2 + ε t.7ε t 1.9ν t 1 (1 [ εt y t = 1.4y t 1.62y t 2
More informationExercises Tutorial at ICASSP 2016 Learning Nonlinear Dynamical Models Using Particle Filters
Exercises Tutorial at ICASSP 216 Learning Nonlinear Dynamical Models Using Particle Filters Andreas Svensson, Johan Dahlin and Thomas B. Schön March 18, 216 Good luck! 1 [Bootstrap particle filter for
More informationLikelihood-free Markov chain Monte
arxiv:1001.2058v1 [stat.me] 13 Jan 2010 Chapter 1 Likelihood-free Markov chain Monte Carlo 1.1 Introduction Scott A. Sisson and Yanan Fan In Bayesian inference, the posterior distribution for parameters
More informationEstimation of Operational Risk Capital Charge under Parameter Uncertainty
Estimation of Operational Risk Capital Charge under Parameter Uncertainty Pavel V. Shevchenko Principal Research Scientist, CSIRO Mathematical and Information Sciences, Sydney, Locked Bag 17, North Ryde,
More informationMCMC for big data. Geir Storvik. BigInsight lunch - May Geir Storvik MCMC for big data BigInsight lunch - May / 17
MCMC for big data Geir Storvik BigInsight lunch - May 2 2018 Geir Storvik MCMC for big data BigInsight lunch - May 2 2018 1 / 17 Outline Why ordinary MCMC is not scalable Different approaches for making
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 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 informationPractical Bayesian Quantile Regression. Keming Yu University of Plymouth, UK
Practical Bayesian Quantile Regression Keming Yu University of Plymouth, UK (kyu@plymouth.ac.uk) A brief summary of some recent work of us (Keming Yu, Rana Moyeed and Julian Stander). Summary We develops
More informationan introduction to bayesian inference
with an application to network analysis http://jakehofman.com january 13, 2010 motivation would like models that: provide predictive and explanatory power are complex enough to describe observed phenomena
More informationCTDL-Positive Stable Frailty Model
CTDL-Positive Stable Frailty Model M. Blagojevic 1, G. MacKenzie 2 1 Department of Mathematics, Keele University, Staffordshire ST5 5BG,UK and 2 Centre of Biostatistics, University of Limerick, Ireland
More informationApproximate Bayesian computation for spatial extremes via open-faced sandwich adjustment
Approximate Bayesian computation for spatial extremes via open-faced sandwich adjustment Ben Shaby SAMSI August 3, 2010 Ben Shaby (SAMSI) OFS adjustment August 3, 2010 1 / 29 Outline 1 Introduction 2 Spatial
More informationAnswers and expectations
Answers and expectations For a function f(x) and distribution P(x), the expectation of f with respect to P is The expectation is the average of f, when x is drawn from the probability distribution P E
More informationProbabilistic Graphical Models Lecture 17: Markov chain Monte Carlo
Probabilistic Graphical Models Lecture 17: Markov chain Monte Carlo Andrew Gordon Wilson www.cs.cmu.edu/~andrewgw Carnegie Mellon University March 18, 2015 1 / 45 Resources and Attribution Image credits,
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 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 information27 : Distributed Monte Carlo Markov Chain. 1 Recap of MCMC and Naive Parallel Gibbs Sampling
10-708: Probabilistic Graphical Models 10-708, Spring 2014 27 : Distributed Monte Carlo Markov Chain Lecturer: Eric P. Xing Scribes: Pengtao Xie, Khoa Luu In this scribe, we are going to review the Parallel
More informationInexact approximations for doubly and triply intractable problems
Inexact approximations for doubly and triply intractable problems March 27th, 2014 Markov random fields Interacting objects Markov random fields (MRFs) are used for modelling (often large numbers of) interacting
More informationeqr094: 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 informationData Analysis I. Dr Martin Hendry, Dept of Physics and Astronomy University of Glasgow, UK. 10 lectures, beginning October 2006
Astronomical p( y x, I) p( x, I) p ( x y, I) = p( y, I) Data Analysis I Dr Martin Hendry, Dept of Physics and Astronomy University of Glasgow, UK 10 lectures, beginning October 2006 4. Monte Carlo Methods
More informationSupplement to A Hierarchical Approach for Fitting Curves to Response Time Measurements
Supplement to A Hierarchical Approach for Fitting Curves to Response Time Measurements Jeffrey N. Rouder Francis Tuerlinckx Paul L. Speckman Jun Lu & Pablo Gomez May 4 008 1 The Weibull regression model
More informationModeling conditional distributions with mixture models: Applications in finance and financial decision-making
Modeling conditional distributions with mixture models: Applications in finance and financial decision-making John Geweke University of Iowa, USA Journal of Applied Econometrics Invited Lecture Università
More informationMCMC and Gibbs Sampling. Kayhan Batmanghelich
MCMC and Gibbs Sampling Kayhan Batmanghelich 1 Approaches to inference l Exact inference algorithms l l l The elimination algorithm Message-passing algorithm (sum-product, belief propagation) The junction
More informationWeb Appendix for Hierarchical Adaptive Regression Kernels for Regression with Functional Predictors by D. B. Woodard, C. Crainiceanu, and D.
Web Appendix for Hierarchical Adaptive Regression Kernels for Regression with Functional Predictors by D. B. Woodard, C. Crainiceanu, and D. Ruppert A. EMPIRICAL ESTIMATE OF THE KERNEL MIXTURE Here we
More informationModeling conditional distributions with mixture models: Theory and Inference
Modeling conditional distributions with mixture models: Theory and Inference John Geweke University of Iowa, USA Journal of Applied Econometrics Invited Lecture Università di Venezia Italia June 2, 2005
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 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 informationControlled sequential Monte Carlo
Controlled sequential Monte Carlo Jeremy Heng, Department of Statistics, Harvard University Joint work with Adrian Bishop (UTS, CSIRO), George Deligiannidis & Arnaud Doucet (Oxford) Bayesian Computation
More informationAdvances and Applications in Perfect Sampling
and Applications in Perfect Sampling Ph.D. Dissertation Defense Ulrike Schneider advisor: Jem Corcoran May 8, 2003 Department of Applied Mathematics University of Colorado Outline Introduction (1) MCMC
More informationLecture 6: Markov Chain Monte Carlo
Lecture 6: Markov Chain Monte Carlo D. Jason Koskinen koskinen@nbi.ku.dk Photo by Howard Jackman University of Copenhagen Advanced Methods in Applied Statistics Feb - Apr 2016 Niels Bohr Institute 2 Outline
More informationAn ABC interpretation of the multiple auxiliary variable method
School of Mathematical and Physical Sciences Department of Mathematics and Statistics Preprint MPS-2016-07 27 April 2016 An ABC interpretation of the multiple auxiliary variable method by Dennis Prangle
More information