Title

Standard Errors For Em Estimation

Keywords

Asymptotic variance-covariance matrix; EM algorithm; Numerical differentiation; Observed information; Precision; SEM algorithm; Slow convergence

Abstract

The EM algorithm is a popular method for computing maximum likelihood estimates. One of its drawbacks is that it does not produce standard errors as a by-product. We consider obtaining standard errors by numerical differentiation. Two approaches are considered. The first differentiates the Fisher score vector to yield the Hessian of the log-likelihood. The second differentiates the EM operator and uses an identity that relates its derivative to the Hessian of the log-likelihood. The well-known SEM algorithm uses the second approach. We consider three additional algorithms: one that uses the first approach and two that use the second. We evaluate the complexity and precision of these three and the SEM algorithm in seven examples. The first is a single-parameter example used to give insight. The others are three examples in each of two areas of EM application: Poisson mixture models and the estimation of covariance from incomplete data. The examples show that there are algorithms that are much simpler and more accurate than the SEM algorithm. Hopefully their simplicity will increase the availability of standard error estimates in EM applications. It is shown that, as previously conjectured, a symmetry diagnostic can accurately estimate errors arising from numerical differentiation. Some issues related to the speed of the EM algorithm and algorithms that differentiate the EM operator are identified.

Publication Date

1-1-2000

Publication Title

Journal of the Royal Statistical Society. Series B: Statistical Methodology

Volume

62

Issue

2

Number of Pages

257-270

Document Type

Article

Personal Identifier

scopus

DOI Link

https://doi.org/10.1111/1467-9868.00230

Socpus ID

0034365803 (Scopus)

Source API URL

https://api.elsevier.com/content/abstract/scopus_id/0034365803

This document is currently not available here.

Share

COinS