Functional data analysis with covariate-dependent mean and covariance structures

Functional data analysis has emerged as a powerful tool in response to the ever-increasing resources and efforts devoted to collecting information about response curves or anything that varies over a continuum. However, limited progress has been made with regard to linking the covariance structures of response curves to external covariates, as most functional models assume a common covariance structure. We propose a new functional regression model with covariate-dependent mean and covariance structures. Particularly, by allowing variances of random scores to be covariate-dependent, we identify eigenfunctions for each individual from the set of eigenfunctions that govern the variation patterns across all individuals, resulting in high interpretability and prediction power. We further propose a new penalized quasi-likelihood procedure that combines regularization and B-spline smoothing for model selection and estimation and establish the convergence rate and asymptotic normality of the proposed estimators. The utility of the developed method is demonstrated via simulations, as well as an analysis of the Avon Longitudinal Study of Parents and Children concerning parental effects on the growth curves of their offspring, which yields biologically interesting results.

Learning delay dynamics for multivariate stochastic processes, with application to the prediction of the growth rate of COVID-19 cases in the United States

Delay differential equations form the underpinning of many complex dynamical systems. The forward problem of solving random differential equations with delay has received increasing attention in recent years. Motivated by the challenge to predict the COVID-19 caseload trajectories for individual states in the U.S., we target here the inverse problem. Given a sample of observed random trajectories obeying an unknown random differential equation model with delay, we use a functional data analysis framework to learn the model parameters that govern the underlying dynamics from the data. We show the existence and uniqueness of the analytical solutions of the population delay random differential equation model when one has discrete time delays in the functional concurrent regression model and also for a second scenario where one has a delay continuum or distributed delay. The latter involves a functional linear regression model with history index. The derivative of the process of interest is modeled using the process itself as predictor and also other functional predictors with predictor-specific delayed impacts. This dynamics learning approach is shown to be well suited to model the growth rate of COVID-19 for the states that are part of the U.S., by pooling information from the individual states, using the case process and concurrently observed economic and mobility data as predictors.

Adaptive smoothing spline estimator for the function-on-function linear regression model

In this paper, we propose an adaptive smoothing spline (AdaSS) estimator for the function-on-function linear regression model where each value of the response, at any domain point, depends on the full trajectory of the predictor. The AdaSS estimator is obtained by the optimization of an objective function with two spatially adaptive penalties, based on initial estimates of the partial derivatives of the regression coefficient function. This allows the proposed estimator to adapt more easily to the true coefficient function over regions of large curvature and not to be undersmoothed over the remaining part of the domain. A novel evolutionary algorithm is developed ad hoc to obtain the optimization tuning parameters. Extensive Monte Carlo simulations have been carried out to compare the AdaSS estimator with competitors that have already appeared in the literature before. The results show that our proposal mostly outperforms the competitor in terms of estimation and prediction accuracy. Lastly, those advantages are illustrated also in two real-data benchmark examples. The AdaSS estimator is implemented in the package , openly available online on CRAN.

A functional proportional hazard cure rate model for interval-censored data

Existing survival models involving functional covariates typically rely on the Cox proportional hazards structure and the assumption of right censorship. Motivated by the aim of predicting the time of conversion to Alzheimer's disease from sparse biomarker trajectories in patients with mild cognitive impairment, we propose a functional mixture cure rate model with both functional and scalar covariates for interval censoring and sparsely sampled functional data. To estimate the nonparametric coefficient function that depicts the effect of the shape of the trajectories on the survival outcome and cure probability, we utilize the functional principal component analysis to extract the functional features from the sparsely and irregularly sampled trajectories. To obtain parameter estimates from the mixture cure rate model with interval censoring, we apply the expectation-maximization algorithm based on Poisson data augmentation. The estimation accuracy of our method is assessed via a simulation study and we apply our model on Alzheimer's disease Neuroimaging Initiative data set.

Functional regression on the manifold with contamination

We propose a new method for functional nonparametric regression with a predictor that resides on a finite-dimensional manifold, but is observable only in an infinite-dimensional space. Contamination of the predictor due to discrete or noisy measurements is also accounted for. By using functional local linear manifold smoothing, the proposed estimator enjoys a polynomial rate of convergence that adapts to the intrinsic manifold dimension and the contamination level. This is in contrast to the logarithmic convergence rate in the literature of functional nonparametric regression. We also observe a phase transition phenomenon related to the interplay between the manifold dimension and the contamination level. We demonstrate via simulated and real data examples that the proposed method has favourable numerical performance relative to existing commonly used methods.

Methods for scalar-on-function regression

Recent years have seen an explosion of activity in the field of functional data analysis (FDA), in which curves, spectra, images, etc. are considered as basic functional data units. A central problem in FDA is how to fit regression models with scalar responses and functional data points as predictors. We review some of the main approaches to this problem, categorizing the basic model types as linear, nonlinear and nonparametric. We discuss publicly available software packages, and illustrate some of the procedures by application to a functional magnetic resonance imaging dataset.

Photocatalytic degradation of furfural in aqueous solution by N-doped titanium dioxide nanoparticles

The photocatalytic degradation of furfural in aqueous solution was investigated using N-doped titanium dioxide nanoparticles under sunlight and ultraviolet radiation (N-TiO₂/Sun and N-TiO₂/UV) in a lab-scale batch photoreactor. The N-TiO₂ nanoparticles prepared using a sol-gel method were characterized using XRD, X-ray photoelectron spectroscopy (XPS), and SEM analyses. Using HPLC to monitor the furfural concentration, the effect of catalyst dosage, contact time, initial solution pH, initial furfural concentration, and sunlight or ultraviolet radiation on the degradation efficiency was studied. The efficiency of furfural removal was found to increase with increased reaction time, nanoparticle loading, and pH for both processes, whereas the efficiency decreased with increased furfural concentration. The maximum removal efficiencies for the N-TiO₂/UV and N-TiO₂/Sun processes were 97 and 78 %, respectively, whereas the mean removal efficiencies were 80.71 ± 2.08 % and 62.85 ± 2.41 %, respectively. In general, the degradation and elimination rate of furfural using the N-TiO₂/UV process was higher than that using the N-TiO₂/Sun process.

Longitudinal Functional Models with Structured Penalties

This article addresses estimation in regression models for longitudinally-collected functional covariates (time-varying predictor curves) with a longitudinal scaler outcome. The framework consists of estimating a time-varying coefficient function that is modeled as a linear combination of time-invariant functions with time-varying coefficients. The model uses extrinsic information to inform the structure of the penalty, while the estimation procedure exploits the equivalence between penalized least squares estimation and a linear mixed model representation. The process is empirically evaluated with several simulations and it is applied to analyze the neurocognitive impairment of HIV patients and its association with longitudinally-collected magnetic resonance spectroscopy (MRS) curves.