Tutorials › Multiple Imputation by Chained Equations (MICE) › Single Imputation and Missing-Value Uncertainty

Why Single Imputation Fails · Tutorial 35 of 497

Single Imputation and Missing-Value Uncertainty

See why one filled-in value per missing observation cannot, by itself, represent the uncertainty that remains in the missing data.

Beginner 9 min read

What You'll Learn

  • Distinguish a completed dataset from the uncertainty about the missing values used to create it.
  • Calculate how fixed and randomly drawn imputations can produce different estimates.
  • Use the variance decomposition to understand what a single completed-data analysis leaves out.
  • Recognize assumptions and warning signs that affect single-imputation results.
  • Compare single imputation with complete-case analysis and multiple imputation.

One completed dataset is not the same as complete information

Replacing a missing value makes a dataset easier to analyze, but it does not reveal what that value actually was. A completed dataset contains a particular substitute; the uncertainty that the substitute might be wrong is not automatically carried into an ordinary analysis of that dataset. This distinction matters because standard errors and confidence intervals are meant to reflect uncertainty, not just the arithmetic consequences of the values currently in the data table.

Single imputation means filling each missing entry once and then analyzing the resulting completed dataset as if those entries had been observed. The replacement may be a fixed value, such as a mean or a prediction, or one random draw from a model for the missing value. These choices differ in how they select the replacement, but one completed dataset still represents only one particular completion.

Definition: Let \(D_{\mathrm{obs}}\) denote the observed data and \(Y_{\mathrm{mis}}\) the unobserved values. An imputation assigns values to \(Y_{\mathrm{mis}}\), producing a completed dataset. Single imputation uses one such assignment for the analysis. Missing-value uncertainty is the uncertainty about the unobserved values that remains after conditioning on the observed data and the assumptions used to model them.

An imputation is not an observation. It is a model-based substitute. A fixed imputation, such as a conditional mean, may summarize the model's expected missing value, but it does not show the range of plausible values around that expectation. A random imputation selects one plausible value, but one draw does not show how the analysis would change across other plausible draws.

For inference to be meaningful, the imputation model must be appropriate for the variable and its relationships with observed information. The assumptions about why data are missing also matter: an imputation based on observed predictors typically needs a missingness assumption such as MAR, along with a correctly specified conditional model. Under MNAR, the distribution of missing values may depend on values that are themselves unobserved, so additional assumptions are required. These assumptions cannot be established merely by completing the dataset.

How one imputation can conceal alternatives

Consider an estimate \(\hat{\theta}\), such as a sample mean or a treatment-effect estimate. If a missing value could plausibly take several values, then the estimate may differ across completed datasets. An analysis of only one completion calculates its estimate and usually its standard error conditional on that particular set of filled-in values. Unless the analysis separately accounts for uncertainty across plausible imputations, that standard error does not capture all the uncertainty created by missing data.

The idea can be expressed with the law of total variance. In a joint probability model, let \(\theta\) be the target parameter. Conditioning on the observed data and then considering the unknown values gives:

$$ \operatorname{Var}(\theta\mid D_{\mathrm{obs}}) = E\!\left[\operatorname{Var}(\theta\mid D_{\mathrm{obs}},Y_{\mathrm{mis}})\mid D_{\mathrm{obs}}\right] + \operatorname{Var}\!\left(E[\theta\mid D_{\mathrm{obs}},Y_{\mathrm{mis}}]\mid D_{\mathrm{obs}}\right). $$

The first term describes average uncertainty remaining when a particular set of missing values is specified. The second describes how the parameter estimate or its conditional expectation varies as the missing values vary over plausible completions. The identity follows by adding and subtracting the conditional expectation:

$$ \theta-E(\theta\mid D_{\mathrm{obs}}) = \{\theta-E(\theta\mid D_{\mathrm{obs}},Y_{\mathrm{mis}})\} + \{E(\theta\mid D_{\mathrm{obs}},Y_{\mathrm{mis}})-E(\theta\mid D_{\mathrm{obs}})\}. $$

