Why Blip-Function Choices Matter
A structural nested mean model describes how a treatment changes an outcome through causal contrasts called blips. Choosing a blip function means deciding how those contrasts vary with treatment amount, treatment time, and the history available before treatment. That choice is substantive: it determines which causal patterns the model can represent and which it rules out.
A compact model can be useful when its restrictions reflect the target scientific question. But a simple equation is not automatically a good choice. A constant treatment coefficient, for example, says that the treatment blip is the same at every modeled occasion and history unless the model includes terms that allow it to vary. The aim is not to make the model as flexible as possible; it is to choose a structure that represents the intended contrast without imposing unsupported restrictions.
This definition makes the reference value and outcome scale explicit. For binary treatment, 0 is commonly the untreated reference and the contrast of interest is treatment versus no treatment. For a continuous dose, the blip at dose \(a_t\) is relative to dose 0. A different reference can be chosen, but the function must then be defined relative to it. In an additive mean model, the blip is measured in outcome units.
Functional-form choice is separate from identification. A causal interpretation still requires assumptions such as consistency, sequential exchangeability given the measured history, and positivity for the treatment levels being compared. Even when those conditions hold, a misspecified blip function may fail to represent the causal contrasts of interest. Conversely, a well-chosen function does not repair a failure of identification assumptions.
Time Structures: One Effect or Several?
When treatment can occur at multiple times, the model must say how the blip varies across those times. A constant effect structure uses one common parameter. A time-specific structure permits separate blips at each occasion. A trend structure lies between those choices: it allows effects to vary, but only according to a specified pattern.
The common-effect form says that a unit of treatment at each modeled time has the same blip, on the chosen scale and conditional on the included history. The time-specific form makes no such equality restriction across visits. The trend form assumes a particular, smooth change in the blip across the chosen time index. Equally spaced visits may justify using the visit number as a linear index; irregularly spaced visits may call for actual elapsed time instead.
Worked Example: Comparing Visit-Specific Blips With a Time Trend
Suppose, for illustration, that the causal blips of interest for a binary treatment at visits 1, 2, and 3 are 0.5, 1.1, and 1.8 outcome units, respectively. These hypothetical values are used to compare model restrictions; they are not estimates from an observed study.
A time-specific form represents the three values directly:
A linear trend fitted through the first and third values has slope \((1.8-0.5)/(3-1)=0.65\) outcome units per visit. Its prediction at visit 2 is \(0.5+0.65(2-1)=1.15\), not 1.1. Thus, a linear time-trend blip cannot represent these visit-specific values exactly. It predicts 0.5, 1.15, and 1.8, so it is close to, but does not exactly reproduce, the stated values.
A least-squares linear trend across all three values also does not reproduce them exactly. The mean visit is \(\bar{t}=2\), and the mean blip is \((0.5+1.1+1.8)/3=1.1333\). The slope is
The fitted intercept is \(1.1333-0.6500(2)=-0.1667\). The fitted blips are therefore 0.4833, 1.1333, and 1.7833 outcome units. The trend is a compact approximation, while the visit-specific form exactly preserves the three stipulated values. Which restriction is appropriate depends on whether the scientific target supports a common trend or requires separate visit contrasts.
Choosing the Treatment-Response Shape
Treatment need not be binary. For a continuous dose, a linear blip assumes that each additional unit of dose changes the causal contrast by the same amount throughout the modeled range. A nonlinear blip can express curvature, thresholds, or diminishing changes. The treatment scale and its units must be stated, because coefficients depend on them.
For any dose-response form with a zero-dose reference, the blip should generally satisfy \(\gamma_t(h_t,0;\psi)=0\). A model with an intercept that violates this condition would assign a nonzero contrast to the reference dose unless it were redefined accordingly. This is a useful structural check before interpreting coefficients.
Worked Example: A Curved Blip for Continuous Dose
Suppose dose \(a\) is measured in units of 10 mg/day, and the hypothetical mean-scale blip is
The coefficient 2 has units of outcome units per dose unit, and the coefficient \(-0.25\) has units of outcome units per squared dose unit. At zero dose, \(\gamma(0)=2(0)-0.25(0)^2=0\), so the reference condition is satisfied.
At \(a=1\), the blip is \(2(1)-0.25(1)^2=1.75\) outcome units. At \(a=2\), it is \(2(2)-0.25(2)^2=4-1=3\) outcome units. At \(a=4\), it is \(2(4)-0.25(4)^2=8-4=4\) outcome units.
A linear model \(\gamma(a)=2a\) would instead predict 2, 4, and 8 outcome units at these doses. The difference is not a minor change in notation: the quadratic form represents diminishing increases in the blip as dose rises, while the linear form does not. Over the range from 0 to 4 dose units, the derivative is \(2-0.5a\), which is nonnegative and falls from 2 to 0. These calculations describe the implications of the specified form; they do not establish that this shape is correct for a particular application.
Matching the Time Index to Treatment Timing
“Time-varying” can refer to different scientific structures. A calendar-time model lets the blip depend on the visit or elapsed study time. A duration model lets it depend on how long the person has been exposed. These are not interchangeable when people start treatment at different times. The index should follow the causal contrast: ask whether the intended change is by study occasion, by time since initiation, or by another pre-treatment feature of the history.
For example, a model indexed only by visit number can assign different blips to visits 1 and 3 even for people with different treatment start dates. A duration-indexed model instead assigns blips according to the time elapsed since initiation. If that duration is determined from treatment history, it must be defined using information available before the treatment whose blip is being modeled.
Worked Example: Calendar Time Versus Time Since Initiation
Imagine a hypothetical setting in which people can initiate treatment at different visits and continue treatment. The intended contrast concerns the blip associated with current treatment, and subject-matter reasoning suggests that this blip falls with time since initiation. Suppose the stipulated blips at durations \(d=1,2,3\) are 2, 1, and 0.5 outcome units.
A duration-specific form can represent these values directly. One compact alternative is an exponential decay form:
Substitution gives \(\gamma(1)=2(1/2)^0=2\), \(\gamma(2)=2(1/2)^1=1\), and \(\gamma(3)=2(1/2)^2=0.5\) outcome units. It reproduces the stipulated pattern exactly. In contrast, a linear function of duration fitted through the first and third values has slope \((0.5-2)/(3-1)=-0.75\) outcome units per duration step, and predicts \(2-0.75=1.25\) at duration 2 rather than 1.
If the scientific question instead concerns how blips differ by study visit, a calendar-time index would be the relevant choice. Using visit number when the target is duration since initiation can mix different exposure durations together; using duration when the target is a calendar-time contrast can obscure visit-specific changes. The example specifies a model for the current-treatment blip, not the total effect of a sustained treatment strategy.
A Practical Specification Process
A disciplined choice starts with the contrast, not the formula. In particular, state the outcome scale, reference treatment, time index, and the kinds of variation the target requires. Then decide which restrictions are defensible: a common treatment effect, a smooth trend, distinct visit effects, or a nonlinear dose curve. If effect modification is relevant, include it only in a form tied to the modifier and contrast of interest; the same logic applies to time and dose.
Name the outcome scale and the treatment level used as the reference. For a dose model, state the dose units and range.
Decide whether the blip should vary by visit, elapsed calendar time, time since initiation, or another feature of the pre-treatment history.
Use a common coefficient, time-specific parameters, a trend, or a nonlinear dose function only if its restrictions are interpretable for the target contrast.
Substitute meaningful treatment levels and times into the function. Check the reference value, units, and whether the resulting blips match the causal pattern the model is meant to express.
These choices matter because G-estimation uses the specified blip model in its estimating conditions. Under the identification assumptions, a correctly specified blip function can support the intended causal interpretation. If the function is misspecified, the fitted parameters need not equal the causal contrasts the analyst intended to describe. A flexible model is not automatically preferable, either: it may require more information than the data provide and may estimate the target less precisely. The central task is to balance scientific alignment and a defensible level of restriction.
Common Mistake
A common mistake is to call any model containing a time term “time-varying enough.” For example, a linear time trend does allow blips to change, but it also requires the changes to follow a straight line in the chosen time index. It cannot generally reproduce arbitrary visit-specific effects, and it may use the wrong index if the scientific target is time since treatment initiation rather than calendar time.
The correction is to interpret the restrictions explicitly. Compare the implied blips at relevant treatment times and levels, as in the worked examples. Then choose a shared effect, trend, or separate effects according to the target contrast and scientific rationale—not merely because one form is convenient to write down.
Check Your Understanding
For each question, focus on the restrictions implied by the chosen blip function.
- What does a common additive blip across visits assume that a time-specific blip does not?
- A linear time trend predicts blips of 0.4 and 1.6 at visits 1 and 3. What does it predict at visit 2, and what restriction produced that prediction?
- Why should a dose-response blip relative to zero dose ordinarily equal zero when the dose is zero?
- When treatment initiation dates vary, what scientific distinction separates a calendar-time index from a time-since-initiation index?
- Why does choosing a plausible blip function not, by itself, establish a causal interpretation?