Tutorials › Marginal Structural Models › Fixed Truncation Thresholds and Asymptotic Bias

Weight Stabilization and Truncation · Tutorial 151 of 501

Fixed Truncation Thresholds and Asymptotic Bias

Derive the probability limit of a fixed-threshold clipped IPW estimator and see how clipping can permanently change the causal estimand.

Advanced 12 min read

What You'll Learn

  • Derive the probability limit of a fixed-cap Hájek estimator for a causal mean.
  • Express its limiting target as a reweighted average of conditional potential-outcome means.
  • Characterize asymptotic bias using a covariance identity.
  • Work out examples where clipping creates bias, creates no bias, or reverses its direction.
  • Distinguish fixed-threshold asymptotic bias from finite-sample variance reduction.
  • Identify practical diagnostics and assumptions that matter when interpreting clipped estimates.

Why a Fixed Threshold Has an Asymptotic Consequence

Clipping large inverse-probability weights is often used to reduce the influence of observations with very small treatment probabilities. That can reduce variability, but it also changes the weights applied to the data. If the clipping threshold stays fixed as the sample grows, the change does not disappear: the estimator can converge to a different population quantity from the intended causal parameter.

The key distinction is between finite-sample error and <strong asymptotic bias. A noisy estimate may become more precise with a larger sample, but a fixed clipping rule can leave the estimator centered on the wrong limiting value. We derive that value for a simple binary-treatment mean estimator, then use it to identify when the discrepancy occurs.

Setup and Fixed-Cap Estimator

Definition: Let \(A\in\{0,1\}\) be treatment, \(X\) baseline covariates, and \(Y\) the observed outcome. Let \(Y^a\) be the potential outcome under treatment \(a\), and define \(p_a(X)=P(A=a\mid X)\). For a fixed cap \(c>0\), the clipped inverse-probability weight for arm \(a\) is \(w_{a,c}(X)=\min\{1/p_a(X),c\}\). The clipped Hájek estimator of the mean potential outcome in arm \(a\) is \(\hat{\mu}_{a,c}=\{\sum_{i=1}^n I(A_i=a)w_{a,c}(X_i)Y_i\}/\{\sum_{i=1}^n I(A_i=a)w_{a,c}(X_i)\}\).

The target causal mean is \(\mu_a=E[Y^a]\); the average treatment effect is \(\tau=E[Y^1-Y^0]=\mu_1-\mu_0\). To identify these means from observed data, assume consistency (\(Y=Y^A\)), conditional exchangeability (\(Y^a\) is independent of \(A\) given \(X\), for each \(a\)), and positivity (\(p_a(X)>0\) almost surely for each arm). Assume observations are independent and identically distributed, \(E[|Y|]<\infty\), and the propensity probabilities in the estimator are the true probabilities. These conditions suffice for the probability limits below; estimated weights require additional analysis.

Write \(m_a(X)=E[Y^a\mid X]\). Under consistency and exchangeability, the conditional mean of the observed outcome among people with \(A=a\), given \(X\), is \(m_a(X)\). The clipped weight is bounded by \(c\), so the weighted outcome has finite absolute expectation under the stated integrability condition.

Probability Limit of the Clipped Mean

Theorem: Under the assumptions above, for fixed \(c>0\), \(\hat{\mu}_{a,c}\xrightarrow{p}\mu_{a,c}\), where \(\mu_{a,c}=\dfrac{E[q_{a,c}(X)m_a(X)]}{E[q_{a,c}(X)]}\) and \(q_{a,c}(X)=\min\{1,c\,p_a(X)\}\). The denominator is strictly positive. In general, \(\mu_{a,c}\) need not equal the causal mean \(\mu_a\).

Proof. For one observation, the denominator contribution has expectation \(E[I(A=a)w_{a,c}(X)]=E[p_a(X)w_{a,c}(X)]=E[\min\{1,c p_a(X)\}]\). This is positive because positivity gives \(p_a(X)>0\) almost surely and \(c>0\). The contribution is bounded by \(c\), so the law of large numbers applies.

