A new approach to particle based smoothed marginal MAP

We present here a new method of finding the MAP state estimator from the weighted particles representation of marginal smoother distribution. This is in contrast to the usual practice, where the particle with the highest weight is selected as the MAP, although the latter is not necessarily the most probable state estimate. The method developed here uses only particles with corresponding filtering and smoothing weights. We apply this estimator for finding the unknown initial state of a dynamical system and addressing the parameter estimation problem

Statistical Inference for Partially Observed Markov Processes via the R Package pomp

Partially observed Markov process (POMP) models, also known as hidden Markov models or state space models, are ubiquitous tools for time series analysis. The R package pomp provides a very flexible framework for Monte Carlo statistical investigations using nonlinear, non-Gaussian POMP models. A range of modern statistical methods for POMP models have been implemented in this framework including sequential Monte Carlo, iterated filtering, particle Markov chain Monte Carlo, approximate Bayesian computation, maximum synthetic likelihood estimation, nonlinear forecasting, and trajectory matching. In this paper, we demonstrate the application of these methodologies using some simple toy problems. We also illustrate the specification of more complex POMP models, using a nonlinear epidemiological model with a discrete population, seasonality, and extra-demographic stochasticity. We discuss the specification of user-defined models and the development of additional methods within the programming environment provided by pomp.Comment: In press at the Journal of Statistical Software. A version of this paper is provided at the pomp package website: http://kingaa.github.io/pom