bssm: Bayesian Inference of Non-linear and Non-Gaussian State Space Models in R

Abstract

We present an R package bssm for Bayesian non-linear/non-Gaussian state space modelling. Unlike the existing packages, bssm allows for easy-to-use approximate inference based on Gaussian approximations such as the Laplace approximation and the extended Kalman filter. The package accommodates also discretely observed latent diffusion processes. The inference is based on fully automatic, adaptive Markov chain Monte Carlo (MCMC) on the hyperparameters, with optional importance sampling post-correction to eliminate any approximation bias. The package implements also a direct pseudo-marginal MCMC and a delayed acceptance pseudo-marginal MCMC using intermediate approximations. The package offers an easy-to-use interface to define models with linear-Gaussian state dynamics with non-Gaussian observation models, and has an Rcpp interface for specifying custom non-linear and diffusion models.

References in corpus (8)

• The No-U-Turn Sampler: Adaptively Setting Path Lengths in Hamiltonian Monte Carlo
• Programming with models: writing statistical algorithms for general model structures with NIMBLE
• Statistical Inference for Partially Observed Markov Processes via the R Package pomp
• KFAS: Exponential Family State Space Models in R
• Mixture Hidden Markov Models for Sequence Data: The seqHMM Package in R
• Importance sampling type estimators based on approximate marginal MCMC
• Importance sampling correction versus standard averages of reversible MCMCs in terms of the asymptotic variance
• Sequential Monte Carlo Methods in the nimble R Package

Cited by in corpus (1)

• ssMousetrack: Analysing computerized tracking data via Bayesian state-space models in {R}