Tutorials › Marginal Structural Models › TMLE Versus Untargeted Sequential Regression

Longitudinal Targeted Maximum Likelihood · Tutorial 397 of 501

TMLE Versus Untargeted Sequential Regression

See how longitudinal TMLE calibrates sequential outcome regressions to the efficient score, and why an untargeted plug-in estimate need not satisfy that calibration.

Advanced 14 min read

What You'll Learn

  • Define the untargeted sequential-regression plug-in estimator for a longitudinal intervention.
  • Derive the longitudinal efficient influence function and its clever covariates.
  • Work out how sequential targeting enforces empirical score equations.
  • Distinguish finite-sample score calibration from consistency and efficiency.
  • Diagnose unstable updates, poor positivity, and overinterpretation of targeting.

Why Compare These Estimators?

Both longitudinal TMLE and untargeted sequential regression begin with models for the outcome process and use them to estimate a counterfactual mean. Their difference is not simply that one uses regression and the other does not: both use sequential outcome regressions. The distinction is that TMLE updates those regressions so their residuals satisfy estimating equations linked to the efficient influence function. The untargeted plug-in estimator does not generally impose those equations.

That calibration matters for inference and for the estimator’s response to errors in nuisance fits. It does not guarantee that the outcome models or treatment mechanism are correct, nor does it repair a positivity failure. The comparison below uses a fixed, finite treatment horizon and makes the assumptions explicit.

Setup and the Untargeted Plug-In Estimate

Let \(O=(W,A_0,S_1,A_1,\ldots,S_K,A_K,Y)\) denote one observed record. Here \(W\) is baseline information; \(S_t\) is the measured history just before treatment \(A_t\); and \(Y\) is the terminal outcome. The history \(S_t\) contains \(W\) and the observed treatment and covariate history available before \(A_t\). Let \(g_t(a\mid s)\) specify the intervention’s distribution of treatment at time \(t\), and let \(\pi_t(a\mid s)=P(A_t=a\mid S_t=s)\) be the observed treatment mechanism.

For a fixed regime \(g=(g_0,\ldots,g_K)\), define \(Q_{K+1}=Y\). Working backward, define \(Q_t(s,a)\) as the conditional expectation of the next intervention-averaged outcome, given the current history and treatment. Write \(\bar Q_t(s)\) for its average under the intervention. This recursion defines the target \(\psi=E[\bar Q_0(W)]\).

$$ \bar Q_{K+1}=Y,\qquad Q_t(s,a)=E[\bar Q_{t+1}(S_{t+1})\mid S_t=s,A_t=a],\qquad \bar Q_t(s)=\sum_a g_t(a\mid s)Q_t(s,a), \qquad \psi=E[\bar Q_0(W)]. $$

For continuous treatments, the sum is replaced by an integral. Given fitted regressions \(\hat Q_t\), untargeted sequential regression applies the same backward recursion using those fits and reports \(P_n\hat{\bar Q}_0(W)\), where \(P_n\) denotes the sample average. In other words, it predicts each person’s outcome under the intervention and averages those predictions. It does not, by construction, adjust the fitted regressions to make the efficient-score equations hold.

Assumptions for a causal interpretation: Consistency links observed outcomes to the outcomes under the treatment actually received. Sequential exchangeability requires no unmeasured confounding of each treatment, conditional on its recorded history. Positivity requires that whenever the intervention assigns positive probability to a treatment at a history occurring under that intervention, the observed treatment mechanism also assigns that treatment positive probability. Identification additionally requires well-defined interventions and the measured histories to be correctly ordered. These assumptions identify the intervention mean; they do not ensure that fitted regressions estimate it accurately.

The Efficient-Score Calibration

For clarity, define the cumulative clever covariate \(H_t\) as the product of intervention probabilities divided by observed treatment probabilities through time \(t\). The longitudinal efficient influence function (EIF) for this intervention mean has one weighted residual for each treatment time, plus a baseline term. Its residual compares the next intervention-averaged prediction with the regression for the treatment actually observed.

$$ H_t=\prod_{j=0}^{t}\frac{g_j(A_j\mid S_j)}{\pi_j(A_j\mid S_j)},\qquad D^*(O)=\sum_{t=0}^{K}H_t\{\bar Q_{t+1}(S_{t+1})-Q_t(S_t,A_t)\} +\bar Q_0(W)-\psi . $$

At the final time, the term \(\bar Q_{K+1}(S_{K+1})\) means \(Y\). For a deterministic regime, \(g_j(A_j\mid S_j)\) is zero unless the observed treatment agrees with the regime. The EIF has mean zero at the true distribution under the identification assumptions and the usual regularity conditions for the statistical model.

Theorem: Suppose the treatment probabilities used in \(H_t\) are positive wherever the intervention assigns positive probability, and the sequential outcome regressions and intervention are defined as above. If the fitted regressions are updated so that, at every time \(t\), \(P_n[H_t\{\hat{\bar Q}_{t+1}(S_{t+1})-\hat Q_t(S_t,A_t)\}]=0\), and the estimate is \(\hat\psi=P_n\hat{\bar Q}_0(W)\), then the empirical mean of the fitted EIF is zero.

