Tutorials › Pattern-Mixture and Selection Models › Building a Trial Model With Current-Outcome Dropout

Diggle-Kenward and Related Models · Tutorial 112 of 499

Building a Trial Model With Current-Outcome Dropout

Learn to build and interpret a joint outcome and current-outcome dropout model for a multi-visit randomized trial.

Advanced 12 min read

What You'll Learn

  • Specify a longitudinal outcome equation with treatment, visit timing, baseline covariates, and within-person variation.
  • Define a monotone dropout history and a visit-specific dropout hazard.
  • Include current and lagged outcomes in a selection equation without confusing their timing.
  • Calculate a dropout probability from a logistic selection model.
  • Recognize assumptions, diagnostics, and limitations of a current-outcome selection model.

Why Build the Outcome and Dropout Models Together?

In a multi-visit randomized trial, outcomes are scheduled over time, but some participants stop providing measurements. If dropout depends on the outcome at the visit where measurement is missed, the selection model must refer to that current outcome even though it is unavailable for a participant who drops out. A useful specification therefore has two connected parts: a model for the outcome trajectory and a model for dropout conditional on the outcome history.

The goal here is to make the model’s ingredients explicit: the treatment assignment, the visit schedule, baseline covariates, observed and latent outcomes, and the sequence of dropout opportunities. A current-outcome selection model is an assumption about how the joint data arise; it does not turn an unobserved outcome into an observed value.

Define the Trial Data and the Dropout History

Suppose participant \(i\) is randomized to treatment \(A_i\), where \(A_i=1\) denotes the active treatment and \(A_i=0\) the control. Let \(Y_{i0}\) be a measured baseline outcome, and let \(Y_{ij}\) be the outcome scheduled at follow-up visit \(j\), for \(j=1,\ldots,J\). The scheduled time of visit \(j\), measured in a specified unit such as weeks, is \(t_j\). Let \(X_i\) be a vector of baseline covariates other than \(Y_{i0}\).

Let \(R_{ij}=1\) if the outcome at visit \(j\) is observed and \(R_{ij}=0\) otherwise. Assume monotone dropout: after a participant stops providing outcomes, no later outcome is observed. Define \(R_{i0}=1\), and define \(D_{ij}=1\) when a participant who was still providing outcomes through visit \(j-1\) drops out before providing \(Y_{ij}\). Thus, among participants with \(R_{i,j-1}=1\), \(D_{ij}=1-R_{ij}\). A participant who already dropped out is not at risk of a new dropout event at later visits.

Definition: In a current-outcome dropout model, the dropout probability at visit \(j\), conditional on being at risk, is allowed to depend on \(Y_{ij}\), the outcome scheduled at that same visit. For a participant who drops out before providing that measurement, \(Y_{ij}\) is latent in the joint model.

This setup describes permanent dropout, not every missed visit. If participants can miss one visit and return later, a single monotone dropout indicator is not enough; the observation process must represent that intermittent missingness separately. The outcome and selection equations should also respect temporal order: baseline covariates and treatment assignment precede follow-up, earlier outcomes precede later visits, and the current outcome in the selection equation is the outcome scheduled at the dropout opportunity.

Specify the Outcome Trajectory

One useful continuous-outcome specification is a linear mixed model with a participant-specific random intercept. It allows average outcome levels to vary by treatment, time, and baseline characteristics, while accounting for repeated measurements from the same participant:

$$ Y_{ij}=\beta_0+\beta_A A_i+\beta_t t_j+\beta_{At}A_i t_j+\gamma^\mathsf{T}X_i+\gamma_0Y_{i0}+b_i+\varepsilon_{ij}, \qquad b_i\sim N(0,\tau^2),\quad \varepsilon_{ij}\sim N(0,\sigma^2). $$

