ISSN:
0945-3245
Keywords:
AMS(MOS): 65D10
;
65F35
;
CR: G.1.2
;
G.1.3
;
G.3
Source:
Springer Online Journal Archives 1860-2000
Topics:
Mathematics
Notes:
Summary We propose a fast Monte-Carlo algorithm for calculating reliable estimates of the trace of the influence matrixA τ involved in regularization of linear equations or data smoothing problems, where τ is the regularization or smoothing parameter. This general algorithm is simply as follows: i) generaten pseudo-random valuesw 1, ...,w n , from the standard normal distribution (wheren is the number of data points) and letw=(w 1, ...,w n ) T , ii) compute the residual vectorw−A τ w, iii) take the ‘normalized” inner-product (w T (w−A τ w))/(w T w) as an approximation to (1/n)tr(I−A τ). We show, both by theoretical bounds and by numerical simulations on some typical problems, that the expected relative precision of these estimates is very good whenn is large enough, and that they can be used in practice for the minimization with respect to τ of the well known Generalized Cross-Validation (GCV) function. This permits the use of the GCV method for choosing τ in any particular large-scale application, with only a similar amount of work as the standard residual method. Numerical applications of this procedure to optimal spline smoothing in one or two dimensions show its efficiency.
Type of Medium:
Electronic Resource
URL:
http://dx.doi.org/10.1007/BF01395775
Permalink