Correction to: Alternatives to the EM algorithm for ML estimation of location, scatter matrix, and degree of freedom of the Student t distribution

In this paper, we consider maximum likelihood estimations of the degree of freedom parameter $\nu$, the location parameter $\mu$ and the scatter matrix $\Sigma$ of the multivariate Student-$t$ distribution. In particular, we are interested in estimating the degree of freedom parameter $\nu$ that determines the tails of the corresponding probability density function and was rarely considered in detail in the literature so far. We prove that under certain assumptions a minimizer of the negative log-likelihood function exists, where we have to take special care of the case $\nu \rightarrow \infty$, for which the Student-$t$ distribution approaches the Gaussian distribution. As alternatives to the classical EM algorithm we propose three other algorithms which cannot be interpreted as EM algorithm. For fixed $\nu$, the first algorithm is an accelerated EM algorithm known from the literature. However, since we do not fix $\nu$, we cannot apply standard convergence results for the EM algorithm. The other two algorithms differ from this algorithm in the iteration step for $\nu$. We show how the objective function behaves for the different updates of $\nu$ and prove for all three algorithms that it decreases in each iteration step. We compare the algorithms as well as some accelerated versions by numerical simulation and apply one of them for estimating the degree of freedom parameter in images corrupted by Student-$t$ noise.

[1]  A. G. Greenhill,et al.  Handbook of Mathematical Functions with Formulas, Graphs, , 1971 .

[2]  Yousef Saad,et al.  Two classes of multisecant methods for nonlinear acceleration , 2009, Numer. Linear Algebra Appl..

[3]  Pradipta Maji,et al.  Spatially Constrained Student’s t-Distribution Based Mixture Model for Robust Image Segmentation , 2018, Journal of Mathematical Imaging and Vision.

[4]  M. Kendall A NEW MEASURE OF RANK CORRELATION , 1938 .

[5]  D. Rubin,et al.  The ECME algorithm: A simple extension of EM and ECM with faster monotone convergence , 1994 .

[6]  Tieyong Zeng,et al.  Variational Approach for Restoring Blurred Images with Cauchy Noise , 2015, SIAM J. Imaging Sci..

[7]  Yiqiu Dong,et al.  Cauchy Noise Removal by Nonconvex ADMM with Convergence Guarantees , 2017, Journal of Scientific Computing.

[8]  Jean-Michel Morel,et al.  A Nonlocal Bayesian Image Denoising Algorithm , 2013, SIAM J. Imaging Sci..

[9]  Jean-François Aujol,et al.  Estimation of the Noise Level Function Based on a Nonparametric Detection of Homogeneous Image Regions , 2015, SIAM J. Imaging Sci..

[10]  Xiao-Li Meng,et al.  The EM Algorithm—an Old Folk‐song Sung to a Fast New Tune , 1997 .

[11]  M. Kendall The treatment of ties in ranking problems. , 1945, Biometrika.

[12]  David E. Tyler,et al.  A curious likelihood identity for the multivariate t-distribution , 1994 .

[13]  J. Pesquet,et al.  Wavelet thresholding for some classes of non–Gaussian noise , 2002 .

[14]  Nikolas P. Galatsanos,et al.  Robust Image Segmentation with Mixtures of Student's t-Distributions , 2007, 2007 IEEE International Conference on Image Processing.

[15]  Donald G. M. Anderson Iterative Procedures for Nonlinear Integral Equations , 1965, JACM.

[16]  Geoffrey J. McLachlan,et al.  Robust Cluster Analysis via Mixtures of Multivariate t-Distributions , 1998, SSPR/SPR.

[17]  Serena Morigi,et al.  Whiteness Constraints in a Unified Variational Framework for Image Restoration , 2018, Journal of Mathematical Imaging and Vision.

[18]  Gabriele Steidl,et al.  Nonlocal Myriad Filters for Cauchy Noise Removal , 2017, Journal of Mathematical Imaging and Vision.

[19]  Milton Abramowitz,et al.  Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables , 1964 .

[20]  G. McLachlan,et al.  The EM algorithm and extensions , 1996 .

[21]  R. Varadhan,et al.  Simple and Globally Convergent Methods for Accelerating the Convergence of Any EM Algorithm , 2008 .

[22]  Ravi Varadhan,et al.  Damped Anderson Acceleration With Restarts and Monotonicity Control for Accelerating EM and EM-like Algorithms , 2018, Journal of Computational and Graphical Statistics.

[23]  Q. M. Jonathan Wu,et al.  Robust Student's-t Mixture Model With Spatial Constraints and Its Application in Medical Image Segmentation , 2012, IEEE Transactions on Medical Imaging.

[24]  Jian Zheng,et al.  Robust Non-Rigid Point Set Registration Using Student's-t Mixture Model , 2014, PloS one.

[25]  Kaare Brandt Petersen,et al.  The Matrix Cookbook , 2006 .

[26]  Gabriele Steidl,et al.  Multivariate Myriad Filters Based on Parameter Estimation of Student-t Distributions , 2018, SIAM J. Imaging Sci..

[27]  Geoffrey J. McLachlan,et al.  Robust mixture modelling using the t distribution , 2000, Stat. Comput..

[28]  Ting-Zhu Huang,et al.  Total variation with overlapping group sparsity for deblurring images under Cauchy noise , 2019, Appl. Math. Comput..

[29]  D. Rubin,et al.  ML ESTIMATION OF THE t DISTRIBUTION USING EM AND ITS EXTENSIONS, ECM AND ECME , 1999 .

[30]  Jeremy MG Taylor,et al.  Robust Statistical Modeling Using the t Distribution , 1989 .

[31]  Guan Gui,et al.  A Convex Constraint Variational Method for Restoring Blurred Images in the Presence of Alpha-Stable Noises , 2018, Sensors.

[32]  David M. Miller,et al.  Handbook of Mathematical Functions With Formulas, Graphs and Mathematical Tables (National Bureau of Standards Applied Mathematics Series No. 55) , 1965 .

[33]  Aristidis Likas,et al.  The mixtures of Student's t-distributions as a robust framework for rigid registration , 2009, Image Vis. Comput..

[34]  Benjamin Schrauwen,et al.  The student-t mixture as a natural image patch prior with application to image compression , 2014, J. Mach. Learn. Res..