Estimating Heterogeneous Treatment Effects by Combining Weak Instruments and Observational Data

Accurately predicting conditional average treatment effects (CATEs) is crucial in personalized medicine and digital platform analytics. Since the treatments of interest often cannot be directly randomized, observational data is leveraged to learn CATEs, but this approach can incur significant bias from unobserved confounding. One strategy to overcome these limitations is to leverage instrumental variables (IVs) as latent quasi-experiments, such as randomized intent-to-treat assignments or randomized product recommendations. This approach, on the other hand, can suffer from low compliance, $\textit{i.e.}$, IV weakness. Some subgroups may even exhibit zero compliance, meaning we cannot instrument for their CATEs at all. In this paper, we develop a novel approach to combine IV and observational data to enable reliable CATE estimation in the presence of unobserved confounding in the observational data and low compliance in the IV data, including no compliance for some subgroups. We propose a two-stage framework that first learns $\textit{biased}$ CATEs from the observational data, and then applies a compliance-weighted correction using IV data, effectively leveraging IV strength variability across covariates. We characterize the convergence rates of our method and validate its effectiveness through a simulation study. Additionally, we demonstrate its utility with real data by analyzing the heterogeneous effects of 401(k) plan participation on wealth.

Paper

Similar papers

Peer review

Reviewer UxMs6/10 · confidence 3/52024-07-03

Summary

This paper introduces a robust two-stage framework leveraging observational and instrumental variable data to predict conditional average treatment effects (CATEs), addressing biases from unobserved confounders and low compliance.

Strengths

- This paper studies leveraging observational data and encouragement data with low compliance to predict conditional average treatment effects (CATEs) accurately. - This paper proposes a two-stage framework that first learns biased CATEs from observational data and subsequently applies a compliance-weighted correction using weak IVs. - This paper demonstrates its utility on real data 401(k) participation on financial wealth.

Weaknesses

**About the identification of CATE.** - The relevance assumption in Assumption 1 cannot ensure $\gamma(x)>0$. - The identification of CATE (Eq.(3) on line 121) would be violated when $P(A^E(1)|x) < P(A^E(0)|x)$ for some $x$. In such cases, even if Assumptions 1 and 2 hold, and $\mathbb{E}[A^E \mid Z^E=1, X^E=x] - \mathbb{E}[A^E \mid Z^E=0, X^E=x] < 0$,  $\gamma(x)$ would be zero and Eq.(3) with $\gamma(x)=0$ no longer holds. - In lines 39-41, the paper provides examples: “certain users on digital platforms may disregard recommendations either altogether or of certain undesired content, and on mobile health platforms certain participants may ignore recommendations (e.g., taking 250 steps per hour) due to time constraints or disinterest.” However, in Eq.(3), the paper does not consider the possibility that $\mathbb{E}[A^E \mid Z^E=1, X^E=x] - \mathbb{E}[A^E \mid Z^E=0, X^E=x] = 0$. This raises concerns about the soundness of the paper regarding the assumptions and theorems. - The paper seems to implicitly assume the monotonicity of instrumental variables. It is necessary for the authors to clearly state the assumptions required for the theory and provide a complete proof of identifiability (Eq.(3)). **I will re-evaluate the soundness of this paper according to the responses of the authors.** **About the related work.** In the presence of unmeasured confounding, there are numerous methods that use proxy variables to estimate heterogeneous treatment effects, including VAE-based methods [Causal effect inference with deep latent-variable models] and negative control methods [A selective review of negative control methods in epidemiology]. **Types:** Line 97 has an extra space. Line 103 should be revised to 'b(x), i.e.'.

Questions

- Does the method proposed in this paper require the additive noise assumption? - Do the observed variables X and latent U exhibit distribution shifts between Observational data and Experimental data? Are the distributions of X and U consistent across both datasets? - Is it feasible to set all Z=1 in the experimental data? Given that observational data effectively corresponds to the part where Z=0.

Rating

6

Confidence

3

Soundness

3

Presentation

3

Contribution

3

Limitations

N/A

