Activity Number:
|
303
|
Type:
|
Topic Contributed
|
Date/Time:
|
Wednesday, August 14, 2002 : 10:30 AM to 12:20 PM
|
Sponsor:
|
Section on Statistical Computing*
|
Abstract - #301001 |
Title:
|
Monte Carlo Standard Errors for Markov Chain Monte Carlo EM
|
Author(s):
|
Galin Jones*+
|
Affiliation(s):
|
University of Minnesota
|
Address:
|
347 Ford Hall 224 Church St S.E., Minneapolis, Minnesota, 55455,
|
Keywords:
|
Markov chain ; Monte Carlo ; EM algorithm ; regenerative simulation ; asymptotic standard errors ; central limit theorem
|
Abstract:
|
The Monte Carlo EM (MCEM) algorithm is a popular method for maximizing intractable likelihoods. A key ingredient of successful implementation of MCEM is controlling the Monte Carlo sample size at each iteration. Here the focus will be on two methods of using regenerative simulation to assess the relevant Monte Carlo error, when the sampling mechanism in each E-step is Markov chain Monte Carlo (MCMC). The first method is to assess the error in the Q-function so that we can be assured that the ascent property enjoyed by deterministic EM also holds (with high probability) for MCEM. The second method is an extension of Booth and Hobert's (1999) method, which assesses the Monte Carlo error of the parameter estimates at each iteration. Conditions on the underlying MCMC algorithm are derived that guarantee valid standard errors are obtained in both cases.
|