Bayesian Generalized Two-way ANOVA Modeling for Functional Data Using INLA
Journal article, 2019

Functional analysis of variance (ANOVA) modeling has been proved particularly useful to investigate the dynamic changes of functional data according to certain categorical factors and their interactions. The current existing methods often encounter difficulties when the functional data are high-dimensional, non-Gaussian, and/or exhibit certain shape characteristics that vary with spatial locations. In this paper, we investigate the functional two-way ANOVA models from a novel Bayesian perspective. A class of highly flexible Gaussian Markov random fields (GMRF) are taken as priors on the functions in the model, which allows us to model various types of functional effects, such as (discrete or continuous) temporal effects and (point-level or areal) spatial effects. The resulting posterior distributions are obtained by an efficient computational tool based on integrated nested Laplace approximations (INLA) Rue, Martino and Chopin (2009). We then employ the excursion method introduced by Bolin and Lindgren (2015) to build simultaneous credible intervals of functional effects and test their significance from a Bayesian point of view. A simulation study and multiple data examples are presented to demonstrate the merits of our method.


Yu (Ryan) Yue

Baruch College

David Bolin

Chalmers, Mathematical Sciences, Applied Mathematics and Statistics

Håvard Rue

King Abdullah University of Science and Technology (KAUST)

Xiao-Feng Wang

Cleveland State University

Statistica Sinica

1017-0405 (ISSN)

Vol. 29 2 741-767

Subject Categories

Probability Theory and Statistics



More information

Latest update