The computations are based on the approximate gradient of the restricted log likelihood. The outer Hessian w.r.t. the log smoothing parameters is computed by central differencing of this gradient, and the standard deviations of the smoothing parameters, or of their inverses (variances), are obtained from it using the delta method.