Reviewer Gwhh6/10 · confidence 3/52024-07-04

Summary

The paper presents a novel method for estimating Conditional Average Treatment Effects (CATEs) by integrating weak instrumental variables (IV) and observational data. This method addresses the challenges of unobserved confounding and low compliance often encountered in causal inference studies. The proposed framework involves a two-stage process: initially, biased CATEs are estimated using observational data, and subsequently, a compliance-weighted correction is applied using IV data. This correction leverages the variability in IV strength across covariates to improve the accuracy of CATE estimation. The method's efficacy is validated through simulations and real-world applications, such as assessing the impact of 401(k) participation on financial wealth.

Strengths

The paper introduces a two-stage framework designed to estimate CATEs by effectively combining observational data with weak instrumental variables. The proposed method is adept at handling unobserved confounding in observational data and low compliance in IV data, including scenarios with zero compliance for some subgroups. The effectiveness of the method is demonstrated through extensive simulations and real-world applications, such as evaluating the effect of 401(k) participation on financial wealth.

Weaknesses

How to understand assumption 2? Is it correct about $Y^{E}(1) - Y^{E}(0) \perp\!\!\!\perp C \mid X^{E}$? What is IV datasets? How could there be IV in the experimental data? IV has a correlation with treatment, but according to Lemma 1, this instrumental propensity does not have any correlation with treatment. In the experimental data, was the treatment not randomly assigned? Experimental data are usually small samples or difficult to collect, which makes the proposed method difficult to apply.

Questions

See weakness

Rating

6

Confidence

3

Soundness

3

Presentation

3

Contribution

3

Limitations

Yes

Reviewer xHgN6/10 · confidence 4/52024-07-07

Summary

The authors propose an approach for conditional average treatment effect estimation using instrumental variables, extending existing work in IVA to settings with weak instruments (i.e., low treatment compliance in some population subgroups). In particular, the approach leverages a two-stage estimation setup: first, a biased CATE estimate is computed from the observational data, then corrected via reweighting based on compliance to obtain a final estimate. Overall, the paper provides good exposition and intuition for a novel approach, and backs it up with convincing theoretical analysis. The theoretical results are fairly insightful, but can be improved with a little more clarification (details below). The empirical results are not quite as well-motivated, and could also be improved with some clarification of the motivation and differentiation from past work. In particular, I felt that better justification of why these datasets/evaluations are the correct ones to test the proposed approach are needed for me to appreciate the paper’s contributions. Furthermore, I’m a little bit unsure if there are sufficient comparisons to baselines in inference under weak IV in the paper as-is (citations below).

Strengths

1. The exposition is clearly written and provides good intuition on IVA (though I am writing from the perspective of someone that has working knowledge of IVA already). 2. The proposed approach is intuitive and well-motivated, with good theoretical properties. Weak instruments are an inevitable problem in instrumental variable methods, so proposals to reduce their downstream impact are a salient area of research. 3. The paper itself is quite clearly written and mostly easy to understand.

Weaknesses

