Why Locate Disagreement?
A single summary of disagreement across candidate missing-not-at-random (MNAR) models can hide where their assumptions matter. Models may give similar predictions early in follow-up but diverge later, or agree overall while producing quite different predictions for one prespecified subgroup. For a treatment comparison, the arms can each show substantial model-specific movement while the treatment contrast remains comparatively stable. Locating these patterns helps focus sensitivity interpretation on the follow-up times and populations where conclusions are most affected.
The task here is not to decide which candidate is best. It is to compare their predictions for the same estimand at specified times and in specified subgroups. Earlier material on model envelopes for sensitivity conclusions describes an overall summary; the time- and subgroup-specific summaries below show where that spread occurs.
Make the Comparisons Like-for-Like
These summaries are meaningful only if the candidate models are being asked the same question. Before calculating them, specify the outcome and its units, the follow-up times, treatment coding, subgroup definitions, and target population. Each \(\hat{\mu}_{k,a,g,t}\) must refer to that same combination. If models use different covariate distributions for standardization, for example, their predictions can differ because of both the MNAR assumptions and the target populations. Use common standardization when the goal is to isolate model sensitivity.
The subgroup should be defined before inspecting where models disagree, and the comparison should use the same subgroup definition in every model. A subgroup-specific mean can be estimated directly for that subgroup or standardized to a specified covariate distribution within it; state which interpretation is intended. Subgroups with little observed support can yield unstable estimates and extrapolations. A large spread in such a subgroup is a reason for caution, not by itself evidence that the subgroup truly behaves differently.
The models should also use a common outcome scale and comparable treatment estimands. For instance, a difference in outcome means is in outcome units, whereas a ratio is unitless; their spreads cannot be compared directly without a justified common scale. If a contrast is interpreted causally, its causal interpretation additionally requires the design and identification assumptions appropriate to the study, such as valid treatment assignment in a randomized trial. The calculation of model spread alone does not establish causality.
Prediction Spread and Treatment-Contrast Spread
At each time and subgroup, first inspect treatment-specific prediction spreads \(S_{g,t,a}\). Then inspect the contrast spread \(C_{g,t}\) if the scientific question concerns the difference between treatments. These are different summaries. Prediction spread describes how much the candidate models move a treatment-specific mean; contrast spread describes how much they move the difference between two means.
To justify the bound, write \(x_k=\hat{\mu}_{k,1,g,t}\) and \(y_k=\hat{\mu}_{k,0,g,t}\). For any two candidate models \(i,j\), $$ |(x_i-y_i)-(x_j-y_j)|=|(x_i-x_j)-(y_i-y_j)|\leq |x_i-x_j|+|y_i-y_j|. $$ The first absolute difference on the right is no greater than \(\max_k x_k-\min_k x_k=S_{g,t,1}\); the second is no greater than \(S_{g,t,0}\). Taking the largest left-hand difference over all \(i,j\) gives \(C_{g,t}\leq S_{g,t,1}+S_{g,t,0}\). The inequality explains why arm-specific results should accompany a treatment-contrast summary when interpretation depends on both.
Worked Example: Finding a Time-Specific Prediction Spread
Suppose three hypothetical MNAR models predict a mean pain score at week 8 of 12.4, 10.9, and 8.8 points for the same treatment and prespecified subgroup. Assume all three predictions refer to the same score scale, subgroup, and target population. Suppose, for this illustration, that a 2-point difference is the prespecified threshold for practical importance.
The model-specific prediction spread is $$ S_{g,8,a}=\max(12.4,10.9,8.8)-\min(12.4,10.9,8.8)=12.4-8.8=3.6\text{ points}. $$
The spread exceeds the stated 2-point threshold by \(3.6-2=1.6\) points. Thus the candidate MNAR assumptions have practically important consequences for this week-8 mean under the chosen threshold. This calculation does not establish which prediction is correct, and it does not measure sampling uncertainty.
Worked Example: Arm Predictions Can Disagree More Than Their Contrast
At week 12 in a hypothetical subgroup, three candidate models give the following mean outcomes, in points:
| Model | Treatment 1 | Treatment 0 | Contrast, treatment 1 minus treatment 0 |
|---|---|---|---|
| 1 | 18 | 14 | 4 |
| 2 | 16 | 13 | 3 |
| 3 | 14 | 12 | 2 |
For treatment 1, the prediction spread is \(18-14=4\) points. For treatment 0, it is \(14-12=2\) points. The model-specific contrasts are calculated as treatment 1 minus treatment 0: \(18-14=4\), \(16-13=3\), and \(14-12=2\) points. Therefore, $$ C_{g,12}=\max(4,3,2)-\min(4,3,2)=4-2=2\text{ points}. $$
Both arms' means move down across these models, but they move in a partly similar way. Consequently, the contrast spread is 2 points, smaller than the 4-point spread for treatment 1. If the scientific target is the treatment contrast, the 2-point spread is directly relevant; the arm-specific spreads still matter for describing the predicted outcomes and checking for cancellation.
Locate the Largest Practical Differences
A useful analysis lays out the summaries across the prespecified times and subgroups, rather than reporting only the largest number. A table or heat map can display \(C_{g,t}\), with arm-specific \(S_{g,t,a}\) available alongside it. Interpret the largest values against a practical threshold on the relevant scale, not merely by their rank. If outcomes use different scales across analyses, raw spreads cannot be compared as though they had a common meaning.
The following hypothetical comparison uses a prespecified 2.5-point threshold for practically important disagreement in the treatment contrast. The values are model-specific contrasts in outcome points; the spread in each cell is the maximum minus the minimum across the three candidate models.
Worked Example: Identifying a Time-by-Subgroup Hotspot
Suppose investigators prespecify two subgroups and two follow-up times. The three hypothetical model-specific treatment contrasts are:
| Subgroup | Time | Contrasts across models (points) | Contrast spread (points) |
|---|---|---|---|
| Younger participants | Week 4 | 3, 2.5, 2 | 1 |
| Younger participants | Week 12 | 4, 3, 1 | 3 |
| Older participants | Week 4 | 2, 1.8, 1.5 | 0.5 |
| Older participants | Week 12 | 3, 0, -2 | 5 |
For younger participants at week 12, \(C=4-1=3\) points. For older participants at week 12, \(C=3-(-2)=5\) points. The week-4 spreads are \(3-2=1\) point and \(2-1.5=0.5\) points, respectively. Thus the largest spread is 5 points for older participants at week 12. It exceeds the stipulated 2.5-point threshold, while the younger-participant week-4 spread does not.
This pattern identifies where model assumptions have the largest apparent practical impact among these prespecified comparisons. It does not show that the treatment effect truly differs by age or that the older-participant estimate is reliable: subgroup support, uncertainty, and the estimand still need scrutiny.
Uncertainty, Diagnostics, and Limitations
Across-model spread and sampling uncertainty answer different questions. \(C_{g,t}\) describes variation in point estimates across the chosen MNAR assumptions. A model-specific standard error or interval describes uncertainty conditional on one model and its assumptions. Neither substitutes for the other. A small across-model spread does not guarantee precise estimates, and a large spread does not quantify the probability that any model is correct.
Where uncertainty is important, show model-specific intervals or uncertainty summaries alongside the cross-model spread. If estimates are obtained by resampling, refit each candidate model within each resample and calculate the same time- and subgroup-specific summaries; this reflects estimation variability under the specified resampling procedure. For Bayesian fits, model-specific posterior summaries can be shown, but they remain conditional on each model and its prior. Do not treat the range of point estimates as a confidence interval or as a posterior interval.
Several diagnostics help reveal when a comparison may be misleading:
- Check the estimand. Confirm that each model predicts the same outcome, time, treatment contrast, subgroup, and target population.
- Inspect support. Small subgroups, few observed outcomes at late visits, or sparse treatment-by-subgroup cells can make estimates unstable and MNAR extrapolation especially influential.
- Review arm-specific predictions. Large shifts in both arms may cancel in the contrast; cancellation should be visible rather than mistaken for model agreement everywhere.
- Check sensitivity to the candidate set. The range only summarizes models included in the comparison. Adding or removing a plausible candidate can change its endpoints.
- Keep scales and thresholds explicit. A numerical spread has no automatic clinical meaning. State the outcome units and the basis for any practical threshold.
- Avoid post hoc searching. Choosing subgroups or times because they show the largest disagreement exaggerates how confirmatory that finding is. Label exploratory searches as exploratory and account for the multiplicity of comparisons in interpretation.
These summaries complement, rather than replace, checks of fit to observed data and reasoned assessment of the extrapolation assumptions. The observed data generally cannot determine the unobserved outcomes after dropout without additional assumptions. Therefore, a time-specific disagreement pattern is evidence of sensitivity to the candidate assumptions, not a test that identifies the true missingness mechanism.
Common Mistake and Practical Alternatives
A common mistake is to see a small contrast spread and conclude that all candidate models make similar predictions. That conclusion can be false: treatment-specific means may each have large spreads that cancel when subtracted. The correction is to report the contrast spread for the treatment question and the arm-specific prediction spreads for context. Conversely, large arm-specific spreads do not automatically imply a large treatment-contrast spread.
A single overall envelope is concise but can conceal a localized time or subgroup sensitivity. A time- and subgroup-specific display is more informative when those dimensions were scientifically specified, though it can become noisy when many sparse cells are examined. Model-fit criteria or observed-data diagnostics address other aspects of candidate models; they do not replace comparing the estimand-specific predictions. Model averaging can produce a combined summary under its own weighting choices, but it may obscure the extent and location of disagreement among assumptions. The appropriate display depends on the question, and no spread measure ranks candidates by truth.
Check Your Understanding
Answer each question using the definitions above. Brief solutions follow the questions.
- Three models predict means of 7.2, 6.1, and 5.8 units for the same outcome and subgroup. What is the prediction spread?
- Two treatment contrasts from three models are 1, 0, and -2 units. Calculate their spread and state what it summarizes.
- Can each treatment arm have a prediction spread of 4 units while the treatment-contrast spread is 0? Explain.
- Why is it problematic to define subgroups only after finding the time and subgroup with the largest model spread?
Solutions. (1) The spread is \(7.2-5.8=1.4\) units. (2) It is \(1-(-2)=3\) units; it summarizes the range of model-specific treatment contrasts, not sampling uncertainty. (3) Yes. If every model shifts both arm means by the same amount, the arm-specific predictions can vary while their differences stay constant; for example, arm means of \((8,6)\), \((6,4)\), and \((4,2)\) have prediction spreads of 4 in both arms but contrasts equal to 2 in every model. (4) Selecting the subgroup after examining the spreads makes the result exploratory and can exaggerate the apparent importance of the selected cell; prespecification supports a more interpretable comparison.