Spatial modeling with R‐INLA: A review

Coming up with Bayesian models for spatial data is easy, but performing inference with them can be challenging. Writing fast inference code for a complex spatial model with realistically‐sized datasets from scratch is time‐consuming, and if changes are made to the model, there is little guarantee th...

Celý popis

Uložené v:
Podrobná bibliografia
Vydané v:Wiley interdisciplinary reviews. Computational statistics Ročník 10; číslo 6; s. e1443 - n/a
Hlavní autori: Bakka, Haakon, Rue, Håvard, Fuglstad, Geir‐Arne, Riebler, Andrea, Bolin, David, Illian, Janine, Krainski, Elias, Simpson, Daniel, Lindgren, Finn
Médium: Journal Article
Jazyk:English
Vydavateľské údaje: Hoboken, USA John Wiley & Sons, Inc 01.11.2018
Wiley Subscription Services, Inc
Predmet:
ISSN:1939-5108, 1939-0068, 1939-0068
On-line prístup:Získať plný text
Tagy: Pridať tag
Žiadne tagy, Buďte prvý, kto otaguje tento záznam!
Popis
Shrnutí:Coming up with Bayesian models for spatial data is easy, but performing inference with them can be challenging. Writing fast inference code for a complex spatial model with realistically‐sized datasets from scratch is time‐consuming, and if changes are made to the model, there is little guarantee that the code performs well. The key advantages of R‐INLA are the ease with which complex models can be created and modified, without the need to write complex code, and the speed at which inference can be done even for spatial problems with hundreds of thousands of observations. R‐INLA handles latent Gaussian models, where fixed effects, structured and unstructured Gaussian random effects are combined linearly in a linear predictor, and the elements of the linear predictor are observed through one or more likelihoods. The structured random effects can be both standard areal model such as the Besag and the BYM models, and geostatistical models from a subset of the Matérn Gaussian random fields. In this review, we discuss the large success of spatial modeling with R‐INLA and the types of spatial models that can be fitted, we give an overview of recent developments for areal models, and we give an overview of the stochastic partial differential equation (SPDE) approach and some of the ways it can be extended beyond the assumptions of isotropy and separability. In particular, we describe how slight changes to the SPDE approach leads to straight‐forward approaches for nonstationary spatial models and nonseparable space–time models. This article is categorized under: Statistical and Graphical Methods of Data Analysis > Bayesian Methods and Theory Statistical Models > Bayesian Models Data: Types and Structure > Massive Data In spatial statistics, an important problem is how to represent spatial models in a way that is computationally efficient, accurate, and convenient to use. Models in R‐INLA focus on sparse precision (inverse covariance) matrices to compute inference quickly. Hence, our implementations of spatial models focus on how to represent the spatial field in such a way that the precision matrix for the "representation" is very sparse. This graphic shows a representation of a Norwegian fjord with a mesh, from which basis functions are built in the finite element method. We use sums of these basis functions to represent the spatial field. This representation has many advantages, but requires some mathematical effort to understand and to set up.
Bibliografia:ObjectType-Article-1
SourceType-Scholarly Journals-1
ObjectType-Feature-2
content type line 14
ISSN:1939-5108
1939-0068
1939-0068
DOI:10.1002/wics.1443