Bayesian Analysis

Bayesian functional {ANOVA} modeling using Gaussian process prior distributions

Cari G. Kaufman and Stephan R. Sain

Full-text: Open access


Functional analysis of variance (ANOVA) models partition a functional response according to the main effects and interactions of various factors. This article develops a general framework for functional ANOVA modeling from a Bayesian viewpoint, assigning Gaussian process prior distributions to each batch of functional effects. We discuss the choices to be made in specifying such a model, advocating the treatment of levels within a given factor as dependent but exchangeable quantities, and we suggest weakly informative prior distributions for higher level parameters that may be appropriate in many situations. We discuss computationally efficient strategies for posterior sampling using Markov Chain Monte Carlo algorithms, and we emphasize useful graphical summaries based on the posterior distribution of model-based analogues of traditional ANOVA decompositions of variance. We illustrate this process of model specification, posterior sampling, and graphical posterior summaries in two examples. The first considers the effect of geographic region on the temperature profiles at weather stations in Canada. The second example examines sources of variability in the output of regional climate models from a designed experiment.

Article information

Bayesian Anal., Volume 5, Number 1 (2010), 123-149.

First available in Project Euclid: 22 June 2012

Permanent link to this document

Digital Object Identifier

Mathematical Reviews number (MathSciNet)

Zentralblatt MATH identifier

Analysis of variance Climate models Functional data Variance components


Kaufman, Cari G.; Sain, Stephan R. Bayesian functional {ANOVA} modeling using Gaussian process prior distributions. Bayesian Anal. 5 (2010), no. 1, 123--149. doi:10.1214/10-BA505.

Export citation