Why Transform an Outcome?
With treatment at multiple times, an observed outcome reflects the history of treatments actually received. A structural nested mean model represents treatment effects through a sequence of blips: causal increments associated with treatment at each time, after later treatment effects have been removed. Blip-down transformations use those increments to construct an outcome with the modeled treatment effects subtracted.
The resulting quantity is useful because it aims to remove the contribution of treatment from the outcome while retaining its other variation. But a transformed value is not automatically the outcome that the same person would have had under no treatment. It is an algebraically constructed pseudo-outcome; its causal interpretation depends on the model and identification assumptions.
The histories \(\bar A_k=(A_1,\ldots,A_k)\) and \(\bar L_k=(L_1,\ldots,L_k)\) include treatment and covariate information through time \(k\). Thus, the blip at time \(k\) can depend on the treatment and covariate history up to that time. The blip is a conditional mean contrast, not necessarily a fixed effect shared by everyone or an individual-level treatment effect.
The Backward Transformation
Suppose there are \(K\) treatment times. Start with the observed outcome, then remove the modeled blip at the latest time. From that result, remove the blip at the preceding time, and continue backward. Write \(H_{K+1}(\psi)=Y\), and define the successive transformed outcomes as follows:
The order matters conceptually. The blip at time \(k\) describes the effect of that treatment on an outcome from which the modeled effects of later treatments have already been removed. It is therefore applied to the outcome produced at the preceding backward step, not used as a substitute for the entire transformation.
Worked Example: Blipping Down Two Treatment Times
Consider one hypothetical participant with observed outcome \(Y=50\). Suppose a specified blip model, evaluated at this participant’s observed histories, gives a time-2 blip of 7 outcome units and a time-1 blip of 5 outcome units. These values are illustrative model outputs, not estimates from a real study.
Start at the last treatment time. The outcome after removing the time-2 blip is \(H_2=50-7=43\) outcome units.
Next remove the time-1 blip: \(H_1=43-5=38\) outcome units.
For this participant, the transformed pseudo-outcome is 38. This calculation alone does not establish that this participant’s unobserved outcome under no treatment would be 38. Under the model and identification assumptions, the population mean of transformed outcomes can represent the reference-treatment mean; a single participant’s transformed value is not that population mean.
Repeated substitution gives a compact expression for the final transformation:
This equation makes the arithmetic transparent: each modeled blip is subtracted once. It does not, by itself, establish a causal interpretation. That interpretation requires the functions to describe the relevant causal mean increments and requires conditions that let those increments be identified from observed data.
Why Subtracting the Blips Can Remove Treatment Effects
The central idea is that a correctly specified blip accounts for the conditional mean change associated with the current treatment, given the history, on an outcome already adjusted for later treatment effects. Subtracting that increment should remove the corresponding treatment-related difference in conditional means.
More formally, let \(H_{k+1}(\psi)\) be the outcome after later blips have been removed. At the true model value \(\psi_0\), the structural mean model represents the conditional mean contrast between setting \(A_k=a\) and setting \(A_k=0\), for the same prior history, by the blip function. In schematic notation, this contrast is:
When the observed blip is subtracted, the conditional mean contrast between treatment \(a\) and reference treatment \(0\) is removed:
This is a conditional-mean statement, not a claim that the transformed values are identical across people or treatment groups. Individuals can differ in their covariates and in outcome variation not explained by the blip. Applying the argument backward at each treatment time motivates the final transformation toward the reference-treatment outcome.
The assumptions have distinct roles. Consistency links an observed outcome to the potential outcome under the treatment history actually received. Sequential exchangeability requires that, conditional on the measured history available before each treatment, treatment assignment carries no further information about the relevant potential outcomes. Positivity requires that each treatment value being compared has positive probability in histories for which that contrast is to be interpreted. Correct specification of the blip functions is also needed for the stated transformation to remove the modeled contrasts. If these conditions fail, the algebra still produces numbers, but the intended causal interpretation may not hold.
Worked Transformations and Their Interpretation
Worked Example: A Participant’s Transformed Value Is Not a Mean
In a hypothetical two-time setting, suppose a participant has observed outcome \(Y=62\). The fitted blips evaluated at that participant’s histories are 8 outcome units at time 2 and 6 outcome units at time 1.
The backward calculation is \(H_2=62-8=54\) outcome units, followed by \(H_1=54-6=48\) outcome units. Equivalently, \(H_1=62-(8+6)=48\) outcome units.
The value 48 describes this participant’s transformed pseudo-outcome under the specified blip model. It is not the population’s reference-treatment mean, and it does not prove that the participant’s individual reference-treatment counterfactual equals 48. To discuss a reference-treatment mean, one must consider transformed outcomes over the relevant population and rely on the model and identification assumptions.
Worked Example: Averaging Transformed Outcomes
Suppose, purely hypothetically, three participants have observed outcomes and evaluated blips as shown below. All values are in outcome units.
| Participant | \(Y\) | Time-2 blip | Time-1 blip | Transformed value |
|---|---|---|---|---|
| A | 50 | 7 | 5 | 38 |
| B | 58 | 4 | 6 | 48 |
| C | 45 | 5 | 2 | 38 |
For participant A, \(50-7-5=38\). For participant B, \(58-4-6=48\). For participant C, \(45-5-2=38\). The sample mean of the transformed values is \[ \bar H_1=\frac{38+48+38}{3}=\frac{124}{3}=41.33 \] outcome units, rounded to two decimal places.
This calculation gives a sample average of constructed values. It does not by arithmetic alone establish a population causal mean. If the blip model and identification assumptions are appropriate, the population mean of transformed outcomes represents the mean under the reference-treatment history; a sample average is an empirical summary intended to inform that target.
Worked Example: Checking the Backward Order
For another hypothetical participant, let \(Y=40\), with an evaluated time-2 blip of 3 and a time-1 blip of 9 outcome units. The backward transformation first gives \(H_2=40-3=37\), then \(H_1=37-9=28\) outcome units.
The cumulative formula checks the result: \(H_1=40-(3+9)=28\) outcome units. Subtracting both values gives the same arithmetic total whichever order the subtraction is performed in, but the model interpretation is backward: the time-1 blip is defined for an outcome whose later-treatment effect has already been removed. Simply adding treatment effects without specifying what outcome each blip refers to can therefore misstate the causal model even when the arithmetic happens to agree.
What the Transformation Does—and Does Not—Do
Blip-down transformations matter because they turn a model for sequential treatment effects into an outcome whose mean is intended to correspond to a reference treatment history. This offers a way to express causal assumptions in terms of how outcome means change with treatment at each time. It also makes explicit why later treatment effects must be addressed before earlier ones.
The transformation does not erase all differences among people. It subtracts the modeled conditional mean increments; residual variation, covariate-related variation, and model error remain. Nor does it automatically solve confounding: the blip values have a causal interpretation only under suitable identification conditions. In particular, if important time-varying confounders are not measured, subtracting a fitted blip does not repair that omission.
A frequent mistake is to treat a transformed value as the individual’s actual no-treatment outcome. The calculation is tempting: it starts from that person’s observed outcome and subtracts effects, so the result can look like a recovered counterfactual. But a mean blip is a conditional mean contrast, not generally that person’s realized individual effect. The correction is to call the result a pseudo-outcome and distinguish its value for one record from the population mean it may represent under the structural model and identification assumptions.
Check Your Understanding
Answer each question by keeping the algebraic transformation separate from its causal interpretation.
- With \(Y=70\), a time-2 blip of 8, and a time-1 blip of 4, what is the fully transformed pseudo-outcome?
- Why does the blip-down transformation proceed from the last treatment time toward the first?
- What does a correct blip function represent: an individual treatment effect or a conditional mean causal increment?
- If one participant’s transformed value is 32 outcome units, what additional quantity and assumptions are needed before interpreting a reference-treatment mean?
- Name two identification conditions required for a causal interpretation of the transformed outcome, and briefly describe their roles.