example1_lgss: State estimation in a linear Gaussian state space model

Description Usage Details Value Note Author(s) References Examples

View source: R/example1-lgss.R


Minimal working example of state estimation in a linear Gaussian state space model using Kalman filtering and a fully-adapted particle filter. The code estimates the bias and mean squared error (compared with the Kalman estimate) while varying the number of particles in the particle filter.




The Kalman filter is a standard implementation without an input. The particle filter is fully adapted (i.e. takes the current observation into account when proposing new particles and computing the weights).


Returns a plot with the generated observations y and the difference in the state estimates obtained by the Kalman filter (the optimal solution) and the particle filter (with 20 particles). Furthermore, the function returns plots of the estimated bias and mean squared error of the state estimate obtained using the particle filter (while varying the number of particles) and the Kalman estimates.

The function returns a list with the elements:


See Section 3.2 in the reference for more details.


Johan Dahlin <uni (at) johandahlin.com.nospam>


Dahlin, J. & Sch<c3><b6>n, T. B. "Getting started with particle Metropolis-Hastings for inference in nonlinear dynamical models." pre-print, arXiv:1511.01707, 2017.



pmhtutorial documentation built on Oct. 7, 2017, 5:02 p.m.