Tutorials › Marginal Structural Models › Variance Reduction from Bounded Weights

Weight Stabilization and Truncation · Tutorial 149 of 501

Variance Reduction from Bounded Weights

Learn why bounded weights limit contribution sizes, how to derive a variance bound, and under which assumptions clipping actually lowers the variance of a simple weighted mean.

Advanced 11 min read

What You'll Learn

  • Derive a finite-sample variance bound for a population-normalized weighted mean with bounded weights and bounded outcomes.
  • Separate a bound on individual contributions from a claim that an estimator’s variance must decrease.
  • Compare raw and clipped weights using an exact variance formula under independence and a mean-zero outcome.
  • Calculate the variance reduction in a two-point weight example.
  • Identify how outcome–weight association can undermine a simple variance comparison.
  • Distinguish variance control from bias control and from guarantees for a sample-normalized ratio estimator.

Why a Bound Can Control Variability

A large weight can make one observation dominate a weighted estimating equation. Bounding weights limits how large that observation’s weighted contribution can be, which can make variability more manageable. But “bounded weights reduce variance” is not a universal theorem: clipping changes the contributions, and whether variance falls depends on the estimator, the outcome, and the assumptions connecting weights to outcomes.

To make the claim precise, start with independent and identically distributed observations \((W_i,Y_i)\), where \(W_i\geq 0\) is a weight and \(Y_i\) is an outcome. For a chosen weight function \(a(W)\), define the population-normalized weighted mean

$$ \hat{\mu}_a=\frac{1}{n}\sum_{i=1}^n \frac{a(W_i)}{E[a(W)]}Y_i, $$

assuming \(0<E[a(W)]<\infty\). This form uses the population mean weight in its denominator, rather than estimating that denominator from the sample. It is useful for isolating the variability of weighted contributions. Its target is \(E[a(W)Y]/E[a(W)]\); it equals the ordinary outcome mean \(E[Y]\) under conditions such as independence of \(W\) and \(Y\).

Definition: Upper truncation at \(B>0\) replaces a nonnegative weight \(W\) with \(W_B=\min(W,B)\). The corresponding normalized weight is \(A_B=W_B/E[W_B]\), provided \(0<E[W_B]<\infty\). Thus \(E[A_B]=1\) and \(0\leq A_B\leq B/E[W_B]\).

The upper bound on \(A_B\) is a bound on the normalized weight, not necessarily \(B\) itself. The normalization matters: even if \(W_B\leq B\), dividing by a small \(E[W_B]\) can make \(A_B\) substantially larger than \(B\).

A Variance Bound for Bounded Contributions

The first result makes the contribution-size argument explicit. It does not require independence between weights and outcomes, but it does require a bounded outcome and a positive, finite normalizing constant.

Theorem: Suppose \((W_i,Y_i)\) are independent and identically distributed, \(W_i\geq0\), \(0<E[W_B]<\infty\), and \(|Y_i|\leq K\) almost surely for some finite \(K\). For \(\hat{\mu}_B=n^{-1}\sum_{i=1}^n W_{B,i}Y_i/E[W_B]\), \[ \operatorname{Var}(\hat{\mu}_B)\leq \frac{K^2B^2}{n\{E[W_B]\}^2}. \]

Proof. Let \(Z_i=W_{B,i}Y_i/E[W_B]\). Since \(0\leq W_{B,i}\leq B\) and \(|Y_i|\leq K\), we have \(|Z_i|\leq KB/E[W_B]\) almost surely. Independence gives \(\operatorname{Var}(\hat{\mu}_B)=\operatorname{Var}(Z_1)/n\). For any square-integrable \(Z_1\), \(\operatorname{Var}(Z_1)=E[Z_1^2]-E[Z_1]^2\leq E[Z_1^2]\). The almost-sure bound therefore implies \(E[Z_1^2]\leq K^2B^2/\{E[W_B]\}^2\), proving the stated inequality. \(\square\)

