Smoothing functional data with a Bayesian hierarchical model

A Bayesian hierarchical model is developed for smoothing functional data. Functional data, with basic data unit being function evaluations (e.g. curves or surfaces) over a continuum, have been frequently encountered in nowadays. While many functional data analysis tools are now available, the issue of simultaneous smoothing is less emphasized. Some methods treat functional data as fully observed while ignoring the measurement noise, others perform smoothing to each functional observation independently thus fail to borrow strength across replications from the same stochastic process. In this paper, we propose a Bayesian hierarchical model to smooth all functional observations simultaneously. The proposed method relies on priors with data-driven hierarchical parameters, which automatically determine the amount of smoothness. It also provides simultaneous estimates for the mean function and covariance. Case studies of simulated and real data demonstrate that this Bayesian method produces more accurate signal estimates and smooth covariance estimate