Here \(b_i\) is a random intercept, and the residuals \(\varepsilon_{ij}\) are assumed independent across visits and independent of \(b_i\), conditional on the model covariates. The treatment-by-time term allows the average treatment contrast to change over time. The coefficient \(\beta_t\) is the control-arm change in mean outcome per unit of time; the active-arm change per unit of time is \(\beta_t+\beta_{At}\). If the trajectory is not plausibly linear in \(t_j\), visit indicators or a suitable time basis can replace the linear time terms. The time coding and units must be stated so coefficients can be interpreted.

This is a model for the outcome process, including outcomes that may be missing. It does not say that every participant’s observed trajectory follows the population mean exactly: \(b_i\) and \(\varepsilon_{ij}\) represent between-participant and within-participant variation. The distributional assumptions should be chosen to suit the outcome; the displayed Gaussian form is an example for a continuous outcome, not a universal requirement.

Worked Example: Outcome Mean at a Scheduled Visit

Suppose the outcome is measured in biomarker units, time is measured in weeks, and the fitted mean equation is \(10+1.0A_i-0.10t_j-0.05A_it_j+0.20Y_{i0}+b_i\). For a hypothetical active-arm participant at week 12, let the centered baseline outcome be \(Y_{i0}=3\) units and the estimated participant random intercept be \(b_i=0.4\) units. Assume the Gaussian outcome model above.

The conditional mean, given \(b_i\) and the covariates, is

$$ 10+1.0(1)-0.10(12)-0.05(1)(12)+0.20(3)+0.4 =10+1.0-1.2-0.6+0.6+0.4 =10.2\text{ units}. $$

Thus, under this hypothetical fitted model, the conditional mean at week 12 is 10.2 biomarker units. This calculation is conditional on the participant effect \(b_i=0.4\); it is not the population mean for all active-arm participants. For the population mean conditional only on treatment and baseline covariates, the random intercept has mean zero and its \(0.4\) contribution would not be included.

Specify Dropout Conditional on Current and Earlier Outcomes

A logistic hazard gives the probability of dropout at visit \(j\), conditional on having remained in follow-up through visit \(j-1\). To make outcome effects comparable across visits, define \(Z_{ij}=(Y_{ij}-c)/s_Y\), where \(c\) and \(s_Y>0\) are fixed centering and scaling constants in outcome units. For the lagged value, set \(Z_{i,j-1}=(Y_{i,j-1}-c)/s_Y\) when \(j\geq2\); at the first follow-up visit, use the baseline outcome \(Y_{i0}\) in place of \(Y_{i,j-1}\). Write that lagged standardized value as \(Z_{i,\mathrm{lag}(j)}\).

$$ \operatorname{logit}\{P(D_{ij}=1\mid R_{i,j-1}=1,\mathbf{Y}_i,A_i,X_i)\} =\alpha_0+\alpha_AA_i+\alpha_t t_j+\alpha_C Z_{ij} +\alpha_L Z_{i,\mathrm{lag}(j)}+\delta^\mathsf{T}X_i. $$

Here \(\mathbf{Y}_i=(Y_{i1},\ldots,Y_{iJ})\) denotes the full outcome history, including latent outcomes. The logit of a probability \(p\) is \(\log\{p/(1-p)\}\). A positive \(\alpha_C\) means that, holding the other terms fixed, higher current outcomes correspond to higher dropout odds. A positive \(\alpha_L\) describes an association with the previous outcome after accounting for the current outcome. These are conditional model associations, not causal effects of changing an outcome.

Including \(t_j\) allows the baseline dropout odds to vary over the scheduled follow-up times. Baseline covariates can adjust dropout probabilities for prognostic characteristics, while treatment can account for different dropout odds by randomized arm. The model can be expanded when justified—for example, with treatment-by-time or treatment-by-current-outcome interactions—but each added term requires adequate information and a clear interpretation. At the first follow-up, the lagged value is baseline; later, it is the immediately preceding follow-up outcome.

Worked Example: Calculating a Current-Outcome Dropout Probability

