Tutorials › Propensity Score Matching and Instrumental Variables › Balance Diagnostics With Missing Covariates

Balance Diagnostics · Tutorial 281 of 699

Balance Diagnostics With Missing Covariates

Learn to distinguish imbalance in whether a covariate is observed from imbalance in its values among observed units, and to calculate and interpret both diagnostics.

Intermediate 10 min read

What You'll Learn

  • Define missingness indicators and calculate their standardized mean differences across treatment groups.
  • Assess covariate balance among units with observed values, keeping the comparison population explicit.
  • Calculate standardized differences for binary and continuous covariates.
  • Interpret missingness and observed-value diagnostics jointly without treating either as a test of a missing-data mechanism.
  • Recognize how complete-case restriction can hide missingness imbalance or change the population being compared.
  • Identify practical diagnostics and limitations when covariate data are incomplete.

Why Missingness Needs Its Own Balance Diagnostic

A covariate can look balanced among people whose values are recorded even when the chance of having a recorded value differs substantially between treatment groups. Conversely, treatment groups may have similar missingness rates while the observed covariate values differ. These are distinct features of the data, so neither should be used as a substitute for the other.

In propensity-score analyses, balance diagnostics describe the treatment groups after the chosen matching or weighting procedure. With incomplete covariates, the diagnostic should make clear both whether a value is missing and, where the value is observed, how the observed values compare. These are descriptive checks of the analytic data, not tests of whether missingness is ignorable and not proof that treatment assignment is unconfounded.

Definition: For covariate \(X_j\), define \(R_j=1\) if its value is observed and \(R_j=0\) if it is missing; equivalently, define the missingness indicator \(M_j=1-R_j\). A missingness-balance diagnostic compares \(M_j\) between treatment groups. An observed-value diagnostic compares \(X_j\) between treatment groups among units with \(R_j=1\). State the population and weights used for each comparison.

Let \(A=1\) denote the treated group and \(A=0\) the comparison group. For a diagnostic on the full analysis cohort, estimate the missingness proportion in group \(a\) as \(p_{M,a}=P(M_j=1\mid A=a)\), or its sample or weighted analogue. A common standardized mean difference (SMD) for a binary indicator is

$$ \operatorname{SMD}_{M_j} = \frac{p_{M,1}-p_{M,0}} {\sqrt{\{p_{M,1}(1-p_{M,1})+p_{M,0}(1-p_{M,0})\}/2}}. $$

The sign indicates which group has the greater missingness proportion; the magnitude expresses the difference in units of a pooled binary standard deviation. If a covariate is observed for only some units, its observed-value SMD is computed using those with \(R_j=1\) in each group. For a continuous covariate, this is the difference in observed-group means divided by the chosen pooled standard deviation. For a binary covariate, use the group proportions among observed units and the corresponding pooled binary standard deviation. Use the same standardization convention as in the other balance summaries, and state any weighting and denominator convention.

Calculating and Interpreting the Two Diagnostics

The binary SMD formula follows directly from the variance of a Bernoulli variable. If an indicator has probability \(p\) of equaling 1, its variance is \(p(1-p)\). Averaging the two group variances and taking the square root gives a common scale for the group difference:

$$ \operatorname{Var}(M_j\mid A=a)=p_{M,a}(1-p_{M,a}), \qquad s_{\mathrm{pooled}} = \sqrt{\frac{\operatorname{Var}(M_j\mid A=1)+\operatorname{Var}(M_j\mid A=0)}{2}}. $$

Dividing \(p_{M,1}-p_{M,0}\) by this scale yields the binary-indicator SMD. The SMD is undefined if the pooled variance is zero, as can happen when the indicator is constant in both groups. In that case, report the group proportions and explain that the standardized measure cannot be calculated; do not replace the undefined value with zero.

Worked Example: Comparing Missingness Rates

Suppose a hypothetical matched cohort contains 100 treated and 100 comparison units. A baseline covariate is missing for 40 treated units and 5 comparison units. Calculate the missingness SMD, defining \(M=1\) for a missing value.

The group missingness proportions are \(p_{M,1}=40/100=0.40\) and \(p_{M,0}=5/100=0.05\). The pooled binary standardizer is

$$ \sqrt{\frac{0.40(0.60)+0.05(0.95)}{2}} = \sqrt{\frac{0.2400+0.0475}{2}} = \sqrt{0.14375} \approx 0.37914. $$

Therefore,

