Multiple Imputation Using Gaussian Copulas

Abstract

Missing observations are pervasive throughout empirical research, especially in the social sciences. Despite multiple approaches to dealing adequately with missing data, many scholars still fail to address this vital issue. In this paper, we present a simple-to-use method for generating multiple imputations using a Gaussian copula. The Gaussian copula for multiple imputation (Hoff, 2007) allows scholars to attain estimation results that have good coverage and small bias. The use of copulas to model the dependence among variables will enable researchers to construct valid joint distributions of the data, even without knowledge of the actual underlying marginal distributions. Multiple imputations are then generated by drawing observations from the resulting posterior joint distribution and replacing the missing values. Using simulated and observational data from published social science research, we compare imputation via Gaussian copulas with two other widely used imputation methods: MICE and Amelia II. Our results suggest that the Gaussian copula approach has a slightly smaller bias, higher coverage rates, and narrower confidence intervals compared to the other methods. This is especially true when the variables with missing data are not normally distributed. These results, combined with theoretical guarantees and ease-of-use suggest that the approach examined provides an attractive alternative for applied researchers undertaking multiple imputations.

An observational study using eye tracking to assess resident and senior anesthetists’ situation awareness and visual perception in postpartum hemorrhage high fidelity simulation

“…The variables extracted were the number of fixations (NBF) and the fixation duration mean (FDM) for each region of interest. Missing data were computed using MCMC method [19,20]. The attentional heat maps were extracted during key steps of care.…”

Coalitions and conflict: A longitudinal analysis of men's politics

“…Missing values (total 8.5%) for status (2%) and physical formidability (11.5%) were imputed using a Bayesian copula approach (Hoff, 2018). See (Hollenbach et al, 2018) for an overview and comparison with other multiple imputation methods.…”

Why Ordinal Variables Can (Almost) Always Be Treated as Continuous Variables: Clarifying Assumptions of Robust Continuous and Ordinal Factor Analysis Estimation Methods

“…If substantial deviations of MCAR are suspected, and a misspecified factor model should be estimated, imputation based approaches might be preferable (Gottschall et al, 2012;Jia and Wu, 2019). The use of sufficiently complex imputation models, such as the Gaussian copula model (Hollenbach et al, 2018), mixture models (Murray and Reiter, 2016), or latent class models (Vermunt et al, 2008;Si and Reiter, 2013) are advantageous to minimize possible distributional misspecifications for MAR data. Appropriate imputation models can also treat specific deviations from MAR (missing not at random; MNAR; Harel and Schafer, 2009;Jung et al, 2011;Kano and Takai, 2011;Zhang and Reiser, 2015;Bartolucci et al, 2018;Kuha et al, 2018;Pohl and Becker, 2020).…”