Tutorials › Pattern-Mixture and Selection Models › Comparing Dropout Models With Different Outcome Histories

Comparing Alternative MNAR Models · Tutorial 178 of 499

Comparing Dropout Models With Different Outcome Histories

Learn to compare outcome-history specifications in dropout models while separating dependence on observed history from dependence on an unobserved outcome.

Advanced 10 min read

What You'll Learn

  • Define latest-outcome, full-history, and summary-history dropout predictors.
  • Distinguish dependence on observed past outcomes from dependence on a missing current outcome.
  • Calculate and compare dropout probabilities under alternative history specifications.
  • Identify when a summary is an exact restriction of a full-history model.
  • Check history-model assumptions and recognize their main diagnostic and interpretation limits.

Why the Outcome History Matters

A dropout model describes how the chance of leaving a longitudinal study changes with information available at a visit. Even when candidate models use the same link and include the same current outcome dependence, they can give different dropout probabilities because they represent past outcomes differently. One model may use only the most recent observed outcome; another may allow each past outcome to contribute separately; a third may use a summary such as the past mean.

These choices are consequential: a model that retains only the latest value cannot distinguish two participants with the same latest value but very different earlier trajectories. A full-history model can represent that distinction, but estimates more parameters and may be unstable when histories are long or sparse. A summary is more compact, but it treats histories as equivalent whenever they have the same summary.

Definition: For participant \(i\) at visit \(t\), let \(D_{it}=1\) indicate dropout before the outcome at visit \(t\) is observed, among participants still under follow-up just before that visit. Let \(Y_{it}\) be the outcome at visit \(t\), \(X_i\) denote baseline covariates, and \(H_{i,t-1}=(Y_{i1},\ldots,Y_{i,t-1})\) denote the observed outcome history before visit \(t\). A latest-outcome model uses \(Y_{i,t-1}\); a full-history model uses the components of \(H_{i,t-1}\); a summary-history model uses a specified function \(S(H_{i,t-1})\). A selection model may also include the current outcome \(Y_{it}\), which is unobserved for someone who drops out before that outcome is recorded.

For clarity, consider logistic dropout models with a common visit-specific intercept and baseline covariate component. In all three specifications below, the coefficient \(\delta\) on \(Y_{it}\) is included as the same form of direct outcome dependence. The models differ in how they use the prior observed history:

$$ \begin{aligned} \operatorname{logit}P(D_{it}=1\mid D_{i,t-1}=0,X_i,H_{i,t-1},Y_{it}) &=\alpha_t+\gamma^\mathsf{T}X_i+\delta Y_{it}+\beta_LY_{i,t-1},\\ \operatorname{logit}P(D_{it}=1\mid D_{i,t-1}=0,X_i,H_{i,t-1},Y_{it}) &=\alpha_t+\gamma^\mathsf{T}X_i+\delta Y_{it}+\sum_{j=1}^{t-1}\beta_jY_{ij},\\ \operatorname{logit}P(D_{it}=1\mid D_{i,t-1}=0,X_i,H_{i,t-1},Y_{it}) &=\alpha_t+\gamma^\mathsf{T}X_i+\delta Y_{it}+\theta S(H_{i,t-1}). \end{aligned} $$

Here \(\operatorname{logit}(p)=\log\{p/(1-p)\}\), so a linear predictor \(\eta\) gives dropout probability \(1/\{1+\exp(-\eta)\}\). The dependence on \(Y_{it}\) makes the selection model MNAR when that outcome is missing at dropout: the dropout probability then depends on an unobserved value. Dependence only on observed past outcomes does not, by itself, create that feature. This distinction matters even though all three candidates are being compared as MNAR models: the common current-outcome term supplies their direct MNAR dependence, while their history terms encode different observed-history relationships.

What the Three History Specifications Assume

The latest-outcome specification assumes that, once \(Y_{i,t-1}\), \(X_i\), and the current outcome term are accounted for, earlier outcomes do not add further information about dropout. This can be a useful working approximation when recent status is expected to dominate, but it is not a consequence of having longitudinal data.

The full-history specification allows each prior visit to have its own coefficient. It can distinguish timing and trajectory shape, but requires more parameters as the number of visits increases. The summary specification imposes a particular compression. For example, if \(S(H_{i,t-1})\) is the unweighted mean of the prior outcomes, every prior visit has the same coefficient per unit of outcome in the linear predictor. Choosing a summary is therefore a substantive modeling restriction, not merely a computational shortcut.

