A Kronecker-based Covariance Specification for Spatially Continuous Multivariate Data
Abstract We propose a covariance specification for modeling spatially continuous multivariate data. This model is based on a reformulation of Kronecker’s product of covariance matrices for Gaussian random fields. We illustrate the case with the Matérn function used for specifying marginal covariances. The structure holds for other choices of covariance functions with parameters varying in their usual domains, which makes the estimation process more accessible. The reduced computational time and flexible generalization for increasing number of variables, make it an attractive alternative for modelling spatially continuous data. Theoretical results for the likelihood function and the derivatives of the covariance matrix are presented. The proposed model is fitted to the literature’s soil250 dataset, and adequacy measures, forecast errors and estimation times are compared with the ones obtained based on classical models. Furthermore, the model is fitted to the classic meuse dataset to illustrate the model’s flexibility in a four-variate analysis. A simulation study is performed considering different parametric scenarios to evaluate the asymptotic properties of the maximum likelihood estimators. The satisfactory results, its simpler structure and the reduced estimation time make the proposed model a candidate approach for multivariate analysis of spatial data.