Tutorials › Marginal Structural Models › Modeling Dropout Driven by Prior Outcomes

Censoring and Informative Dropout · Tutorial 231 of 501

Modeling Dropout Driven by Prior Outcomes

Learn to include prior health outcomes in discrete-time dropout models and understand the conditional-independence assumption that makes censoring adjustment defensible.

Advanced 11 min read

What You'll Learn

  • Define the risk set and retention hazard when dropout may depend on prior outcomes.
  • Specify a discrete-time censoring model using observed lagged health measures.
  • Explain the sequential conditional-independence assumption behind censoring adjustment.
  • Derive why inverse retention probabilities recover full-data history expectations under stated assumptions.
  • Calculate retention probabilities and weights, and identify positivity and timing failures.
  • Distinguish observed-history-dependent dropout from dependence on unobserved outcomes.

Why Prior Outcomes Belong in a Dropout Model

In a longitudinal study, a participant’s decision or ability to remain under observation may depend on how their health has been changing. For example, participants whose symptoms were severe at the previous visit may be more likely to miss the next visit. If the censoring model ignores that observed symptom history, it can assign the same retention probability to people whose follow-up prospects are quite different.

The key modeling question is therefore not simply whether dropout is associated with the outcome. It is whether the variables used to model retention include the outcome information that was observed before each retention decision. A lagged outcome can be a valid predictor even when the outcome at the visit being missed is unavailable.

Definition: Let \(R_t=1\) if a participant remains observed through visit \(t\), and \(R_t=0\) otherwise, with \(R_0=1\). Dropout is monotone: once \(R_t=0\), all later \(R_s=0\). Let \(Y_t\) denote the outcome at visit \(t\) that would be defined in the absence of dropout, and let \(H_{t-1}\) be the history available before the retention decision at \(t\), including baseline covariates and observed outcomes through \(t-1\). For someone with \(R_{t-1}=1\), the conditional retention probability is \(q_t(H_{t-1})=P(R_t=1\mid R_{t-1}=1,H_{t-1})\).

A discrete-time model for \(q_t\) can include the most recent outcome, a summary of earlier outcomes, and other predictors known before the decision. For example, a pooled logistic model might specify \(\operatorname{logit}\{q_t(H_{t-1})\}=\alpha_t+\beta^\top X+\gamma Y_{t-1}\), where \(X\) is a baseline covariate vector and \(\alpha_t\) allows the baseline retention rate to vary by visit. The form should reflect the study’s observation schedule and substantive knowledge; the equation is a model, not a guarantee that the true retention probability has that form.

The timing rule is essential. If the decision to attend visit \(t\) occurs before \(Y_t\) is measured, \(Y_t\) cannot be a predictor in the model for retention at \(t\). It is unavailable for participants who do not attend. By contrast, \(Y_{t-1}\) can be included if it was recorded before that decision. A model that uses information measured after dropout, or treats an unobserved outcome as though it were observed, does not define a usable observed-history censoring model.

Conditional Independence and the Weighting Rationale

The relevant identifying assumption is sequential conditional independence: among people still observed at the start of a visit, retention at that visit is independent of the future full-data outcomes once the measured past is given. More formally, for each \(t\), assume \((Y_t,\ldots,Y_K)\perp R_t\mid R_{t-1}=1,H_{t-1}\). This assumption says that, within levels of the recorded history, those retained at \(t\) represent those who would have been retained or lost with respect to their future outcomes. Including a prior outcome is important when it helps make this conditional statement plausible.

Also require positivity: for histories occurring in the target population, \(q_t(H_{t-1})>0\). To estimate a population quantity using inverse retention probabilities, the relevant conditional probabilities must be correctly specified or consistently estimated, and the target expectations must be finite. Positivity is not secured merely by fitting a logistic model: fitted probabilities close to zero indicate that some histories have little support among those retained.

Theorem: Suppose the follow-up is monotone; the sequential conditional-independence assumption holds at every visit through \(t\); \(q_j(H_{j-1})>0\) on the histories of interest; and \(g(H_t)\) is integrable. Then \(E\left[\dfrac{R_t g(H_t)}{\prod_{j=1}^{t}q_j(H_{j-1})}\right]=E[g(H_t)]\). Here \(H_t\) is the full-data history through visit \(t\), which is observed whenever \(R_t=1\).

