Published by the Foundation for Open Access Statistics
Editors-in-chief: Bettina Grün, Torsten Hothorn, Edzer Pebesma, Achim Zeileis    ISSN 1548-7660; CODEN JSSOBK
On Circulant Embedding for Gaussian Random Fields in R | Davies | Journal of Statistical Software
Authors: Tilman M. Davies, David Bryant
Title: On Circulant Embedding for Gaussian Random Fields in R
Abstract: The high-dimensionality typically associated with discretized approximations to Gaussian random fields is a considerable hinderance to computationally efficient methods for their simulation. Many direct approaches require spectral decompositions of the associated covariance matrix and so are unable to complete the solving process in a timely fashion, if at all. However under certain conditions, we may construct block-circulant versions of the covariance matrix at hand thereby allowing access to fast-Fourier methods to perform the required operations with impressive speed. We demonstrate how circulant embedding and subsequent simulation can be performed directly in the R language. The approach is currently implemented in C for the R package RandomFields, and used in the recently released package lgcp. Motivated by applications dealing with spatial point processes we restrict attention to stationary Gaussian fields on R2, where sparsity of the covariance matrix cannot necessarily be assumed.

Page views:: 2221. Submitted: 2012-01-16. Published: 2013-11-13.
Paper: On Circulant Embedding for Gaussian Random Fields in R     Download PDF (Downloads: 2927)
v55i09.R: R example and replication code Download (Downloads: 284; 21KB)

DOI: 10.18637/jss.v055.i09

This work is licensed under the licenses
Paper: Creative Commons Attribution 3.0 Unported License
Code: GNU General Public License (at least one of version 2 or version 3) or a GPL-compatible license.