Tutorials › Regression Analysis (Undergraduate) › Multiple Imputation as Repeated Plausible Completion

Missing Data in Regression Analyses · Tutorial 855 of 991

Multiple Imputation as Repeated Plausible Completion

See how repeated plausible completions lead to a regression estimate that reflects both variation within completed-data analyses and differences across imputations.

Intermediate 10 min read

What You'll Learn

  • Explain why multiple imputation creates several completed data sets rather than one replacement.
  • Describe how each completed data set is used in the same regression analysis.
  • Calculate a simple average of regression estimates across imputations.
  • Distinguish sampling uncertainty within a completed-data analysis from missing-data uncertainty across imputations.
  • Recognize conditions and limitations relevant to interpreting multiple-imputation results.

Why Complete the Data More Than Once?

When values are missing, a regression cannot use those values as if they were observed. Replacing each missing value with one fixed guess produces a single completed data set, but it hides uncertainty about what the missing values might have been. Regression imputation with residual variation takes a step toward representing that uncertainty. Multiple imputation extends the idea by creating several completed data sets, each with a different plausible set of missing values.

The goal is not to choose the completion that looks most convenient or to treat the several data sets as if they were new independent samples. They are alternative representations of the same incomplete sample. The analyst fits the same regression to each one, then combines the estimates and their uncertainty. Variation among the estimates carries information about how much the missing values matter to the result.

Definition: In multiple imputation, each missing value is replaced by a draw from a predictive distribution informed by the observed data and an imputation model. Repeating this process creates \(m\) completed data sets. The same intended analysis is carried out on every completed data set, and the resulting estimates are combined to account for both uncertainty within each analysis and uncertainty about the missing values.

The Basic Sequence

Suppose the scientific question concerns a regression of an outcome \(Y\) on a predictor \(X\). For example, \(Y\) might be a blood-pressure measurement and \(X\) might be hours of exercise per week. The intended regression specification is fixed before examining the results—for a simple linear regression, it could be \(Y_i=\beta_0+\beta_1X_i+\varepsilon_i\), where \(\beta_1\) is the population slope and \(\varepsilon_i\) represents deviations from the regression mean.

If some \(X\) or \(Y\) values are missing, multiple imputation follows this general sequence:

1
Specify an imputation model.
Use a model for the missing values that represents their plausible variation given relevant observed information. For valid use, its assumptions and specification need to be reasonable for the data and the analysis.
2
Generate several completions.
Draw plausible values for the missing entries repeatedly, retaining the observed entries. Each repetition produces one completed data set.
3
Fit the same regression in each completion.
Keep the outcome, predictors, functional form, and target coefficient the same across all \(m\) completed data sets.
4
Combine estimates and uncertainty.
Summarize the coefficient estimates across completions and account for both the regression uncertainty within each fit and the variation caused by different plausible missing values.

The imputations should be stochastic: repeated completions need not assign exactly the same values to missing entries. The draws should also reflect uncertainty in estimating the imputation model, rather than merely inserting predictions from one fixed fitted model. Otherwise the completed data sets can look more certain than the information warrants.

What the Two Sources of Uncertainty Mean

Each completed-data regression has ordinary sampling uncertainty: even if that particular completed data set were known, a coefficient estimate would vary from sample to sample. This is represented by the estimate's standard error within that regression. There is also missing-data uncertainty: different plausible values for the unobserved entries can lead to different fitted coefficients.

The distinction can be expressed using the law of total variance. Let \(Z_{\text{mis}}\) denote the missing values and let \(Q\) be a regression quantity of interest, such as a slope. Conditional on the observed data, uncertainty about \(Q\) can be separated into uncertainty remaining after the missing values are specified and uncertainty arising because those values are not known:

$$ \operatorname{Var}(Q\mid\text{observed data}) = E\!\left[\operatorname{Var}(Q\mid Z_{\text{mis}},\text{observed data})\mid\text{observed data}\right] + \operatorname{Var}\!\left(E[Q\mid Z_{\text{mis}},\text{observed data}]\mid\text{observed data}\right). $$

The first term corresponds to average within-completion uncertainty; the second corresponds to between-completion variation. Multiple imputation aims to represent both. This explains why averaging estimates alone does not finish the uncertainty calculation: the estimates' variation across imputations and their within-regression standard errors also matter. The exact pooling calculations and inferential procedures are separate details; the essential point here is that neither source of uncertainty should be discarded.

Worked Example: Averaging a Slope Across Plausible Completions

Consider a hypothetical study of a continuous health score regressed on weekly exercise hours. Exercise hours are missing for some participants. After generating \(m=3\) plausible completed data sets and fitting the same linear regression to each, the estimated slopes are \(1.8\), \(2.1\), and \(2.4\) health-score points per additional exercise hour.

A simple combined point estimate is the arithmetic mean of the three fitted slopes:

$$ \bar{Q}=\frac{1}{m}\sum_{j=1}^{m}Q_j =\frac{1.8+2.1+2.4}{3} =\frac{6.3}{3} =2.1. $$

