We consider the problem of consistently estimating the conditional distribution of a functional data object given covariates in a general space, assuming that and 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 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 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 is an element of a separable Hilbert space, and when is an element of the Banach space of continuous functions. The latter results imply that sets that specify path properties of , 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.
View on arXiv