IDEAS home Printed from https://ideas.repec.org/a/eee/csdana/v91y2015icp4-19.html
   My bibliography  Save this article

Convergent stochastic Expectation Maximization algorithm with efficient sampling in high dimension. Application to deformable template model estimation

Author

Listed:
  • Allassonnière, Stéphanie
  • Kuhn, Estelle

Abstract

Estimation in the deformable template model is a big challenge in image analysis. The issue is to estimate an atlas of a population. This atlas contains a template and the corresponding geometrical variability of the observed shapes. The goal is to propose an accurate estimation algorithm with low computational cost and with theoretical guaranties of relevance. This becomes very demanding when dealing with high dimensional data, which is particularly the case of medical images. The use of an optimized Monte Carlo Markov Chain method for a stochastic Expectation Maximization algorithm, is proposed to estimate the model parameters by maximizing the likelihood. A new Anisotropic Metropolis Adjusted Langevin Algorithm is used as transition in the MCMC method. First it is proven that this new sampler leads to a geometrically uniformly ergodic Markov chain. Furthermore, it is proven also that under mild conditions, the estimated parameters converge almost surely and are asymptotically Gaussian distributed. The methodology developed is then tested on handwritten digits and some 2D and 3D medical images for the deformable model estimation. More widely, the proposed algorithm can be used for a large range of models in many fields of applications such as pharmacology or genetic. The technical proofs are detailed in an appendix.11The appendix is available as supplementary material (see Appendix A).

