Back To Index Previous Article Next Article Full Text

Statistica Sinica 29 (2019), 2205-2227

A MATRIX-FREE METHOD FOR SPATIAL-
TEMPORAL GAUSSIAN STATE-SPACE MODELS
Debashis Mondal and Chunxiao Wang
Oregon State University

Abstract: This study develops a scalable matrix-free h-likelihood method for spatial- temporal Gaussian state-space models. The state vectors are constructed in such a way that they follow spatial-temporal Gaussian autoregressions that are consistent with the conditional formulation of auto-normal spatial fields. The proposed h-likelihood method provides the same inferences as those obtained from the Kalman filter and residual maximum likelihood analyses. However, for data from a large number of spatial sites, our method has significant computational advantages. Furthermore, we describe inferences in small time steps and indicate how the proposed method can be adapted to other complex spatial-temporal dynamical models based on stochastic partial differential equations. The proposed method applies to data with regularly or irregularly sampled spatial locations. Lastly, we illustrate our method by means of a simulation study and a data example on atmospheric concentrations of total nitrate across eastern North America.

Key words and phrases: Advection-diffusion, conditional autoregression, discrete cosine transform, Gaussian Markov random field, H-likelihood, incomplete Cholesky, Kalman filter, Lanczos algorithm, Residual likelihood, stochastic partial differential equation, trust region.

Back To Index Previous Article Next Article Full Text