Proof. For a participant retained through \(t-1\), condition on the history \(H_{t-1}\). By the definition of \(q_t\) and sequential conditional independence, the conditional expectation of \(R_t g(H_t)\) is \(q_t(H_{t-1})\) times the conditional expectation of \(g(H_t)\). Dividing by \(q_t(H_{t-1})\) therefore removes the selection at visit \(t\). Repeating this argument backward, one visit at a time, cancels each factor \(q_j(H_{j-1})\) in the product and gives \(E[g(H_t)]\). Participants with \(R_t=0\) contribute zero to the weighted expression; histories needed to evaluate its nonzero terms are observed. Integrability and positivity permit the conditional expectations and divisions used in this argument. \(\square\)

For \(g(H_t)=Y_t\), the identity explains how weighting retained participants can recover the full-data mean at visit \(t\). The lagged outcome is not added to redefine that mean. It helps model who remains observed, so that the weighted observed histories can represent the histories that dropout would otherwise remove.

Building and Interpreting the Model

A practical data set for a discrete-time retention model has one row per participant-visit at which the participant is at risk of dropout. The response indicates whether the participant remains observed at the next visit; the predictors are measured before that retention decision. A model may include visit indicators, baseline covariates, the last observed outcome, and a small number of justified functions of the outcome history. The exact specification depends on the measurement schedule and the hypothesized dropout process.

The choice of lagged-outcome terms should reflect how prior health could drive follow-up. If current health status is predictive, the last value may be appropriate. If deterioration or improvement matters, a recent change such as \(Y_{t-1}-Y_{t-2}\) may carry information beyond the level. If the relationship is not plausibly linear on the log-odds scale, categories, splines, or interactions with visit can be considered. These choices should be based on available pre-decision information, not on whichever specification yields the most convenient weights.

Worked Example: A Retention Model With a Lagged Symptom Score

Suppose higher symptom scores predict lower retention. A fitted model is \(\operatorname{logit}\{q_t\}=1.5-0.8Z+0.4X\), where \(Z=1\) indicates a high symptom score at the previous visit and \(X=1\) indicates a baseline risk factor. For \(X=0\) and a low prior score (\(Z=0\)), the linear predictor is \(1.5\), so \(q_t=1/(1+e^{-1.5})\approx0.8176\).

For the same baseline group with a high prior score (\(Z=1\)), the linear predictor is \(1.5-0.8=0.7\), giving \(q_t=1/(1+e^{-0.7})\approx0.6682\). Thus the model predicts lower retention after a high symptom score. For \(X=1,Z=1\), the predictor is \(1.5-0.8+0.4=1.1\), and retention is \(1/(1+e^{-1.1})\approx0.7503\).

These are conditional retention probabilities, not probabilities of having a good health outcome. The lagged score changes the modeled probability of remaining observed. If it was measured before the visit-\(t\) retention decision, it is available for this use even when \(Y_t\) is missing after dropout.

Worked Example: Multiplying Retention Probabilities Across Visits

Consider a participant with a low outcome at each of two observed visits. Suppose the fitted conditional retention probabilities are \(0.8808\) for remaining through visit 1 and \(0.8320\) for remaining through visit 2, given the history available at each decision. The estimated probability of remaining through both visits is the product: \(0.8808\times0.8320=0.7328256\), which rounds to \(0.7328\).

The corresponding unstabilized inverse retention weight through visit 2 is \(1/0.7328256\approx1.3646\). Using the unrounded probabilities from a fitted model gives a product of about \(0.732839\) and a weight of about \(1.3646\) as well. The product uses conditional probabilities at successive decisions; it is not a single visit-2 probability applied twice. This calculation illustrates the weight construction, while the earlier result explains why such weighting can recover full-data history expectations under its assumptions.

Worked Example: Why Conditioning on Prior Health Matters