Thus the combined slope estimate is \(2.1\) health-score points per additional weekly exercise hour. It summarizes the three plausible completed-data analyses; it is not the slope from a data set in which every missing exercise value was replaced by the same fixed guess. The interpretation remains an estimated association, not evidence that increasing exercise causes the score to change.

The same regression specification is important. If one completion uses a different predictor set or functional form, the resulting estimates no longer answer exactly the same question. Averages across unlike analyses can be difficult to interpret, and their variation cannot be attributed simply to uncertainty about missing values.

Worked Example: Reading Between-Imputation Variation

Use the three hypothetical slope estimates from the health-score example: \(1.8\), \(2.1\), and \(2.4\) points per exercise hour. Their mean is \(2.1\). To describe how much the estimates vary across completions, calculate their sample variance:

$$ B=\frac{(1.8-2.1)^2+(2.1-2.1)^2+(2.4-2.1)^2}{3-1} =\frac{0.09+0+0.09}{2} =0.09. $$

The corresponding standard deviation across estimates is \(\sqrt{0.09}=0.30\) points per exercise hour. This is a descriptive measure of between-imputation variation, not a complete standard error for the combined estimate.

Suppose the three regression standard errors are \(0.40\), \(0.35\), and \(0.45\) points per exercise hour. Their average squared value is

$$ \frac{0.40^2+0.35^2+0.45^2}{3} =\frac{0.1600+0.1225+0.2025}{3} =0.1617. $$

These within-fit standard errors describe sampling uncertainty conditional on each completed data set. The dispersion of the slopes describes an additional source of uncertainty. The values \(0.30\) and \(0.1617\) have different meanings and cannot be substituted for one another; a full pooled uncertainty calculation must combine the two components using the appropriate rules.

Conditions and Limits of Interpretation

Multiple imputation is a way to propagate uncertainty under a model; it does not make missingness harmless by itself. A common working condition is that missingness is missing at random (MAR) conditional on information included in the analysis of missingness. In this setting, the probability that an entry is missing may depend on observed information, but not on the unobserved value after conditioning on that information. If missingness still depends on the unobserved value after conditioning, a MAR-based analysis may not capture the relevant uncertainty.

The imputation model also needs to be credible. It should preserve important relationships needed for the regression analysis, and its distributional assumptions should be suitable for the variables being imputed. In addition, the draws need to represent plausible variability rather than repeated copies of a single prediction. A large number of completed data sets cannot repair a systematically inappropriate model or a missingness assumption that does not fit the setting.

A practical consequence is that multiple imputation is not a method for recovering the unknowable true values of every missing observation. The individual imputed values are plausible draws, not observed facts. The target is a regression analysis whose estimate and uncertainty appropriately reflect the available data and the assumptions used to complete it.

Worked Example: Interpreting Several Imputed Regressions

A hypothetical clinic dataset records a symptom score and a treatment-duration measure for 120 patients; treatment duration is missing for some patients. An analyst creates five completed data sets, fits the same regression of symptom score on treatment duration in each, and obtains slopes of \(-0.6\), \(-0.4\), \(-0.5\), \(-0.7\), and \(-0.3\) symptom-score points per additional treatment week.

The combined point estimate is

$$ \bar{Q}=\frac{-0.6-0.4-0.5-0.7-0.3}{5} =\frac{-2.5}{5} =-0.5 $$

symptom-score points per additional treatment week. The negative sign indicates that the fitted association is toward lower symptom scores at longer treatment durations. The estimates differ across completions, so a report should not treat \(-0.5\) as if it came from a single, fully observed data set. Its uncertainty should reflect both the within-regression standard errors and the spread of these five estimates. The association is not, on its own, a causal effect.

Why Repeated Plausible Completion Matters

A single completed data set can make the analysis appear more precise than it is because it acts as though the inserted values were known. Multiple imputation instead asks how the regression would look across several plausible completions of the same incomplete sample. When the resulting estimates are similar, the missing values have relatively little influence on that coefficient under the imputation model. When they differ substantially, the uncertainty about the missing entries is consequential and should be visible in the final inference.

A common mistake is to fit the regression in every imputed data set, average the coefficient estimates, and then report the standard error from just one of those fits. That standard error reflects uncertainty within only that completion. Correction: use a pooling procedure that incorporates within-completion sampling uncertainty and between-completion variation. The average coefficient is a useful combined point estimate, but it does not by itself provide a complete account of uncertainty.

Key takeaway: Multiple imputation creates several plausible versions of the incomplete data, applies the same regression to each, and combines the results. Its purpose is to carry uncertainty about the missing values into the regression inference, alongside ordinary sampling uncertainty—not to declare the imputed values observed or to guarantee that a missing-data assumption is correct.

Check Your Understanding

Answer these questions using the distinctions between completed data sets and sources of uncertainty.

  1. Why can one completed data set understate uncertainty about a regression coefficient?
  2. What should remain the same across the regressions fitted to multiple imputed data sets?
  3. In your own words, distinguish within-completion uncertainty from between-completion variation.
  4. Three completed-data regressions give slopes of \(0.9\), \(1.2\), and \(1.5\). Calculate their simple average and state what it represents.
  5. Does creating many imputations guarantee valid conclusions if the imputation model or missingness assumption is inappropriate? Explain.