Functional ANOVA for Modelling Temperature Profiles in Ecuador

The experiments obtained from global and regional climate models have recently become popular and are widely used in the Statistical literature to analyze and understand the evolution of natural phenomena such as climate change in a limited area, using discretized versions of physical processes that they are modeled by systems of differential equations. This methodology is appropriate for the analysis of time series of climatic variables; therefore, the present work is motivated by the need to compare sources of variability in the projections made by the computational models used to model the climate. We propose to use a functional data modeling through an ANOVA, assigning prior distributions of Gaussian processes in each batch of the functional effects. Some alternatives to specify the model will be discussed, and computationally efficient strategies will be considered for the simulation of the samples and estimation of the posterior distribution using a Markov Chain Monte Carlo algorithm for the first case and a Laplace approximation for the second case. The methodology will be illustrated through two cases: the effect of the geographical region on the temperature profiles in the meteorological stations of Ecuador and the sources of variability in the output of regional climate models.

[1]  Hong Kong,et al.  Stochastic modeling , 2021, 2016 IEEE Electrical Design of Advanced Packaging and Systems (EDAPS).

[2]  P. Jones,et al.  Version 4 of the CRU TS monthly high-resolution gridded multivariate climate dataset , 2020, Scientific Data.

[3]  Finn Lindgren,et al.  Advanced Spatial Modeling with Stochastic Partial Differential Equations Using R and INLA , 2018 .

[4]  Yura Kim,et al.  Spatial analysis of future East Asian seasonal temperature using two regional climate model simulations , 2016, Asia-Pacific Journal of Atmospheric Sciences.

[5]  Hans-Georg Müller,et al.  Functional Data Analysis , 2016 .

[6]  H. Rue,et al.  Scaling intrinsic Gaussian Markov random field priors in spatial modelling , 2014 .

[7]  Jin-Ting Zhang,et al.  Analysis of Variance for Functional Data , 2013 .

[8]  R. Furrer,et al.  Conjugate distributions in hierarchical Bayesian ANOVA for computational efficiency and assessments of both practical and statistical significance , 2013, 1303.3390.

[9]  John K Kruschke,et al.  Bayesian data analysis. , 2010, Wiley interdisciplinary reviews. Cognitive science.

[10]  S. Sain,et al.  Bayesian functional ANOVA modeling using Gaussian process prior distributions , 2010 .

[11]  H. Rue,et al.  Approximate Bayesian inference for latent Gaussian models by using integrated nested Laplace approximations , 2009 .

[12]  Leonhard Held,et al.  Gaussian Markov Random Fields: Theory and Applications , 2005 .

[13]  A. Gelman Discussion of "Analysis of variance--why it is more important than ever" by A. Gelman , 2005, math/0508530.

[14]  J. Faraway Regression analysis for a functional response , 1997 .

[15]  Saba Infante,et al.  Estimation of Ordinary Differential Equations Solutions with Gaussian Processes and Polynomial Chaos Expansion , 2021, Information and Communication Technologies.

[16]  F. Lutscher Spatial Variation , 2019, Interdisciplinary Applied Mathematics.

[17]  H. Rue,et al.  Statistica Sinica Preprint No : SS-2016-0055 . R 3 Title Bayesian Generalized Two-way ANOVA Modeling for Functional Data Using INLA , 2017 .

[18]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[19]  J. Norris Appendix: probability and measure , 1997 .

[20]  Gregory Piatetsky-Shapiro,et al.  Advances in Knowledge Discovery and Data Mining , 2004, Lecture Notes in Computer Science.

[21]  Julian J. Faraway,et al.  An F test for linear models with functional responses , 2004 .