期刊
JOURNAL OF THE AMERICAN STATISTICAL ASSOCIATION
卷 103, 期 484, 页码 1545-1555出版社
AMER STATISTICAL ASSOC
DOI: 10.1198/016214508000000959
关键词
Compactly supported correlation function; Covariance estimation; Estimating equation; Gaussian process
Maximum likelihood is an attractive method of estimating covariance parameters in spatial models based oil Gaussian processes. But calculating, the likelihood can he computationally infeasible for large data sets, requiring O(m(3)) calculations for a data set with a observations. This article proposes the method of covariance tapering to approximate the likelihood in this setting. In this approach. covariance matrixes are tapered. or multiplied element wise by a sparse correlation matrix. The resulting matrixes can their be manipulated using efficient sparse matrix algorithms. We propose two approximations to the Gaussian likelihood using tapering. One of the approximations simply replaces the model covariance with a tapered version, whereas the other is motivated by the theory of unbiased estimating equations. Focusing on the particular case of the Matern class of covariance functions, we give conditions under which estimators maximizing the tapering approximations are. like the maximum likelihood estimator, Strongly consistent. Moreover, we show in a simulation study that the tapering estimators can have sampling densities quite similar to that of the maximum likelihood estimator even when the degree of tapering is severe. We illustrate the accuracy and computational gains of the tapering methods in an analysis of yearly total precipitation anomalies at weather stations in the United States.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据