Let \(H=1\) indicate poor prior health, with \(P(H=1)=P(H=0)=0.5\). Suppose the future outcome \(Y\) has conditional means \(E[Y\mid H=1]=0.8\) and \(E[Y\mid H=0]=0.2\). Retention probabilities are \(q(1)=0.5\) and \(q(0)=0.9\). The full-data mean is \(E[Y]=0.5(0.8)+0.5(0.2)=0.5\).

Among retained participants, the unnormalized masses for \(H=1\) and \(H=0\) are \(0.5(0.5)=0.25\) and \(0.5(0.9)=0.45\). The retained sample’s mean outcome is therefore \(\{0.25(0.8)+0.45(0.2)\}/(0.25+0.45)=0.29/0.70\approx0.4143\), below the full-data mean.

Weight by the inverse probability conditional on \(H\). The weighted numerator is \(E[RY/q(H)]=0.5(0.5)(0.8)/0.5+0.5(0.9)(0.2)/0.9=0.4+0.1=0.5\). The weighted denominator is \(E[R/q(H)]=0.5(0.5)/0.5+0.5(0.9)/0.9=1\). Their ratio is \(0.5\), recovering the full-data mean. The calculation depends on modeling retention using the prior health group; an unconditional retention probability does not remove the differential selection shown here.

Model Checks, Alternatives, and Limits

The model should be checked for problems that undermine useful adjustment. Review whether each predictor was actually measured before the relevant retention decision, whether the risk set and visit indicators are coded correctly, and whether predicted probabilities are plausible across important histories. Very small probabilities imply large inverse weights and weak empirical support for some retained histories. Sparse combinations of visit, prior outcome, and covariates can also make fitted probabilities unstable. These are warnings about the data and model, not proof that a particular assumption has failed or succeeded.

A model using baseline covariates alone is simpler, but it can leave dropout associated with prior health within baseline strata. A richer observed-history model can better represent that association, though unnecessary terms can increase variance and estimation instability. Outcome regression or imputation are alternative ways to address missing outcomes, but they require their own assumptions and model choices. They do not make the censoring model’s conditional-independence assumption automatically true.

No observed-data diagnostic can establish that dropout is independent of unobserved future outcomes conditional on observed history. If people with the same recorded history still differ in retention because of an unmeasured impending health decline, the assumption may fail. This is dependence on information unavailable at the decision point, not a problem solved by inserting the missing \(Y_t\) into the retention model. In that setting, weighting alone cannot identify the full-data distribution without additional assumptions or information.

  • Common mistake: “The outcome predicts dropout, so I should include the outcome at the visit that was missed.” That outcome is unavailable for dropouts and may occur after the retention decision.
  • Correction: Include outcomes and other variables observed before the decision, such as \(Y_{t-1}\). State the resulting conditional-independence assumption explicitly and assess whether the measured history plausibly captures the factors linking dropout to future outcomes.
  • Another pitfall: Treating a good fit among retained participants as evidence that conditional independence holds. Model fit assesses the observed retention process, not its dependence on unobserved future outcomes.

The central limitation is thus substantive as well as statistical: conditional modeling can address dropout driven by recorded prior outcomes only to the extent that those outcomes and other observed history capture the relevant selection. It cannot recover information about dropout mechanisms that remain unmeasured.

Check Your Understanding

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

  1. A participant’s symptom score at visit \(t\) is recorded only after the participant attends. Can it predict retention at visit \(t\)?
  2. Why might a model containing only baseline covariates be inadequate when dropout depends on prior health?
  3. In the two-visit example, what is the probability of remaining observed through both visits and the resulting inverse weight?
  4. Does a well-fitting retention model prove that dropout is conditionally independent of future outcomes?

Solutions. (1) No, not if attendance is the prerequisite for measuring the score or if the decision precedes measurement. Use information available before the decision, such as the prior visit’s score. (2) Participants with the same baseline covariates may have different prior health and different retention probabilities; omitting that history can leave selection associated with future outcomes. (3) The product is \(0.8808(0.8320)=0.7328256\), or \(0.7328\) rounded; its reciprocal is about \(1.3646\). (4) No. Observed-data fit cannot verify independence from outcomes that are unobserved after dropout.

Clinical Biostats