For the numerator, consistency and conditional exchangeability give \(E[Y\mid A=a,X]=m_a(X)\). Iterated expectation therefore yields \(E[I(A=a)w_{a,c}(X)Y]=E[p_a(X)w_{a,c}(X)m_a(X)]=E[q_{a,c}(X)m_a(X)]\). The absolute expectation is finite: \(w_{a,c}\leq c\) and \(E[|Y|]<\infty\). Thus the law of large numbers applies to the numerator as well. Since the denominator average converges in probability to a positive constant, the ratio converges to the ratio of these expectations by the continuous mapping theorem. This is the stated limit. \(\square\)

The factor \(q_{a,c}(X)\) shows what clipping does to the covariate distribution represented in the weighted arm. Without clipping, \(w=1/p_a(X)\), so \(p_aw=1\) and all covariate values receive equal weight in the target mean. With a fixed cap, values with \(p_a(X)<1/c\) instead receive factor \(c p_a(X)<1\). They remain represented under strict positivity, but are downweighted relative to the original target population.

Worked Example: Clipping Shifts the Treated Mean Upward

Let \(X\) be binary with \(P(X=0)=P(X=1)=1/2\), and suppose \(p_1(0)=0.1\) and \(p_1(1)=0.5\). Let the treated potential outcome be \(Y^1=0\) when \(X=0\) and \(Y^1=10\) when \(X=1\). The target treated mean is \(\mu_1=(0+10)/2=5\). Use cap \(c=3\).

At \(X=0\), the raw treated weight is \(1/0.1=10\), so the clipped weight is 3. At \(X=1\), the raw weight is \(1/0.5=2\), so the clipped weight remains 2. The corresponding factors in the limit are \(q_{1,c}(0)=0.1(3)=0.3\) and \(q_{1,c}(1)=0.5(2)=1\). Hence \(\mu_{1,c}=\{(1/2)(0.3)(0)+(1/2)(1)(10)\}/\{(1/2)(0.3)+(1/2)(1)\}=5/0.65=100/13\approx7.69\).

The fixed-cap limit exceeds the target by \(100/13-5=35/13\approx2.69\). The less common treated group at \(X=0\), whose treated potential outcome is lower, has been relatively downweighted. More observations estimate this altered average more precisely; they do not restore the original mean.

When Does the Limiting Value Equal the Target?

Result: If \(E[|m_a(X)|]<\infty\), then \(\mu_{a,c}-\mu_a=\dfrac{\operatorname{Cov}(q_{a,c}(X),m_a(X))}{E[q_{a,c}(X)]}\). Consequently, the clipped limit equals the target exactly when this covariance is zero. In particular, it equals the target if \(m_a(X)\) is constant almost surely, or if clipping does not change the weights almost surely.

Derivation. Since \(\mu_a=E[m_a(X)]\), subtract this quantity from the ratio defining \(\mu_{a,c}\): \[ \mu_{a,c}-\mu_a =\frac{E[q_{a,c}m_a]-E[q_{a,c}]E[m_a]}{E[q_{a,c}]} =\frac{\operatorname{Cov}(q_{a,c},m_a)}{E[q_{a,c}]}. \] The denominator is positive, so the difference is zero exactly when the covariance is zero. If \(m_a\) is constant, its covariance with any integrable \(q_{a,c}\) is zero. If no weights are changed, then \(q_{a,c}=1\) almost surely, which also makes the covariance zero. \(\square\)

For the average treatment effect estimated as \(\hat{\mu}_{1,c}-\hat{\mu}_{0,c}\), the limiting discrepancy is the difference between the two arm-specific discrepancies: \[ (\mu_{1,c}-\mu_{0,c})-\tau =\frac{\operatorname{Cov}(q_{1,c},m_1)}{E[q_{1,c}]} -\frac{\operatorname{Cov}(q_{0,c},m_0)}{E[q_{0,c}]}. \] Thus the direction of bias in a contrast cannot be inferred from clipping alone. Each arm can be shifted differently, and the shifts can reinforce or offset one another.

Worked Example: Clipping Without Asymptotic Bias

Keep the two covariate groups and treatment probabilities from the previous example, and retain \(c=3\). Now let \(Y^1=4\) for both values of \(X\). The target mean is \(\mu_1=4\). The clipped denominator expectation is \((1/2)(0.3)+(1/2)(1)=0.65\).