$$ \operatorname{SMD}_{M} = \frac{0.40-0.05}{\sqrt{0.14375}} \approx \frac{0.35}{0.379144} \approx 0.9231. $$

The treated group has a higher missingness proportion, and the difference is about 0.923 pooled standard deviations. This is a substantial descriptive discrepancy that should be made visible in the balance report. It does not establish why the values are missing, nor does it determine the appropriate missing-data method.

The observed-value comparison answers a different question: among units for whom \(X_j\) is recorded, how similar are the recorded values? The denominator must use the observed values in each group, while the report should also show how many observations contribute. If matching or weighting is used, calculate the group summaries with the analysis weights and the course’s chosen weighted-SMD convention. The full-cohort missingness comparison should use those same analysis weights, rather than silently switching to an unweighted diagnostic.

Worked Example: A Binary Covariate Among Observed Units

In a hypothetical analysis, a binary baseline covariate is observed for some units in each treatment group. Among observed units, its proportion is 0.08 in the treated group and 0.10 in the comparison group. Calculate the observed-value SMD using the pooled binary standard deviation.

The two Bernoulli variances are \(0.08(1-0.08)=0.0736\) and \(0.10(1-0.10)=0.0900\). Thus the pooled standardizer is

$$ \sqrt{\frac{0.0736+0.0900}{2}} = \sqrt{0.0818} \approx 0.2860. $$

The standardized difference is

$$ \operatorname{SMD}_{X\mid R=1} = \frac{0.08-0.10}{\sqrt{0.0818}} \approx \frac{-0.0200}{0.286007} \approx -0.06993. $$

Among units with observed covariate values, the treated-group proportion is 0.02 lower, or about 0.06993 pooled standard deviations lower. This small observed-value difference does not erase the need to report the missingness SMD: the calculation says nothing about whether the groups differ in \(P(R=1\mid A)\).

Worked Example: A Continuous Covariate With Incomplete Values

Consider hypothetical observed values of a continuous baseline covariate. There are 80 observed treated units with mean 52 and sample standard deviation 10, and 90 observed comparison units with mean 48 and sample standard deviation 12. Calculate the observed-value SMD using the pooled sample standard deviation. The covariate’s missingness rate must be reported separately; it is not provided here.

With \(n_1=80\), \(n_0=90\), \(s_1=10\), and \(s_0=12\), the pooled sample standard deviation is

$$ s_{\mathrm{pooled}} = \sqrt{\frac{(80-1)(10^2)+(90-1)(12^2)}{80+90-2}} = \sqrt{\frac{7900+12816}{168}} = \sqrt{123.3095} \approx 11.1045. $$

The observed-value SMD is

$$ \operatorname{SMD}_{X\mid R=1} = \frac{52-48}{11.1045} \approx 0.3602. $$

Among the observed units, the treated-group mean is 4 units higher, or about 0.3602 pooled standard deviations higher. The covariate’s measurement units are those of \(X\); the SMD itself is unitless. Because the example provides no missing-value counts, it cannot support a conclusion about missingness balance. A complete report would include the missingness rates and their SMD as well.

Read Missingness and Observed Values Together

A useful balance table gives separate rows or clearly labeled columns for the missingness indicator and the observed covariate values. For every covariate, show the treatment-group missingness proportions, their SMD, the observed-value summaries, their SMD, and the observed counts. For a binary variable, report proportions; for a continuous variable, report means and a spread summary. If weights are used, identify them and provide weighted summaries consistent with the analysis.

Missingness diagnosticObserved-value diagnosticInterpretation to consider
Substantial difference in missingnessSmall observed-value SMDRecorded values may look similar, but the groups differ in which units have values recorded.
Small difference in missingnessSubstantial observed-value SMDAvailability is similar, but the observed covariate distributions remain different.
Substantial differences in bothSubstantial differences in bothBoth availability and the observed values warrant attention; neither diagnostic explains the other.
Small differences in bothSmall differences in bothThese diagnostics show less evidence of imbalance on these summaries, but do not establish overall exchangeability.

One way to organize the diagnostic workflow is:

1
Define the analysis cohort and weights.
Use the cohort and treatment coding relevant to the matching or weighting analysis. Keep the target population and weighting scheme explicit.
2
Summarize availability.
For each covariate, calculate the observed and missing counts or proportions in each treatment group, then compare the missingness indicators.
3
Compare recorded values.
Among units with \(R_j=1\), calculate group summaries and an observed-value SMD. Report observed counts so readers can see the comparison’s basis.
4
Interpret the pair of results.
Describe availability imbalance separately from value imbalance. Do not infer a missing-data mechanism or causal validity from either SMD alone.

