Track papers across labs, topics, and trends.
Search papers, labs, and topics across Lattice.
This paper introduces a novel sequential whitening approach to conformal prediction for spatially dependent data, addressing the breakdown of exchangeability assumptions in standard split conformal methods. By conditioning on calibration residuals sequentially, the method effectively removes predictable spatial variation from both calibration and target errors, leading to improved efficiency and stability. The proposed technique achieves exact finite-sample coverage under correct covariance assumptions and any spatial design, and demonstrates superior performance in simulations and a real-world PM2.5 application.
Spatially dependent data breaks standard conformal prediction; this method achieves exact coverage and oracle efficiency by sequentially whitening residuals, even with nonrandom sampling.
Split conformal prediction uses prediction errors on held-out (calibration) data to determine how wide the prediction intervals should be. It guarantees distribution-free finite-sample coverage when these errors and the error at the target site are exchangeable. This assumption may fail under spatial dependence and nonrandom sampling geometry. Existing spatial methods use fitting residuals to remove the predictable part of spatial variation from calibration and target errors. However, the spatial variation that only the calibration residuals can predict remains in both the target and calibration errors, reducing the efficiency and stability of the interval. We address this by additionally conditioning on the calibration residuals sequentially, which scales to large networks through nearest-neighbour approximations. Under a correct working covariance and an elliptical residual law, the resulting interval has exact finite-sample coverage under any spatial design, and under further conditions it is asymptotically oracle efficient. We also bound coverage loss under covariance misspecification and develop a diagnostic that identifies regions at risk of undercoverage. In simulated data, our method produces narrower and more stable intervals than global and localized state-of-the-art alternatives. In a national PM2.5 application, it produces narrower intervals within the network and identifies regions at risk of coverage failure.