Proof. Form the fitted EIF by replacing \(Q_t\), \(\bar Q_t\), and \(\psi\) in the display above by their fitted versions and \(\hat\psi\). Averaging that expression gives

$$ P_n\hat D^* =\sum_{t=0}^{K}P_n[H_t\{\hat{\bar Q}_{t+1}(S_{t+1})-\hat Q_t(S_t,A_t)\}] +P_n\hat{\bar Q}_0(W)-\hat\psi. $$

Each summand in the sum is zero by the stated calibration condition. The last two terms cancel because \(\hat\psi=P_n\hat{\bar Q}_0(W)\). Therefore \(P_n\hat D^*=0\). This is a finite-sample algebraic identity, not a claim that \(\hat D^*\) equals the true EIF. \(\square\)

A standard targeting step obtains these equations by fitting a fluctuation submodel for each sequential regression. For a binary outcome, a common choice is a logistic fluctuation with the current regression as an offset and the relevant clever covariate as a covariate. The score equation for its fluctuation parameter is the weighted residual equation above. The regressions are updated backward so each step uses the already updated next-time prediction. The exact fluctuation depends on the outcome support; it must preserve that support.

Worked Comparisons

Worked Example: One-Time Treatment and a Score-Driven Update

Suppose \(A\) is randomized with \(P(A=1)=0.5\), the intervention sets \(A=1\), and an initial regression predicts \(Q(A)=0.4\) for either treatment. Four observations are \((A,Y)=(1,1),(1,0),(0,0),(0,1)\). The untargeted plug-in estimate is \(0.4\), since it averages the initial predictions under \(A=1\).

Here the clever covariate is \(H=I(A=1)/0.5=2I(A=1)\). Use an identity-link fluctuation \(Q_\epsilon(A)=0.4+\epsilon H\). Its empirical score equation is \(P_n[H\{Y-Q_\epsilon(A)\}]=0\). First calculate the score at \(\epsilon=0\):

$$ P_n[H(Y-0.4)] =\frac{2(1-0.4)+2(0-0.4)+0+0}{4} =\frac{0.4}{4}=0.1. $$

Also, \(P_n(H^2)=\{4+4+0+0\}/4=2\). Solving the linear score equation gives \(\hat\epsilon=0.1/2=0.05\). For \(A=1\), the updated prediction is \(0.4+0.05(2)=0.5\), so the targeted plug-in estimate is \(0.5\). It equals the observed treated mean, \((1+0)/2=0.5\), and the treated residuals now sum to zero. The example shows the calibration mechanism; it does not imply that TMLE must equal the treated mean in every design or model.

Worked Example: Two Treatment Times and Cumulative Weights

Consider two treatment times, a deterministic intervention setting both treatments to 1, and an observed treatment probability of 0.5 at each time for every history. For a person who receives \(A_0=A_1=1\), the cumulative clever covariates are \(H_0=1/0.5=2\) and \(H_1=(1/0.5)(1/0.5)=4\). If either observed treatment differs from 1, the corresponding deterministic-regime numerator is zero, so the relevant cumulative covariate is zero.

Suppose, for one record on the intervention path, the current fitted values are \(Q_0(S_0,A_0)=0.45\), \(\bar Q_1(S_1)=0.55\), \(Q_1(S_1,A_1)=0.50\), and \(Y=0.70\). The two weighted residual contributions to the EIF are

$$ H_0\{\bar Q_1(S_1)-Q_0(S_0,A_0)\}=2(0.55-0.45)=0.20, \qquad H_1\{Y-Q_1(S_1,A_1)\}=4(0.70-0.50)=0.80. $$

Their sum is \(1.00\), before adding the baseline term. An untargeted regression fit has no general requirement that these contributions average to zero across the sample. A targeted fit instead updates the sequential regressions to solve the time-specific weighted equations. This one-record calculation illustrates the weights and residuals; the empirical calibration is imposed across all records, not separately on each one.

Worked Example: A Large Clever Covariate Is a Warning

Suppose at a history the intervention assigns treatment 1 with probability \(g_t(1\mid s)=1\), but the observed probability is only \(\pi_t(1\mid s)=0.05\). For a person receiving treatment 1, that time’s likelihood ratio is \(1/0.05=20\). If the next-outcome residual is \(0.30\), its contribution to the weighted score at that time is \(20(0.30)=6\), much larger than the unweighted residual.

This is not evidence that the residual itself is 6, nor is it a correction for absent treatment support. It shows how a few observations can dominate the targeting equation when positivity is weak. Inspect the distribution of cumulative \(H_t\), the largest values, and whether updates are driven by only a small number of records. Truncating or otherwise stabilizing weights changes the estimating procedure and may introduce bias; it is not a free repair.

What Calibration Does—and Does Not—Establish

Under suitable regularity conditions, empirical EIF calibration supports an asymptotically linear estimator. A precise way to see its role is to state the conditions needed for that conclusion rather than treating the finite-sample score equation as a guarantee of efficiency.