I'd love to see the following points addressed to clear up any misunderstandings on my part: 1. My biggest criticism is about the motivation and coverage of the empirical results. Coverage-wise, while I appreciate the comparison with a vanilla LATE estimator ($\tau^E$ if I followed correctly?), additional comparisons to baselines in learning with weak instruments would help strengthen the paper, such as the ones cited by the authors as closely related [1, 2] — it’s not clear to me the empirical lift provided by the proposed approach compared to these methods. For the rebuttal phase, in lieu of new results, perhaps a precise explanation of how the authors’ proposed approach shares similarities and differs with [1, 2] would be most helpful. 2. The analysis of the 401k dataset results is slightly imprecise — I’m not sure I buy that $\hat{\tau}(x)$ closely tracks with $\tau^E$; I don’t have a prior for what’s “close enough.” Could the authors clarify more precisely (1) the motivation behind experiments on the 401k dataset (beyond extension to real-world data), (2) and why the results show proof-of-concept for the proposed approach? 3. I have a couple of overarching concerns about the theoretical results. As a first-order comment, how do the generalization bounds of the proposed approach compare to similar bounds for IVA/what are the trade-offs compared to past bounds? As a second-order comment, if I know that $\tau^O$ will be biased (i.e., due to unobserved confounding), why is it imperative for estimation error w.r.t. $\tau^O$ to be low as well (maybe this is so that, given a good estimate of $\theta$ — we get a good estimate of $b(x)$, and therefore correct for the bias)? Nits (points that would improve the paper in my opinion, but are not urgent): 1. There are a few critical derivations where the clarity could be improved somewhat. It took me quite some time to follow how Eq. 4 was derived — while it’s probably okay to reserve most details to the Appendix, some intuition about which terms are being substituted where would be helpful, and making it clear in Eq. 4 that you’re taking the squared difference (example-wise) of the pseudo-outcome of Eq. 3 (as fitted on the intention-to-treat dataset) and a reweighted version of $b(x) + \tau^O(x)$ as fitted on the observational dataset (since it is equal to $\tau(x)$ by definition). 2. I notice that $\hat{\tau}^E$ and $\hat{\tau}^O$ are fitted on different subsets of the data — somewhat reminiscent of cross-fitting based estimators (e.g., [3]). Out of curiosity, is it necessary to fit the two estimators on different data splits for the theoretical guarantees to hold? [1] Abadie, A., Gu, J., & Shen, S. (2024). Instrumental variable estimation with first-stage heterogeneity. Journal of Econometrics, 240(2), 105425. \ [2] Coussens, S., & Spiess, J. (2021). Improving inference from simple instruments through compliance estimation. arXiv preprint arXiv:2108.03726. \ [3] Kennedy, E. H. (2023). Towards optimal doubly robust estimation of heterogeneous causal effects. Electronic Journal of Statistics, 17(2), 3008-3049.

Questions

I think my largest concerns were expressed above in the Weaknesses. I’d be happy to raise my score with a thorough and precise response to especially 1-2. Addressing the following might help strengthen the paper even more, but I don’t consider them as high-priority. 1. I’m slightly confused about the interpretation of Lemma 1 — to me, it seems like we have simply chosen to use IPW to estimate the numerator in Eq. 3. Is this understanding correct? If so, why not consider alternative estimators for the numerator (e.g., doubly-robust methods such as [1])? This is not a huge issue — just making sure I parsed the equation correctly. 2. Did the authors consider/evaluate the sensitivity of the proposed approach to different instantiations of the base learners (e.g., something besides RF/T-learner)? This is not a dealbreaker but would be a nice result to have. 3. Re: “The weighting scheme in Equation 4 creates a weighted distribution…” (L172) — what is this a weighted distribution of? The paper goes on to claim that the difference between the weighted distribution and the target distribution creates a transfer learning problem, but I’m having a bit of trouble understanding what these distributions are. Are these distributions defined over the covariates? 4. Re: Assumption 3-4 (realizability of $b(x)$) — could the authors expand on why this assumption is a reasonable one to make, or point to some works that have made similar assumptions? I think I’ve seen versions of the other assumptions, so I buy those, but I’m not sure if assuming $b(x)$ is linear in the representation $\phi$ is too limiting. Similarly, in Sec. 4.2, there’s an assumption that $\tau^O$ and $\tau$ — one of which is biased — have some shared representation. Is there a more concrete reason for why this is reasonable? [1] Kennedy, E. H. (2023). Towards optimal doubly robust estimation of heterogeneous causal effects. Electronic Journal of Statistics, 17(2), 3008-3049.

Rating

6

Confidence

4

Soundness

3

Presentation

3

Contribution

3

Limitations

The authors address their limitations in the Appendix. As a tiny suggestion, I think it is important to note that the lack of unobserved confounding is statistically unverifiable yet necessary for many causal inference approaches. Overall, I find the discussion of the limitations to be complete and well-written. As a very, very minor nitpick, I do believe it important to acknowledge the limitations of an approach more up-front (often in the conclusion) — it does not detract from my appreciation of the method and relegating discussion of limitations to the Appendix does not feel quite right to me.