Suggested Citation

  • Allassonnière, Stéphanie & Kuhn, Estelle, 2015. "Convergent stochastic Expectation Maximization algorithm with efficient sampling in high dimension. Application to deformable template model estimation," Computational Statistics & Data Analysis, Elsevier, vol. 91(C), pages 4-19.
    • Handle: RePEc:eee:csdana:v:91:y:2015:i:c:p:4-19
    DOI: 10.1016/j.csda.2015.04.011
    as

    Download full text from publisher

    File URL: http://www.sciencedirect.com/science/article/pii/S0167947315001164
    Download Restriction: Full text for ScienceDirect subscribers only.
    File URL: https://libkey.io/10.1016/j.csda.2015.04.011?utm_source=ideas
    LibKey link: if access is restricted and if your library uses this service, LibKey will redirect you to where you can use your library subscription to access this item
    ---><---

    As the access to this document is restricted, you may want to search for a different version of it.

    References listed on IDEAS

    as
    1. C. A. Glasbey & K. V. Mardia, 2001. "A penalized likelihood approach to image warping," Journal of the Royal Statistical Society Series B, Royal Statistical Society, vol. 63(3), pages 465-492.
    2. O. Stramer & R. L. Tweedie, 1999. "Langevin-Type Models II: Self-Targeting Candidates for MCMC Algorithms," Methodology and Computing in Applied Probability, Springer, vol. 1(3), pages 307-328, October.
    3. Jarner, Søren Fiig & Hansen, Ernst, 2000. "Geometric ergodicity of Metropolis algorithms," Stochastic Processes and their Applications, Elsevier, vol. 85(2), pages 341-361, February.
    4. Su, J. & Srivastava, A. & Huffer, F.W., 2013. "Detection, classification and estimation of individual shapes in 2D and 3D point clouds," Computational Statistics & Data Analysis, Elsevier, vol. 58(C), pages 227-241.
    5. O. Stramer & R. L. Tweedie, 1999. "Langevin-Type Models I: Diffusions with Given Stationary Distributions and their Discretizations," Methodology and Computing in Applied Probability, Springer, vol. 1(3), pages 283-306, October.
    6. Mark Girolami & Ben Calderhead, 2011. "Riemann manifold Langevin and Hamiltonian Monte Carlo methods," Journal of the Royal Statistical Society Series B, Royal Statistical Society, vol. 73(2), pages 123-214, March.
    7. S. Allassonnière & Y. Amit & A. Trouvé, 2007. "Towards a coherent statistical framework for dense deformable template estimation," Journal of the Royal Statistical Society Series B, Royal Statistical Society, vol. 69(1), pages 3-29, February.
    8. Kuhn, E. & Lavielle, M., 2005. "Maximum likelihood estimation in nonlinear mixed effects models," Computational Statistics & Data Analysis, Elsevier, vol. 49(4), pages 1020-1038, June.
    Full references (including those not matched with items on IDEAS)

    Citations

    Citations are extracted by the CitEc Project, subscribe to its RSS feed for this item.
    as


    Cited by:

    1. Maire, Florian & Moulines, Eric & Lefebvre, Sidonie, 2017. "Online EM for functional data," Computational Statistics & Data Analysis, Elsevier, vol. 111(C), pages 27-47.

    Most related items

    These are the items that most often cite the same works as this one and are cited by the same works as this one.
    1. Dalalyan, Arnak S. & Karagulyan, Avetik, 2019. "User-friendly guarantees for the Langevin Monte Carlo with inaccurate gradient," Stochastic Processes and their Applications, Elsevier, vol. 129(12), pages 5278-5311.
    2. Fort, G. & Moulines, E., 2003. "Polynomial ergodicity of Markov transition kernels," Stochastic Processes and their Applications, Elsevier, vol. 103(1), pages 57-99, January.
    3. Arnak S. Dalalyan, 2017. "Theoretical guarantees for approximate sampling from smooth and log-concave densities," Journal of the Royal Statistical Society Series B, Royal Statistical Society, vol. 79(3), pages 651-676, June.
    4. O. F. Christensen & J. Møller & R. P. Waagepetersen, 2001. "Geometric Ergodicity of Metropolis-Hastings Algorithms for Conditional Simulation in Generalized Linear Mixed Models," Methodology and Computing in Applied Probability, Springer, vol. 3(3), pages 309-327, September.
    5. G. O. Roberts & O. Stramer, 2002. "Langevin Diffusions and Metropolis-Hastings Algorithms," Methodology and Computing in Applied Probability, Springer, vol. 4(4), pages 337-357, December.
    6. Fort, Gersende & Moulines, Eric, 2000. "V-Subgeometric ergodicity for a Hastings-Metropolis algorithm," Statistics & Probability Letters, Elsevier, vol. 49(4), pages 401-410, October.
    7. Yves F. Atchadé, 2006. "An Adaptive Version for the Metropolis Adjusted Langevin Algorithm with a Truncated Drift," Methodology and Computing in Applied Probability, Springer, vol. 8(2), pages 235-254, June.
    8. Samuel Livingstone, 2021. "Geometric Ergodicity of the Random Walk Metropolis with Position-Dependent Proposal Covariance," Mathematics, MDPI, vol. 9(4), pages 1-14, February.
    9. Ranjan, Rakesh & Sen, Rijji & Upadhyay, Satyanshu K., 2021. "Bayes analysis of some important lifetime models using MCMC based approaches when the observations are left truncated and right censored," Reliability Engineering and System Safety, Elsevier, vol. 214(C).
    10. Bournakis, Ioannis & Tsionas, Mike G., 2023. "A Non-Parametric Estimation of Productivity with Idiosyncratic and Aggregate Shocks: The Role of Research and Development (R&D) and Corporate Tax," MPRA Paper 118100, University Library of Munich, Germany.
    11. Chen, Zhongfei & Wanke, Peter & Tsionas, Mike G., 2018. "Assessing the strategic fit of potential M&As in Chinese banking: A novel Bayesian stochastic frontier approach," Economic Modelling, Elsevier, vol. 73(C), pages 254-263.
    12. Atkinson, Scott E. & Tsionas, Mike G., 2021. "Generalized estimation of productivity with multiple bad outputs: The importance of materials balance constraints," European Journal of Operational Research, Elsevier, vol. 292(3), pages 1165-1186.
    13. Maire, Florian & Moulines, Eric & Lefebvre, Sidonie, 2017. "Online EM for functional data," Computational Statistics & Data Analysis, Elsevier, vol. 111(C), pages 27-47.
    14. Caroline Khan & Mike G. Tsionas, 2021. "Constraints in models of production and cost via slack-based measures," Empirical Economics, Springer, vol. 61(6), pages 3347-3374, December.
    15. Olivier Collier & Arnak S, Dalalyan, 2013. "Curve registration by Nonparametric goodness-of-fit Testing," Working Papers 2013-33, Center for Research in Economics and Statistics.
    16. Liu, Jia & Maheu, John M & Song, Yong, 2023. "Identification and Forecasting of Bull and Bear Markets using Multivariate Returns," MPRA Paper 119515, University Library of Munich, Germany.
    17. Dimitrakopoulos, Stefanos & Tsionas, Mike, 2019. "Ordinal-response GARCH models for transaction data: A forecasting exercise," International Journal of Forecasting, Elsevier, vol. 35(4), pages 1273-1287.
    18. Vanhatalo, Jarno & Veneranta, Lari & Hudd, Richard, 2012. "Species distribution modeling with Gaussian processes: A case study with the youngest stages of sea spawning whitefish (Coregonus lavaretus L. s.l.) larvae," Ecological Modelling, Elsevier, vol. 228(C), pages 49-58.
    19. Stephen G. Hall & Heather D. Gibson & G. S. Tavlas & Mike G. Tsionas, 2020. "A Monte Carlo Study of Time Varying Coefficient (TVC) Estimation," Computational Economics, Springer;Society for Computational Economics, vol. 56(1), pages 115-130, June.
    20. Will Penny & Biswa Sengupta, 2016. "Annealed Importance Sampling for Neural Mass Models," PLOS Computational Biology, Public Library of Science, vol. 12(3), pages 1-25, March.

    Corrections

    All material on this site has been provided by the respective publishers and authors. You can help correct errors and omissions. When requesting a correction, please mention this item's handle: RePEc:eee:csdana:v:91:y:2015:i:c:p:4-19. See general information about how to correct material in RePEc.

    If you have authored this item and are not yet registered with RePEc, we encourage you to do it here. This allows to link your profile to this item. It also allows you to accept potential citations to this item that we are uncertain about.

    If CitEc recognized a bibliographic reference but did not link an item in RePEc to it, you can help with this form .

    If you know of missing items citing this one, you can help us creating those links by adding the relevant references in the same way as above, for each refering item. If you are a registered author of this item, you may also want to check the "citations" tab in your RePEc Author Service profile, as there may be some citations waiting for confirmation.

    For technical questions regarding this item, or to correct its authors, title, abstract, bibliographic or download information, contact: Catherine Liu (email available below). General contact details of provider: http://www.elsevier.com/locate/csda .

    Please note that corrections may take a couple of weeks to filter through the various RePEc services.

    IDEAS is a RePEc service. RePEc uses bibliographic data supplied by the respective publishers.