Covariate Shift Corrected Conditional Randomization Test

Conditional independence tests are crucial across various disciplines in determining the independence of an outcome variable $Y$ from a treatment variable $X$, conditioning on a set of confounders $Z$. The Conditional Randomization Test (CRT) offers a powerful framework for such testing by assuming known distributions of $X \mid Z$; it controls the Type-I error exactly, allowing for the use of flexible, black-box test statistics. In practice, testing for conditional independence often involves using data from a source population to draw conclusions about a target population. This can be challenging due to covariate shift -- differences in the distribution of $X$, $Z$, and surrogate variables, which can affect the conditional distribution of $Y \mid X, Z$ -- rendering traditional CRT approaches invalid. To address this issue, we propose a novel Covariate Shift Corrected Pearson Chi-squared Conditional Randomization (csPCR) test. This test adapts to covariate shifts by integrating importance weights and employing the control variates method to reduce variance in the test statistics and thus enhance power. Theoretically, we establish that the csPCR test controls the Type-I error asymptotically. Empirically, through simulation studies, we demonstrate that our method not only maintains control over Type-I errors but also exhibits superior power, confirming its efficacy and practical utility in real-world scenarios where covariate shifts are prevalent. Finally, we apply our methodology to a real-world dataset to assess the impact of a COVID-19 treatment on the 90-day mortality rate among patients.

Paper

Similar papers

Peer review

Reviewer REEU7/10 · confidence 4/52024-07-04

Summary

The paper introduces a novel approach to addressing covariate shift in conditional independence tests. By leveraging importance weights and the control variates method,the paper proposes Covariate Shift Corrected Pearson Chi-squared Conditional Randomization (csPCR) test, which maintains asymptotic Type-I error control. In addition, a power enhancement is proposed, which reduces the variance in the test statistic and thus improve the power towards alternatives. The methodology is validated both theoretically and empirically, demonstrating superior performance in simulation studies and practical application to a COVID-19 treatment dataset, highlighting its efficacy in real-world scenarios.

Strengths

Originality The paper presents original contributions by introducing a new approach to conditional independence testing under covariate shift. This method combines the idea of conditional randomization test, importance weighting, the control variates technique. Quality The work is of high quality, demonstrated through rigorous statistical theories on the validity and power, and comprehensive empirical validations. Clarity The paper is clearly written and well-organized. Each section has a clear message and readers can quickly grasp the main idea. Significance I consider this work as a significant contribution to the field of conditional independence testing and covariate shift. The method is easy-to-implement, and is more powerful than existing methods as shown in numerical experiments.

Weaknesses

The paper does not leave serious technical weaknesses, but I am concerned about the limited applicability of this method in real-world problems. Based on my understanding, one would test for conditional independence directly on the target distribution if the outcomes in the target population were easy to collect. So csPCR is useful only when the outcomes are costly to collect. This questions several application examples mentioned in the paper. First, in the college admission example (from Line 37-53), why don't the economists just collect the college admission results from the target population and then conduct the vanilla PCR? Typically, the target outcomes have to be collected for followup analysis after the conditional independence test. Second, Section 5 introduces a real-world application of this method on Covid-19 pandemic data, where the source and the target data set are segmented based on time. Given such a sequential data collecting process, it is unclear whether the iid assumption still holds in the source and target distribution. See Questions for more details. Despite of the concerns above, I do believe that this method is useful for certain genetic and biomedical problems, where the outcomes are difficult to collect, but other covariates (such as SNPs) are relatively easy. Therefore, I suggest the authors clarify the application examples a bit by discussing more reasonable applications with some references.

Questions

1. Page 5, Algorithm 1: Besides the variance issues considered in Algorithm 2, it is clear that the choice of test statistic in Eq (4) also affects the power. Can you comment on how to choose the test statistic? In the numerical experiment, the test statistic is specified to be the product of $Y$ and $X$, which makes sense as it measures the "correlation" between Y and X. However, $Z$ is not included in the test statistic, which seems to decrease its power. Have you tried any other test statistics? Do difference choices affect the power significantly? For instance, one alternative is to use $(Y - \beta_y Z)(X - \beta_x Z)$, where $\beta_y and \beta_x$ are regression coefficients for $Y\sim Z$ and $X\sim Z$ on an independent data set from source distribution. 2. Page 4, Line 170-177: The paper argues that they used PCR as the baseline test because PCR is more powerful than the vanilla CRT. But on the other hand, CRT is also better than PCR in the sense that it has finite-sample type I error control, rather than asymptotic type I error control. It's worth mentioning this point as well. While I notice that the simulations in the paper shows satisfactory type I error control, the type I errors may get inflated if, for example, the data have heavy tails. 3. Page 9, Section 9: Relating to the weaknesses, have you run any tests to check the iid assumption on the Covid-19 data set?

