Maximum likelihood smoothing estimation in state-space models: An incomplete-information based approach

Abstract

This paper revisits classical works of Rauch (1963, et al. 1965) and develops a novel method for maximum likelihood (ML) smoothing estimation from incomplete information/data of stochastic state-space systems. Score function and conditional observed information matrices of incomplete data are introduced and their distributional identities are established. Using these identities, the ML smoother xk ns =xk f(xk,xk+1 ns, y0:nθ), k≤ n-1, is presented. The result shows that the ML smoother gives an estimate of state xk with more adherence of loglikehood having less standard errors than that of the ML state estimator xk=xk f(xk,y0:kθ), with xn ns=xn. Recursive estimation is given in terms of an EM-gradient-particle algorithm which extends the work of Lange for ML smoothing estimation. The algorithm has an explicit iteration update which lacks in (Ramadan) EM-algorithm for smoothing. A sequential Monte Carlo method is developed for valuation of the score function and observed information matrices. A recursive equation for the covariance matrix of estimation error is developed to calculate the standard errors. In the case of linear systems, the method shows that the Rauch-Tung-Striebel (RTS) smoother is a fully efficient smoothing state-estimator whose covariance matrix coincides with the Cram\'er-Rao lower bound, the inverse of expected information matrix. Furthermore, the RTS smoother coincides with the Kalman filter having less covariance matrix. Numerical studies are performed, confirming the accuracy of the main results.

0

Turn this paper into a lesson

ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…