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)]\).
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.
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.
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.
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
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\):
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
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.
Derivation. Average the assumed representation and use the calibrated equation:
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.
- 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?
- In the first worked example, why does the targeted prediction at \(A=1\) increase by \(0.10\), even though \(\hat\epsilon=0.05\)?
- Does \(P_n\hat D^*=0\) by itself establish that the estimator is asymptotically efficient? State the missing requirement.
- 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.