Why Tilt Treatment Odds?
A deterministic regime asks what would happen if treatment were set to a specified value, while a stochastic intervention specifies a probability distribution for treatment. An incremental intervention offers a particular stochastic alternative: it changes the odds of treatment relative to the observed treatment mechanism without requiring treatment to be assigned with probability one. The incremental factor describes the size and direction of the shift.
This construction is especially useful when deterministic regimes would demand treatment levels rarely observed for some histories. The intervention preserves the support of the observed treatment mechanism, while shifting its probabilities. Its targeting covariate is a cumulative product of intervention-to-observed treatment density ratios. Because the intervention itself is defined using the propensity score, however, that covariate alone does not describe every contribution to the efficient influence function.
Definition and Treatment Density
Let \(H_t\) denote the measured history immediately before the binary treatment \(A_t\in\{0,1\}\) at time \(t\), and let \(g_t(a\mid h)=P(A_t=a\mid H_t=h)\) be the observed treatment mechanism. For a fixed incremental factor \(\delta>0\), define the intervention distribution by multiplying the odds of treatment by \(\delta\).
The intervention odds equal \(\delta\) times the observed odds whenever \(0<g_t(h)<1\): \[ \frac{q_{\delta,t}(1\mid h)}{q_{\delta,t}(0\mid h)} =\delta\frac{g_t(h)}{1-g_t(h)}. \] Thus \(\delta=1\) leaves the mechanism unchanged, \(\delta>1\) increases treatment odds, and \(0<\delta<1\) decreases them. For finite positive \(\delta\), the intervention probabilities sum to one because their numerators sum to \(d_{\delta,t}(h)\).
Worked Example: A Single Decision and Its Density Ratio
Suppose \(P(A=1\mid H=h)=g(h)=0.20\), and choose \(\delta=3\). The denominator is \(d=1+(3-1)(0.20)=1.40\). Therefore \(q_3(1\mid h)=3(0.20)/1.40=3/7\) and \(q_3(0\mid h)=0.80/1.40=4/7\).
The observed odds are \(0.20/0.80=1/4\). The intervention odds are \((3/7)/(4/7)=3/4\), exactly three times the observed odds. The intervention-to-observed density ratio for a treated person is \((3/7)/(1/5)=15/7\). For an untreated person it is \((4/7)/(4/5)=5/7\). Averaging these ratios under the observed treatment distribution gives \((1/5)(15/7)+(4/5)(5/7)=3/7+4/7=1\), as a normalized density ratio should.
Longitudinal Density Ratios and Targeting Covariates
For decisions at times \(t=0,\ldots,K\), the intervention distribution for the treatment sequence, conditional on the evolving histories, is the product of its conditional probabilities. Relative to the observed treatment sequence, its density ratio is therefore the product of the corresponding one-time ratios.
This covariate is the analogue of the clever covariate for a stochastic intervention. In a sequential targeting step, it weights the residual at that time according to how much more or less likely the observed treatment path is under the intervention than under the observed mechanism. It is not the probability of the observed treatment path under the intervention: it is the ratio of that probability to its observed-mechanism counterpart.
Worked Example: Two Treatment Decisions
Consider two decisions and let \(\delta=2\). At time 0, suppose \(g_0(H_0)=0.20\) and the person receives \(A_0=1\). Then \(d_0=1+(2-1)(0.20)=1.20\), so the time-0 ratio is \(2/1.20=5/3\).
At time 1, suppose the new history has \(g_1(H_1)=0.50\) and the person receives \(A_1=0\). Then \(d_1=1+(2-1)(0.50)=1.50\), so the time-1 ratio is \(1/1.50=2/3\). The cumulative targeting covariate through time 1 is \(C_{\delta,1}=(5/3)(2/3)=10/9\).
For comparison, the intervention probabilities are \(q_2(1\mid H_0)=0.40/1.20=1/3\) and \(q_2(0\mid H_1)=0.50/1.50=1/3\). The observed probabilities of the realized sequence, conditional on its histories, are \(0.20\) and \(0.50\), whose product is \(0.10\). The intervention sequence probability is \((1/3)(1/3)=1/9\). Their ratio is \((1/9)/(0.10)=10/9\), agreeing with the cumulative covariate.
Proof. By the chain rule for conditional probabilities, the observed conditional probability of a treatment sequence \(a_0,\ldots,a_K\), given the histories generated along that sequence, is \(\prod_{t=0}^{K}g_t(a_t\mid h_t)\). Under the intervention it is \(\prod_{t=0}^{K}q_{\delta,t}(a_t\mid h_t)\). Where the observed probability is positive, dividing these expressions gives \[ \frac{\prod_{t=0}^{K}q_{\delta,t}(a_t\mid h_t)} {\prod_{t=0}^{K}g_t(a_t\mid h_t)} =\prod_{t=0}^{K}\frac{q_{\delta,t}(a_t\mid h_t)}{g_t(a_t\mid h_t)}. \] The partial product through \(t\) is the corresponding density ratio for the sequence up to time \(t\). \(\square\)
Worked Example: Three Decisions and a Bounded Ratio
Let \(\delta=0.5\). Suppose the three observed treatment probabilities are \(g_0=0.40\), \(g_1=0.10\), and \(g_2=0.70\), and the realized treatments are \((A_0,A_1,A_2)=(0,1,1)\). The denominators are \(d_0=1-0.5(0.40)=0.80\), \(d_1=1-0.5(0.10)=0.95\), and \(d_2=1-0.5(0.70)=0.65\).
The three ratio factors are \(1/0.80=1.25\), \(0.5/0.95=10/19\), and \(0.5/0.65=10/13\). Hence \(C_{\delta,2}=1.25(10/19)(10/13)=1250/2470\), approximately \(0.506\). The intervention has reduced treatment odds, so this particular sequence, which includes two treated decisions, receives a cumulative ratio below one.
For fixed \(\delta\), each binary factor is bounded: it is \(1/d\) for \(A_t=0\) and \(\delta/d\) for \(A_t=1\), where \(d=1+(\delta-1)g\). Since \(g\in[0,1]\), these factors are bounded above by \(\max(1,\delta,1/\delta)\). With \(K+1\) fixed decisions, the cumulative ratio is bounded above by \(\max(1,\delta,1/\delta)^{K+1}\).
Why the Propensity-Score Dependence Matters
For a fixed stochastic intervention, the cumulative ratio weights the outcome-regression residuals in the usual longitudinal targeting construction. An incremental intervention differs in one important respect: \(q_{\delta,t}\) is defined using \(g_t\). Changing the treatment mechanism therefore changes both the observed-data distribution and the intervention being evaluated. Consequently, the cumulative ratio is the outcome-residual targeting covariate, but it is not by itself the entire influence-function contribution for this parameter.
Derivation. Fix a history \(h\) and write \(g=g(h)\). The intervention mean conditional on \(h\) is \(V(h)=q_\delta(1\mid h)Q(1,h)+q_\delta(0\mid h)Q(0,h)\). Differentiating the intervention probabilities with respect to \(g\) gives \[ \frac{\partial q_\delta(1\mid h)}{\partial g}=\frac{\delta}{d(h)^2}, \qquad \frac{\partial q_\delta(0\mid h)}{\partial g}=-\frac{\delta}{d(h)^2}. \] Thus the derivative of \(V(h)\) from changing \(g\) is \(\{\delta/d(h)^2\}\{Q(1,h)-Q(0,h)\}\dot g(h)\). A score for the Bernoulli treatment mechanism along this change is \(\{A-g(h)\}\dot g(h)/\{g(h)(1-g(h))\}\). The conditional covariance of \(\{A-g(h)\}\) with that score is \(\dot g(h)\), because \(E[(A-g)^2\mid H=h]=g(h)\{1-g(h)\}\). Multiplying \(A-g(h)\) by \(\{\delta/d(h)^2\}\{Q(1,h)-Q(0,h)\}\) therefore reproduces the required derivative. The outcome-residual term follows from the conditional mean derivative: its multiplier is \(q_\delta(A\mid H)/g(A\mid H)\). The centered baseline term accounts for variation in the distribution of \(H\). \(\square\)
In a longitudinal problem, the same distinction applies at each decision: the outcome-residual covariate is the cumulative ratio through that decision, while dependence of the intervention on each \(g_t\) also contributes to the derivative of the target. The relevant contrast at time \(t\) is between the downstream intervention values after setting \(A_t=1\) and \(A_t=0\); histories at time \(t\) are distributed under the earlier intervention. A targeting procedure must account for the intervention’s dependence on the treatment mechanisms rather than treating the estimated \(q_{\delta,t}\) as fixed without justification.
Worked Example: The Extra Single-Time Contribution
Continue with \(g=0.20\) and \(\delta=3\), so \(d=1.40\). Suppose \(Q(1,h)=0.60\) and \(Q(0,h)=0.30\). The intervention mean is \(V(h)=(3/7)(0.60)+(4/7)(0.30)=1.8/7+1.2/7=3/7\).
The coefficient on \(A-g(h)\) in the treatment-mechanism contribution is \(\delta\{Q(1,h)-Q(0,h)\}/d^2=3(0.30)/(1.40^2)=0.90/1.96=45/98\), approximately \(0.4592\). If \(A=1\), this contribution is \((45/98)(1-0.20)=(45/98)(0.80)=18/49\), approximately \(0.3673\). If \(A=0\), it is \((45/98)(-0.20)=-9/98\), approximately \(-0.0918\).
Its conditional expectation is zero, since \(0.20(18/49)+0.80(-9/98)=18/245-36/490=0\). The calculation shows why the contribution is centered at the observed propensity score and why it does not replace the outcome-residual term.
Assumptions, Diagnostics, and Limitations
For a causal interpretation, the usual longitudinal requirements remain relevant: consistency of observed outcomes with the treatment actually received, sequential exchangeability given the measured histories, and adequate measurement of confounders at each decision. The treatment and outcome histories must be ordered so that \(H_t\) is available before \(A_t\). For standard regularity arguments, one commonly assumes a fixed finite number of decisions, suitable integrability of the outcome, and treatment probabilities bounded away from zero and one. The incremental intervention is defined for any finite \(\delta>0\).
Incremental interventions preserve the observed treatment support: if \(g_t(1\mid h)=0\), then \(q_{\delta,t}(1\mid h)=0\); if \(g_t(1\mid h)=1\), then \(q_{\delta,t}(0\mid h)=0\). This can avoid the unsupported treatment assignments required by some deterministic regimes. But it does not create information about treatment effects for histories at which one treatment is never observed. Nor does boundedness of the density-ratio factors guarantee precise estimation: long products, noisy nuisance estimates, and rare histories can still cause unstable finite-sample behavior.
Useful diagnostics include inspecting the fitted \(g_t(H_t)\) distributions at every decision, checking that computed \(q_{\delta,t}(1\mid H_t)\) lies in \([0,1]\), and checking the mean and upper tail of each cumulative ratio. Verify that the ratios are computed as \(q/g\), not \(g/q\), and that the product includes only decisions up to the time of the residual being targeted. Also check whether the fitted mechanism and the intervention used in targeting are updated consistently.
A common mistake is to interpret the incremental factor as an addition to the treatment probability, or to use \(\delta g\) as the intervention probability. Neither operation preserves the intended odds tilt or necessarily gives a probability. The correction is to normalize the tilted probabilities by \(1+(\delta-1)g\). Another mistake is to claim that the cumulative ratio alone gives the full efficient influence function. It gives the outcome-residual multiplier; because \(q_\delta\) depends on \(g\), there is also a treatment-mechanism derivative contribution.
Compared with deterministic treatment regimes, incremental interventions make smaller, graded shifts and retain the observed mechanism’s support. Compared with a stochastic intervention specified independently of \(g\), they adapt to the observed propensity score, which is their defining feature and also adds a derivative term to inference and targeting. They do not automatically solve confounding, model misspecification, or sparse-data problems; the factor \(\delta\) defines a different causal target, rather than a correction for bias.
Check Your Understanding
Try each question, then check the solutions below.
- If \(g(h)=0.60\) and \(\delta=0.5\), calculate both intervention probabilities and verify the odds relation.
- For the same values, calculate \(q_\delta(1\mid h)/g(1\mid h)\) and \(q_\delta(0\mid h)/g(0\mid h)\).
- With two decisions and ratios \(3/2\) and \(2/3\), calculate the cumulative targeting covariate. What does it represent?
- Why is the longitudinal targeting covariate not, by itself, the full influence-function contribution for an incremental intervention?
Solutions. 1. The denominator is \(d=1+(-0.5)(0.60)=0.70\). Thus \(q(1)=0.5(0.60)/0.70=3/7\) and \(q(0)=0.40/0.70=4/7\). Their odds are \(3/4\); the observed odds are \(0.60/0.40=3/2\), and half of these is \(3/4\). 2. The treated ratio is \((3/7)/(3/5)=5/7\), and the untreated ratio is \((4/7)/(2/5)=10/7\). Their observed-treatment average is \((3/5)(5/7)+(2/5)(10/7)=1\). 3. The product is \((3/2)(2/3)=1\). It is the intervention-to-observed density ratio for the treatment sequence through the second decision. 4. The intervention probabilities depend on the propensity score, so changing that mechanism also changes the target intervention. The influence function therefore includes a treatment-mechanism derivative in addition to the outcome-residual term weighted by the cumulative ratio.