Hostname: page-component-cd9895bd7-jkksz Total loading time: 0 Render date: 2024-12-23T17:19:33.156Z Has data issue: false hasContentIssue false

Coupling a stochastic approximation version of EMwith an MCMC procedure

Published online by Cambridge University Press:  15 September 2004

Estelle Kuhn
Affiliation:
Université Paris Sud, Bât. 425, 91400 Orsay, France; [email protected].
Marc Lavielle
Affiliation:
Université René Descartes and Université Paris Sud, Bât. 425, 91400 Orsay, France; [email protected].
Get access

Abstract

The stochastic approximation version of EM (SAEM) proposed by Delyon et al. (1999) isa powerful alternative to EM when the E-step is intractable. Convergence ofSAEM toward a maximum of the observed likelihood is established whenthe unobserved data are simulated at each iteration under the conditionaldistribution. We show that this very restrictive assumption can be weakened. Indeed, the results of Benveniste et al. for stochastic approximationwith Markovian perturbations are used to establish the convergenceof SAEM when it is coupled with a Markov chain Monte-Carloprocedure. This result is very useful for many practical applications. Applications to the convolution model and the change-points model are presented to illustrate the proposed method.

Type
Research Article
Copyright
© EDP Sciences, SMAI, 2004

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

A. Benveniste, M. Métivier and P. Priouret, Adaptive algorithms and stochastic approximations. Springer-Verlag, Berlin (1990). Translated from the French by Stephen S. Wilson.
Brandière, O. and Duflo, M., Les algorithmes stochastiques contournent-ils les pièges ? C. R. Acad. Sci. Paris Ser. I Math. 321 (1995) 335338.
Chen, H.F., Lei, G. and Gao, A.J., Convergence and robustness of the Robbins-Monro algorithm truncated at randomly varying bounds. Stochastic Process. Appl. 27 (1988) 217231. CrossRef
Concordet, D. and Nunez, O.G., A simulated pseudo-maximum likelihood estimator for nonlinear mixed models. Comput. Statist. Data Anal. 39 (2002) 187201. CrossRef
Delyon, B., Lavielle, M. and Moulines, E., Convergence of a stochastic approximation version of the EM algorithm. Ann. Statist. 27 (1999) 94128.
Dempster, A.P., Laird, N.M. and Rubin, D.B., Maximum likelihood from incomplete data via the EM algorithm. J. R. Stat. Soc. Ser. B 39 (1977) 138.
Gu, M.G. and Kong, F.H., A stochastic approximation algorithm with Markov chain Monte-Carlo method for incomplete data estimation problems. Proc. Natl. Acad. Sci. USA 95 (1998) 72707274 (electronic). CrossRef
Gu, M.G. and Zhu, H.-T., Maximum likelihood estimation for spatial models by Markov chain Monte Carlo stochastic approximation. J. R. Stat. Soc. Ser. B 63 (2001) 339355. CrossRef
Lange, K., A gradient algorithm locally equivalent to the EM algorithm. J. R. Stat. Soc. Ser. B 57 (1995) 425437.
Lavielle, M. and Lebarbier, E., An application of MCMC methods to the multiple change-points problem. Signal Processing 81 (2001) 3953. CrossRef
Lavielle, M. and Moulines, E., A simulated annealing version of the EM algorithm for non-Gaussian deconvolution. Statist. Comput. 7 (1997) 229236. CrossRef
Meng, X.-L. and Rubin, D.B., Maximum likelihood estimation via the ECM algorithm: a general framework. Biometrika 80 (1993) 267278. CrossRef
Mengersen, K.L. and Tweedie, R.L., Rates of convergence of the Hastings and Metropolis algorithms. Ann. Statist. 24 (1996) 101121.
S.P. Meyn and R.L. Tweedie, Markov chains and stochastic stability, Springer-Verlag London Ltd., London. Comm. Control Engrg. Ser. (1993).
Jeff Wu, C.-F., On the convergence properties of the EM algorithm. Ann. Statist. 11 (1983) 95103.
Yao, J.-F., On recursive estimation in incomplete data models. Statistics 34 (2000) 2751 (English). CrossRef