These specifications can be nested. A latest-outcome model is a restriction of the full-history model when all coefficients on earlier outcomes are zero. A summary model is a restriction of the full-history model when the full-history coefficient vector has the specified weights of that summary. This gives a precise way to understand what information a compressed history retains.

$$ \sum_{j=1}^{t-1}\beta_jY_{ij} =\theta\left(\sum_{j=1}^{t-1}w_jY_{ij}\right) \quad\text{if}\quad \beta_j=\theta w_j\ \text{for every }j. $$

Thus, for the weighted summary \(S(H)=\sum_jw_jY_{ij}\), the summary and full-history linear predictors agree exactly under those coefficient restrictions. If the restrictions fail, the summary cannot reproduce the full-history predictor for every possible history. This algebraic fact does not establish that either restriction is scientifically appropriate.

Worked Example: Calculating Dropout Probabilities

Suppose a hypothetical participant has prior outcomes \(Y_{i1}=2\), \(Y_{i2}=4\), and \(Y_{i3}=5\), and current outcome \(Y_{i4}=6\). For illustration, use the following candidate linear predictors, with no baseline covariates:

Latest outcome: \(\eta_L=-1+0.5Y_{i4}+0.3Y_{i3}\). Full history: \(\eta_F=-1+0.5Y_{i4}+0.1Y_{i1}+0.2Y_{i2}+0.3Y_{i3}\). Mean summary: \(\eta_S=-1+0.5Y_{i4}+0.3\{(Y_{i1}+Y_{i2}+Y_{i3})/3\}\). These hypothetical coefficients are illustrative, not estimates from data.

For the latest-outcome model, \(\eta_L=-1+0.5(6)+0.3(5)=-1+3+1.5=3.5\). Therefore,

$$ P_L(D_{i4}=1)=\frac{1}{1+\exp(-3.5)}\approx0.9707. $$

For the full-history model, \(\eta_F=-1+0.5(6)+0.1(2)+0.2(4)+0.3(5)=-1+3+0.2+0.8+1.5=4.5\). Hence,

$$ P_F(D_{i4}=1)=\frac{1}{1+\exp(-4.5)}\approx0.9890. $$

The prior mean is \((2+4+5)/3=11/3\approx3.667\). Thus \(\eta_S=-1+3+0.3(3.667)\approx3.100\), and \(P_S(D_{i4}=1)=1/\{1+\exp(-3.100)\}\approx0.9569\). The three predictions differ because the history terms differ. The result is not evidence that one model is better: model coefficients would ordinarily be estimated under their respective specifications, and this arithmetic only illustrates how each specification uses a given history.

When Summaries Hide Trajectory Differences

A direct way to see what is lost by compression is to compare histories with the same summary. Suppose the two hypothetical histories are \((2,4,5)\) and \((4,2,5)\). Both have mean \(11/3\) and latest value 5. Consequently, the latest-outcome and mean-summary predictors in the preceding example assign the same dropout probabilities to both histories, provided the current outcome and other covariates are also the same.

The full-history predictor can distinguish them because the coefficients vary by visit. For \((2,4,5)\), the history contribution is \(0.1(2)+0.2(4)+0.3(5)=2.5\). For \((4,2,5)\), it is \(0.1(4)+0.2(2)+0.3(5)=2.3\). With the shared contribution \(-1+0.5(6)=2\), the respective linear predictors are \(4.5\) and \(4.3\). The corresponding probabilities are approximately \(0.9890\) and \(1/\{1+\exp(-4.3)\}\approx0.9866\). The latest and mean models cannot express this difference, because their inputs are identical for these two histories.

Worked Example: Checking Whether a Summary Can Reproduce a Full-History Model

At a visit with three observed past outcomes, suppose the full-history contribution is \(0.2Y_{i1}+0.2Y_{i2}+0.2Y_{i3}\). Can a model using the prior mean reproduce this contribution for every history?

Let \(S=(Y_{i1}+Y_{i2}+Y_{i3})/3\). Then the full-history contribution equals \(0.2(Y_{i1}+Y_{i2}+Y_{i3})=0.6S\). A mean-summary term with coefficient \(\theta=0.6\) therefore reproduces the full-history contribution exactly for every possible set of three outcomes.