The two bracketed terms have zero conditional cross-product in expectation, because the first has conditional expectation zero given \(D_{\mathrm{obs}}\) and \(Y_{\mathrm{mis}}\). Squaring, taking expectations given \(D_{\mathrm{obs}}\), and expanding therefore yields the variance decomposition above. A single completed dataset supplies just one set of \(Y_{\mathrm{mis}}\); it does not, by itself, estimate the variation represented by the second term.

This decomposition is a way to understand the issue, not a claim that every frequentist analysis uses this exact parameter-variance formula. In practice, the key distinction remains: uncertainty calculated conditional on one completion is not the same as uncertainty that also reflects alternatives for the missing values.

Worked Example: A Fixed Substitute for a Missing Measurement

Suppose a hypothetical dataset contains measurements of \(4\), \(6\), and \(8\) units, plus one missing measurement. For a simple illustration, assume a predictive model gives equal weight to possible values of \(2\) and \(10\) units. Its conditional mean is:

$$ \frac{2+10}{2}=6 \text{ units}. $$

If the missing entry is replaced by this fixed mean, the completed-data average is:

$$ \frac{4+6+8+6}{4}=\frac{24}{4}=6 \text{ units}. $$

But the model's two possible values produce different completed-data averages. Substituting \(2\) gives \((4+6+8+2)/4=20/4=5\) units. Substituting \(10\) gives \((4+6+8+10)/4=28/4=7\) units. Thus, the fixed replacement gives an average of \(6\) units, while these two plausible completions give averages of \(5\) and \(7\) units. The fixed value is not evidence that the missing measurement was exactly \(6\); it hides the fact that the estimate depends on the unobserved value.

Why one random draw is not enough

Random imputation can represent uncertainty in the value being filled in better than a deterministic prediction: the imputed value is allowed to vary according to a specified predictive distribution. Yet after drawing once, the analyst still has only one completed dataset. The draw's randomness is not the same as an analysis that incorporates the spread of results across possible draws. In addition, a random draw from a model with estimated parameters may fail to account for uncertainty in those parameters unless that uncertainty is included in the draw procedure.

Worked Example: A Single Random Draw for a Missing Outcome

In a hypothetical group of six patients, five recorded outcome values are \(10\), \(12\), \(14\), \(16\), and \(18\) units. Suppose the conditional predictive model for the missing sixth value assigns probability \(0.5\) to \(8\) units and probability \(0.5\) to \(20\) units. The observed values sum to \(70\) units.

If the one random draw is \(8\), the completed-data mean is:

$$ \frac{70+8}{6}=\frac{78}{6}=13 \text{ units}. $$

If instead the draw is \(20\), the mean is:

$$ \frac{70+20}{6}=\frac{90}{6}=15 \text{ units}. $$

The model-implied expected missing value is \((8+20)/2=14\) units, which would yield a mean of \((70+14)/6=14\) units if inserted as a fixed value. That number is between the two random-draw results, but it does not represent the variation between them. A single random completion reports one of the possible means, \(13\) or \(15\) units, and does not reveal the other possibility or the model's uncertainty about which is appropriate.

Worked Example: Missing Binary Responses

Imagine a hypothetical survey of 10 participants about whether they improved after a program. Six of the eight observed responses indicate improvement; two responses are missing. For illustration, assume each missing response has conditional probability \(0.5\) of indicating improvement, independently of the other, given the information in the model.

The number of improvements among the two missing responses can be \(0\), \(1\), or \(2\). Under the stated assumptions, those possibilities have probabilities \(0.25\), \(0.50\), and \(0.25\), respectively. The completed-data proportions are therefore:

  • Zero additional improvements: \(6/10=0.60\), or \(60\%\).
  • One additional improvement: \(7/10=0.70\), or \(70\%\).
  • Two additional improvements: \(8/10=0.80\), or \(80\%\).

