Tutorials › Marginal Structural Models › Depression Study Dropout After Symptom Changes

Case Studies · Tutorial 487 of 501

Depression Study Dropout After Symptom Changes

Learn to estimate marginal remission under antidepressant strategies when retention depends on observed symptoms and treatment history, and see why complete-case estimates may differ.

Advanced 13 min read

What You'll Learn

  • Define the treatment and retention histories needed for a depression dropout analysis.
  • State the assumptions for identifying remission under a sustained antidepressant strategy.
  • Derive the censoring-only identity and the combined treatment-and-censoring weighted risk.
  • Calculate inverse-probability weights across multiple follow-up intervals.
  • Compare weighted remission estimates with complete-case estimates.
  • Recognize diagnostics and limitations when symptom-related dropout is common.

Why Symptom-Related Dropout Changes the Estimate

In a depression study, participants may stop attending after symptoms worsen, improve, or change in response to treatment. The remission rate among participants who complete follow-up then describes a selected group, not necessarily the full population that began the study. Weighting can correct this selection when dropout is explained by the measured history, but the weights must account for treatment assignment as well as retention.

Let \(L_t\) denote symptom severity measured at time \(t\), \(A_t\) the treatment received, and \(Y\) the indicator of remission at the end of follow-up, with \(Y=1\) for remission. Let \(R_t=1\) mean that the participant is still observed at time \(t\), with \(R_0=1\). Write \(H_t\) for the observed history available before treatment at time \(t\), including baseline covariates and prior symptoms and treatments. The history used to model retention after treatment at time \(t\) includes \(H_t\) and \(A_t\).

Definition: A treatment strategy \(g\) specifies the treatment \(g_t(H_t)\) to be received at each time \(t\), possibly as a function of observed history. The potential outcome \(Y^g\) is the remission outcome that would occur if the participant followed strategy \(g\) throughout follow-up. The marginal remission risk under \(g\) is \(E[Y^g]\).

For example, \(g\) might specify antidepressant treatment at every visit, while a comparison strategy might specify no antidepressant treatment at every visit. The target \(E[Y^g]\) concerns the full baseline population under that strategy; it is not the remission proportion among the people who happened to follow it and remain observed.

Assumptions and the Two Sources of Weighting

The analysis requires consistency: a participant who follows \(g\) has observed outcome \(Y=Y^g\). It also requires sequential treatment exchangeability: conditional on the measured history \(H_t\), treatment at time \(t\) is independent of the potential outcomes under the strategies being compared. In addition, sequential independent censoring requires that, among those still observed, retention at the next interval be independent of \(Y^g\) conditional on the measured history and treatment received. This is an outcome-specific condition; it does not assert independence from every unobserved future quantity.

Define the treatment probability \(\pi_t(a\mid H_t)=P(A_t=a\mid R_t=1,H_t)\) and the retention probability \(q_t(H_t,A_t)=P(R_{t+1}=1\mid R_t=1,H_t,A_t)\). Both probabilities must be positive for the histories that can occur under \(g\): \(\pi_t(g_t(H_t)\mid H_t)>0\) and \(q_t(H_t,g_t(H_t))>0\). These are treatment and retention positivity conditions. They rule out histories at which the strategy is never followed or retention is impossible.

For a participant who is observed through the final outcome and whose observed treatment path follows \(g\), the combined inverse-probability weight is

$$ W^g= \frac{I(\bar A=\bar g)\,R_K} {\displaystyle\prod_{t=0}^{K-1} \pi_t\!\left(g_t(H_t)\mid H_t\right) q_t\!\left(H_t,g_t(H_t)\right)}. $$

Here \(I(\bar A=\bar g)\) indicates that treatment followed the strategy at every interval, and \(R_K\) indicates that the final outcome was observed. Histories and probabilities in this expression are evaluated along the participant’s observed path. The treatment factors correct for the fact that only some participants follow \(g\); the retention factors correct for the fact that only some of those participants remain observed.

Two Weighting Identities

First consider a population actually generated under intervention \(g\), rather than observational data in which participants may receive other treatments. Write \(R_K^g\) for retention under that intervention and \(q_t^g\) for the retention probability along its history. Under the stated sequential independent-censoring condition for \(Y^g\), the censoring-only identity is

Result: If \(Y^g\) is integrable, retention has positive probability at each relevant history, and \(Y^g\) is conditionally independent of each next retention indicator given retention through time \(t\) and the intervention history, then
$$ E\left[\frac{R_K^gY^g}{\prod_{t=0}^{K-1}q_t^g}\right]=E[Y^g]. $$

To see the logic, at an interval with history fixed, the conditional expectation of the retention factor is one. In particular, when \(R_t^g=1\), independent censoring for this outcome gives