Consider a hypothetical dropout equation \(\operatorname{logit}(p_{ij})=-1.40+0.25A_i+0.10j+0.50Z_{ij}-0.20Z_{i,\mathrm{lag}(j)}+0.15F_i\), where \(p_{ij}\) is the dropout probability at visit \(j\), \(F_i=1\) for a specified baseline group and 0 otherwise, and visit number \(j\) is used as the time scale. Suppose a participant is in the active arm, is at visit 2, has \(Z_{ij}=-1.0\), \(Z_{i,\mathrm{lag}(j)}=-0.5\), and \(F_i=1\).

First substitute these values:

$$ \eta_{ij}=-1.40+0.25(1)+0.10(2)+0.50(-1.0)-0.20(-0.5)+0.15(1)=-1.20. $$

Convert the log odds to a probability:

$$ p_{ij}=\frac{\exp(-1.20)}{1+\exp(-1.20)} =\frac{0.3010}{1+0.3010} \approx0.2315. $$

The modeled probability of dropout at visit 2, conditional on still being at risk and on the specified covariates and outcomes, is 23.15%. This is a model-based conditional probability, not the observed dropout proportion for the trial. The calculation assumes the stated logistic form and uses the current outcome value as part of the model even though that value would be unavailable if this participant actually dropped out before measurement.

How the Two Equations Form a Joint Model

The outcome equation supplies a distribution for the full outcome trajectory; the dropout equation supplies conditional probabilities for the response history given that trajectory. A convenient conceptual factorization is the outcome density multiplied by the visit-specific dropout probabilities for participants still at risk. For example, a participant observed at visit 1 and dropping out before visit 2 contributes an outcome density for the observed \(Y_{i1}\), a probability of remaining observed at visit 1, and a visit-2 dropout probability that depends on \(Y_{i2}\) as well as the history. The \(Y_{i2}\) term is latent in this dropout probability.

The reason is the chain rule for conditional probabilities: the model specifies the distribution of outcomes and then specifies dropout at each eligible visit conditional on the outcomes and covariates. For a dropout event at visit \(j\), the selection factor is \(P(D_{ij}=1\mid R_{i,j-1}=1,\mathbf{Y}_i,A_i,X_i)\), not a probability conditioned only on outcomes already observed. Because \(Y_{ij}\) is missing for that event, calculations based on observed data must account for its uncertainty under the joint model. This is a justification for the model structure, not a recipe for replacing a missing current outcome with a guessed value.

Worked Example: Reading a Dropout History Correctly

Suppose a hypothetical participant provides outcomes at visits 1 and then drops out before visit 2, with no later outcomes. The response history is \(R_{i1}=1,R_{i2}=0,R_{i3}=0\), and dropout occurs at visit 2. Suppose the visit-2 equation is the one in the preceding example, but this participant is in control (\(A_i=0\)), has \(Z_{i2}=-0.8\), \(Z_{i1}=0.4\), and \(F_i=0\).

The visit-2 logit is \(-1.40+0.25(0)+0.10(2)+0.50(-0.8)-0.20(0.4)+0.15(0)=-1.68\). Therefore,

$$ P(D_{i2}=1\mid R_{i1}=1,\mathbf{Y}_i,A_i,X_i) =\frac{\exp(-1.68)}{1+\exp(-1.68)} =\frac{0.1864}{1.1864} \approx0.1571. $$

Under the specified model and these hypothetical outcome values, the conditional visit-2 dropout probability is 15.71%. The participant’s observed history tells us that \(Y_{i1}\) is available but \(Y_{i2}\) is not. Thus \(Z_{i2}=-0.8\) is illustrative of a possible latent value, not a value that could be read from this participant’s observed record. At visit 3 the participant is no longer at risk, so the model does not assign a new dropout hazard there.

Assumptions, Diagnostics, and Alternatives

A coherent specification requires several assumptions. The dropout process must be monotone for the hazard definition above; randomization and baseline covariates must be correctly recorded; the outcome distribution and time structure must represent the repeated measurements adequately; and the logistic equation must capture the chosen dependence on current outcome, earlier outcome, treatment, time, and baseline covariates. There must also be meaningful support for the covariate and outcome combinations used to interpret the model. In particular, a current-outcome association in the dropout equation is difficult to learn from observed data alone because the relevant outcome is missing precisely for some dropout events.