The expected number of improvements among the two missing responses is \(2(0.5)=1\), so the expected completed-data proportion is \(7/10=0.70\), or \(70\%\). If a single random imputation happens to assign one improvement, the completed dataset reports \(70\%\). That result alone does not communicate that, under this illustrative model, completed proportions of \(60\%\) and \(80\%\) also remain possible. The probabilities here depend on the explicitly stated conditional model; they are not facts learned from the two missing responses.

Assumptions, diagnostics, and alternatives

A single-imputation result depends on both the imputation rule and its assumptions. A fixed substitute can be especially misleading when the missing values differ systematically from the observed values after accounting for the predictors used. A random substitute can also be misleading if the predictive distribution is poorly specified, relevant predictors are omitted, or uncertainty about model parameters is ignored. No random number generator can correct a misspecified model or identify an untestable missingness mechanism.

Useful checks include comparing the observed data with the values the imputation model predicts, examining whether predictions are plausible on the variable's scale, and checking whether conclusions change under reasonable alternative imputation assumptions. For random imputations, rerunning the draw and noting how estimates change can reveal sensitivity to the particular draw. These checks do not prove the assumptions are right; they can expose instability or implausible imputations.

Single imputation also has practical limitations. A completed dataset may make later analyses convenient, but reporting an ordinary standard error or confidence interval from it as if every value had been observed can understate uncertainty. The point estimate can be distorted as well: the expected substitute need not preserve the target relationship or estimate when missingness and outcome structure are more complex. Neither a fixed substitute nor a random draw is automatically unbiased.

ApproachWhat it doesKey limitation or distinction
Complete-case analysisAnalyzes records with no missing values for the variables in the analysis.Uses fewer records and can be biased unless its assumptions justify excluding incomplete cases.
Single imputationFills each missing entry once and analyzes one completed dataset.One completion generally does not carry forward the uncertainty across plausible missing values.
Multiple imputationCreates and analyzes multiple completed datasets, then combines results using an appropriate pooling procedure.Still relies on a suitable imputation model and missingness assumptions; it is not a cure for model misspecification.

The central comparison is not that single imputation is always worse for every descriptive purpose. A filled dataset can be useful for a clearly limited task, and a fixed value may be a transparent convention when its purpose is not inferential. The problem arises when one imputed dataset is treated as though it contains the same information as the unobserved complete data, particularly when uncertainty statements are then based only on that one completion.

Common mistake

  • Wrong reasoning: “I used a random draw, so the uncertainty from missing values is included in my standard error.”
  • Correction: A random draw reflects one realization from the imputation model. An analysis of that single completed dataset still cannot show how estimates vary over other plausible draws. Accounting for missing-value uncertainty requires an inferential method that includes that variation, not merely a randomized replacement.
Key takeaway: A completed dataset contains the values chosen for one completion, not proof that those values are known. Fixed imputations suppress alternatives by construction; one random imputation selects one alternative. Interpret uncertainty from a single completed-data analysis accordingly.

Check Your Understanding

Answer each question, then compare with the brief solutions below.

  1. A missing measurement is replaced by the conditional mean predicted by a model. Does this show that the unobserved measurement equals that mean? Why or why not?
  2. In the six-patient example, what is the completed-data mean if the missing value is \(20\) units? How does it compare with the mean if the value is \(8\) units?
  3. In the binary-response example, how many improvements are expected among the two missing responses under the stated model, and what completed-data proportion does that imply?
  4. Why does one random imputation not, by itself, account for variation across plausible missing values?

Solutions. (1) No. The conditional mean is a model-based expected value; it does not establish the actual unobserved value, which may differ from it. (2) The mean is \((70+20)/6=90/6=15\) units, compared with \((70+8)/6=78/6=13\) units, a difference of \(2\) units. (3) The expected count is \(2(0.5)=1\) improvement, giving \((6+1)/10=0.70\), or \(70\%\). (4) One draw selects one completion. It provides no estimate of how the result would vary across other plausible completions unless that variation is separately incorporated.

Clinical Biostats