This routine estimates SCAM coefficients given log smoothing parameters using the Newton-Raphson method.
The estimation of the smoothing parameters by the GCV/UBRE score optimization is outer to the model fitting. Routine
gcv.ubre_grad evaluates the first derivatives of the smoothness selection scores with respect to the
log smoothing parameters. Routine bfgs_gcv.ubre estimates the smoothing parameters using the BFGS method.
The function is not normally called directly, but rather service routines for scam.
scam.fit(G, sp, maxit=200, maxHalf.fit=40, devtol.fit=1e-8, steptol.fit=1e-8,
gamma=1, start=NULL, etastart=NULL, mustart=NULL, env=env)A list of items needed to fit a SCAM.
The vector of smoothing parameters.
Maximum iterations in the Newton-Raphson procedure.
If a step of the Newton-Raphson optimization method leads to a worse penalized deviance, then the step length of the model coefficients is halved. This is the number of halvings to try before giving up.
A positive scalar giving the tolerance at which the scaled distance between two successive penalized deviances is considered close enough to zero to terminate the algorithm.
A positive scalar giving the tolerance at which the scaled distance between two successive iterates is considered close enough to zero to terminate the algorithm.
This constant allows to inflate the model degrees of freedom in the GCV or UBRE/AIC score.
Initial values for the model coefficients.
Initial values for the linear predictor.
Initial values for the expected values.
Get the enviroment for the model coefficients, their derivatives and the smoothing parameter.
The routine applies step halving to any step that increases the penalized deviance substantially.
Pya, N. and Wood, S.N. (2015) Shape constrained additive models. Statistics and Computing, 25(3), 543-559
Pya, N. (2010) Additive models with shape constraints. PhD thesis. University of Bath. Department of Mathematical Sciences
Wood, S.N. (2008) Fast stable direct fitting and smoothness selection for generalized additive models. Journal of the Royal Statistical Society (B) 70(3):495-518
Wood, S.N. (2011) Fast stable restricted maximum likelihood and marginal likelihood estimation of semiparametric generalized linear models. Journal of the Royal Statistical Society (B) 73(1):3-36