Functional Horseshoe Smoothing for Functional Trend Estimation
share
Due to developments in instruments and computers, functional observations are increasingly popular. However, effective methodologies for flexibly estimating the underlying trends with valid uncertainty quantification for a sequence of functional data (e.g. functional time series) are still scarce. In this work, we develop a locally adaptive smoothing method, called functional horseshoe smoothing, by introducing a shrinkage prior to the general order of differences of functional variables. This allows us to capture abrupt changes by taking advantage of the shrinkage capability and also to assess uncertainty by Bayesian inference. The fully Bayesian framework also allows the selection of the number of basis functions via the posterior predictive loss. Also, by taking advantage of the nature of functional data, this method is able to handle heterogeneously observed data without data augmentation. We show the theoretical properties of the proposed prior distribution and the posterior mean, and finally demonstrate them through simulation studies and applications to a real-world dataset.
READ FULL TEXT