$$ E\left[\left.\frac{R_{t+1}^g}{q_t^g}Y^g\right| R_t^g=1,H_t^g,Y^g\right] = Y^g\frac{P(R_{t+1}^g=1\mid R_t^g=1,H_t^g,Y^g)}{q_t^g} =Y^g. $$

Applying this equality successively from the last interval back to the first removes each retention indicator and its probability factor, leaving \(E[Y^g]\). The argument uses independence for \(Y^g\), not for an arbitrary full-data quantity.

In observational data, this identity alone does not identify \(E[Y^g]\): a participant who received a different treatment path does not reveal \(Y^g\) by consistency. Treatment weighting is also necessary. Under consistency, sequential treatment exchangeability, sequential independent censoring for \(Y^g\), and both positivity conditions, the observed-data identification result is

Theorem: For an integrable outcome, the combined treatment-and-retention weight satisfies
$$ E\left[ \frac{I(\bar A=\bar g)R_KY} {\displaystyle\prod_{t=0}^{K-1} \pi_t(g_t(H_t)\mid H_t)q_t(H_t,g_t(H_t))} \right] =E[Y^g]. $$

A direct justification follows the sequential data-generating process. At each interval, the indicator that treatment matches \(g_t\), divided by its conditional probability, replaces the observational treatment draw with the strategy’s specified treatment. Among those retained, the next retention indicator divided by its conditional probability removes the selection due to dropout. Iterating these conditional expectation steps over \(t=0,\ldots,K-1\) leaves the distribution of histories and outcomes under \(g\), and hence \(E[Y^g]\). The treatment step uses exchangeability to equate the outcome distribution among those who received the specified treatment, conditional on history, with the corresponding potential-outcome distribution. The retention step uses the outcome-specific censoring condition. Neither step licenses replacing the strategy-specific potential outcome by an observed outcome for someone who followed another treatment path. \(\square\)

For independent participants, a simple inverse-probability weighted estimator is \(n^{-1}\sum_{i=1}^n W_i^gY_i\). A common finite-sample alternative is the Hájek ratio \(\sum_iW_i^gY_i/\sum_iW_i^g\), which normalizes by the observed sum of weights. Both target the same risk under the identification assumptions; the ratio is not exactly unbiased in finite samples. Consistency of these sample estimates additionally requires suitable sampling and moment conditions, such as independent sampling and finite expectations for the weighted terms.

Worked Examples

Worked Example: Dropout After Symptom Worsening

Suppose that under one antidepressant strategy, half the population has severe symptoms at an intermediate visit and half has mild symptoms. The eventual remission probabilities are 0.20 for the severe group and 0.80 for the mild group. Severe symptoms are followed by retention probability 0.25; mild symptoms by retention probability 0.75. Assume retention is independent of eventual remission conditional on symptom group.

The target remission risk is the population average:

\(E[Y^g]=0.50(0.20)+0.50(0.80)=0.10+0.40=0.50.\)

Among retained participants, the severe group contributes \(0.50(0.25)=0.125\) of the starting population and the mild group contributes \(0.50(0.75)=0.375\). Thus the retained sample has severe fraction \(0.125/(0.125+0.375)=0.25\) and mild fraction \(0.375/0.500=0.75\). Its complete-case remission risk is \(0.25(0.20)+0.75(0.80)=0.05+0.60=0.65\), not 0.50.

Inverse retention weighting recovers the target. In the starting population, the weighted remission contribution is \(0.50(0.25)(0.20)/0.25+0.50(0.75)(0.80)/0.75=0.10+0.40=0.50\). The weighted denominator is \(0.50(0.25)/0.25+0.50(0.75)/0.75=0.50+0.50=1\). Their ratio is therefore 0.50. The complete-case distortion occurs because symptom severity predicts both retention and remission.

Worked Example: Combining Treatment and Retention Probabilities

At two treatment intervals, suppose a participant’s probability of receiving the specified strategy treatment is 0.60 at the first history and 0.50 at the second. Given that treatment path, the probabilities of remaining observed through the next intervals are 0.80 and 0.50. The participant follows the specified treatments, remains observed through the final visit, and is in remission.

The treatment-path probability is \(0.60(0.50)=0.30\). Conditional on that path, the retention probability is \(0.80(0.50)=0.40\). The joint probability represented in the denominator is \(0.30(0.40)=0.12\), so the combined weight is \(1/0.12=8.333\ldots\). This participant contributes \(8.333\ldots\) to the weighted remission total and to the weighted denominator, before any normalization by the cohort size.

Using only the retention weight would give \(1/0.40=2.5\), which fails to correct for the selective probability of following the treatment path. Using only the treatment weight \(1/0.30=3.333\ldots\) would fail to correct for dropout. The combined weight addresses both selection steps, provided the probabilities are correctly specified and the assumptions hold.