This is an absolute bound, not yet a comparison with the untruncated estimator. If outcomes are unbounded, bounded weights alone do not bound \(W_BY\): an extreme outcome can still produce an extreme contribution. A variance bound then needs additional moment assumptions, for example \(E[Y^2]<\infty\), and its value depends on the outcome distribution as well as the weight bound.

Worked Example: Calculating a Contribution Bound

Suppose \(W\) equals \(0.5\) with probability \(0.8\) and \(3\) with probability \(0.2\). Truncate at \(B=1\), and suppose \(|Y|\leq2\). First calculate the normalizer: \(E[W_B]=0.8(0.5)+0.2(1)=0.6\). The normalized truncated weights are \(0.5/0.6=5/6\) and \(1/0.6=5/3\). Each weighted contribution \(W_BY/E[W_B]\) has absolute value at most \(2(5/3)=10/3\).

The theorem consequently gives \(\operatorname{Var}(\hat{\mu}_B)\leq (1/n)(10/3)^2=100/(9n)\). The bound is valid even if \(W\) and \(Y\) are associated, provided the stated bounded-outcome and sampling conditions hold. It may be conservative: it uses the largest possible contribution rather than the actual distribution of contributions.

When Clipping Gives a Direct Variance Comparison

A clean comparison is possible in a special but instructive setting. Suppose \(W\geq0\), \(0<E[W]<\infty\), \(Y\) is independent of \(W\), and \(E[Y]=0\), \(E[Y^2]=\sigma^2<\infty\). Normalize either the raw or truncated weights by their population means. Independence and the zero outcome mean make both resulting estimators unbiased for zero, so their variances can be compared directly.

Theorem: Under those assumptions, for any nonnegative weight function \(a(W)\) with \(0<E[a(W)]<\infty\) and \(E[a(W)^2]<\infty\), \[ \operatorname{Var}(\hat{\mu}_a)=\frac{\sigma^2}{n}\frac{E[a(W)^2]}{\{E[a(W)]\}^2}. \] Thus upper truncation has lower variance than raw weights precisely when \[ \frac{E[W_B^2]}{\{E[W_B]\}^2}<\frac{E[W^2]}{\{E[W]\}^2}. \]

Proof. Independence gives \(E[a(W)Y]=E[a(W)]E[Y]=0\), so each summand \(Z_i=a(W_i)Y_i/E[a(W)]\) has mean zero. Again by independence, \(E[Z_i^2]=E[a(W)^2]E[Y^2]/\{E[a(W)]\}^2\). The summands are independent, so the variance of their average is their common variance divided by \(n\). Substituting \(E[Y^2]=\sigma^2\) gives the formula. Applying it first with \(a(W)=W_B\) and then with \(a(W)=W\) yields the comparison. \(\square\)

The theorem gives a condition for the comparison; it does not assert that the inequality holds for every distribution and every cutoff. In this setting, the normalized second moment \(E[a(W)^2]/E[a(W)]^2\) determines the variance. A cutoff is helpful when it reduces that quantity. The conditions also explain why this is not a general causal-effect result: independence and a mean-zero outcome are strong simplifications, and clipping can alter a weighted target when outcomes depend on weights.

Worked Example: A Strict Variance Reduction

Let \(W=0.5\) with probability \(0.8\) and \(W=3\) with probability \(0.2\), as above. Let \(Y\) be independent of \(W\), with \(P(Y=1)=P(Y=-1)=0.5\). Then \(E[Y]=0\) and \(\sigma^2=1\). For the raw weights, \(E[W]=0.8(0.5)+0.2(3)=1\), and \(E[W^2]=0.8(0.25)+0.2(9)=0.2+1.8=2\). Therefore \(\operatorname{Var}(\hat{\mu}_{W})=2/n\).

Now truncate at \(B=1\). We have \(E[W_B]=0.6\) and \(E[W_B^2]=0.8(0.25)+0.2(1)=0.4\). The variance is \(\operatorname{Var}(\hat{\mu}_{W_B})=(1/n)(0.4/0.6^2)=(10/9)/n\). The ratio of the truncated to raw variance is \((10/9)/2=5/9\). In this example, truncation reduces variance by \(1-5/9=4/9\), or about \(44.4\%\). The calculation verifies the comparison using the theorem’s assumptions and the actual normalized second moments.