Rating

7

Confidence

4

Soundness

4

Presentation

4

Contribution

3

Limitations

The papers discuss about the limitations of their method under model misspecifications in Section 4, pointing out that the powerful enhancement can disappear when there is a full nonlinear component in the model. The authors have addressed the limitations adequately.

Reviewer r2377/10 · confidence 3/52024-07-11

Summary

This paper proposes a new variation of CRT to be applied in the presence of covariate shifts. The paper presents a method and one extension for each, with higher power. Then, the authors present some needed theoretical results and finish with experiments.

Strengths

- The paper proposes a theoretically correct approach for testing CI under covariate shift; - The paper presents their algorithms and theoretical results; - The paper had convincing experiments;

Weaknesses

- It would be interesting if the authors could estimate the full density ratio in their simulations (possibly in a high-d scenario) and then compare their results with the IS approach; - The paper does not conduct any experiment where Type-I error control is shown in a real dataset, e.g., like for example in [1,2]. References [1] Pogodin, R., Schrab, A., Li, Y., Sutherland, D. J., & Gretton, A. (2024). Practical Kernel Tests of Conditional Independence. arXiv preprint arXiv:2402.13196. [2] Maia Polo, Felipe, Yuekai Sun, and Moulinath Banerjee. "Conditional independence testing under misspecified inductive biases." Advances in Neural Information Processing Systems 36 (2023): 58577-58612.

Questions

- How can we understand the role of the effective sample sizes (weight imbalance) in the effectiveness of your method? - Could you elaborate more on the surrogate variables? I am familiar with the literature on CI testing, but I haven't seen the presence of these variables in other works.

Rating

7

Confidence

3

Soundness

3

Presentation

3

Contribution

3

Limitations

The authors comment on limitations.

Reviewer r2372024-08-09

Thank you for your reply! - W1: Thank you for your new experiment! (this suggestion is related to the ESS question because in high-d we expect the ESS to be lower) - W2: I think the idea is good. When shuffling Y, you will need to do that within each value of Z, correct? If your Z is continuous, you would probably need to use some binning. - Q1: My question was related to the paper you mentioned here by Martino et al. In that paper, they argue that the ESS can be defined as a variance ratio and the implications of that are clear in the IS literature. I was wondering if you could extract some more meaningful relationships here as well. - Q2: Thank you for the detailed explanation! I would probably try to input some introduction on the surrogate variable in the abstract since you mention it but it might not be clear you are going to use them throughout the paper (you are working in this specific setup, which is not the same as the original CRT setup). I have increased my score.

Authorsrebuttal2024-08-10

We express sincere thanks to your insightful comments and positive recognition on our work. $\textbf{W1}$: Thanks for pointing out that an increasing dimension tends to cause higher variational importance weights and lower ESS. We will follow this to design our additional simulation on power v.s. ESS (changing with the variance of importance weights). $\textbf{W2}$: Thanks! When shuffling $Y$, we think both (i) marginal permutation and (ii) conditioning on (the full set or a subset of) Z will cause Y and X to be independent conditional on Z. Thus, we will try both setups and you are correct that we can use binning if Z is continuous and relatively high-dimensional. $\textbf{Q1}$: Thanks for the question! We think the effective sample size of our csPCR estimator (without control variate) can be defined in a similar way as their paper, i.e., the original source sample size * the ratio between the variances of (i) and (ii), where (i) stands for the unweighted indicator on target (which one would obtain if they could observe and use the same amount of labeled sample on target for PCR), and (ii) stands for the weighted indicator variable on source. This ratio can be shown to be smaller or equal to 1, with the “=” holding only when there is no covariate shift between the source and target. Thus, covariate shift will cause loss of the effective sample size to our method and larger degree of covariate shift tends to induce larger variance of the weighted estimator on source and larger loss of the ESS. This interpretation is quite analog to that in IS literature. Further, we can write down the ESS for our power enhanced (PE) estimator and apply our Theorem 2 to show that it is larger than the ESS of csPCR without PE. $\textbf{Q2}$: Thanks for the suggestion and we will highlight more on the specific setup with surrogate in the abstract and introduction section!

