A STOCHASTIC ESTIMATOR OF THE TRACE OF THE INFLUENCE MATRIX FOR LAPLACIAN SMOOTHING SPLINES
A STOCHASTIC ESTIMATOR OF THE TRACE OF THE INFLUENCE MATRIX FOR LAPLACIAN SMOOTHING SPLINES
复制标题
DOI:
10.1080/03610919008812866
复制
发表时间:
1990-01-01
影响因子:
0.9
通讯作者:
HUTCHINSON, MF
中科院分区:
文献类型:
--
作者:
HUTCHINSON, MF
An unbiased stochastic estimator of tr(I–A), where A is the influence matrix associated with the calculation of Laplacian smoothing splines, is described. The estimator is similar to one recently developed by Girard but satisfies a minimum variance criterion and does not require the simulation of a standard normal variable. It uses instead simulations of the discrete random variable which takes the values 1, -1 each with probability 1/2. Bounds on the variance of the estimator, similar to those established by Girard, are obtained using elementary methods. The estimator can be used to approximately minimize generalised cross validation (GCV) when using discretized iterative methods for fitting Laplacian smoothing splines to very large data sets. Simulated examples show that the estimated trace values, using either the estimator presented here or the estimator of Girard, perform almost as well as the exact values when applied to the minimization of GCV for n as small as a few hundred, where n is the number of data points.