Particle Markov Chain Monte Carlo Methods in Marine Biogeochemistry. Lawrence Murray and Emlyn Jones CSIRO Mathematics, Informatics and Statistics
|
|
- Collin Wells
- 5 years ago
- Views:
Transcription
1 Particle Markov Chain Monte Carlo Methods in Marine Biogeochemistry Lawrence Murray and Emlyn Jones CSIRO Mathematics, Informatics and Statistics
2 Outline Case studies in marine biogeochemistry and generic challenges. Conventional state-space models Collapsed state-space models (to deal with the absence of a closed-form transition density) The particle marginal Metropolis-Hastings (PMMH) sampler in this context. Results and implications. Acknowledgements: Eddy Campbell (CSIRO), John Parslow (CSIRO), Nugzar Margvelashvili (CSIRO), Noel Cressie (Ohio State University). Lawrence Murray PMCMC in marine biogeochemistry : 2 of 4
3 Marine biogeochemical model Sea Surface Phytoplankton Growth P (Phytoplankton) N (DIN) Zooplankton Grazing Remineralisation Z (Zooplankton) D (Detritus) Zooplankton Mortality Messy Feeding Base of Mixed Layer Mixing Sinking Lawrence Murray PMCMC in marine biogeochemistry : 3 of 4
4 Issues, issues, issues... These models tend to have long memory and don t mix well. Available observations may be sparse. The model may be strongly nonlinear,......it might even be chaotic. The transition density is unlikely to have a closed form. Lawrence Murray PMCMC in marine biogeochemistry : 4 of 4
5 Conventional state-space model U U 2... U T X X X 2... X T Y Y 2... Y T Lawrence Murray PMCMC in marine biogeochemistry : 5 of 4
6 Conventional state-space model Sampling with sequential Monte Carlo (auxiliary particle filter)... For particle i at time t, extend x i t by drawing xi t q(x t ) and weight with: wt i = p(y t x i t)p(x i t x i t ) q(x i wt. i t) Lawrence Murray PMCMC in marine biogeochemistry : 6 of 4
7 Conventional state-space model Sampling with sequential Monte Carlo (auxiliary particle filter)... For particle i at time t, extend x i t by drawing xi t q(x t ) and weight with: wt i = p(y t x i t)p(x i t x i t ) q(x i wt. i t) We don t have that! Lawrence Murray PMCMC in marine biogeochemistry : 7 of 4
8 Conventional state-space model U U 2... U T X X X 2... X T Y Y 2... Y T Lawrence Murray PMCMC in marine biogeochemistry : 8 of 4
9 Collapsed state-space model... U U 2 U T X Y Y 2... Y T Lawrence Murray PMCMC in marine biogeochemistry : 9 of 4
10 Collapsed state-space model Sampling with sequential Monte Carlo (auxiliary particle filter)... For particle i at time t, extend u i :t by drawing ui t q(u t ) and weight with: wt i = p(y t u i :t, x i )p(u i t u i :t, x i ) q(u i wt. i t) Lawrence Murray PMCMC in marine biogeochemistry : of 4
11 Collapsed state-space model Sampling with sequential Monte Carlo (auxiliary particle filter)... For particle i at time t, extend u i :t by drawing ui t q(u t ) and weight with: wt i = p(y t u i :t, x i )p(u i t) q(u i wt. i t) We do have that! Lawrence Murray PMCMC in marine biogeochemistry : of 4
12 What should q( ) be? PF: bootstrap PF: as PF plus one-step single-pilot lookahead and resample. MUPF: UKF run offline, use time marginals p N (u t y :t ) as proposals. MUPF: as MUPF plus one-step single-pilot lookahead and resample. CUPF: UKF conditioned on each particle, so use p N (u t u :t, y :t ). CUPF: as CUPF plus UKF lookahead and resample. (UKF = Unscented Kalman Filter) Lawrence Murray PMCMC in marine biogeochemistry : 2 of 4
13 Particle marginal Metropolis-Hastings The target (posterior) density is π(u :T, x, θ y :T ), factorised as either: or π (u :T, x θ, y :T )π 2 (θ y :T ) π (u :T x, θ, y :T )π 2 (x, θ y :T ). In either case π is targeted with an auxiliary particle filter, π 2 with Metropolis-Hastings [Andrieu et al., 2, Pitt et al., 2]. The second factorisation should be preferred when prior information over x is scarce, so that importance sampling of it will be degenerate. Lawrence Murray PMCMC in marine biogeochemistry : 3 of 4
14 PZ model Simple phytoplankton-zooplankton model [Jones et al., 2]. Lotka-Volterra with stochastic growth rate and quadratic mortality term. dp dt dz dt α t = α t P cp Z = ecp Z m l Z m q Z 2 N (µ, σ). P is observed with noise. Simulated data used here. Lawrence Murray PMCMC in marine biogeochemistry : 4 of 4
15 PZ model: state estimate Z Prior Posterior Observed Truth P α t t Lawrence Murray PMCMC in marine biogeochemistry : 5 of 4
16 NPZD model The model features nine noise terms, ξ i for i =,..., 9, each coupled to a univariate autoregressive process B i. B i (t + t) = B i (t) ( t/τ P ) + (µ i + P DF σ i ξ i ) t/τ P, where t is a discrete time step (one day), µ i a parameter to be estimated, P DF a common diversity factor to be estimated, σ i a prescribed scaling factor, and τ P a common characteristic time scale, also prescribed. Lawrence Murray PMCMC in marine biogeochemistry : 6 of 4
17 NPZD model A multiplicative temperature correction T c is applied to all rate processes, for which a Q formulation for dependence on temperature, T, is used: T c = Q (T T ref)/, where T ref is a reference temperature, and Q a prescribed constant. Lawrence Murray PMCMC in marine biogeochemistry : 7 of 4
18 NPZD model The zooplankton grazing rate (gr) is dependent on the phytoplankton concentration (zooplankton functional response): gr = T c I Z A υ ( + A υ ), () where υ is a given power. Lawrence Murray PMCMC in marine biogeochemistry : 8 of 4
19 NPZD model The relative availability of phytoplankton, A, is A = Cl Z P I Z ; I Z is the maximum zooplankton ingestion rate (mg P per mg Z per day); Cl Z is the maximum clearance rate (volume in m 3 swept clear per mg Z per day). For υ =, this takes the form of a Type-2 functional response (standard rectangular hyperbola) [Holling, 966], and for υ > a Type-3 sigmoid functional response. Lawrence Murray PMCMC in marine biogeochemistry : 9 of 4
20 NPZD model A quadratic formulation for zooplankton mortality is adopted after Steele [976] and Steele and Henderson [992]: m = T c m Q Z, where the quadratic mortality rate m Q has units of d (mgzm 3 ). Lawrence Murray PMCMC in marine biogeochemistry : 2 of 4
21 NPZD model The detrital remineralization rate is dependent only on temperature: r = T c r D, where r D prescribes the remineralisation rate at a reference temperature. Lawrence Murray PMCMC in marine biogeochemistry : 2 of 4
22 NPZD model The phytoplankton specific growth rate, g, depends on temperature, T, available light or irradiance, E, and dissolved inorganic nutrient, N: g = T c g max h E h N /(h E + h N ). Lawrence Murray PMCMC in marine biogeochemistry : 22 of 4
23 NPZD model The light-limitation factor is given by h E = exp( α λ max E/g max ), where α is the initial slope of the photosynthesis versus irradiance curve (mg C mg Chla mol photon m 2 ), and λ max is the maximum Chla : C ratio (mg Chla mg C ). Lawrence Murray PMCMC in marine biogeochemistry : 23 of 4
24 NPZD model E is the mean photosynthetic available radiation (PAR) in the mixed layer and is given by E = E.( exp( Kz))/Kz, where E is the mean daily photosynthetically available radiation (PAR) just below the air-sea interface and Kz is given by Kz = (K W + a Ch Chla) MLD. (2) K W is attenuation due to the seawater and a Ch. Lawrence Murray PMCMC in marine biogeochemistry : 24 of 4
25 NPZD model The nutrient-limitation factor is given by h N = N (g max T c/a N ) + N, where a N is the maximum specific affinity for nitrogen uptake (d mg N m 3 ). Lawrence Murray PMCMC in marine biogeochemistry : 25 of 4
26 NPZD model The phytoplankton N : C ratio, χ, predicted by the model is given by χ = χmin h E + χ max h N h E + h N, where χ min and χ max are the prescribed minimum and maximum N : C ratios (mg N mg C ). Lawrence Murray PMCMC in marine biogeochemistry : 26 of 4
27 NPZD model The equations governing interactions between the remaining state variables {P, Z, D, N} are: dp = g P gr Z + κ dt MLD (BCP P ) dz = E Z gr Z m Z dt dd dt dn dt = ( E Z ) f D gr Z + m Z r D S D κ (BCD D) MLD = g P + ( E Z ) ( f D ) gr Z + r D + κ (BCN N). MLD D MLD + Lawrence Murray PMCMC in marine biogeochemistry : 27 of 4
28 State estimation (observed) 35 3 DIN (µ g N l ) ln(chla (µ g Chla l )) Prior 95% Posterior 95% Observations Prior median Posterior median Lawrence Murray PMCMC in marine biogeochemistry : 28 of 4
29 State estimation (unobserved) 6 ln(p(µ g N l ) ln(z(µ g N l ) ln(d(µ g N l ) Prior 95% Posterior 95% Prior median Posterior median Lawrence Murray PMCMC in marine biogeochemistry : 29 of 4
30 Parameter estimation ln(k ) W ln(a ) Ch s D ln(f D ) ln(pdf) ln(zdf) ln(µ g max) ln(µ λ max) ln(µ ) RN ln(µ ) an 2 3 ln(µ ) IZ ln(µ ) ClZ ln(µ ) EZ 3 2 ln(µ ) rd ln(µ ) mq Posterior Prior Lawrence Murray PMCMC in marine biogeochemistry : 3 of 4
31 State forecast 25 6 ) DIN(µ g N l )) ln(p(µ g N l /74 2/75 9/75 4 7/74 2/75 9/ )) ln(z(µ g N l /74 2/75 9/75 )) ln(d(µ g N l /74 2/75 9/75 )) ln(chla(µ g N l /74 2/75 9/75 Prior 95% Posterior 95% Forecast 95% Observations Prior median Posterior median Forecast median Lawrence Murray PMCMC in marine biogeochemistry : 3 of 4
32 PZ model: convergence.2. PF MUPF CUPF PF MUPF CUPF.2. PF MUPF CUPF PF MUPF CUPF.8.8 R p Step Step Convergence rates of Markov chains for the PZ case study, (left) particle-matched, and (right) compute-matched. Each line shows the evolution of the ˆR p statistic of Brooks and Gelman [998] for a particular method. Lawrence Murray PMCMC in marine biogeochemistry : 32 of 4
33 NPZD model: convergence PF MUPF CUPF PF MUPF CUPF PF MUPF CUPF PF MUPF CUPF Step Step Convergence rates of Markov chains for the NPZD case study, (left) particle-matched, and (right) compute-matched. Each line shows the evolution of the ˆR p statistic of Brooks and Gelman [998] for a particular method. Lawrence Murray PMCMC in marine biogeochemistry : 33 of 4
34 PZ model: acceptance rates σ PF MUPF CUPF σ PF MUPF CUPF µ µ µ Lawrence Murray PMCMC in marine biogeochemistry : 34 of 4
35 PZ model: acceptance rates, compute-matched σ PF MUPF CUPF σ PF MUPF CUPF µ µ µ Lawrence Murray PMCMC in marine biogeochemistry : 35 of 4
36 NPZD model: acceptance rates.8.8 Cumulative density.6.4 Cumulative density.6.4 PF.2 MUPF CUPF PF MUPF CUPF Conditional acceptance rate (CAR) PF.2 MUPF CUPF PF MUPF CUPF Conditional acceptance rate (CAR) Empirical cdfs of the acceptance rates for each method on the NPZD case study, (left) particle-matched, and (right) compute-matched. Lawrence Murray PMCMC in marine biogeochemistry : 36 of 4
37 NPZD model: parameters vs acceptance rates K W a Ch S D f D PDF ZDF µ g max µ λ max µ RN µ an µ IZ µ ClZ µ EZ µ rd µ mq Lawrence Murray PMCMC in marine biogeochemistry : 37 of 4
38 Lessons for PMMH Initialisation is extremely important, sampler will fail if initialised in a region where the acceptance rate is low. A particularly informative prior distribution may bias the posterior into regions where the acceptance rate is low. Is the rule-of-thumb 23% acceptance rate appropriate in this context? Lawrence Murray PMCMC in marine biogeochemistry : 38 of 4
39 Summary Environmental models such as these in marine biogeochemistry pose significant challenges to Monte Carlo methodology. One particular challenge is the absence of a closed-form transition density. The collapsed state-space model facilitates one means of improving sampler performance without needing this. The particle marginal Metropolis-Hastings sampler can have significant mixing and acceptance rate problems, initialisation and careful consideration of the prior distribution are important. Lawrence Murray PMCMC in marine biogeochemistry : 39 of 4
40 References C. Andrieu, A. Doucet, and R. Holenstein. Particle Markov chain Monte Carlo methods. Journal of the Royal Statistical Society Series B, 72:269 32, 2. S. P. Brooks and A. Gelman. General methods for monitoring convergence of iterative simulations. Journal of Computational and Graphical Statistics, 7: , 998. C. S. Holling. The functional response of predators to prey density and its role in mimicry and population regulation. Memoirs of the Entomology Society of Canada, 45:6, 966. E. Jones, J. Parslow, and L. M. Murray. A Bayesian approach to state and parameter estimation in a phytoplankton-zooplankton model. Australian Meteorological and Oceanographic Journal, 59(SP):7 6, 2. M. K. Pitt, R. S. Silva, P. Giordani, and R. Kohn. Auxiliary particle filtering within adaptive Metropolis-Hastings sampling. 2. URL J. Steele. Role of predation in ecosystem models. Marine Biology, 35():9, 976. ISSN J. Steele and E. Henderson. The role of predation in plankton models. Journal of Plankton Research, 4():57 72, 992. Lawrence Murray PMCMC in marine biogeochemistry : 4 of 4
41 CSIRO Mathematics, Informatics and Statistics Lawrence Murray Phone: lawrence.murray@csiro.au Web: Contact Us Phone: or enquiries@csiro.au Web:
arxiv: v1 [stat.ap] 7 Nov 2012
Bayesian Learning and Predictability in a Stochastic Nonlinear Dynamical Model John Parslow a, Noel Cressie b, Edward P. Campbell c, Emlyn Jones a, Lawrence Murray c arxiv:.77v [stat.ap] 7 Nov a CSIRO
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 informationECE276A: Sensing & Estimation in Robotics Lecture 10: Gaussian Mixture and Particle Filtering
ECE276A: Sensing & Estimation in Robotics Lecture 10: Gaussian Mixture and Particle Filtering Lecturer: Nikolay Atanasov: natanasov@ucsd.edu Teaching Assistants: Siwei Guo: s9guo@eng.ucsd.edu Anwesan Pal:
More informationBiogeochemical modelling and data assimilation: status in Australia
Biogeochemical modelling and data assimilation: status in Australia Richard Matear, Andrew Lenton, Matt Chamberlain, Mathieu Mongin, Emlyn Jones, Mark Baird www.cmar.csiro.au/staff/oke/ Biogeochemical
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 informationMH I. Metropolis-Hastings (MH) algorithm is the most popular method of getting dependent samples from a probability distribution
MH I Metropolis-Hastings (MH) algorithm is the most popular method of getting dependent samples from a probability distribution a lot of Bayesian mehods rely on the use of MH algorithm and it s famous
More informationAn Introduction to Particle Filtering
An Introduction to Particle Filtering Author: Lisa Turner Supervisor: Dr. Christopher Sherlock 10th May 2013 Abstract This report introduces the ideas behind particle filters, looking at the Kalman filter
More informationSequential Monte Carlo Methods (for DSGE Models)
Sequential Monte Carlo Methods (for DSGE Models) Frank Schorfheide University of Pennsylvania, PIER, CEPR, and NBER October 23, 2017 Some References These lectures use material from our joint work: Tempered
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 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 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 informationF 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 informationEstimation of positive sum-to-one constrained parameters with ensemble-based Kalman filters: application to an ocean ecosystem model
Estimation of positive sum-to-one constrained parameters with ensemble-based Kalman filters: application to an ocean ecosystem model Ehouarn Simon 1, Annette Samuelsen 1, Laurent Bertino 1 Dany Dumont
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 informationAn introduction to Sequential Monte Carlo
An introduction to Sequential Monte Carlo Thang Bui Jes Frellsen Department of Engineering University of Cambridge Research and Communication Club 6 February 2014 1 Sequential Monte Carlo (SMC) methods
More informationfor CESM Jessica Luo, Matthew Long, Keith Lindsay, Mike Levy NCAR Climate and Global Dynamics OMWG / BGC Working Group Meeting, Jan 12, 2018
Constructing a sizestructured plankton model for CESM Jessica Luo, Matthew Long, Keith Lindsay, Mike Levy NCAR Climate and Global Dynamics OMWG / BGC Working Group Meeting, Jan 12, 2018 Challenge: predicting
More informationThe Unscented Particle Filter
The Unscented Particle Filter Rudolph van der Merwe (OGI) Nando de Freitas (UC Bereley) Arnaud Doucet (Cambridge University) Eric Wan (OGI) Outline Optimal Estimation & Filtering Optimal Recursive Bayesian
More informationAuxiliary Particle Methods
Auxiliary Particle Methods Perspectives & Applications Adam M. Johansen 1 adam.johansen@bristol.ac.uk Oxford University Man Institute 29th May 2008 1 Collaborators include: Arnaud Doucet, Nick Whiteley
More informationarxiv: v1 [stat.co] 1 Jun 2015
arxiv:1506.00570v1 [stat.co] 1 Jun 2015 Towards automatic calibration of the number of state particles within the SMC 2 algorithm N. Chopin J. Ridgway M. Gerber O. Papaspiliopoulos CREST-ENSAE, Malakoff,
More informationParticle Metropolis-adjusted Langevin algorithms
Particle Metropolis-adjusted Langevin algorithms Christopher Nemeth, Chris Sherlock and Paul Fearnhead arxiv:1412.7299v3 [stat.me] 27 May 2016 Department of Mathematics and Statistics, Lancaster University,
More informationStatistical Machine Learning Lecture 8: Markov Chain Monte Carlo Sampling
1 / 27 Statistical Machine Learning Lecture 8: Markov Chain Monte Carlo Sampling Melih Kandemir Özyeğin University, İstanbul, Turkey 2 / 27 Monte Carlo Integration The big question : Evaluate E p(z) [f(z)]
More informationMonte Carlo Approximation of Monte Carlo Filters
Monte Carlo Approximation of Monte Carlo Filters Adam M. Johansen et al. Collaborators Include: Arnaud Doucet, Axel Finke, Anthony Lee, Nick Whiteley 7th January 2014 Context & Outline Filtering in State-Space
More informationSensor Fusion: Particle Filter
Sensor Fusion: Particle Filter By: Gordana Stojceska stojcesk@in.tum.de Outline Motivation Applications Fundamentals Tracking People Advantages and disadvantages Summary June 05 JASS '05, St.Petersburg,
More informationThe Hierarchical Particle Filter
and Arnaud Doucet http://go.warwick.ac.uk/amjohansen/talks MCMSki V Lenzerheide 7th January 2016 Context & Outline Filtering in State-Space Models: SIR Particle Filters [GSS93] Block-Sampling Particle
More informationMarkov Chain Monte Carlo Methods for Stochastic
Markov Chain Monte Carlo Methods for Stochastic Optimization i John R. Birge The University of Chicago Booth School of Business Joint work with Nicholas Polson, Chicago Booth. JRBirge U Florida, Nov 2013
More informationBlind Equalization via Particle Filtering
Blind Equalization via Particle Filtering Yuki Yoshida, Kazunori Hayashi, Hideaki Sakai Department of System Science, Graduate School of Informatics, Kyoto University Historical Remarks A sequential Monte
More informationData-driven Particle Filters for Particle Markov Chain Monte Carlo
Data-driven Particle Filters for Particle Markov Chain Monte Carlo Patrick Leung, Catherine S. Forbes, Gael M. Martin and Brendan McCabe September 13, 2016 Abstract This paper proposes new automated proposal
More informationBayesian statistical data assimilation for ecosystem models using Markov Chain Monte Carlo
Available online at www.sciencedirect.com Journal of Marine Systems 68 (2007) 439 456 www.elsevier.com/locate/jmarsys Bayesian statistical data assimilation for ecosystem models using Markov Chain Monte
More informationBayesian parameter inference for stochastic biochemical network models using particle MCMC
Bayesian parameter inference for stochastic biochemical network models using particle MCMC Andrew Golightly Darren J. Wilkinson May 16, 2011 Abstract Computational systems biology is concerned with the
More informationThe Kalman Filter ImPr Talk
The Kalman Filter ImPr Talk Ged Ridgway Centre for Medical Image Computing November, 2006 Outline What is the Kalman Filter? State Space Models Kalman Filter Overview Bayesian Updating of Estimates Kalman
More informationPackage RcppSMC. March 18, 2018
Type Package Title Rcpp Bindings for Sequential Monte Carlo Version 0.2.1 Date 2018-03-18 Package RcppSMC March 18, 2018 Author Dirk Eddelbuettel, Adam M. Johansen and Leah F. South Maintainer Dirk Eddelbuettel
More informationOn some properties of Markov chain Monte Carlo simulation methods based on the particle filter
On some properties of Markov chain Monte Carlo simulation methods based on the particle filter Michael K. Pitt Economics Department University of Warwick m.pitt@warwick.ac.uk Ralph S. Silva School of Economics
More informationLearning Static Parameters in Stochastic Processes
Learning Static Parameters in Stochastic Processes Bharath Ramsundar December 14, 2012 1 Introduction Consider a Markovian stochastic process X T evolving (perhaps nonlinearly) over time variable T. We
More informationLecture 2: From Linear Regression to Kalman Filter and Beyond
Lecture 2: From Linear Regression to Kalman Filter and Beyond Department of Biomedical Engineering and Computational Science Aalto University January 26, 2012 Contents 1 Batch and Recursive Estimation
More informationAn efficient stochastic approximation EM algorithm using conditional particle filters
An efficient stochastic approximation EM algorithm using conditional particle filters Fredrik Lindsten Linköping University Post Print N.B.: When citing this work, cite the original article. Original Publication:
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 informationA new iterated filtering algorithm
A new iterated filtering algorithm Edward Ionides University of Michigan, Ann Arbor ionides@umich.edu Statistics and Nonlinear Dynamics in Biology and Medicine Thursday July 31, 2014 Overview 1 Introduction
More informationIntroduction to Particle Filters for Data Assimilation
Introduction to Particle Filters for Data Assimilation Mike Dowd Dept of Mathematics & Statistics (and Dept of Oceanography Dalhousie University, Halifax, Canada STATMOS Summer School in Data Assimila5on,
More informationA Note on Auxiliary Particle Filters
A Note on Auxiliary Particle Filters Adam M. Johansen a,, Arnaud Doucet b a Department of Mathematics, University of Bristol, UK b Departments of Statistics & Computer Science, University of British Columbia,
More informationIntroduction to Machine Learning
Introduction to Machine Learning Brown University CSCI 1950-F, Spring 2012 Prof. Erik Sudderth Lecture 25: Markov Chain Monte Carlo (MCMC) Course Review and Advanced Topics Many figures courtesy Kevin
More informationAn Brief Overview of Particle Filtering
1 An Brief Overview of Particle Filtering Adam M. Johansen a.m.johansen@warwick.ac.uk www2.warwick.ac.uk/fac/sci/statistics/staff/academic/johansen/talks/ May 11th, 2010 Warwick University Centre for Systems
More informationL09. PARTICLE FILTERING. NA568 Mobile Robotics: Methods & Algorithms
L09. PARTICLE FILTERING NA568 Mobile Robotics: Methods & Algorithms Particle Filters Different approach to state estimation Instead of parametric description of state (and uncertainty), use a set of state
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 spatial hierarchical modeling for temperature extremes
Bayesian spatial hierarchical modeling for temperature extremes Indriati Bisono Dr. Andrew Robinson Dr. Aloke Phatak Mathematics and Statistics Department The University of Melbourne Maths, Informatics
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 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 informationA statistical overview and perspectives on data assimilation for marine biogeochemical models
Special Issue Paper Received: 6 October 2013, Revised: 22 January 2014, Accepted: 24 January 2014, Published online in Wiley Online Library (wileyonlinelibrary.com) DOI: 10.1002/env.2264 A statistical
More informationSequential Monte Carlo Methods for Bayesian Computation
Sequential Monte Carlo Methods for Bayesian Computation A. Doucet Kyoto Sept. 2012 A. Doucet (MLSS Sept. 2012) Sept. 2012 1 / 136 Motivating Example 1: Generic Bayesian Model Let X be a vector parameter
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 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 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 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 Methods for Stochastic Optimization
Markov Chain Monte Carlo Methods for Stochastic Optimization John R. Birge The University of Chicago Booth School of Business Joint work with Nicholas Polson, Chicago Booth. JRBirge U of Toronto, MIE,
More informationPseudo-marginal Metropolis-Hastings: a simple explanation and (partial) review of theory
Pseudo-arginal Metropolis-Hastings: a siple explanation and (partial) review of theory Chris Sherlock Motivation Iagine a stochastic process V which arises fro soe distribution with density p(v θ ). Iagine
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 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 informationParticle Filtering Approaches for Dynamic Stochastic Optimization
Particle Filtering Approaches for Dynamic Stochastic Optimization John R. Birge The University of Chicago Booth School of Business Joint work with Nicholas Polson, Chicago Booth. JRBirge I-Sim Workshop,
More informationIntroduction to Probability and Statistics (Continued)
Introduction to Probability and Statistics (Continued) Prof. icholas Zabaras Center for Informatics and Computational Science https://cics.nd.edu/ University of otre Dame otre Dame, Indiana, USA Email:
More informationFlexible Particle Markov chain Monte Carlo methods with an application to a factor stochastic volatility model
Flexible Particle Markov chain Monte Carlo methods with an application to a factor stochastic volatility model arxiv:1401.1667v4 [stat.co] 23 Jan 2018 Eduardo F. Mendes School of Applied Mathematics Fundação
More informationLearning of state-space models with highly informative observations: a tempered Sequential Monte Carlo solution
Learning of state-space models with highly informative observations: a tempered Sequential Monte Carlo solution Andreas Svensson, Thomas B. Schön, and Fredrik Lindsten Department of Information Technology,
More informationAn intro to particle methods for parameter inference in state-space models
An intro to particle methods for parameter inference in state-space models PhD course FMS020F NAMS002 Statistical inference for partially observed stochastic processes, Lund University http://goo.gl/sx8vu9
More informationAn introduction to particle filters
An introduction to particle filters Andreas Svensson Department of Information Technology Uppsala University June 10, 2014 June 10, 2014, 1 / 16 Andreas Svensson - An introduction to particle filters Outline
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 informationComputer Vision Group Prof. Daniel Cremers. 11. Sampling Methods
Prof. Daniel Cremers 11. Sampling Methods Sampling Methods Sampling Methods are widely used in Computer Science as an approximation of a deterministic algorithm to represent uncertainty without a parametric
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 informationDownloaded from:
Camacho, A; Kucharski, AJ; Funk, S; Breman, J; Piot, P; Edmunds, WJ (2014) Potential for large outbreaks of Ebola virus disease. Epidemics, 9. pp. 70-8. ISSN 1755-4365 DOI: https://doi.org/10.1016/j.epidem.2014.09.003
More informationBayesian room-acoustic modal analysis
Bayesian room-acoustic modal analysis Wesley Henderson a) Jonathan Botts b) Ning Xiang c) Graduate Program in Architectural Acoustics, School of Architecture, Rensselaer Polytechnic Institute, Troy, New
More informationBayesian Modeling of Conditional Distributions
Bayesian Modeling of Conditional Distributions John Geweke University of Iowa Indiana University Department of Economics February 27, 2007 Outline Motivation Model description Methods of inference Earnings
More informationLecture 6: Multiple Model Filtering, Particle Filtering and Other Approximations
Lecture 6: Multiple Model Filtering, Particle Filtering and Other Approximations Department of Biomedical Engineering and Computational Science Aalto University April 28, 2010 Contents 1 Multiple Model
More informationData-driven Particle Filters for Particle Markov Chain Monte Carlo
Data-driven Particle Filters for Particle Markov Chain Monte Carlo Patrick Leung, Catherine S. Forbes, Gael M. Martin and Brendan McCabe December 21, 2016 Abstract This paper proposes new automated proposal
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 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 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 informationIntroduction to Mobile Robotics Bayes Filter Particle Filter and Monte Carlo Localization
Introduction to Mobile Robotics Bayes Filter Particle Filter and Monte Carlo Localization Wolfram Burgard, Cyrill Stachniss, Maren Bennewitz, Kai Arras 1 Motivation Recall: Discrete filter Discretize the
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 informationCheng Soon Ong & Christian Walder. Canberra February June 2018
Cheng Soon Ong & Christian Walder Research Group and College of Engineering and Computer Science Canberra February June 2018 Outlines Overview Introduction Linear Algebra Probability Linear Regression
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 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 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 information28 : Approximate Inference - Distributed MCMC
10-708: Probabilistic Graphical Models, Spring 2015 28 : Approximate Inference - Distributed MCMC Lecturer: Avinava Dubey Scribes: Hakim Sidahmed, Aman Gupta 1 Introduction For many interesting problems,
More informationParticle MCMC for Bayesian Microwave Control
Particle MCMC for Bayesian Microwave Control P. Minvielle 1, A. Todeschini 2, F. Caron 3, P. Del Moral 4, 1 CEA-CESTA, 33114 Le Barp, France 2 INRIA Bordeaux Sud-Ouest, 351, cours de la Liberation, 33405
More information1 / 31 Identification of nonlinear dynamic systems (guest lectures) Brussels, Belgium, June 8, 2015.
Outline Part 4 Nonlinear system identification using sequential Monte Carlo methods Part 4 Identification strategy 2 (Data augmentation) Aim: Show how SMC can be used to implement identification strategy
More informationAN EFFICIENT TWO-STAGE SAMPLING METHOD IN PARTICLE FILTER. Qi Cheng and Pascal Bondon. CNRS UMR 8506, Université Paris XI, France.
AN EFFICIENT TWO-STAGE SAMPLING METHOD IN PARTICLE FILTER Qi Cheng and Pascal Bondon CNRS UMR 8506, Université Paris XI, France. August 27, 2011 Abstract We present a modified bootstrap filter to draw
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 informationarxiv: v1 [stat.co] 2 Nov 2017
Binary Bouncy Particle Sampler arxiv:1711.922v1 [stat.co] 2 Nov 217 Ari Pakman Department of Statistics Center for Theoretical Neuroscience Grossman Center for the Statistics of Mind Columbia University
More informationBayesian Computations for DSGE Models
Bayesian Computations for DSGE Models Frank Schorfheide University of Pennsylvania, PIER, CEPR, and NBER October 23, 2017 This Lecture is Based on Bayesian Estimation of DSGE Models Edward P. Herbst &
More informationComputer Vision Group Prof. Daniel Cremers. 14. Sampling Methods
Prof. Daniel Cremers 14. Sampling Methods Sampling Methods Sampling Methods are widely used in Computer Science as an approximation of a deterministic algorithm to represent uncertainty without a parametric
More informationMicrobial Grazers Lab
Microbial Grazers Lab Objective: Measure the rate at which bacteria are consumed by predators. Overview Size based food webs Microbial loop concepts acterial predators Methods to assess microbial grazing
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 informationBiogeochemistry of trace elements and isotopes in the Indian Ocean
Biogeochemistry of trace elements and isotopes in the Indian Ocean Sunil Kumar Singh Geosciences Division Physical Research Laboratory Ahmedabad 380009 Ministry of Earth Sciences Government of India 2
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 informationWater-vegetation interactions in a sea of water: Plankton in mesoscale turbulence
Water-vegetation interactions in a sea of water: Plankton in mesoscale turbulence A. Provenzale ISAC-CNR, Torino and CIMA, Savona, Italy Work done with: Annalisa Bracco, Claudia Pasquero Adrian Martin,
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 information16 : Markov Chain Monte Carlo (MCMC)
10-708: Probabilistic Graphical Models 10-708, Spring 2014 16 : Markov Chain Monte Carlo MCMC Lecturer: Matthew Gormley Scribes: Yining Wang, Renato Negrinho 1 Sampling from low-dimensional distributions
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 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 informationEstimation and Inference on Dynamic Panel Data Models with Stochastic Volatility
Estimation and Inference on Dynamic Panel Data Models with Stochastic Volatility Wen Xu Department of Economics & Oxford-Man Institute University of Oxford (Preliminary, Comments Welcome) Theme y it =
More informationSequential Monte Carlo Samplers for Applications in High Dimensions
Sequential Monte Carlo Samplers for Applications in High Dimensions Alexandros Beskos National University of Singapore KAUST, 26th February 2014 Joint work with: Dan Crisan, Ajay Jasra, Nik Kantas, Alex
More informationGaussian Mixture Model
Case Study : Document Retrieval MAP EM, Latent Dirichlet Allocation, Gibbs Sampling Machine Learning/Statistics for Big Data CSE599C/STAT59, University of Washington Emily Fox 0 Emily Fox February 5 th,
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 information