Why Missing-Data Assumptions Matter for Exposure-Response Decisions
An exposure-response analysis uses measured or modeled exposure to predict a clinical outcome. If some outcomes are missing because participants discontinue, miss visits, or become too unwell to provide measurements, the observed participants may not represent everyone whose outcomes matter to the decision. A prediction based only on observed outcomes can therefore depend on an implicit assumption about what would have happened after dropout.
This tutorial focuses on sensitivity of decision-relevant exposure-response predictions to those assumptions. The central comparison is between missing at random (MAR), where missingness can be explained by information already observed, and informative dropout beyond observed information, where missingness also depends on unobserved outcomes. The latter is often described as missing not at random (MNAR). Neither label is a diagnosis that can usually be established from the observed data alone.
A decision target must be stated before comparing assumptions. For example, it may be the mean week-12 response predicted at a specified exposure in a target population, or the probability that this mean exceeds a development threshold. An exposure-response prediction describes an association unless the design and analysis also support a causal interpretation. Missing-data sensitivity analysis addresses uncertainty about unobserved outcomes; it does not by itself resolve confounding, exposure measurement error, or misspecification of the exposure-response relationship.
What MAR Allows—and What It Does Not
Under MAR, participants with the same observed exposure, covariates, and response history provide information about outcomes for participants who dropped out with that same observed history. A likelihood analysis can use this information without explicitly modeling the dropout process, provided the outcome model is adequate and the usual ignorability conditions hold. In particular, the parameters governing the outcome model and missingness mechanism must be distinct.
The factorization explains the point. Write \(Y_{\mathrm{obs}}\) and \(Y_{\mathrm{mis}}\) for observed and missing outcomes. Under MAR, the probability of the observed missingness pattern, conditional on outcomes and covariates, depends on observed outcomes but not on \(Y_{\mathrm{mis}}\). If the outcome parameter \(\theta\) and missingness parameter \(\psi\) are distinct, the joint likelihood can be written as:
The second factor contains no \(\theta\), so it does not change likelihood-based estimation of \(\theta\). This is why an ignorable likelihood analysis need not fit a dropout model under MAR and parameter distinctness. This does not mean the outcome model can be careless: it must represent the relevant exposure, covariate, and observed response-history relationships. Any missingness model used by a chosen method, such as inverse-probability weighting, must also be adequate. The covariates needed for prediction must be observed for the participants to whom the analysis applies, and there must be sufficient overlap in observed exposure and history patterns.
MAR is an assumption about dependence, not a claim that dropout is unrelated to prognosis. If participants with poor observed responses are more likely to discontinue, that can still be MAR when the analysis conditions appropriately on those observed responses. In contrast, if participants discontinue because of an unobserved worsening that is not captured by their recorded history, MAR may fail.
Worked Example: Predicting Missing Outcomes From Observed History
Hypothetical data: At a specified exposure, 40 participants are intended to contribute a week-12 response measured in response units. Thirty have an observed week-12 response, with mean 14 units. The other 10 dropped out. Their recorded week-4 responses place six in a history group for which the fitted MAR outcome model predicts a week-12 mean of 13 units, and four in a group for which it predicts 5 units. Assume the outcome model is appropriate for these histories and exposure.
Under this MAR model, the predicted total response among the 10 missing participants is \(6(13)+4(5)=98\) response-unit equivalents, so their predicted mean is \(98/10=9.8\) units. The predicted mean for all 40 participants is:
The complete-case mean is 14 units, whereas the MAR-based all-participant prediction is 12.95 units. The difference arises because the missing participants' observed histories lead the model to predict lower week-12 responses. This calculation illustrates the assumption: within the relevant observed histories, missing participants are assigned the outcome distribution represented by the fitted model. It does not prove that their unobserved outcomes actually follow that distribution.
Representing Informative Dropout With Sensitivity Parameters
When outcomes after dropout may be worse or better than MAR predicts, a sensitivity analysis can shift the predicted missing-outcome distribution. A simple pattern-mixture specification adds a parameter \(\delta\) to the MAR-predicted mean for missing outcomes. In a group with a common shift:
Here \(\delta\) is measured in the outcome's units. A negative value represents outcomes lower than the MAR prediction; a positive value represents higher outcomes. The data generally cannot estimate \(\delta\) from the observed outcomes alone, so analysts examine a range justified by clinical knowledge, external information, or explicit decision stress tests. The pattern-mixture approach makes the departure transparent, but a constant shift may be too simple if the departure varies by exposure, response history, or time.
A tipping point is the value of a sensitivity parameter at which a specified decision changes. It is a useful summary of decision robustness, not evidence that the tipping-point value is the true missing-data mechanism.
Worked Example: Finding the Tipping Point for a Response Threshold
Hypothetical data: At a target exposure, 60 of 80 participants have observed week-12 responses averaging 12 units. For the 20 participants with missing outcomes, the MAR model predicts a mean of 8 units. The development criterion is a mean response of at least 10.5 units. Let \(\delta\) be a common shift, in response units, applied to the missing-outcome predictions.
The all-participant mean as a function of the shift is:
At MAR, \(\delta=0\), the predicted mean is 11 units, which exceeds the 10.5-unit criterion. To find the tipping point, set the predicted mean equal to the criterion:
Thus, a shift of \(-2\) units brings the prediction exactly to the threshold; a more negative shift puts it below the threshold. For example, at \(\delta=-3\), the missing-outcome mean is \(8-3=5\) units and the all-participant mean is \((720+20(5))/80=820/80=10.25\) units. The conclusion changes under that scenario. The tipping point does not tell us whether a \(-2\)-unit shift is plausible; that judgment requires substantive justification.
Comparing Assumptions and Checking the Analysis
A useful analysis reports the decision prediction under a primary assumption and under alternatives that represent plausible departures. For example, a likelihood analysis under MAR can be compared with one or more MNAR pattern-mixture shifts. A selection model instead specifies dropout probability conditional on observed information and the potentially unobserved outcome. These approaches parameterize the missingness assumption differently; neither makes MNAR departures identifiable merely by fitting the observed data. Complete-case analysis is a simpler alternative, but it estimates the target mean without bias only under additional conditions, such as missingness being independent of the outcome within the analyzed groups. It is not a general substitute for a missing-data assumption.
Diagnostics can reveal weaknesses in the observed-data analysis, but cannot verify MAR against MNAR using observed outcomes alone. Check whether dropout rates vary with exposure, baseline covariates, and recorded response history; whether predicted outcomes are supported by participants with similar observed histories; and whether model fit, convergence, and uncertainty are acceptable. Sparse histories or exposure ranges with few observed outcomes signal extrapolation. A highly influential dropout pattern or unstable predictions across reasonable outcome models also warrants caution.
For the primary outcome model, examine whether its exposure-response form and relevant covariate and history effects describe the observed data adequately. For methods that explicitly model dropout, examine that missingness model as well. Do not require a dropout model solely because the likelihood analysis assumes MAR: under MAR with parameter distinctness, an adequate likelihood for the outcome can be ignorable without one. Equally, do not mistake adequate fit among observed outcomes for evidence that the unobserved outcomes follow MAR predictions.
Specify the population, exposure range, outcome time, summary, and decision threshold before comparing assumptions.
State the observed covariates and response history used, evaluate model adequacy, and identify where predictions rely on sparse data or extrapolation.
Explain why MAR is a reasonable primary analysis, and specify any alternative shifts or dropout models and their clinical interpretation.
Report the estimates and uncertainty under each scenario, identify any decision change, and describe the assumptions at the tipping point.
The practical limitation is fundamental: observed data alone typically cannot distinguish MAR from MNAR. Sensitivity results are therefore conditional statements: if the missing outcomes differ from MAR predictions by the stated amount, the decision prediction is as reported. They should be presented with the assumptions, units, and uncertainty, rather than as a single assumption-free answer.
Worked Example: Different Dropout Departures Can Change an Exposure-Response Contrast
Hypothetical data: Consider two exposure groups, with 50 participants in each. In the lower-exposure group, 40 week-12 outcomes are observed with mean 10 response units; MAR predicts a mean of 9 units for the 10 missing outcomes. In the higher-exposure group, 30 outcomes are observed with mean 14 units; MAR predicts a mean of 12 units for the 20 missing outcomes. Treat exposure as a two-level grouping for this illustration.
Under MAR, the predicted means are \( (40(10)+10(9))/50=490/50=9.8\) units at lower exposure and \( (30(14)+20(12))/50=660/50=13.2\) units at higher exposure. The higher-minus-lower contrast is \(13.2-9.8=3.4\) response units.
Now consider an informative-dropout scenario in which missing outcomes are 1 unit below MAR at lower exposure and 3 units below MAR at higher exposure. The adjusted missing-outcome means are 8 and 9 units, respectively. The resulting group means are:
The contrast is now \(12.0-9.6=2.4\) response units, rather than 3.4 units under MAR. If the decision requires a higher-exposure mean of at least 12.5 units, the MAR prediction meets the criterion and this MNAR scenario does not. The calculation shows why the departure may need to vary by exposure: a common shift would not represent this particular scenario. The two-level contrast is illustrative and is not, by itself, a continuous exposure-response model.
Common Pitfall
Wrong reasoning: “Dropout is related to prognosis, so a missingness model must be fitted for a valid MAR likelihood analysis.” Prognosis can affect dropout through observed response history, while MAR still holds conditional on that history. Under MAR and parameter distinctness, the likelihood for the outcome parameters can be ignorable without modeling dropout.
Correction: Adequately model the outcome relationships relevant to the target prediction, and ensure the required covariates and histories are observed with sufficient overlap. Assess a missingness model when the chosen method uses one. Then evaluate how the decision changes under explicit departures from MAR, while recognizing that observed-data diagnostics cannot establish which missingness assumption is true.
Check Your Understanding
Answer each question, then compare your reasoning with the solutions below.
- In a longitudinal analysis, dropout depends on a participant's recorded week-4 response but not on later unobserved responses after conditioning on week-4 response and covariates. Is this compatible with MAR?
- Why can a likelihood analysis under MAR omit an explicit dropout model?
- If a missing-outcome mean is predicted as 15 mg/L under MAR and the sensitivity shift is \(-2\) mg/L, what mean is used for that missing-outcome pattern?
- Can a good fit to observed outcomes prove that the data are MAR? Explain briefly.
- In a decision analysis, the MAR prediction exceeds the target but a clinically plausible MNAR scenario falls below it. What should be communicated?
Solutions: (1) Yes. Dropout may depend on observed history under MAR, provided it does not additionally depend on the unobserved outcomes conditional on that history and other observed data. (2) Under MAR and distinct outcome and missingness parameters, the observed-data likelihood factors into an outcome likelihood and a missingness factor that contains no outcome parameter; the latter can be omitted for estimating the outcome parameters. (3) \(15-2=13\) mg/L. (4) No. Fit to observed outcomes assesses the specified outcome model among observed data, but cannot reveal whether missing outcomes differ from its MAR predictions. (5) Report both predictions, the assumptions and units for each, the threshold, and the fact that the decision depends on the missing-data assumption; explain the plausibility basis for the alternative scenario.