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

DOI10.1016/J.CSDA.2015.04.011MaRDI QIDQ1663188zbMATH OpenOpenAlexFDO

Authors Stéphanie Allassonnière, Estelle Kuhn

Publication date 21 August 2018

Published in Computational Statistics and Data Analysis (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1207.5938

zbMATH Keywords

maximum likelihood estimationMCMCstochastic EM algorithmhigh dimensiondeformable templategeometric variabilityanisotropic MALAmissing variable

Mathematics Subject Classification ID

Computational methods for problems pertaining to statistics (62-08) Applications of statistics to biology and medical sciences; meta analysis (62P10)

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 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. We propose to use an optimized Monte Carlo Markov Chain method into a stochastic Expectation Maximization algorithm in order to estimate the model parameters by maximizing the likelihood. In this paper, we present a new Anisotropic Metropolis Adjusted Langevin Algorithm which we use as transition in the MCMC method. We first prove that this new sampler leads to a geometrically uniformly ergodic Markov chain. We prove 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.

Recommendations

• Construction of Bayesian deformable models via a stochastic approximation algorithm: a convergence study
• A new class of stochastic EM algorithms. Escaping local maxima and handling intractable sampling
• Towards a Coherent Statistical Framework for Dense Deformable Template Estimation
• Convergence of a stochastic approximation version of the EM algorithm
• scientific article; zbMATH DE number 6458272
• Stochastic algorithm for Bayesian mixture effect template estimation
• scientific article; zbMATH DE number 788229

Cites work

• scientific article; zbMATH DE number 3567782 (Why is no real title available?)
• scientific article; zbMATH DE number 840151 (Why is no real title available?)
• scientific article; zbMATH DE number 6458272 (Why is no real title available?)
• A penalized likelihood approach to image warping. (With discussion).
• An adaptive approach to Langevin MCMC
• An adaptive version for the Metropolis adjusted Langevin algorithm with a truncated drift
• Computational anatomy: an emerging discipline
• Construction of Bayesian deformable models via a stochastic approximation algorithm: a convergence study
• Convergence of a stochastic approximation version of the EM algorithm
• Coupling a stochastic approximation version of EM with an MCMC procedure
• Detection, classification and estimation of individual shapes in 2D and 3D point clouds
• Exponential convergence of Langevin distributions and their discrete approximations
• Geometric ergodicity of Metropolis algorithms
• Langevin-type models II: Self-targeting candidates for MCMC algorithms
• Langevin-type models. I: Diffusions with given stationary distributions and their discretizations
• Markov chains and stochastic stability
• Maximum likelihood estimation in nonlinear mixed effects models
• On the consistency of Fréchet means in deformable models for curve and image analysis
• On the ergodicity properties of some adaptive MCMC algorithms
• Riemann manifold Langevin and Hamiltonian Monte Carlo methods. With discussion and authors' reply
• Stability of Stochastic Approximation under Verifiable Conditions
• Stochastic algorithm for Bayesian mixture effect template estimation
• Towards a Coherent Statistical Framework for Dense Deformable Template Estimation

Cited in

(8)

• Estimating blood vessel areas in ultrasound images using a deformable template model
• Construction of Bayesian deformable models via a stochastic approximation algorithm: a convergence study
• Stochastic algorithm for Bayesian mixture effect template estimation
• Template shape estimation: correcting an asymptotic bias
• Online EM for functional data
• On the consistency of Fréchet means in deformable models for curve and image analysis
• Bayesian mixed effect atlas estimation with a diffeomorphic deformation model
• Statistical models for deformable templates in image and shape analysis

Describes a project that uses

Uses Software

• LDDMM

This page was built for publication: Convergent stochastic expectation maximization algorithm with efficient sampling in high dimension. Application to deformable template model estimation

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q1663188)