Result: Let \(D^*\) be the true EIF, with \(E[D^*]=0\) and finite variance. Suppose the estimator satisfies \(P_n\hat D^*=0\), and write \(\hat D^*=D^*-(\hat\psi-\psi)+R_n\), where \(P_nR_n=o_p(n^{-1/2})\). Then \(\hat\psi-\psi=P_nD^*+o_p(n^{-1/2})\). If observations are independent and identically distributed, the central limit theorem then gives an asymptotically normal estimator with variance \(\operatorname{Var}(D^*)/n\).

Derivation. Average the assumed representation and use the calibrated equation:

$$ 0=P_n\hat D^* =P_nD^*-(\hat\psi-\psi)+P_nR_n. $$

Rearranging gives \(\hat\psi-\psi=P_nD^*+P_nR_n\). The assumed bound on the remainder yields the stated expansion. Since the observations are independent and identically distributed and \(D^*\) has finite variance, the classical central limit theorem applies to \(P_nD^*\). The crucial qualification is the remainder assumption: empirical score calibration alone does not prove it. In practice, controlling the remainder requires adequate nuisance estimation, positivity, and appropriate conditions for the estimation procedure, including cross-fitting where needed. \(\square\)

An untargeted plug-in estimator can still be consistent, and it can sometimes have good large-sample properties. For example, if its sequential regressions consistently estimate the required conditional means under the identification assumptions, the plug-in recursion can estimate the intervention mean. But there is no general finite-sample reason for its weighted residuals to satisfy the EIF equations. Conversely, targeting an incorrect set of nuisance fits does not automatically make the target estimate consistent or efficient. The earlier material on product-rate conditions gives the relevant type of control needed for asymptotic linearity.

Diagnostics, Alternatives, and Limitations

The main practical comparison is between standardization through fitted outcome regressions and score-calibrated standardization. The untargeted estimate is direct and easy to inspect as an average of predicted outcomes. TMLE adds a principled update tied to the EIF and, for bounded outcomes with an appropriate fluctuation, can keep predictions within the outcome range. Neither feature removes the need to examine the data and assumptions. In particular, the update may be unstable when cumulative clever covariates are extreme or the fluctuation approaches a boundary.

Useful checks include confirming that the intervention and observed treatment probabilities use the same treatment coding, verifying the time order in every \(H_t\), checking the range of each updated prediction, and calculating the empirical weighted residual at every time point. For a converged targeting fit, those residual averages should be close to zero (up to numerical tolerance). Also compare initial and updated predictions: a very large change may signal weak positivity, separation, or a poorly behaved fluctuation rather than successful learning from abundant information.

  • Mistake: “TMLE is efficient because its fitted EIF averages to zero.” That is not enough. The fitted EIF must also approximate the true EIF and the remaining bias must be negligible at the required rate.
  • Mistake: “An untargeted estimate is always biased.” Not so. Correctly estimated sequential regressions can yield a consistent plug-in estimate. The distinction is the absence of automatic efficient-score calibration.
  • Mistake: “Targeting fixes positivity.” It does not create observations for treatment-history combinations that are absent or nearly absent. Extreme clever covariates can make the update unstable.
  • Mistake: “A bounded estimate is necessarily well calibrated.” Respecting the outcome bounds is useful, but it does not establish consistency, adequate uncertainty quantification, or a small bias.

Longitudinal inverse-probability weighting is another standard alternative, but it estimates the intervention mean through treatment weights rather than through this sequential-regression plug-in construction. The contrast here is narrower: targeting adds EIF-based calibration to sequential regression, while an untargeted plug-in estimate reports the initial recursion. The choice between methods depends on nuisance estimation, support, stability, and the inferential goal—not on a universal ranking.

Check Your Understanding

Answer each question, then compare with the solutions.

  1. For a deterministic intervention that sets \(A_0=A_1=1\), with observed probabilities \(0.4\) and \(0.5\), what is \(H_1\) for a person who receives both assigned treatments?
  2. In the first worked example, why does the targeted prediction at \(A=1\) increase by \(0.10\), even though \(\hat\epsilon=0.05\)?
  3. Does \(P_n\hat D^*=0\) by itself establish that the estimator is asymptotically efficient? State the missing requirement.
  4. At a history with intervention probability 1 and observed treatment probability 0.02, what is the treatment likelihood ratio for a person receiving that treatment? What practical concern does it raise?

Solutions. (1) \(H_1=(1/0.4)(1/0.5)=2.5\times2=5\). (2) The fluctuation is \(Q_\epsilon=0.4+\epsilon H\), and \(H=1/0.5=2\) for \(A=1\), so the change is \(0.05\times2=0.10\). (3) No. The fitted EIF must approximate the true EIF and the remainder in the asymptotic expansion must be \(o_p(n^{-1/2})\), in addition to suitable regularity conditions. (4) The ratio is \(1/0.02=50\). Such a large value can make the score update dominated by a small number of observations, indicating weak positivity and potential instability.

Clinical Biostats