STOCHASTIC APPROXIMATION OF SCORE FUNCTIONS FOR GAUSSIAN PROCESSES
STOCHASTIC APPROXIMATION OF SCORE FUNCTIONS FOR GAUSSIAN PROCESSES
复制标题
高斯过程分数函数的随机逼近
DOI:
--
复制
发表时间:
2013
期刊:
影响因子:
--
通讯作者:
M. Anitescu
中科院分区:
文献类型:
--
作者:
M. Stein;Jie Chen;M. Anitescu
We discuss the statistical properties of a recently introduced unbiased stochastic approximation to the score equations for maximum likelihood calculation for Gaussian processes. Under certain conditions, including bounded condition number of the covariance matrix, the approach achieves O(n) storage and nearly O(n) computational eort per optimization step, where n is the number of data sites. Here, we prove that if the condition number of the covariance matrix is bounded, then the approximate score equations are nearly optimal in a well-dened sense. Therefore not only is the approximation efcient to compute, but it also has comparable statistical properties to the exact maximum likelihood estimates. We discuss a modication of the stochastic approximation in which design elements of the stochastic terms mimic patterns from a 2 n factorial design. We prove these designs are always at least as good as the unstructured design, and we demonstrate through simulation that they can produce a substantial improvement over random designs. Our ndings are validated by numerical experiments on simulated datasets of up to 1 million observations. We apply the approach to t a space-time