Reviewer MwvS5/10 · confidence 3/52024-07-12

Summary

This paper introduces the Covariate Shift Corrected Pearson Chi-squared Conditional Randomization (csPCR) test, designed to address covariate shift in conditional independence testing. The csPCR method incorporates importance weights and employs the control variates method to enhance test power and reduce variance in the statistical analysis. Theoretical contributions demonstrate that csPCR controls Type-I error asymptotically. Empirical validations through simulation studies and a real-world application assessing COVID-19 treatment effectiveness showcase the method's effectiveness.

Strengths

1. **Theoretical contribution:** The paper makes a theoretical contribution by addressing the covariate shift in conditional independence testing. 2. **Methodological Innovation:** Introduction of the Covariate Shift Corrected Pearson Chi-squared Conditional Randomization (csPCR) test that incorporates importance weights and control variates method to manage covariate shifts effectively. 3. **Empirical Validation:** Extensive simulation studies and real-world application (assessment of COVID-19 treatment on 90-day mortality) demonstrate the practical efficacy and superior power of the proposed csPCR test over traditional methods.

Weaknesses

1. **Dependence on Accurate Estimations of Density Ratio:** The csPCR test's performance is critically dependent on the accurate estimation of density ratios. This dependence could pose significant challenges in practical scenarios characterized by limited, noisy, or high-dimensional data, where reliable density ratio estimation becomes inherently difficult. 2. **Clarity on the Advantage Over Resampling Methods:** The paper does not explicitly clarify the advantages of the csPCR test over simpler resampling-based methods that also utilize estimated density ratios. While the csPCR integrates density ratios directly into the test statistic calculation and employs variance reduction techniques, it remains crucial for the authors to demonstrate why these features offer substantial improvements over traditional methods that adjust the sample weights during resampling. The discussion should address whether the complexity of csPCR provides tangible benefits, such as improved error rates or robustness in more varied practical applications, compared to potentially simpler resampling approaches.

Questions

- Could the authors provide a detailed comparison between the csPCR test and traditional resampling methods that also use estimated density ratios? This would further strengthen the paper - The paper mentions the use of control variates to reduce variance introduced by importance weights. Could the authors discuss how this approach compares to variance reduction techniques used in resampling methods? - Could the authors elaborate on any theoretical limitations or potential failure modes of the csPCR method? - Are there additional empirical validations planned or underway to test the csPCR method across more varied datasets?

Rating

5

Confidence

3

Soundness

3

Presentation

3

Contribution

2

Limitations

The authors are encouraged to add a separate limitation section in either the main paper or appendix.

Authorsrebuttal2024-08-07

We extend our sincere thanks for your reviewing work and insightful comments. Please see our responses to your comments and questions as below in the rebuttal.

Area Chair oa132024-08-11

Reminder to engage in the discussion

Hello reviewer MwvS. The reviewer-author discussion period is between Aug 7-13 (AoE). Please read the authors’ rebuttal, see if it addresses your questions, and engage in a discussion with the authors. You are strongly encouraged to read the official reviews posted by other reviewers. Some of your questions may have already been answered. **You must acknowledge the authors’ rebuttal and give a reason for why it did or did not address your concerns. If the rebuttal does not change your rating, please also state the reason.** Thank you for your services.

Reviewer G2Ah4/10 · confidence 1/52024-07-15

Summary

This paper addresses the issue of conditional independence under covariate shift. The authors' goal is to test the conditional independence for causal inference in the target data. The authors can use the source data whose distribution is potentially different from that of the target data. For this problem, the authors propose using the covariate-shift adaptation method and develop a method for testing the null hypothesis.

Strengths