The clipped numerator expectation is \((1/2)(0.3)(4)+(1/2)(1)(4)=0.6+2=2.6\). Therefore \(\mu_{1,c}=2.6/0.65=4\), exactly the target. In the covariance formula, \(m_1(X)=4\) is constant, so \(\operatorname{Cov}(q_{1,c},m_1)=0\). Clipping changes the represented covariate distribution here, but that change does not affect the mean because the conditional outcome mean does not vary with \(X\).

Worked Example: Clipping Can Shift the Mean Downward

Use the same covariate probabilities, treatment probabilities, and cap \(c=3\). Set \(Y^1=10\) at \(X=0\) and \(Y^1=0\) at \(X=1\). The target is again \(\mu_1=(10+0)/2=5\). The clipped factors remain \(0.3\) and 1, and the denominator is \(0.65\).

The numerator expectation is \((1/2)(0.3)(10)+(1/2)(1)(0)=1.5\), so \(\mu_{1,c}=1.5/0.65=30/13\approx2.31\). The bias is \(30/13-5=-35/13\approx-2.69\). In this case, the group that is downweighted has the higher treated potential-outcome mean, so the clipped limit moves downward instead. The same clipping rule produced upward bias in the first example and downward bias here.

Interpretation, Diagnostics, and Alternatives

A fixed cap bounds each individual weight, but bounded weights are not by themselves a guarantee of consistency for the original causal target. They ensure that the specified estimator has a well-defined limiting weighted average under the theorem’s conditions. Whether that average is the desired causal mean depends on the relationship between the clipping factors and the conditional outcome means.

Useful diagnostics include reporting the cap and the proportion of weights affected, checking how estimates change across scientifically plausible fixed caps, and examining whether clipping materially changes weighted covariate distributions. These checks reveal sensitivity; they do not prove that a chosen cap removes bias. A small proportion of clipped observations can still matter if those observations represent covariate strata with substantially different conditional outcomes.

The main alternative is the untruncated IPW estimator: under the identification assumptions and suitable integrability, its population weight factor is \(p_a(X)/p_a(X)=1\), recovering \(E[Y^a]\). Its weights may be highly variable when treatment probabilities are small. Outcome standardization instead estimates conditional outcome means and averages them over the target covariate distribution; that approach has its own modeling requirements. Truncation trades some protection from extreme weights for a potentially altered target, rather than making the original estimator uniformly better.

The formula above is for an arm-specific Hájek mean with a clipped unstabilized inverse-probability weight and known propensity. A stabilized weight with a covariate-dependent numerator, a longitudinal product weight, or a different estimating equation can have a different limiting functional. The relevant calculation is always to take the expectation of the actual clipped weighted estimating equation; the simple \(q_{a,c}\) formula should not be transferred mechanically to those settings.

Common Pitfall

Common mistake: “Clipping only affects a few extreme observations, so its effect vanishes as \(n\) grows.”

Correction: With a fixed cap, the same rule continues to alter weights in the population, no matter how large the sample becomes. The law of large numbers makes the empirical clipped ratio converge to its clipped population ratio, \(\mu_{a,c}\), not automatically to \(\mu_a\). Equality requires conditions such as zero covariance between the clipping factor and the conditional potential-outcome mean.

Check Your Understanding

Answer each question, then compare with the solutions.

  1. For the fixed-cap treated mean, what covariate factor determines the limiting target after clipping?
  2. Suppose \(m_1(X)=6\) almost surely and clipping changes some weights. What is the clipped limiting treated mean, and why?
  3. If \(\operatorname{Cov}(q_{1,c},m_1)<0\), what is the direction of bias in the clipped treated mean?
  4. Why does the law of large numbers not establish consistency for the original causal mean in the fixed-cap example?

Solutions. (1) The factor is \(q_{1,c}(X)=\min\{1,c p_1(X)\}\), so the limit is its weighted average of \(m_1(X)\). (2) The limit is 6, because a constant conditional mean has zero covariance with the clipping factor. (3) The bias is negative because its denominator \(E[q_{1,c}]\) is positive. (4) The law of large numbers establishes convergence to the ratio of expectations for the clipped numerator and denominator. That ratio is \(\mu_{1,c}\), which equals the target only under additional conditions such as zero covariance.

Clinical Biostats