New Confidence Regions for Linear Regression Parameters with Stationary-Ergodic Dependent Errors
This study addresses the challenge of constructing valid joint confidence regions for linear regression coefficients when regression errors exhibit unknown serial dependence and are jointly stationary and ergodic with the covariates. The authors propose a novel approach that avoids explicit modeling of the error dependence structure by introducing independent auxiliary samples and applying stochastic smoothing with a decaying bandwidth to both the regression function and second moments. Coupled with data-driven bandwidth selection and mild truncation, this method yields Wald-type confidence regions and simultaneous confidence intervals. It does not rely on long-run variance estimation or parametric assumptions about dependence, achieving coverage probabilities close to nominal levels across diverse dependence structures—including ARMA, ARFIMA, copula-based Markov processes, and fractional Gaussian noise—while producing smaller confidence region volumes than Newey–West HAC and MAC methods. The approach is successfully demonstrated in an analysis of Beijing PM2.5 data.