In contrast, coefficients \(0.1, 0.2, 0.3\) cannot be represented as a coefficient times the unweighted mean: that would require all three full-history coefficients to be equal. Histories with the same mean but different ordering can therefore yield different full-history predictors, as in the preceding calculation. The exact equivalence depends on the specified summary weights and coefficient restrictions, not on the general idea of summarizing.

Assumptions, Comparison, and Failure Modes

A meaningful comparison requires aligned models. Use the same definition of dropout, outcome scale, link, baseline covariates, current-outcome term, visit structure, and analysis population. Otherwise, a change in predictions cannot be attributed specifically to the history specification. The history must also be available as defined for every participant at risk; when prior visits themselves are missing, the rule for forming or modeling the history must be explicit. A model using a mean of observed prior values, for example, is not automatically the same as one using a mean over all scheduled prior visits.

The logistic form assumes a linear relationship between each included predictor and the log odds of dropout, conditional on the other terms. It also assumes the chosen history representation is adequate for the intended comparison. Inference further relies on a correctly specified likelihood and appropriate treatment of repeated observations within participants. Predictions should not be extrapolated to covariate or outcome histories unsupported by the data; estimates of current-outcome dependence are especially consequential because that outcome is missing for dropouts.

Model comparisons should therefore examine what the predictor structure can and cannot express, not just compare a single predicted probability. The following table summarizes the central trade-offs.

SpecificationHistory information usedMain restriction or cost
Latest outcomeMost recent observed outcomeEarlier outcomes add no contribution once the latest value is included
Full historySeparate terms for prior outcomesFlexible, but uses more coefficients and may be imprecise with limited histories
Summary historyA defined function, such as a meanHistories with the same summary are treated alike by that history term

Warning signs include implausibly large or unstable history coefficients, fitted probabilities near 0 or 1 for many participants, and substantial changes in coefficients when a small number of history terms are added or removed. These symptoms can reflect sparse histories, strong correlation among repeated outcomes, or an overly ambitious full-history model. A summary model can also fail in a less conspicuous way: its fit may appear stable while its imposed equivalence of distinct trajectories is scientifically implausible.

Observed-data fit alone cannot settle how dropout depends on the missing current outcome. Different assumptions about that unobserved value can be compatible with the same observed information. The history comparison is useful for exposing one part of the modeling choice, but it does not validate the unobserved-outcome dependence. Interpret conclusions as conditional on the specified selection model, and distinguish observed-history flexibility from assumptions about outcomes after dropout.

Common Mistake

Common mistake: “The dropout model uses prior outcomes, so it is MNAR.” That conclusion does not follow. If dropout depends only on outcomes already observed (and included covariates), the dependence is on observed information. MNAR dependence arises here from including the current outcome when it is unobserved for participants who drop out. Conversely, calling a model “full history” does not mean it uses all information about a participant: it uses only the specified observed outcomes and functional form.

A second pitfall is treating the latest-outcome model as simply a smaller version that must be adequate. It is a restriction: earlier outcomes have zero additional contribution conditional on the latest outcome and other terms. Likewise, a summary is not neutral compression; it makes histories sharing that summary equivalent in the dropout predictor. State these restrictions before interpreting differences among candidate models.

Check Your Understanding

Answer each question, then check the brief solutions.

  1. If dropout at visit \(t\) depends on \(Y_{i,t-1}\) but not on \(Y_{it}\), is that dependence alone on an observed or unobserved outcome?
  2. A full-history predictor has coefficients \(0.4,0.4,0.4\) on three prior outcomes. Express its history contribution using the prior mean.
  3. Two participants have the same latest outcome and the same prior mean, but different outcome ordering. Which of the three history specifications can distinguish them without additional terms?
  4. Give one reason a full-history model may be unstable when there are many visits.
  5. Why must candidate models share their current-outcome term if the aim is to isolate differences in observed-history specification?

Solutions. (1) It is dependence on an observed prior outcome, not by itself dependence on an unobserved current outcome. (2) The contribution is \(0.4(Y_{i1}+Y_{i2}+Y_{i3})=1.2S\), where \(S\) is the three-outcome mean. (3) The full-history specification can distinguish them if its visit-specific coefficients differ; latest-outcome and mean-summary specifications cannot. (4) It estimates many potentially correlated coefficients, which can produce imprecision or instability with limited data. (5) Otherwise prediction differences could arise from changed MNAR assumptions about the missing current outcome, rather than from the history representation being compared.

Clinical Biostats