Abstract
In statistical models involving constrained or missing data, likelihoods containing integrals emerge. In the case of both constrained and missing data, the result is a ratio of integrals, which for multivariate data may defy exact or approximate analytic expression. Seeking maximum‐likelihood estimates in such settings, we propose Monte Carlo approximants for these integrals, and subsequently maximize the resulting approximate likelihood. Iteration of this strategy expedites the maximization, while the Gibbs sampler is useful for the required Monte Carlo generation. As a result, we handle a class of models broader than the customary EM setting without using an EM‐type algorithm. Implementation of the methodology is illustrated in two numerical examples.
Original language | English (US) |
---|---|
Pages (from-to) | 303-311 |
Number of pages | 9 |
Journal | Canadian Journal of Statistics |
Volume | 21 |
Issue number | 3 |
DOIs | |
State | Published - Sep 1993 |
Externally published | Yes |
Keywords
- EM algorithm
- Gibbs sampler
- Monte Carlo approximant