Consider a heteroscedastic regression model Y=m(X) +σ(X)ε, where the functions m and σ are “smooth”, and ε is independent of X. An estimator of the distribution of ε based on non-parametric regression residuals is proposed and its weak convergence is obtained. Applications to prediction intervals and goodness-of-fit tests are discussed.
Akritas, M. G., & Van Keilegom, I. (2001). Non-parametric Estimation of the Residual Distribution. Scandinavian Journal of Statistics : theory and applications, 28(3), 549-567. https://doi.org/10.1111/1467-9469.00254 (Original work published 2001)