This repository contains the R code used for research done in the paper "Flexible Quantile Contour Estimation for Multivariate Functional Data: Beyond Convexity". It also contains the links to the sources of datasets used as an application in this paper.
Nowadays, multivariate functional data are frequently observed in many scientific fields, and the estimation of quantiles of these data is essential in data analysis. Unlike in the univariate setting, quantiles are more challenging to estimate for multivariate data, let alone multivariate functional data. This article proposes a new method to estimate the quantiles for multivariate functional data. The proposed multivariate functional quantile model is a nonparametric, time-varying coefficient model, and basis functions are used for the estimation and prediction. The estimated quantile contours can characterize non-Gaussian and even nonconvex multivariate distributions marginally, and the estimated multivariate quantile function is a continuous function of time for a fixed quantile level. Computationally, the proposed method is efficient for high dimensions, and can handle more than just bivariate functional data. The monotonicity, uniqueness, and consistency of the estimated multivariate quantile function have been established. The proposed method was demonstrated on bivariate and trivariate functional data in the simulation studies, and was applied to study the joint distribution of