Skip to main content

Lund University Publications

LUND UNIVERSITY LIBRARIES

An explicit link between Gaussian fields and Gaussian Markov random fields; The SPDE approach

Lindgren, Finn LU ; Lindström, Johan LU orcid and Rue, Håvard (2010) In Preprints in Mathematical Sciences
Abstract
Continuously indexed Gaussian fields (GFs) is the most important ingredient in spatial statistical modelling and geo-statistics. The specification through the covariance function gives an intuitive interpretation of its properties. On the computational side, GFs are hampered with the "big-n" problem, since the cost of factorising dense matrices is cubic in the dimension. Although the computational power today is all-time-high, this fact seems still to be a computational bottleneck in applications. Along with GFs, there is the class of Gaussian Markov random fields (GMRFs) which are discretely indexed. The Markov property makes the involved precision matrix sparse which enables the use of numerical algorithms for sparse matrices, that for... (More)
Continuously indexed Gaussian fields (GFs) is the most important ingredient in spatial statistical modelling and geo-statistics. The specification through the covariance function gives an intuitive interpretation of its properties. On the computational side, GFs are hampered with the "big-n" problem, since the cost of factorising dense matrices is cubic in the dimension. Although the computational power today is all-time-high, this fact seems still to be a computational bottleneck in applications. Along with GFs, there is the class of Gaussian Markov random fields (GMRFs) which are discretely indexed. The Markov property makes the involved precision matrix sparse which enables the use of numerical algorithms for sparse matrices, that for fields in R^2 only use the square-root of the time required by general algorithms. The specification of a GMRF is through its full conditional distributions but its marginal properties are not transparent in such a parametrisation.



In this paper, we show that using an approximate stochastic weak solution to (linear) stochastic partial differential equations (SPDEs), we can, for some GFs in the Matérn class, provide an

explicit link, for any triangulation of R^d, between GFs and GMRFs. The consequence is that we can take the best from the two worlds and do the modelling using GFs but do the computations using GMRFs. Perhaps more importantly, our approach generalises to other covariance functions generated by SPDEs, including oscillating and non-stationary GFs, as well as GFs on manifolds. We illustrate our approach by analysing global temperature data with a non-stationary model defined on a sphere. (Less)
Please use this url to cite or link to this publication:
author
; and
organization
publishing date
type
Book/Report
publication status
published
subject
keywords
Gaussian Markov random fields, Generalised additive mixed models, Laplace approximation, Sparse matrices, Stochastic partial differential equations, Approximate Bayesian inference, Parallel computing, Structured additive regression models.
in
Preprints in Mathematical Sciences
pages
39 pages
publisher
[Publisher information missing]
report number
3
ISSN
1403-9338
language
English
LU publication?
yes
id
4572c7b3-3b2c-41cb-89f5-1b32d9ea76f3 (old id 1581110)
alternative location
http://www.math.ntnu.no/~hrue/r-inla.org/papers/S5-2010.pdf
date added to LUP
2016-04-04 09:23:46
date last changed
2019-03-25 13:53:53
@techreport{4572c7b3-3b2c-41cb-89f5-1b32d9ea76f3,
  abstract     = {{Continuously indexed Gaussian fields (GFs) is the most important ingredient in spatial statistical modelling and geo-statistics. The specification through the covariance function gives an intuitive interpretation of its properties. On the computational side, GFs are hampered with the "big-n" problem, since the cost of factorising dense matrices is cubic in the dimension. Although the computational power today is all-time-high, this fact seems still to be a computational bottleneck in applications. Along with GFs, there is the class of Gaussian Markov random fields (GMRFs) which are discretely indexed. The Markov property makes the involved precision matrix sparse which enables the use of numerical algorithms for sparse matrices, that for fields in R^2 only use the square-root of the time required by general algorithms. The specification of a GMRF is through its full conditional distributions but its marginal properties are not transparent in such a parametrisation.<br/><br>
<br/><br>
In this paper, we show that using an approximate stochastic weak solution to (linear) stochastic partial differential equations (SPDEs), we can, for some GFs in the Matérn class, provide an<br/><br>
explicit link, for any triangulation of R^d, between GFs and GMRFs. The consequence is that we can take the best from the two worlds and do the modelling using GFs but do the computations using GMRFs. Perhaps more importantly, our approach generalises to other covariance functions generated by SPDEs, including oscillating and non-stationary GFs, as well as GFs on manifolds. We illustrate our approach by analysing global temperature data with a non-stationary model defined on a sphere.}},
  author       = {{Lindgren, Finn and Lindström, Johan and Rue, Håvard}},
  institution  = {{[Publisher information missing]}},
  issn         = {{1403-9338}},
  keywords     = {{Gaussian Markov random fields; Generalised additive mixed models; Laplace approximation; Sparse matrices; Stochastic partial differential equations; Approximate Bayesian inference; Parallel computing; Structured additive regression models.}},
  language     = {{eng}},
  number       = {{3}},
  series       = {{Preprints in Mathematical Sciences}},
  title        = {{An explicit link between Gaussian fields and Gaussian Markov random fields; The SPDE approach}},
  url          = {{https://lup.lub.lu.se/search/files/5312973/1581115.pdf}},
  year         = {{2010}},
}