MONTE CARLO SMOOTHING FOR NON-LINEAR TIME SERIES
Simon Godsill and Arnaud Doucet
University of Cambridge, UK
and
Mike West
ISDS
October 2002,
based on the original
September 2000 version listed (and still available as ISDS
Discussion Paper #00-01). The current version is published in
JASA 2004
We develop methods for performing smoothing computations in general state-space models. The methods rely on a particle representation of the filtering distributions, and their evolution through time using sequential importance sampling and resampling ideas. In particular, novel techniques are presented for generation of sample realizations of historical state sequences. This is carried out in a forward-filtering backward-smoothing procedure which can be viewed as the non-linear, non-Gaussian counterpart of standard Kalman filter-based simulation smoothers in the linear Gaussian case. Convergence in the mean-squared error sense of the smoothed trajectories is proved, showing the validity of our proposed method. The methods are tested in a substantial application for the processing of speech signals represented by a time-varying autoregression and parameterised in terms of time-varying partial correlation coefficients, comparing the results of our algorithm with those from a simple smoother based upon the filtered trajectories.
The manuscript is available in pdf format