Worked Example: Weighted and Complete-Case Contrasts

Consider a small illustrative analysis with four observed participants following each of two strategies. The listed weights are their realized combined treatment-and-retention weights. Remission is coded 1 and no remission 0.

StrategyParticipant weightsRemission outcomesComplete-case riskWeighted risk
Antidepressant strategy \(g_1\)2, 2, 1, 11, 0, 0, 01/42/6 = 1/3
Comparison strategy \(g_0\)1, 1, 2, 21, 1, 1, 03/44/6 = 2/3

For \(g_1\), the complete-case risk is \((1+0+0+0)/4=1/4\). The weighted numerator is \(2(1)+2(0)+1(0)+1(0)=2\), and the sum of weights is \(2+2+1+1=6\), giving the Hájek risk \(2/6=1/3\).

For \(g_0\), the complete-case risk is \((1+1+1+0)/4=3/4\). The weighted numerator is \(1(1)+1(1)+2(1)+2(0)=4\); the weight sum is \(1+1+2+2=6\), giving \(4/6=2/3\). The complete-case risk difference, \(g_1-g_0\), is \(1/4-3/4=-1/2\). The weighted difference is exactly \(1/3-2/3=-1/3\), approximately \(-0.333\). The example shows how weighting changes the contrast when high-weight participants have different outcomes from low-weight participants; it does not, by itself, establish that either estimate is correct for a real study.

Diagnostics, Alternatives, and Limits

The central diagnostic question is whether the observed histories support the required probability models. Check whether the retention model includes recent symptom severity and treatment history when these predict dropout. Compare predicted retention with observed retention across symptom and treatment-history groups, and inspect whether important histories have estimated probabilities near zero. Likewise, assess whether treatment probabilities are plausible across histories. A model that omits symptom changes can leave the very selection mechanism of concern unaddressed.

Complete-case analysis is simple but generally targets the retained, selected population when retention depends on prognosis. Inverse-probability weighting instead targets the baseline population under each strategy if its assumptions hold. Likelihood-based methods and multiple imputation are alternative approaches to missing outcomes, but they also require assumptions about the missingness process and a correctly specified model; they do not remove the need to justify how dropout relates to observed and unobserved information. The introductory material on informative loss to follow-up provides broader context for these approaches.

Weighting cannot recover information absent from the measured history. If dropout depends on unrecorded worsening symptoms even after conditioning on recorded history, sequential independent censoring fails. If few participants with a particular symptom history remain observed, positivity may fail or be practically weak; estimates can then be unstable and highly dependent on a few participants. Model misspecification, measurement error in symptoms, and treatment deviations also threaten the interpretation. A weighted estimate is therefore not automatically unbiased merely because weights were calculated.

  • Common mistake: Treating the retention-weighted outcome average among people who followed the strategy as an estimate for the full population, without treatment weights. The correction is to include both treatment-path and retention probabilities, under the corresponding exchangeability and positivity assumptions.
  • Another pitfall: Assuming that censoring is ignorable because dropout rates are similar overall across treatment groups. Dropout can still depend strongly on recent symptoms within groups; assess retention conditional on the relevant measured history.
  • Rounding check: Calculate contrasts from unrounded risks. In the example, \(1/3-2/3=-1/3\), approximately \(-0.333\); subtracting rounded values can produce a spurious final digit.

The practical conclusion is that symptom-related dropout is a feature of the causal identification problem, not just a loss of sample size. State the target population and strategy, model retention using the observed symptom and treatment history, include treatment weights when estimating observational strategy risks, and report the assumptions that connect the weighted analysis to marginal remission.

Check Your Understanding

Answer each question, then compare with the solutions below.

  1. Why is inverse retention weighting alone insufficient to estimate remission under a strategy in observational data?
  2. In the first worked example, what fraction of retained participants have severe symptoms, and why does that differ from the starting population fraction?
  3. If treatment-path probability is 0.40 and retention probability along that path is 0.50, what is the combined inverse-probability weight for an observed strategy follower?
  4. Give one assumption under which weighting cannot correct symptom-related dropout.

Solutions. (1) Retention weighting addresses selection into observation, but not selection into following the strategy; treatment-path indicators and treatment probabilities are also needed. (2) The fraction is \(0.125/(0.125+0.375)=0.25\), because the severe group has lower retention; it was 0.50 in the starting population. (3) The joint probability is \(0.40(0.50)=0.20\), so the combined weight is \(1/0.20=5\). (4) For example, weighting fails to identify the target if dropout depends on unrecorded symptoms that predict remission, violating sequential independent censoring conditional on the measured history.

Clinical Biostats