Estimating the conditional distribution in functional regression problems

TL;DR

This paper develops methods for consistently estimating the conditional distribution of a functional response given covariates in a functional linear regression framework, applicable in various function spaces.

Contribution

It introduces two estimation approaches based on residual bootstrapping and parametric modeling, with theoretical guarantees under certain conditions.

Findings

Consistent estimation achieved under general conditions.

Methods applicable to Hilbert and Banach spaces.

Validated through simulations and real data analysis.

Abstract

We consider the problem of consistently estimating the conditional distribution $P(Y \in A |X)$ of a functional data object $Y=(Y(t): t\in[0,1])$ given covariates $X$ in a general space, assuming that $Y$ and $X$ are related by a functional linear regression model. Two natural estimation methods are proposed, based on either bootstrapping the estimated model residuals, or fitting functional parametric models to the model residuals and estimating $P(Y \in A |X)$ via simulation. Whether either of these methods lead to consistent estimation depends on the consistency properties of the regression operator estimator, and the space within which $Y$ is viewed. We show that under general consistency conditions on the regression operator estimator, which hold for certain functional principal component based estimators, consistent estimation of the conditional distribution can be achieved, both when $Y$ is an element of a separable Hilbert space, and when $Y$ is an element of the Banach space of continuous functions. The latter results imply that sets $A$ that specify path properties of $Y$, which are of interest in applications, can be considered. The proposed methods are studied in several simulation experiments, and data analyses of electricity price and pollution curves.