What the Comparison Does Not Guarantee

Variance depends on the distribution of the full contribution, not on the weight alone. For example, even a pointwise reduction in some weights need not make the products \(W_iY_i\) more similar to one another. Consider two equally likely paired values: \((W,Y)=(1,10)\) and \((10,1.1)\). The raw products are \(10\) and \(11\), with variance \(0.25\). Truncating at \(B=1\) changes them to \(10\) and \(1.1\), with variance \((10-1.1)^2/4=19.8025\). This calculation concerns the distribution of the contributions, not a comparison of two unbiased estimators for one unchanged target; indeed, the association between \(W\) and \(Y\) is exactly why the simple theorem above does not apply.

There are two further distinctions to keep in view. First, a variance reduction is not automatically an improvement in overall estimation error: truncation can change the weighted target and introduce bias. The consequences for mean squared error are a separate question. Second, the estimator defined here uses the population normalizer \(E[a(W)]\). A sample-normalized ratio, such as \(\sum_i a(W_i)Y_i/\sum_i a(W_i)\), has a random denominator. Its variance is not given exactly by the displayed formula for \(\hat{\mu}_a\).

In practice, inspect more than the maximum weight. Check the largest weighted outcome contributions, the distribution of normalized weights, and whether a small number of observations account for a large share of the weighted sum. Compare results across reasonable cutoffs and consider whether the observations most affected by clipping have systematically different outcomes. These checks can reveal instability or an outcome–weight association that defeats a simple variance argument.

Key takeaway: A bound limits weight size and, with bounded outcomes, gives an explicit bound on each weighted contribution. A variance reduction requires more: for the simple independent, mean-zero setting, it occurs exactly when truncation lowers the normalized second moment of the weights.

Worked Example: Diagnosing a Variance Claim

An analyst reports that clipping weights at \(B=2\) must reduce the variance of a weighted mean because every clipped weight is no larger than its raw value. Is that conclusion justified?

No. The pointwise inequality \(W_B\leq W\) does show that a nonnegative weighted contribution \(W_BY\) is no larger than \(WY\) when \(Y\geq0\). But a variance depends on squared deviations from the contribution’s mean, not just on whether individual values have decreased. Moreover, clipping changes the normalizing constant, and if outcomes are associated with weights, it can change the weighted target. The analyst should calculate or estimate the relevant variance under stated assumptions rather than infer a variance reduction from \(W_B\leq W\) alone.

Check Your Understanding

Answer each question, then compare with the solutions.

  1. If \(W_B\leq B\), what is the upper bound on the normalized weight \(W_B/E[W_B]\)?
  2. For independent mean-zero \(Y\) with variance \(\sigma^2\), write the variance of \(\hat{\mu}_a\).
  3. In the two-point example with \(W=0.5\) with probability \(0.8\), \(W=3\) with probability \(0.2\), and \(B=1\), calculate \(E[W_B]\) and \(E[W_B^2]\).
  4. Why does bounded weight alone not give a bounded weighted contribution when \(Y\) is unbounded?
  5. Does the variance formula for the population-normalized estimator apply exactly to a sample-normalized ratio? Why or why not?

Solutions. (1) The bound is \(B/E[W_B]\), assuming \(E[W_B]>0\). (2) It is \(\sigma^2 E[a(W)^2]/\{nE[a(W)]^2\}\). (3) \(E[W_B]=0.8(0.5)+0.2(1)=0.6\), and \(E[W_B^2]=0.8(0.25)+0.2(1)=0.4\). (4) The product can still be arbitrarily large if the outcome can be arbitrarily large. (5) No. A sample-normalized ratio has a random denominator, while the stated formula uses the fixed population normalizer \(E[a(W)]\).

Clinical Biostats