We consider the problem of conditional independence testing: given a response Y and covariates (X,Z), we test the null hypothesis that Y is independent of X given Z. The conditional randomization test (CRT) was recently proposed as a way to use distributional information about X|Z to exactly (non-asymptotically) control Type-I error using any test statistic in any dimensionality without assuming anything about Y|(X,Z). This flexibility in principle allows one to derive powerful test statistics from complex prediction algorithms while maintaining statistical validity. Yet the direct use of such advanced test statistics in the CRT is prohibitively computationally expensive, especially with multiple testing, due to the CRT's requirement to recompute the test statistic many times on resampled data. We propose the distilled CRT, a novel approach to using state-of-the-art machine learning algorithms in the CRT while drastically reducing the number of times those algorithms need to be run, thereby taking advantage of their power and the CRT's statistical guarantees without suffering the usual computational expense. In addition to distillation, we propose a number of other tricks like screening and recycling computations to further speed up the CRT without sacrificing its high power and exact validity. Indeed, we show in simulations that all our proposals combined lead to a test that has similar power to the most powerful existing CRT implementations but requires orders of magnitude less computation, making it a practical tool even for large data sets. We demonstrate these benefits on a breast cancer dataset by identifying biomarkers related to cancer stage.
翻译:我们考虑的是有条件独立测试的问题:根据Y和Covariates(X,Z)的答复,我们测试了Y独立于X给Z的无效假设。 有条件随机测试(CRT)最近被提议作为一种方法,利用X ⁇ 的分布信息来(非非非同步)控制(非同步)类型I错误,在任何维度中使用任何测试统计,而没有假定Y ⁇ (X,Z)任何数据。这种灵活性原则上允许人们从复杂的预测算法中获取强大的测试统计数据,同时保持统计有效性。然而,在CRT中直接使用这种高级测试统计数据,在计算成本方面代价过高,特别是多次测试。由于CRT要求多次重新计算测试统计数据,我们提出了将X ⁇ (X,X,Z)的分布信息进行精确控制,这是在任何维度上使用最先进的机器学习算法的新办法,同时大幅降低这些算法运行的次数,从而在保持统计有效性和CRT统计保障的同时不承受通常的计算成本。除了估算外,我们提议在C级中进行大规模的计算,我们提议在模拟阶段里进行一个更高的计算。