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.

Bayesian inference for transductive learning of kernel matrix using the Tanner-Wong data augmentation algorithm

Model-based Transductive Learning of the Kernel Matrix

A Kernel Theory of Modern Data Augmentation

Bayesian tensor-on-tensor regression with efficient computation

A Bayesian Approach To Analysing Training Data Attribution In Deep Learning

A Bayesian Matrix Factorization Model for Relational Data

Bayesian Tensor Analysis

Hebbian Learning with Kernel-Based Embedding of Input Data

A Data-Adaptive Prior for Bayesian Learning of Kernels in Operators

Scalable Bayesian Tensor Ring Factorization for Multiway Data Analysis

Bayesian Kernelized Tensor Factorization as Surrogate for Bayesian Optimization

Bayesian Max-margin Multi-Task Learning with Data Augmentation.

Bayesian Tensor Modeling for Image-based Classification of Alzheimer's Disease

The Calculation of Posterior Distributions by Data Augmentation

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

Bayesian Nonparametric Models for Multiway Data Analysis

Convex Point Estimation using Undirected Bayesian Transfer Hierarchies

Data transforming augmentation for heteroscedastic models

Estimation of Distribution Algorithms with Matrix Transpose in Bayesian Learning

Bayesian Method for Learning Graphical Models with Incompletely Categorical Data

Bayesian learning of joint distributions of objects