Reviewer znEs6/10 · confidence 4/52024-07-17

Summary

The paper tackles the problem of estimating CATE when unobserved confounding is present in an observation study but an IV experiments is accessible, though the instrument could be weak. The paper proposes a two stage framework to first learn a biased CATE from observational data and makes a bias correction using complican weighted IV samples. The paper then demonstrate the effectiveness using simulation studies and a real world example on 401k.

Strengths

Originality: There is a lot of work on combining observational data with experimental data to better estimate causal estimands. This paper is the first to consider combining with an IV study with potentially weak instruments. Quality: The paper has clear theoretical results covering two cases and both simulation study and real world data example. Clarity: I found the paper easy to follow, with clear related work, contributions of the paper, motivation, theoretical results and experiments.

Weaknesses

Experiments: 1. It would be better to have some baseline comparisons, for example some debiased CATE estimation methods; 2. The simulation study and the real world example both use an equal size IV/observational data. What would happen if we only have a much smaller sample size IV study? I believe this is more common in real life since observational data is cheaper to get. 3. It would be interesting to see what happens in high-dimensional settings.

Questions

1. Do you have any results on violation of different assumptions? For example the realizibility. 2. Also see weakness on questions on experiments.

Rating

6

Confidence

4

Soundness

3

Presentation

3

Contribution

3

Limitations

The authors list the limitations in the appendix.

Reviewer xHgN2024-08-09

Thanks for the response! I appreciate the detailed answers, and I think most of my misunderstandings have been cleared up. Here's where I stand now: **W1:** Ok, I think I partially follow this — just to double-check, am I correct in saying that we cannot apply such past works because they can assist in targeting $\tau(x)$ (which is equivalent to $\tau^E(x)$ in the setting considered), but one of the core contributions of the work is incorporating $\tau^O(x)$ (or rather, the observational dataset in general) into estimates of $\tau^E(x)$ (i.e., as stated at the top of Section 4)? **W2:** Thanks for the tables — yeah, I can definitely see that the error bars would overlap. This would be great to turn into a figure for a camera-ready/future revision. **W3:** The clarifications make sense and address my concerns. I can see why alternative bounds in the literature would be incomparable, since they operate in a completely different problem setting from the proposed work. On a 2nd glance, the derived bounds indeed depend on both parameters of the observational and experimental datasets, which makes a lot of sense. **Questions:** Thanks for the clarifications — they've cleared up the misunderstandings. Good connections to shared-representation based models as well; I think I was a tiny bit concerned about assuming so such shared structure such that the $\tau - \tau^O$ is *linear* in some $\phi(x)$; the cited shared-representation approaches (TarNET [4] and DragonNet [5]; using citation numbers from rebuttal) still "split" the shared representation into neural-net based heads for each counterfactual distribution, right? Overall, while weakening this assumption would be nice, it's not a dealbreaker. My assessment of the paper has definitely improved after re-evaluation — the detailed responses targeted my concerns and have helped me gain a better understanding + appreciation of this work. I'm upgrading my score to 6 (WA).

Authorsrebuttal2024-08-13

Thank you for reconsidering and raising your score! **W1:** Yes, you are correct. Additionally, the other works target LATE, which is a compliance-weighted average of $\tau(x)$ over the population $\mathcal{X}$ (similar to ATE in observational studies). They use compliance weighting to reduce the variance of the LATE estimator under a homogeneous linear IV model, which we do not assume here. **W2:** We will include that figure in the camera-ready version. **Questions:** You’re right—while the shared representation is theoretically split for different heads, in most applications (e.g., [3]), the part after the split is typically just a linear transformation. We will discuss how to potentially weaken this assumption in the camera-ready version. Please let us know if there are any other questions we can address.

Reviewer Gwhh2024-08-12

Thank you for the rebuttal. My concerns have been addressed, and I wish to maintain my score.

Program Chairsdecision2024-09-25

Decision

Accept (poster)

© 2026 NYSGPT2525 LLC