This paper addresses the issue of conditional independence under covariate shift, which I find to be an intriguing problem. However, I was unable to comprehend the underlying assumptions of this problem, making it difficult to proceed with the reading. I understand the objective of this manuscript to be as follows: 1. The authors aim to test the null hypothesis $H_0: X \perp Y \mid Z$ for the target data (Here, I denote the independence by $\perp$). 2. There is a possibility that $X \perp Y \mid Z$ does not hold simultaneously for both the source and target data. 3. To test the null hypothesis for the target data, the authors employ an algorithm adapted to covariate shift. My primary concern is whether the assumption "there is a possibility that $X \perp Y \mid Z$ does not hold simultaneously for both the source and target data" is justified in the first place. The authors assume that, in order to use the covariate shift adaptation algorithm, - $Y \mid X, Z, V$ is the same in both the source and target data. Under this assumption, is it not the case that the situation where $X \perp Y \mid Z$ holds in one dataset and does not hold in the other cannot exist? $X \perp Y \mid Z$ implies that $p(y, x \mid z) = p(y \mid x, z) p(x \mid z) = p(y \mid z) p(x \mid z)$, meaning that $p(y \mid x, z) = p(y \mid z)$. Given that we assume $Y \mid X, Z, V$ to be the same in both the source and target data, does this not imply that whether $X \perp Y \mid Z$ holds cannot differ between the source and target data? Due to this ambiguity, I was unable to further evaluate the paper. I would like the authors to clarify this point so that I can proceed with the evaluation.

Weaknesses

See above.

Questions

See above.

Rating

4

Confidence

1

Soundness

2

Presentation

2

Contribution

2

Limitations

See above.

Reviewer G2Ah2024-08-11

Re: Rebuttal by Authors

I appreciate the authors' reply, which has deepened my understanding of your contributions. It seems I might have caused some misunderstanding with my previous question. Here is what I intended to ask: Objective of this study: To investigate whether conditional independence holds in the target data. Assumption in this study: Covariate shift = meaning that P(Y | X, Z, V) remains unchanged between the source data and the target data. Under this situation, wouldn't testing only the source data suffice to achieve the objective, given the assumption? I wanted to confirm whether I misunderstood the problem setting.

Authorsrebuttal2024-08-12

Thank you for your reply and clarification! You are correct that our goal is to test $X \perp Y \mid Z$ on the target data, with the assumption that $P(Y \mid X, Z, V)$ remains consistent between the source and target domains. However, even if $P(Y \mid X, Z, V)$ is the same for both the source and target, $P(V \mid X, Z)$ can still differ between them. As we demonstrated in our example in the rebuttal, this difference in $P(V \mid X, Z)$ led to $X \perp Y \mid Z$ holding true on the source data but not on the target. If we only test on the source data, this discrepancy can significantly inflate the Type-I error rate, as illustrated in Figure 1 of our paper. Please Let us know if you have any further concerns.

Reviewer G2Ah2024-08-13

Re: Official Comment by Authors

Thank you for your response. I have gained a deeper understanding of the authors' contributions. However, to be honest, I am not entirely clear about the motivation behind this study. If the authors assume that the conditional distributions are equal between the source and target data, wouldn't it suffice to test only for the source data? I am not convinced about the necessity of concerning with the Type 1 error in test for the target data. In the other words, it is unclear for me why we conduct test for the target data. I also have doubts about the claim that the Type 1 error changes. I think this is due to covariate shift, but if we consider covariates as non-stochastic, wouldn't the Type 1 error remain unchanged? My background is actually in economics. While the authors provide examples from economics, in social science data analysis, it is common to estimate models using only the source data and then conduct counterfactual simulations on the target data under the assumption that the conditional distributions are equal. That is, it is enough to care only about the model on the source data, and we can transfer the results under the assumption that the conditional distribution is invariant. Given these typical analysis procedures, I find the problem setting somewhat difficult to understand. I do not oppose acceptance, but I have some reservations about the problem setting, so I will keep my current score. In the future, since this problem setting may not be intuitively accepted, it might be worthwhile to focus more on justifying the problem setting itself more rather than on the technical issue of controlling the Type 1 error.

Reviewer REEU2024-08-08

Thanks for the clarification. I don't have further concerns.

Program Chairsdecision2024-09-25

Decision

Accept (poster)

© 2026 NYSGPT2525 LLC