Useful checks begin with the model’s construction and observed implications:

  • Confirm that the visit times, time units, treatment coding, and baseline covariate coding match the trial design.
  • Verify that each dropout hazard is defined only for participants still at risk, and that the current and lagged outcomes refer to the intended visits.
  • Compare the fitted outcome trajectories with observed summaries by treatment and visit, while recognizing that observed summaries can be selected by dropout.
  • Inspect fitted dropout probabilities across visits and across plausible covariate and outcome values. Probabilities that are nearly zero or one over ordinary ranges may indicate an implausible specification or sparse information.
  • Check whether the assumed linear time effect is credible; systematic departures may motivate visit-specific effects or a more flexible time basis.

A lagged-outcome dropout model uses earlier outcomes but not the current latent outcome; it is a different selection assumption. A model depending only on observed covariates and observed history is another alternative, and a pattern-mixture model instead organizes assumptions by missingness pattern. These are alternative parameterizations and assumptions, not interchangeable labels for the same model. The appropriate choice depends on the scientific question and on how the unobserved outcomes are to be related to dropout.

The central limitation is that observed records do not directly determine how dropout depends on an outcome that is missing at dropout. Different choices for this dependence can be compatible with similar observed data while implying different distributions for missing outcomes. Therefore, current-outcome coefficients should not be presented as established causal effects or as empirically certain merely because a model has been fitted. Their role is to encode a selection assumption, which should be made explicit and examined for scientific plausibility.

Common Mistakes

  • Using only the last observed outcome while calling the model current-outcome dependent. A term involving \(Y_{i,j-1}\) is lagged; current-outcome dependence requires \(Y_{ij}\).
  • Treating the missing current outcome as observed. In a dropout event, \(Y_{ij}\) enters the joint model as a latent quantity. It is not available just because the selection equation includes it.
  • Applying a dropout hazard after dropout. A participant with \(R_{i,j-1}=0\) is not at risk at visit \(j\) under monotone dropout.
  • Ignoring the visit schedule or units. A time coefficient has no unambiguous interpretation unless time coding and units are stated.
  • Interpreting a selection coefficient causally. A conditional association between outcome level and dropout does not show that changing the outcome would cause a change in dropout.
Key takeaway: A current-outcome dropout model must specify both the outcome trajectory and the dropout probability conditional on the full outcome history. The current outcome at a dropout visit is latent, so the model states an assumption about selection rather than supplying an observed value.

Check Your Understanding

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

  1. In a trial with visits at weeks 4, 12, and 24, what does the coefficient of \(A_i t_j\) represent in the outcome equation?
  2. Using the dropout equation from the worked probability example, calculate the dropout probability for a control participant at visit 1 with current standardized outcome \(0\), lagged standardized outcome \(0\), and \(F_i=0\).
  3. A participant has response history \(1,1,0,0\). At which visit is this participant at risk of dropout for the first time, and should a new hazard be applied at the final visit?
  4. Why can a current-outcome coefficient not be interpreted as a causal effect of the outcome on dropout?

Solutions. (1) The coefficient \(\beta_{At}\) is the difference between treatment arms in mean outcome change per week: the control slope is \(\beta_t\), and the active-arm slope is \(\beta_t+\beta_{At}\). (2) The logit is \(-1.40+0.25(0)+0.10(1)+0.50(0)-0.20(0)+0.15(0)=-1.30\). Thus \(p=\exp(-1.30)/\{1+\exp(-1.30)\}=0.2142\), or 21.42%. (3) The participant is at risk at visit 1 and remains at risk at visit 2 because both outcomes are observed; dropout occurs at visit 3. No new hazard applies at visit 4 because the participant has already dropped out. (4) The coefficient describes a conditional association in a selection model involving an outcome that is unobserved for dropout events. It does not arise from an intervention on the outcome, and its value depends on the model’s assumptions about the joint outcome and dropout process.

Clinical Biostats