zbMATH — the first resource for mathematics

A comparison of GCV and GML for choosing the smoothing parameter in the generalized spline smoothing problem. (English) Zbl 0596.65004
The special spline smoothing model to be considered is given by \(y_ i=f(t_ i)+\epsilon_ i\), \(i=1,...,n\), \(t_ i\in [0,1]\), where \(\epsilon =(\epsilon_ 1,...,\epsilon_ n)\sim N(0,\sigma^ 2I_{n\times n})\), \(\sigma^ 2\) is unknown, and f(\(\cdot)\) is some function in the Sobolev space \(W_ 2^ m[0,1]=\{f:\) \(f,f^ 1,...,f^{(m-1)}\) absolutely continuous, \(f^{(m)}\in L_ 2[0,1]\}\) and the smoothing spline estimate \(f_{n,\lambda}\) of f is the minimizer in \(W^ m_ 2[0,1]\) of \[ 1/n\sum^{n}_{i=1}(f(t_ i)-y_ i)^ 2+\lambda \int^{1}_{0}(f^{(m)}(t))^ 2dt. \] A generalization of the maximum likelihood (GML) estimate for the smoothing parameter \(\lambda\) is obtained, and this estimate is compared with the generalized cross validation (GCV) estimate both analytically and by Monte Carlo methods. The theoretical results are shown to extend to the generalized spline smoothing model, which includes the estimate of functions given noisy values of various integrals of them.
Reviewer: Y.Sun

65D10 Numerical smoothing, curve fitting
65R20 Numerical methods for integral equations
Full Text: DOI