Abstract:The idea of data augmentation arises naturally in missing value problems, as exemplified by the standard ways of filling in missing cells in balanced two-way tables. Thus data augmentation refers to a scheme of augmenting the observed data so as to make it more easy to analyze. This device is used to great advantage by the EM algorithm (Dempster, Laird, and Rubin 1977) in solving maximum likelihood problems. In situations when the likelihood cannot be approximated closely by the normal likelihood, maximum likelihood estimates and the associated standard errors cannot be relied upon to make valid inferential statements. From the Bayesian point of view, one must now calculate the posterior distribution of parameters of interest. If data augmentation can be used in the calculation of the maximum likelihood estimate, then in the same cases one ought to be able to use it in the computation of the posterior distribution. It is the purpose of this article to explain how this can be done. The basic idea is quite simple. The observed data y is augmented by the quantity z, which is referred to as the latent data. It is assumed that if y and z are both known, then the problem is straightforward to analyze, that is, the augmented data posterior p(θ | y, z) can be calculated. But the posterior density that we want is p(θ | y), which may be difficult to calculate directly. If, however, one can generate multiple values of z from the predictive distribution p(z | y) (i.e., multiple imputations of z), thenp(θ | y) can be approximately obtained as the average ofp(θ | y, z) over the imputed z's. However, p(z | y) depends, in turn, onp(θ | y). Hence if p(θ | y) was known, it could be used to calculate p(z | y). This mutual dependency between p(θ | y) and p(z | y) leads to an iterative algorithm to calculate p(θ | y). Analytically, this algorithm is essentially the method of successive substitution for solving an operator fixed point equation. We exploit this fact to prove convergence under mild regularity conditions. Typically, to implement the algorithm, one must be able to sample from two distributions, namely p(θ | y, z) andp(z | θ, y). In many cases, it is straightforward to sample from either distribution. In general, though, either sampling can be difficult, just as either the E or the M step can be difficult to implement in the EM algorithm. For p(θ | y, z) arising from parametric submodels of the multinomial, we develop a primitive but generally applicable way to approximately sample θ. The idea is first to sample from the posterior distribution of the cell probabilities and then to project to the parametric surface that is specified by the submodel, giving more weight to those observations lying closer to the surface. This procedure should cover many of the common models for categorical data. There are several examples given in this article. First, the algorithm is introduced and motivated in the context of a genetic linkage example. Second, we apply this algorithm to an example of inference from incomplete data regarding the correlation coefficient of the bivariate normal distribution. It is seen that the algorithm recovers the bimodal nature of the posterior distribution. Finally, the algorithm is used in the analysis of the traditional latent-class model as applied to data from the General Social Survey.

Three Contributions to Statistical Computing

Forthcoming Papers

Stochastic Expectation Maximization with Variance Reduction.

The ECME Algorithm: A Simple Extension of EM and ECM with Faster Monotone Convergence

The Calculation of Posterior Distributions by Data Augmentation

Topics in Transformation-based Statistical Methods

Efficient data augmentation techniques for some classes of state space models

Using Em to Obtain Asymptotic Variance - Covariance Matrices - the Sem Algorithm

Regularized EM Algorithms: A Unified Framework and Statistical Guarantees

Central Limit Theorems of a Recursive Stochastic Algorithm with Applications to Adaptive Designs

Acceleration of the EM Algorithm

From EM to Data Augmentation: The Emergence of MCMC Bayesian Computation in the 1980s

Fast Mixing of Data Augmentation Algorithms: Bayesian Probit, Logit, and Lasso Regression

ML-ESTIMATION OF THE T-DISTRIBUTION USING EM AND ITS EXTENSIONS, ECM AND ECME

A Tutorial on MM Algorithms

Applying the Expectation-Maximization Algorithm to Multivariate Signal Extraction

Big Learning Expectation Maximization

Convergence analysis of data augmentation algorithms for Bayesian robust multivariate linear regression with incomplete data

Efficient Data Augmentation in Dynamic Models for Binary and Count Data

An Asynchronous Distributed Expectation Maximization Algorithm for Massive Data: The DEM Algorithm

Acceleration of the EM algorithm using the vector epsilon algorithm