Calibration and selection of stochastic palaeoclimate models
|
|
- Erica Thornton
- 5 years ago
- Views:
Transcription
1 Calibration and selection of stochastic palaeoclimate models Michel Crucifix, Jake Carson, Simon Preston, Richard Wilkinson, Jonathan Rougier IRM - Internal Seminar - Uccle 2 April 2014
2 Posing the problem Principles and methods of model calibration and selection Applications
3 2 2.5 Benthic d 18 O (normalised) time (ka) time (ka) L. E. Lisiecki and M. E. Raymo, Paleoceanography, 20, PA
4 External factors? t t + 1 Internal structure / parameters
5
6 Astronomical forcing as quasi-periodic signal Precession : , , , , ,... years Obliquity : 41090, , , , 41811,... years Eccentricity : , , , ,... years from A. Berger and M.F. Loutre, Quaternary Science Reviews, 10,
7 The Langevin model (1 dimension) dx = ( U 0 (x)+f (t))dt + dw derivative of a cimate potential stochastic forcing deterministic forcing Nicolis and Nicolis, Tellus (1981) ; Nicolis, Tellus, (1982)
8 Competing paradigm : Oscillators (or oscillator-like) dx = ( + F (t) y)dt dy = (y y 3 + x) dt
9 X Y V D Forcing LR time (ka) M. Crucifix, How can a glacial inception be predicted?, The Holocene, 21,
10 Forced stochastic van der pol Model dx = ( + F (t) y)dt dy = (y y 3 + x) dt + dw
11 V D model with stochastic perturbation X V Y D time (ka) Figure 3. Same model as Figure 2 but with a small stochastic M. Crucifix, How can a glacial inception be predicted?, The Holocene, 21,
12 Ice age models founded on the hypothesis of a synchronised limit cycle feature a strange non-chaotic attractor Lyapunov stable BUT structurally unstable M. Crucifix, Why could ice ages be unpredictable?, Climate of the Past, 9, Takahito Mitsui and Kazuyuki Aihara, Dynamics between order and chaos in conceptual models of glacial cycles, Climate Dynamics, see also :
13 Posing the problem Identification of structure and parameters (using information emerging from observation and simulations with GCMs, ice sheet models, etc. ) in non-autonomous systems that may be structurally unstable... with sparse-data, and uncertain dating
14 Posing the problem Principles and methods of model calibration and selection Applications
15 / Z One observes Z = z
16 Prior : P ( ) Likelihood : P (z ) Model Evidence : P (Z = z) = R 2 P (z )P ( ) Posterior : P ( z) = P (z )P ( ) P (Z=z)
17 The latent variable as an intermediate step / X (hidden or latent variable ) Z P (z ) = Z X P (z, (X))P (X)dX z is a shorthand for Z = z
18 Dynamical system : many latent variables (F t 1, ) (F t, ) (F t+1, ) X 0 /... / X t 1 / X t / X t+1 /... Z t 1 Z t Z t+1
19 Intractable and unsolved : some thoughts on statistical data assimilation with uncertain static parameters Jonathan Rougier Three years ago, I started a collaboration with Michel Crucifix at the Université catholique de Louvain, on glacial cycles. Our intention was to use simple phenomenological models [1], which would be calibrated against measurements from ice core and benthic datasets, and used to address questions such as Ruddiman s Early Anthropogenic hypothesis (reviewed in [2]). I thought we would use an approach suggested by Liu & West [3], which was a sophisticated treatment of data assimilation with uncertain static parameters (defined below), in which the state vector is augmented with the static parameters. One year later, I was discussing our lack of progress with Prof. Christophe Andrieu, an authority on statistical data assimilation. He was not at all surprised, because in his view the problem was, in general, intractable and unsolved. Regardless, I simply upped the computational sophistication of our methods, including switching to his own recent development, particle Markov chain Monte Carlo (P-MCMC, [4]) and running on a faster computer /rsta Phil. Trans. R. Soc. A 28 May 2013 vol. 371 no
20 stochastic dynamical system implemented numerically as a map (numerical scheme) random vector (X t 1,! t 1 ) deterministic! X t inference on the! i equivalent to inference on the X i
21 State inference : P (X z, ) Joint state-parameter inference : P (,X z) Parameter inference : R X P ( z,x)
22 State inference : P (X z, ) Sequential algorithms Kalman Filter (exact if linear and Gaussian) Extended, Unscented, Ensemble Kalman Filters Sequential Monte Carlo (Particle Filter) Variational (bit by bit) Global approaches Variational (global) Annealing
23 Kalman filters Sequential estimation of state (assimilation of data) The information being propagated is The mean and the variance Their dynamics are estimated using the Linear Tangeant (Extended Kalman Filter) an ensemble of points (Ensemble of Points) the dynamics of a well chosen set of points : the Sigma Points (Unscented Kalman Filter, Julier and Uhlman) EKF and UKF are deterministic and biased
24 Particle Filter X w i,j = p(! i,j )p(z i X i,j )/q(! i,j )... Lik j = P w i,j Lik( ) = P j (Lik j)! i,j q(! i,j ) t j t j+1
25 Joint state-parameter inference : P (,X z) 1. Augmented state approach ( Liu and West, 2001) add =0 2. PMCMC (Andrieux et al. 2010) 3. SMC 2 (Chopin et al. 2012)
26 SMC-2 (Chopin et al. 2012) z (1) (2) x2 x2 x3 x1 x2 x1 x2 x1 t 0 z t 1 t 1 t 0 t 2 (1) Sample (,X 0 ) and propagate forward until filter degenerates (e.g. : at time t 1 ) (2) Resample, perturb leaving the posterior distribution at t 1 invariant, and continue formard
27 SMC-2 (Chopin et al. 2012) z (1) (2) x2 x2 x3 x1 x2 x1 x2 x1 t 0 z t 1 t 1 t 0 t 2 (1) Sample (,X 0 ) and propagate forward until filter degenerates (e.g. : at time t 1 ) (2) Resample, perturb leaving the posterior distribution at t 1 invariant, and continue formard
28 Posing the problem Principles and methods of model calibration and selection Applications
29 application 1 : Liu and West Algorithm M. Crucifix and J. Rougier, European Physics Journal Special Topics, 174,
30 The Satlzman & Maash model Ice Sheets (ζ) [1] B. Saltzman and K. A. Maasch. A first-order global model of late Cenozoic climate. Trans. R. Soc. Edinburgh Earth Sci, 81: , Ocean Circulation (ψ) Carbon Dioxide (theta) = 1 (c 2 )µ 3 k 2 k R 2 R(t) + I = K ( 5 6 ) + F µ (t) + = 1 2 I 3 + F (t) + K = 1 ( )
31 Augmented-state strategy : handle filter degeneracy 2e+13 9e+13 Parameter: a1 Parameter: b1 4e Parameter: b2 Parameter: b3 Parameter: b4 0 3e 04 2e 05 Parameter: b Years (ky, CE) M. Crucifix and J. Rougier, European Physics Journal Special Topics, 174,
32 application 2 : Estimating likelihood: Kalman vs Particle unpublished note
33 linear model linear model dx = ( + F (t) y)dt + 1 dw 1 dy = (y + x)dt + 2 dw 2 with observations z t = x t + t, t N (0, obs)
34 Log lik estimator : linear model as a function of TAU log likelihood KF PF (n=200) PF (n=1500) tau
35 Lik. estimator : linear model i log likelihood KF PF (n=100) log of mean(lik) PF
36 Forced van der pol Model non- linear model dx = ( + F (t) y)dt + 1 dw 1 dy = (y y 3 + x)dt + 2 dw 2 with observations z t = x t + t, t N (0, obs)
37 Log lik estimator : NON linear model as a function of TAU log likelihood UKF UKF PF (n=200) PF (n=1500) tau
38 application 3 : Model Calibration : PMCMC vs SMC2 poster, EGU Vienna, 2013
39 Forced slow-fast model (CR12) Dynamics dx = [ x + 2 (x 2 1) y + + ]dt + 1 dw 1 dy = [ (y y 3 /3) + x) ]dt + 2 dw 2 state Observations determinstic drift = 30 (slow-fast system) astronomical forcing additive di usion ((w 1,w 2 ) =Wiener) z t = x t + t, t N (0, obs)
40 LR04 stack Time [ka] Huybers 04 stack Time [ka] benthic d18o [std] ODP929 Time [ka] benthic d18o [per mil]
41 Calibration on LR04 β 0 β 1 β 2 δ γ π Density Density Density Density Density γ ε σ x σ y σ z Density Density Density Density PMCMC SMC2 prior
42 application 4 : Model selection with SMC2 Jake Carson s thesis : in prep.
43 Bayes Factor = Model Evidence (M1) Model Evidence (M2)
44 Observations = data simulated with the Forced SM91 SS-F Model Evidence (Forced Model) Evidence (Unforced Model) SM T PP
45 Observations = LR04 stack (astronomically tuned) LR04 Model Evidence (Forced Model) Evidence (Unforced Model) SM T PP
46 Observations = non-astronomically tuned data ODP982b Model Evidence (Forced Model) Evidence (Unforced Model) SM T PP12 * -
47 Particle Filter... is an art choice of number of particles, proposal, comparison between algorithm, ofter problem dependent Augmented space approach for parameter estimation may be deceptive. We no longer recommend it. PMCMC was proposed as a solution may have poor mixing : stationary solutions not trusted SMC^2 conceptually more advanced, and gives model evidence as a bonus : model selection is possible [ The practical added value of PF in additive, Gaussian noise problem still needs to be demonstrated ]
Climate sensitivity and dynamical systems
Michel Crucifix, 01 February 2016, p. 1. Climate sensitivity and dynamical systems Michel Crucifix Université catholique de Louvain & Fonds National de la Recherche Scientifique michel.crucifix@uclouvain.be
More informationTime-series analysis for palaeoclimate data
Time-series analysis for palaeoclimate data Michel Crucifix acknowledgements: Guillaume Lenoir, Carlos Almeida (ULouvain), Jonathan Rougier (Bristol), Richard Wilkinson, Jake Carson (Nottingham) Disclaimer
More informationIce age challenges. Michel Crucifix. Université catholique de Louvain & Fonds National de la Recherche Scientifique
Michel Crucifix, Dynamics Days Exeter 2015, p. 1. Ice age challenges Michel Crucifix Université catholique de Louvain & Fonds National de la Recherche Scientifique michel.crucifix@uclouvain.be Dynamics
More informationBayesian model selection for the glacial interglacial cycle
Appl. Statist. (2018) 67, Part 1, pp. 25 54 Bayesian model selection for the glacial interglacial cycle Jake Carson, Imperial College London, UK Michel Crucifix, Earth and Life Institute, Université catholique
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 informationBayesian model selection for the glacial-interglacial cycle
arxiv:1511.03467v1 [stat.ap] 11 Nov 2015 Bayesian model selection for the glacial-interglacial cycle Jake Carson 1, Michel Crucifix 2, Simon Preston 3, Richard Wilkinson 4 1 University of Warwick, Department
More informationPATTERN RECOGNITION AND MACHINE LEARNING CHAPTER 13: SEQUENTIAL DATA
PATTERN RECOGNITION AND MACHINE LEARNING CHAPTER 13: SEQUENTIAL DATA Contents in latter part Linear Dynamical Systems What is different from HMM? Kalman filter Its strength and limitation Particle Filter
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 informationLecture 2: From Linear Regression to Kalman Filter and Beyond
Lecture 2: From Linear Regression to Kalman Filter and Beyond January 18, 2017 Contents 1 Batch and Recursive Estimation 2 Towards Bayesian Filtering 3 Kalman Filter and Bayesian Filtering and Smoothing
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 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 informationLagrangian Data Assimilation and Its Application to Geophysical Fluid Flows
Lagrangian Data Assimilation and Its Application to Geophysical Fluid Flows Laura Slivinski June, 3 Laura Slivinski (Brown University) Lagrangian Data Assimilation June, 3 / 3 Data Assimilation Setup:
More informationA Note on the Particle Filter with Posterior Gaussian Resampling
Tellus (6), 8A, 46 46 Copyright C Blackwell Munksgaard, 6 Printed in Singapore. All rights reserved TELLUS A Note on the Particle Filter with Posterior Gaussian Resampling By X. XIONG 1,I.M.NAVON 1,2 and
More informationKalman filtering and friends: Inference in time series models. Herke van Hoof slides mostly by Michael Rubinstein
Kalman filtering and friends: Inference in time series models Herke van Hoof slides mostly by Michael Rubinstein Problem overview Goal Estimate most probable state at time k using measurement up to time
More information(Palaeo-)climate sensitivity: ideas and definitions from the NPG literature
(Palaeo-)climate sensitivity: ideas and definitions from the NPG literature Michel Crucifix Université catholique de Louvain & Belgian National Fund of Scientific Research Ringberg Grand Challenge Workshop:
More informationStochastic Collocation Methods for Polynomial Chaos: Analysis and Applications
Stochastic Collocation Methods for Polynomial Chaos: Analysis and Applications Dongbin Xiu Department of Mathematics, Purdue University Support: AFOSR FA955-8-1-353 (Computational Math) SF CAREER DMS-64535
More informationOn a Data Assimilation Method coupling Kalman Filtering, MCRE Concept and PGD Model Reduction for Real-Time Updating of Structural Mechanics Model
On a Data Assimilation Method coupling, MCRE Concept and PGD Model Reduction for Real-Time Updating of Structural Mechanics Model 2016 SIAM Conference on Uncertainty Quantification Basile Marchand 1, Ludovic
More informationShort tutorial on data assimilation
Mitglied der Helmholtz-Gemeinschaft Short tutorial on data assimilation 23 June 2015 Wolfgang Kurtz & Harrie-Jan Hendricks Franssen Institute of Bio- and Geosciences IBG-3 (Agrosphere), Forschungszentrum
More informationWhy do we care? Examples. Bayes Rule. What room am I in? Handling uncertainty over time: predicting, estimating, recognizing, learning
Handling uncertainty over time: predicting, estimating, recognizing, learning Chris Atkeson 004 Why do we care? Speech recognition makes use of dependence of words and phonemes across time. Knowing where
More informationRecent Developments in the Theory of Glacial Cycles
Recent Developments in the Theory of Richard McGehee Seminar on the Mathematics of Climate Change School of Mathematics October 6, 010 Hansen, et al, Target atmospheric CO: Where should humanity aim? Open
More informationLecture 6: Bayesian Inference in SDE Models
Lecture 6: Bayesian Inference in SDE Models Bayesian Filtering and Smoothing Point of View Simo Särkkä Aalto University Simo Särkkä (Aalto) Lecture 6: Bayesian Inference in SDEs 1 / 45 Contents 1 SDEs
More informationCorrelation integral likelihood for stochastic differential equations
Correlation integral likelihood for stochastic differential equations Heikki Haario, Janne Hakkarainen, Ramona Maraia, Sebastian Springer Abstract A new approach was recently introduced for the task of
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 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 informationLinear Dynamical Systems
Linear Dynamical Systems Sargur N. srihari@cedar.buffalo.edu Machine Learning Course: http://www.cedar.buffalo.edu/~srihari/cse574/index.html Two Models Described by Same Graph Latent variables Observations
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 informationThe Particle Filter. PD Dr. Rudolph Triebel Computer Vision Group. Machine Learning for Computer Vision
The Particle Filter Non-parametric implementation of Bayes filter Represents the belief (posterior) random state samples. by a set of This representation is approximate. Can represent distributions that
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 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 informationDynamic System Identification using HDMR-Bayesian Technique
Dynamic System Identification using HDMR-Bayesian Technique *Shereena O A 1) and Dr. B N Rao 2) 1), 2) Department of Civil Engineering, IIT Madras, Chennai 600036, Tamil Nadu, India 1) ce14d020@smail.iitm.ac.in
More informationA new Hierarchical Bayes approach to ensemble-variational data assimilation
A new Hierarchical Bayes approach to ensemble-variational data assimilation Michael Tsyrulnikov and Alexander Rakitko HydroMetCenter of Russia College Park, 20 Oct 2014 Michael Tsyrulnikov and Alexander
More informationBayesian Inference and the Symbolic Dynamics of Deterministic Chaos. Christopher C. Strelioff 1,2 Dr. James P. Crutchfield 2
How Random Bayesian Inference and the Symbolic Dynamics of Deterministic Chaos Christopher C. Strelioff 1,2 Dr. James P. Crutchfield 2 1 Center for Complex Systems Research and Department of Physics University
More informationLagrangian Data Assimilation and Manifold Detection for a Point-Vortex Model. David Darmon, AMSC Kayo Ide, AOSC, IPST, CSCAMM, ESSIC
Lagrangian Data Assimilation and Manifold Detection for a Point-Vortex Model David Darmon, AMSC Kayo Ide, AOSC, IPST, CSCAMM, ESSIC Background Data Assimilation Iterative process Forecast Analysis Background
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 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 informationState Estimation using Moving Horizon Estimation and Particle Filtering
State Estimation using Moving Horizon Estimation and Particle Filtering James B. Rawlings Department of Chemical and Biological Engineering UW Math Probability Seminar Spring 2009 Rawlings MHE & PF 1 /
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 informationApproximate Bayesian Computation and Particle Filters
Approximate Bayesian Computation and Particle Filters Dennis Prangle Reading University 5th February 2014 Introduction Talk is mostly a literature review A few comments on my own ongoing research See Jasra
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 information9 Multi-Model State Estimation
Technion Israel Institute of Technology, Department of Electrical Engineering Estimation and Identification in Dynamical Systems (048825) Lecture Notes, Fall 2009, Prof. N. Shimkin 9 Multi-Model State
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 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 informationGaussian Process Approximations of Stochastic Differential Equations
Gaussian Process Approximations of Stochastic Differential Equations Cédric Archambeau Dan Cawford Manfred Opper John Shawe-Taylor May, 2006 1 Introduction Some of the most complex models routinely run
More informationState and Parameter Estimation in Stochastic Dynamical Models
State and Parameter Estimation in Stochastic Dynamical Models Timothy DelSole George Mason University, Fairfax, Va and Center for Ocean-Land-Atmosphere Studies, Calverton, MD June 21, 2011 1 1 collaboration
More informationAutonomous Mobile Robot Design
Autonomous Mobile Robot Design Topic: Extended Kalman Filter Dr. Kostas Alexis (CSE) These slides relied on the lectures from C. Stachniss, J. Sturm and the book Probabilistic Robotics from Thurn et al.
More informationObliquity pacing of the late Pleistocene glacial terminations
Obliquity pacing of the late Pleistocene glacial terminations Peter Huybers 1 & Carl Wunsch 2 1 Woods Hole Oceanographic Institution, Woods Hole, USA 2 Massachusetts Institute of Technology, Cambridge,
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 informationMonte Carlo Methods. Geoff Gordon February 9, 2006
Monte Carlo Methods Geoff Gordon ggordon@cs.cmu.edu February 9, 2006 Numerical integration problem 5 4 3 f(x,y) 2 1 1 0 0.5 0 X 0.5 1 1 0.8 0.6 0.4 Y 0.2 0 0.2 0.4 0.6 0.8 1 x X f(x)dx Used for: function
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 informationAutonomous Navigation for Flying Robots
Computer Vision Group Prof. Daniel Cremers Autonomous Navigation for Flying Robots Lecture 6.2: Kalman Filter Jürgen Sturm Technische Universität München Motivation Bayes filter is a useful tool for state
More informationThe Ensemble Kalman Filter:
p.1 The Ensemble Kalman Filter: Theoretical formulation and practical implementation Geir Evensen Norsk Hydro Research Centre, Bergen, Norway Based on Evensen 23, Ocean Dynamics, Vol 53, No 4 p.2 The Ensemble
More informationA nested sampling particle filter for nonlinear data assimilation
Quarterly Journal of the Royal Meteorological Society Q. J. R. Meteorol. Soc. : 14, July 2 A DOI:.2/qj.224 A nested sampling particle filter for nonlinear data assimilation Ahmed H. Elsheikh a,b *, Ibrahim
More informationSeminar: Data Assimilation
Seminar: Data Assimilation Jonas Latz, Elisabeth Ullmann Chair of Numerical Mathematics (M2) Technical University of Munich Jonas Latz, Elisabeth Ullmann (TUM) Data Assimilation 1 / 28 Prerequisites Bachelor:
More informationObliquity pacing of the late Pleistocene glacial terminations
Obliquity pacing of the late Pleistocene glacial terminations Peter Huybers 1 & Carl Wunsch 2 1 Woods Hole Oceanographic Institution, Woods Hole, USA 2 Massachusetts Institute of Technology, Cambridge,
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 informationWhy do we care? Measurements. Handling uncertainty over time: predicting, estimating, recognizing, learning. Dealing with time
Handling uncertainty over time: predicting, estimating, recognizing, learning Chris Atkeson 2004 Why do we care? Speech recognition makes use of dependence of words and phonemes across time. Knowing where
More informationLearning Energy-Based Models of High-Dimensional Data
Learning Energy-Based Models of High-Dimensional Data Geoffrey Hinton Max Welling Yee-Whye Teh Simon Osindero www.cs.toronto.edu/~hinton/energybasedmodelsweb.htm Discovering causal structure as a goal
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 informationA State Space Model for Wind Forecast Correction
A State Space Model for Wind Forecast Correction Valrie Monbe, Pierre Ailliot 2, and Anne Cuzol 1 1 Lab-STICC, Université Européenne de Bretagne, France (e-mail: valerie.monbet@univ-ubs.fr, anne.cuzol@univ-ubs.fr)
More informationNonlinear and/or Non-normal Filtering. Jesús Fernández-Villaverde University of Pennsylvania
Nonlinear and/or Non-normal Filtering Jesús Fernández-Villaverde University of Pennsylvania 1 Motivation Nonlinear and/or non-gaussian filtering, smoothing, and forecasting (NLGF) problems are pervasive
More informationProbabilistic Graphical Models
Probabilistic Graphical Models Brown University CSCI 2950-P, Spring 2013 Prof. Erik Sudderth Lecture 13: Learning in Gaussian Graphical Models, Non-Gaussian Inference, Monte Carlo Methods Some figures
More informationSandwiching the marginal likelihood using bidirectional Monte Carlo. Roger Grosse
Sandwiching the marginal likelihood using bidirectional Monte Carlo Roger Grosse Ryan Adams Zoubin Ghahramani Introduction When comparing different statistical models, we d like a quantitative criterion
More informationConstrained State Estimation Using the Unscented Kalman Filter
16th Mediterranean Conference on Control and Automation Congress Centre, Ajaccio, France June 25-27, 28 Constrained State Estimation Using the Unscented Kalman Filter Rambabu Kandepu, Lars Imsland and
More informationMethods of Data Assimilation and Comparisons for Lagrangian Data
Methods of Data Assimilation and Comparisons for Lagrangian Data Chris Jones, Warwick and UNC-CH Kayo Ide, UCLA Andrew Stuart, Jochen Voss, Warwick Guillaume Vernieres, UNC-CH Amarjit Budiraja, UNC-CH
More informationEnsemble Kalman Filter
Ensemble Kalman Filter Geir Evensen and Laurent Bertino Hydro Research Centre, Bergen, Norway, Nansen Environmental and Remote Sensing Center, Bergen, Norway The Ensemble Kalman Filter (EnKF) Represents
More informationSequential Bayesian Inference for Dynamic State Space. Model Parameters
Sequential Bayesian Inference for Dynamic State Space Model Parameters Arnab Bhattacharya and Simon Wilson 1. INTRODUCTION Dynamic state-space models [24], consisting of a latent Markov process X 0,X 1,...
More informationPOMP inference via iterated filtering
POMP inference via iterated filtering Edward Ionides University of Michigan, Department of Statistics Lecture 3 at Wharton Statistics Department Thursday 27th April, 2017 Slides are online at http://dept.stat.lsa.umich.edu/~ionides/talks/upenn
More informationRobotics. Mobile Robotics. Marc Toussaint U Stuttgart
Robotics Mobile Robotics State estimation, Bayes filter, odometry, particle filter, Kalman filter, SLAM, joint Bayes filter, EKF SLAM, particle SLAM, graph-based SLAM Marc Toussaint U Stuttgart DARPA Grand
More informationInference for partially observed stochastic dynamic systems: A new algorithm, its theory and applications
Inference for partially observed stochastic dynamic systems: A new algorithm, its theory and applications Edward Ionides Department of Statistics, University of Michigan ionides@umich.edu Statistics Department
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 informationSmoothers: Types and Benchmarks
Smoothers: Types and Benchmarks Patrick N. Raanes Oxford University, NERSC 8th International EnKF Workshop May 27, 2013 Chris Farmer, Irene Moroz Laurent Bertino NERSC Geir Evensen Abstract Talk builds
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 informationGaussian Process Approximations of Stochastic Differential Equations
Gaussian Process Approximations of Stochastic Differential Equations Cédric Archambeau Centre for Computational Statistics and Machine Learning University College London c.archambeau@cs.ucl.ac.uk CSML
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 informationEnKF-based particle filters
EnKF-based particle filters Jana de Wiljes, Sebastian Reich, Wilhelm Stannat, Walter Acevedo June 20, 2017 Filtering Problem Signal dx t = f (X t )dt + 2CdW t Observations dy t = h(x t )dt + R 1/2 dv t.
More information1 / 32 Summer school on Foundations and advances in stochastic filtering (FASF) Barcelona, Spain, June 22, 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 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 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 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 informationSequential Monte Carlo Methods in High Dimensions
Sequential Monte Carlo Methods in High Dimensions Alexandros Beskos Statistical Science, UCL Oxford, 24th September 2012 Joint work with: Dan Crisan, Ajay Jasra, Nik Kantas, Andrew Stuart Imperial College,
More informationSTA 4273H: Statistical Machine Learning
STA 4273H: Statistical Machine Learning Russ Salakhutdinov Department of Statistics! rsalakhu@utstat.toronto.edu! http://www.utstat.utoronto.ca/~rsalakhu/ Sidney Smith Hall, Room 6002 Lecture 11 Project
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 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 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 informationDual Estimation and the Unscented Transformation
Dual Estimation and the Unscented Transformation Eric A. Wan ericwan@ece.ogi.edu Rudolph van der Merwe rudmerwe@ece.ogi.edu Alex T. Nelson atnelson@ece.ogi.edu Oregon Graduate Institute of Science & Technology
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 informationPalaeoclimate dynamics : a voyage through scales
Palaeoclimate dynamics : a voyage through scales Michel Crucifix, Takahito Mitsui, Guillaume Lenoir EGU General Assembly 2015 16 April 2015 EGU2015-3632 1976... other is that the small control eccentricthe
More informationAdvanced Computational Methods in Statistics: Lecture 5 Sequential Monte Carlo/Particle Filtering
Advanced Computational Methods in Statistics: Lecture 5 Sequential Monte Carlo/Particle Filtering Axel Gandy Department of Mathematics Imperial College London http://www2.imperial.ac.uk/~agandy London
More information2D Image Processing (Extended) Kalman and particle filter
2D Image Processing (Extended) Kalman and particle filter Prof. Didier Stricker Dr. Gabriele Bleser Kaiserlautern University http://ags.cs.uni-kl.de/ DFKI Deutsches Forschungszentrum für Künstliche Intelligenz
More informationROBOTICS 01PEEQW. Basilio Bona DAUIN Politecnico di Torino
ROBOTICS 01PEEQW Basilio Bona DAUIN Politecnico di Torino Probabilistic Fundamentals in Robotics Gaussian Filters Course Outline Basic mathematical framework Probabilistic models of mobile robots Mobile
More informationOn the (Nonlinear) Causes of Abrupt Climate Change During the Last Ice Age
On the (Nonlinear) Causes of Abrupt Climate Change During the Last Ice Age J.A. Rial Wave Propagation Lab, University of North Carolina-Chapel Chapel Hill The Astronomical Theory of the Ice Ages Precession
More informationCS491/691: Introduction to Aerial Robotics
CS491/691: Introduction to Aerial Robotics Topic: State Estimation Dr. Kostas Alexis (CSE) World state (or system state) Belief state: Our belief/estimate of the world state World state: Real state of
More informationFundamental Probability and Statistics
Fundamental Probability and Statistics "There are known knowns. These are things we know that we know. There are known unknowns. That is to say, there are things that we know we don't know. But there are
More informationGaussian processes for inference in stochastic differential equations
Gaussian processes for inference in stochastic differential equations Manfred Opper, AI group, TU Berlin November 6, 2017 Manfred Opper, AI group, TU Berlin (TU Berlin) inference in SDE November 6, 2017
More informationSequential Monte Carlo and Particle Filtering. Frank Wood Gatsby, November 2007
Sequential Monte Carlo and Particle Filtering Frank Wood Gatsby, November 2007 Importance Sampling Recall: Let s say that we want to compute some expectation (integral) E p [f] = p(x)f(x)dx and we remember
More informationConvergence of the Ensemble Kalman Filter in Hilbert Space
Convergence of the Ensemble Kalman Filter in Hilbert Space Jan Mandel Center for Computational Mathematics Department of Mathematical and Statistical Sciences University of Colorado Denver Parts based
More informationParticle Filters. Outline
Particle Filters M. Sami Fadali Professor of EE University of Nevada Outline Monte Carlo integration. Particle filter. Importance sampling. Degeneracy Resampling Example. 1 2 Monte Carlo Integration Numerical
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 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 informationState-Space Methods for Inferring Spike Trains from Calcium Imaging
State-Space Methods for Inferring Spike Trains from Calcium Imaging Joshua Vogelstein Johns Hopkins April 23, 2009 Joshua Vogelstein (Johns Hopkins) State-Space Calcium Imaging April 23, 2009 1 / 78 Outline
More information