Table of Contents
Fetching ...

Prediction-Augmented Trees for Reliable Statistical Inference

Vikram Kher, Argyris Oikonomou, Manolis Zampetakis

TL;DR

Prediction-Augmented Trees for Reliable Statistical Inference develops PART, a CART-style adaptive partitioning estimator, and PAQ, a quadrature-based estimator, to safely exploit ML predictions in mean and regression inference with scarce labeled data and abundant unlabeled data. PART provides asymptotic distributions and CI construction while achieving consistently tighter confidence intervals than PPI/PPI++ by local variance reduction; PAQ analyzes PART in a deep-tree limit and proves variance decline $O(N^{-1}+n^{-4})$ under smooth residuals, outperforming the traditional $O(N^{-1}+n^{-1})$ rate. The paper validates these methods on diverse real-world datasets (deforestation, astronomy, wine, housing, proteomics, census) and also develops a PAQ theory in one and higher dimensions, linking ML predictions to principled statistical inference. Together, PART and PAQ offer scalable, reliable tools for discovery that leverage unlabeled data and ML predictors to produce tighter, valid confidence intervals across scientific domains.

Abstract

The remarkable success of machine learning (ML) in predictive tasks has led scientists to incorporate ML predictions as a core component of the scientific discovery pipeline. This was exemplified by the landmark achievement of AlphaFold (Jumper et al. (2021)). In this paper, we study how ML predictions can be safely used in statistical analysis of data towards scientific discovery. In particular, we follow the framework introduced by Angelopoulos et al. (2023). In this framework, we assume access to a small set of $n$ gold-standard labeled samples, a much larger set of $N$ unlabeled samples, and a ML model that can be used to impute the labels of the unlabeled data points. We introduce two new learning-augmented estimators: (1) Prediction-Augmented Residual Tree (PART), and (2) Prediction-Augmented Quadrature (PAQ). Both estimators have significant advantages over existing estimators like PPI and PPI++ introduced by Angelopoulos et al. (2023) and Angelopoulos et al. (2024), respectively. PART is a decision-tree based estimator built using a greedy criterion. We first characterize PART's asymptotic distribution and demonstrate how to construct valid confidence intervals. Then we show that PART outperforms existing methods in real-world datasets from ecology, astronomy, and census reports, among other domains. This leads to estimators with higher confidence, which is the result of using both the gold-standard samples and the machine learning predictions. Finally, we provide a formal proof of the advantage of PART by exploring PAQ, an estimation that arises when considering the limit of PART when the depth its tree grows to infinity. Under appropriate assumptions in the input data we show that the variance of PAQ shrinks at rate of $O(N^{-1} + n^{-4})$, improving significantly on the $O(N^{-1}+n^{-1})$ rate of existing methods.

Prediction-Augmented Trees for Reliable Statistical Inference

TL;DR

Prediction-Augmented Trees for Reliable Statistical Inference develops PART, a CART-style adaptive partitioning estimator, and PAQ, a quadrature-based estimator, to safely exploit ML predictions in mean and regression inference with scarce labeled data and abundant unlabeled data. PART provides asymptotic distributions and CI construction while achieving consistently tighter confidence intervals than PPI/PPI++ by local variance reduction; PAQ analyzes PART in a deep-tree limit and proves variance decline under smooth residuals, outperforming the traditional rate. The paper validates these methods on diverse real-world datasets (deforestation, astronomy, wine, housing, proteomics, census) and also develops a PAQ theory in one and higher dimensions, linking ML predictions to principled statistical inference. Together, PART and PAQ offer scalable, reliable tools for discovery that leverage unlabeled data and ML predictors to produce tighter, valid confidence intervals across scientific domains.

Abstract

The remarkable success of machine learning (ML) in predictive tasks has led scientists to incorporate ML predictions as a core component of the scientific discovery pipeline. This was exemplified by the landmark achievement of AlphaFold (Jumper et al. (2021)). In this paper, we study how ML predictions can be safely used in statistical analysis of data towards scientific discovery. In particular, we follow the framework introduced by Angelopoulos et al. (2023). In this framework, we assume access to a small set of gold-standard labeled samples, a much larger set of unlabeled samples, and a ML model that can be used to impute the labels of the unlabeled data points. We introduce two new learning-augmented estimators: (1) Prediction-Augmented Residual Tree (PART), and (2) Prediction-Augmented Quadrature (PAQ). Both estimators have significant advantages over existing estimators like PPI and PPI++ introduced by Angelopoulos et al. (2023) and Angelopoulos et al. (2024), respectively. PART is a decision-tree based estimator built using a greedy criterion. We first characterize PART's asymptotic distribution and demonstrate how to construct valid confidence intervals. Then we show that PART outperforms existing methods in real-world datasets from ecology, astronomy, and census reports, among other domains. This leads to estimators with higher confidence, which is the result of using both the gold-standard samples and the machine learning predictions. Finally, we provide a formal proof of the advantage of PART by exploring PAQ, an estimation that arises when considering the limit of PART when the depth its tree grows to infinity. Under appropriate assumptions in the input data we show that the variance of PAQ shrinks at rate of , improving significantly on the rate of existing methods.
Paper Structure (51 sections, 16 theorems, 117 equations, 9 figures, 3 algorithms)

This paper contains 51 sections, 16 theorems, 117 equations, 9 figures, 3 algorithms.

Key Result

Lemma 1

Let $\widehat{\mu}_C$ be the estimate from estimator: non-adaptive, and assume $N \gg n$. Then $\widehat{\mu}_C$ is an unbiased estimator for the mean $\mu = \mathbb{E}[Y]$. Moreover, when $\operatorname{Var}(Y - f(X))$ is finite, and as $n \to \infty$,

Figures (9)

  • Figure 1: Theoretical coverage of \ref{['thm:tree-ci-coverage']} (dashed lines) versus empirical coverage fo the PART estimator (solid lines) at a 95% confidence level.
  • Figure 2: Estimating the mean deforestation rate in the Amazon rainforest from satellite imagery Bullock2020 (see \ref{['sec:experiments']} for testing methodology).
  • Figure 3: Estimating the mean fraction of spiral galaxies in the universe using Galaxy Zoo 2 data Willett2013 (see \ref{['sec:experiments']} for testing methodology).
  • Figure 4: Estimating the mean quality of Portugese Vinho Verde wine wine_data (see \ref{['sec:experiments']} for testing methodology).
  • Figure 5: Estimating the average median house price in a census block using California census data housing_data (see \ref{['sec:experiments']} for testing methodology).
  • ...and 4 more figures

Theorems & Definitions (40)

  • Lemma 1: Properties of \ref{['estimator: non-adaptive']}
  • Corollary 2: Weak Variance Reduction
  • proof : Sketch of Proof
  • Example 1: Strict Variance Reduction
  • Remark 1
  • Remark 2
  • Theorem 1: Coverage of \ref{['alg:semi-supervised_tree']}
  • Remark 3
  • Theorem 2
  • Theorem 3: Asymptotic Coverage
  • ...and 30 more