Why Residual Imbalance Matters Through the Outcome
A misspecified treatment model can leave covariates imbalanced after weighting, but imbalance alone does not determine the direction or size of an effect error. The relevant question is how the weighting changes the covariate distribution in relation to the conditional outcome mean. If the outcome mean is constant across the residual imbalance, the imbalance does not shift the weighted mean. If it varies with the covariates, the same imbalance can move the estimated effect.
The connection can be stated exactly for a point-treatment setting. The result is a population calculation: it describes the limit targeted by normalized inverse-probability weighting under the stated assumptions, not every finite-sample realization.
Setup and Exact Bias Identity
Let \(A\in\{0,1\}\) be treatment, \(L\) a vector of baseline covariates, and \(Y\) the observed outcome. Write \(Y^a\) for the potential outcome under treatment \(a\), and define the true propensity probabilities \(e_a(l)=P(A=a\mid L=l)\). Let \(q_a(l)\) be the fitted treatment probabilities used to construct weights, with \(q_0(l)+q_1(l)=1\) and \(q_a(l)>0\) wherever \(e_a(l)>0\). The outcome regression is \(m_a(l)=E[Y\mid A=a,L=l]\).
For identification of the average treatment effect, assume consistency, conditional exchangeability \(Y^a\perp A\mid L\), and positivity \(e_a(L)>0\) almost surely for \(a=0,1\). Also assume the relevant expectations are finite. These conditions give \(m_a(l)=E[Y^a\mid L=l]\), so the target means are \(\mu_a=E[m_a(L)]=E[Y^a]\), and the target effect is \(\tau=\mu_1-\mu_0\).
The population counterpart of the normalized weighted mean in arm \(a\) is
Define the ratio \(r_a(L)=e_a(L)/q_a(L)\). Require \(0<E[r_a(L)]<\infty\) and \(E[r_a(L)|m_a(L)|]<\infty\). Conditional expectation gives the following identity.
Proof. Given \(L\), \(P(A=a\mid L)=e_a(L)\), and \(E[I(A=a)Y\mid L]=e_a(L)m_a(L)\). Therefore
Their ratio is the stated expression for \(\mu_a^{\,q}\). Since \(E[r_am_a]=E[r_a]E[m_a]+\operatorname{Cov}(r_a,m_a)\), subtracting \(\mu_a=E[m_a]\) gives the arm-specific bias identity. Subtracting the two target arm means from the two weighted arm means gives the effect identity. \(\square\)
The same calculation describes covariate imbalance. If \(F\) is the population distribution of \(L\), the weighted distribution in arm \(a\) has density relative to \(F\)
Thus \(\mu_a^{\,q}=\int m_a(l)\,dG_a^q(l)\), while \(\mu_a=\int m_a(l)\,dF(l)\). The weighted mean shifts when this change in distribution places different weight on covariate values with different outcome means. This also shows why a misspecified treatment model does not guarantee imbalance: if \(r_a(l)\) is constant over \(l\), then \(G_a^q=F\).
Worked Example: A Binary Confounder and a Constant Fitted Probability
Let \(L\) be binary with \(P(L=1)=P(L=0)=0.5\). Suppose \(e_1(1)=0.8\), \(e_1(0)=0.2\), but the fitted model uses \(q_1(l)=q_0(l)=0.5\) at both values of \(L\). Let \(m_1(l)=2+3l\) and \(m_0(l)=1+l\).
The target arm means are \(\mu_1=0.5(5)+0.5(2)=3.5\) and \(\mu_0=0.5(2)+0.5(1)=1.5\). Hence \(\tau=2\). For treated participants, \(r_1(1)=0.8/0.5=1.6\) and \(r_1(0)=0.2/0.5=0.4\), so \(E[r_1]=0.5(1.6)+0.5(0.4)=1\). The weighted proportion with \(L=1\) is \(0.5(1.6)/1=0.8\). The weighted treated mean is therefore \(0.8(5)+0.2(2)=4.4\), a bias of \(0.9\).
For controls, \(r_0(1)=0.2/0.5=0.4\) and \(r_0(0)=0.8/0.5=1.6\), again with \(E[r_0]=1\). The weighted proportion with \(L=1\) is \(0.5(0.4)=0.2\). The weighted control mean is \(0.2(2)+0.8(1)=1.2\), a bias of \(-0.3\). The weighted effect is \(4.4-1.2=3.2\), so its bias is \(3.2-2=1.2\), equal to \(0.9-(-0.3)\).
The covariance terms verify the same result. For treatment, the deviations of \(r_1\) from its mean are \(0.6\) at \(L=1\) and \(-0.6\) at \(L=0\); the deviations of \(m_1\) from its mean are \(1.5\) and \(-1.5\). Thus \(\operatorname{Cov}(r_1,m_1)=0.5(0.6)(1.5)+0.5(-0.6)(-1.5)=0.9\). For controls, the corresponding deviations are \(-0.6,0.6\) for \(r_0\) and \(0.5,-0.5\) for \(m_0\), giving \(\operatorname{Cov}(r_0,m_0)=-0.3\). With both denominators equal to 1, the effect bias is \(0.9-(-0.3)=1.2\).
Continuous Covariate and a Bias Bound
The identity applies to continuous covariates as well. It can also be bounded in terms of how much the ratio varies and how strongly it covaries with the outcome regression.
Derivation. The bias identity and the Cauchy–Schwarz inequality give \[ |\mu_a^{\,q}-\mu_a| = \frac{|\operatorname{Cov}(r_a,m_a)|}{E[r_a]} \leq \frac{\sqrt{\operatorname{Var}(r_a)\operatorname{Var}(m_a)}}{E[r_a]}. \] All quantities are finite by assumption, and the positive denominator permits division. \(\square\) This bound is not generally sharp: it controls the magnitude but does not reveal the sign, and the two arm-specific biases can partially cancel or reinforce one another in the effect contrast.
Worked Example: Continuous Confounding With Linear Outcome Regressions
Let \(L\sim\operatorname{Uniform}(-1,1)\), so \(E[L]=0\) and \(\operatorname{Var}(L)=1/3\). Suppose \(e_1(L)=0.5+0.2L\), which lies between \(0.3\) and \(0.7\), while the fitted probabilities are \(q_1(L)=q_0(L)=0.5\). Take \(m_1(L)=2+3L\) and \(m_0(L)=1+L\). The target means are \(\mu_1=2\) and \(\mu_0=1\), so \(\tau=1\).
Here \(r_1(L)=1+0.4L\) and \(r_0(L)=1-0.4L\), with \(E[r_1]=E[r_0]=1\). The treated weighted covariate mean is \(E[r_1L]=E[L]+0.4E[L^2]=0.4/3=2/15\). Thus \(\mu_1^{\,q}=2+3(2/15)=12/5\), and the treated-arm bias is \(12/5-2=2/5\).
For controls, the weighted covariate mean is \(E[r_0L]=E[L]-0.4E[L^2]=-2/15\). Hence \(\mu_0^{\,q}=1-2/15=13/15\), with bias \(13/15-1=-2/15\). The weighted effect is \(12/5-13/15=23/15\), and its bias is \(23/15-1=8/15\). Directly, the covariance terms are \(0.4(3)\operatorname{Var}(L)=2/5\) for treatment and \(-0.4(1)\operatorname{Var}(L)=-2/15\) for controls, yielding the same effect bias \(2/5-(-2/15)=8/15\).
Worked Example: Misspecification Without Residual Imbalance
Let \(P(L=1)=0.3\), and suppose the true treatment probability is \(e_1(L)=0.5\) for both values of \(L\). Use the incorrect fitted probabilities \(q_1(L)=0.25\) and \(q_0(L)=0.75\). Then \(r_1(L)=0.5/0.25=2\) and \(r_0(L)=0.5/0.75=2/3\), each constant across \(L\). Consequently, normalizing the weights leaves the distribution of \(L\) unchanged in either arm.
For example, let \(m_1(L)=4+2L\) and \(m_0(L)=1-L\). The target means are \(\mu_1=4+2(0.3)=4.6\) and \(\mu_0=1-0.3=0.7\), giving \(\tau=3.9\). Because the ratios are constant, \(\mu_1^{\,q}=4.6\) and \(\mu_0^{\,q}=0.7\) as well, so \(\tau^q=3.9\). The fitted probabilities are wrong, but neither the normalized covariate distributions nor these weighted means are shifted. This is why model misspecification alone is not a proof of effect bias.
Longitudinal Interpretation
For a longitudinal regime \(\bar a=(a_0,\ldots,a_K)\), let \(H_t\) denote the observed history immediately before treatment at time \(t\). Write \(e_t(a_t\mid H_t)\) for the true conditional treatment probability and \(q_t(a_t\mid H_t)\) for its fitted counterpart. Under the usual sequential identification conditions, the ratio comparing the true and fitted treatment factors along a regime is
In the population weighted regime sample, the regime-history distribution is reweighted by \(R\) and then normalized. Accordingly, a history-specific outcome mean that is higher in histories receiving relatively large \(R\) contributes more to the weighted mean. This is the longitudinal counterpart of the point-treatment covariance relationship. It does not imply that fitted probabilities being wrong must change the normalized history distribution: if \(R\) is constant over histories under the regime distribution, normalization removes that constant. Nor does history imbalance alone determine outcome bias; its relationship to the conditional outcome mean matters.
This interpretation requires the regime to be supported by the observed treatment process, valid sequential exchangeability and consistency, and finite weighted expectations. With many treatment times, small ratio errors can multiply, and histories with very large or very small product ratios can make estimates unstable. The point-treatment identity is a guide to the mechanism, not a substitute for checking longitudinal support and weight behavior.
Diagnostics, Alternatives, and Limitations
Post-weighting balance diagnostics examine whether the weighted covariate distribution resembles the intended target. They are useful for detecting some forms of residual imbalance, but they do not directly measure the covariance between the weight ratio and the outcome regression. Balance on a few means can miss distributional differences that matter for nonlinear outcome relationships; conversely, measurable imbalance may have little effect when the outcome mean changes little across the imbalanced values. Outcome-blind diagnostics therefore cannot establish that weighted effect bias is absent.
The identity also clarifies the roles of alternative approaches. Inverse-probability weighting relies on the treatment model to construct a weighted population and is vulnerable when its ratios remain associated with outcome-relevant covariates. Outcome-regression standardization instead relies on modeling the conditional outcome mean and averaging predictions over the target covariate distribution; its vulnerability is outcome-model misspecification rather than residual imbalance induced by treatment weights. These approaches make different modeling commitments, and neither is automatically protected from all misspecification.
A frequent mistake is to argue, “the weighted covariates are imbalanced, so the treatment effect must be biased by that amount.” This is wrong: an imbalance is a change in a covariate distribution, not an outcome contrast. The correction is to consider how the change acts on \(m_a(L)\), arm by arm, as the covariance identity does. A second mistake is to infer that an incorrect treatment model necessarily creates imbalance. Constant true-to-fitted ratios within an arm produce no normalized distribution shift, as the third example demonstrates.
Check Your Understanding
Answer each question, then compare with the solutions.
- In the point-treatment setup, what condition on \(r_a(L)\) ensures that normalized weighting leaves the distribution of \(L\) unchanged in arm \(a\)?
- If \(r_a(L)\) varies but \(m_a(L)\) is constant, what is the arm-specific weighted-mean bias?
- Why can a weighted covariate mean be unchanged even when the full weighted covariate distribution differs from its target?
- In a longitudinal regime, what quantity plays the role of the point-treatment ratio \(r_a(L)\)? Why does its being constant over regime histories matter?
Solutions. (1) A constant \(r_a(L)\) almost surely; dividing it by its expectation gives a density ratio of 1. (2) Zero, because \(\operatorname{Cov}(r_a,m_a)=0\) when \(m_a\) is constant. (3) Different distributions can have the same mean; for example, a nonlinear outcome regression may respond to differences in spread or tail probabilities that a mean-balance check does not detect. (4) The product \(R(\bar H)=\prod_t e_t(a_t\mid H_t)/q_t(a_t\mid H_t)\). If it is constant, normalization removes it, leaving the regime-history distribution unchanged.