Diagnostics, Assumptions, and Failure Modes

These comparisons are most interpretable when covariates are baseline variables measured before treatment, have consistent definitions across groups, and are evaluated in the same analytic cohort used for the balance assessment. The comparison also assumes that coding of missing values is correct: a blank, an explicit “unknown” category, and a genuine numerical value must not be accidentally conflated. For weighted analyses, weights should be finite and the weighted summaries should reflect the intended target population. Severe weight concentration or very few observed values can make summaries unstable, so counts and weight concentration deserve attention alongside SMDs.

Observed-value SMDs are conditional comparisons: they compare \(X_j\) among units with \(R_j=1\), not among all units. If observation depends on treatment, on the covariate itself, or on other characteristics, the observed subset can differ from the full group. Thus a small observed-value SMD does not show that the unobserved values would have been balanced. Likewise, similar missingness proportions do not show that the missing values are missing completely at random or missing at random. Balance summaries alone cannot distinguish those mechanisms.

A useful diagnostic review asks whether missingness is concentrated in particular treatment groups, whether only a small number of observed units contribute to the value comparison, and whether the observed-value summaries change markedly under the analysis weights. If missingness patterns involve several covariates, it can also be informative to tabulate common patterns of observed and missing values by treatment group. Such summaries can reveal a problem that separate marginal indicators obscure, although sparse patterns may be difficult to interpret.

Avoid coding a missing continuous value as zero merely to calculate a mean: that changes the covariate’s meaning and can create artificial imbalance. Also avoid reporting only a complete-case balance table. Within a complete-case subset, every included unit has \(R_j=1\), so the missingness indicator is constant and its SMD is undefined; this restriction has discarded precisely the information needed to describe differential availability in the original cohort. A complete-case comparison may still be a chosen analysis, but it does not replace a report of missingness in the cohort before that restriction.

How These Diagnostics Relate to Other Approaches

Missingness balance is a diagnostic, not a missing-data remedy. Complete-case analysis restricts the analysis to units with the required measurements; it can change the target population and does not by itself correct selection associated with missingness. Imputation supplies values under a specified model and assumptions, while a missing-category approach treats missingness as a category for some modeling purposes. Those choices have different implications and should not be selected solely because one balance table looks better. The present diagnostics help describe the data relevant to those choices; they do not establish which method is valid.

Balance assessment also differs from a significance test. An SMD describes the size of a group difference on a standardized scale and does not become evidence of equivalence simply because a test is not statistically significant. Nor does a small SMD guarantee balance in every feature of a covariate distribution. The wider set of balance tools, including distributional and nonlinear checks, remains relevant where appropriate. The role of the missingness diagnostics is to ensure that treatment-group differences in data availability are not hidden by summaries of observed values alone.

Key takeaway: Report two distinct comparisons for an incomplete covariate: balance in its missingness indicator across the full analysis cohort, and balance in its values among units where those values are observed. Interpret each with its population, counts, weighting, and limitations stated.

Check Your Understanding

Answer each question, then compare your reasoning with the solutions below.

  1. In a hypothetical cohort, 30 of 120 treated units and 12 of 120 comparison units have a covariate missing. What are the missingness proportions and the binary missingness SMD?
  2. If the observed-value SMD is close to zero but the missingness SMD is large, what can and cannot be concluded?
  3. Why does a complete-case balance table not diagnose missingness balance in the original cohort?
  4. Do similar missingness rates establish that missingness is completely at random? Explain.

Solutions. (1) The missingness proportions are \(30/120=0.25\) and \(12/120=0.10\). The pooled standardizer is \(\sqrt{\{0.25(0.75)+0.10(0.90)\}/2}=\sqrt{(0.1875+0.0900)/2}=\sqrt{0.13875}\approx0.37249\). Therefore the SMD is \((0.25-0.10)/0.37249\approx0.4027\). The treated group has the higher missingness proportion. (2) The recorded values appear similar by that summary among observed units, but availability differs; nothing in the two SMDs reveals the unobserved values or establishes a missing-data mechanism. (3) In the complete-case subset the indicator is always \(R=1\), so it has no variation and its SMD cannot be calculated. The restriction has removed the missingness information that must be summarized in the original cohort. (4) No. Equal rates describe only the marginal proportions missing in the two treatment groups. They do not show whether missingness depends on measured or unmeasured characteristics, or on the unobserved values themselves.

Clinical Biostats