Abstract:
|
While Bayesian functional mixed models have been shown effective to model functional data with various complex structures, their application to extremely high-dimensional data is limited due to computational challenges involved in posterior sampling. We introduce a new computational framework that enables ultra-fast approximate inference for high-dimensional data in functional form. This framework adopts parsimonious basis to represent functional observations, which facilitates efficient compression and parallel computing in basis space. Instead of performing expensive Markov chain Monte Carlo sampling, we approximate the posterior distribution using variational Bayes and adopt fast iterative algorithms to estimate parameters of the approximate distribution. Our approach provides a multiple testing procedure in basis space, which can be used to identify significant local regions on functional data that are nonzero or different across groups of samples. We perform two simulation studies to assess the accuracy of approximate inference and the computational benefits, and demonstrate applications of the proposed approach by using two real datasets.
|