Smoothing Parameter and Model Selection for General Smooth Models. Wood, S. N., Pya, N., & Säfken, B. Journal of the American Statistical Association, 111(516):1548–1563, 2016.
Smoothing Parameter and Model Selection for General Smooth Models [link]Paper  doi  abstract   bibtex   
This paper discusses a general framework for smoothing parameter estimation for models with regular likelihoods constructed in terms of unknown smooth functions of covariates. Gaussian random effects and parametric terms may also be present. By construction the method is numerically stable and convergent, and enables smoothing parameter uncertainty to be quantified. The latter enables us to fix a well known problem with AIC for such models. The smooth functions are represented by reduced rank spline like smoothers, with associated quadratic penalties measuring function smoothness. Model estimation is by penalized likelihood maximization, where the smoothing parameters controlling the extent of penalization are estimated by Laplace approximate marginal likelihood. The methods cover, for example, generalized additive models for non-exponential family responses (for example beta, ordered categorical, scaled t distribution, negative binomial and Tweedie distributions), generalized additive models for location scale and shape (for example two stage zero inflation models, and Gaussian location-scale models), Cox proportional hazards models and multivariate additive models. The framework reduces the implementation of new model classes to the coding of some standard derivatives of the log likelihood.
@article{Wood2015Smoothing,
 abstract = {This paper discusses a general framework for smoothing parameter estimation

for models with regular likelihoods constructed in terms of unknown smooth

functions of covariates. Gaussian random effects and parametric terms may also

be present. By construction the method is numerically stable and convergent,

and enables smoothing parameter uncertainty to be quantified. The latter

enables us to fix a well known problem with AIC for such models. The smooth

functions are represented by reduced rank spline like smoothers, with

associated quadratic penalties measuring function smoothness. Model estimation

is by penalized likelihood maximization, where the smoothing parameters

controlling the extent of penalization are estimated by Laplace approximate

marginal likelihood. The methods cover, for example, generalized additive

models for non-exponential family responses (for example beta, ordered

categorical, scaled t distribution, negative binomial and Tweedie

distributions), generalized additive models for location scale and shape (for

example two stage zero inflation models, and Gaussian location-scale models),

Cox proportional hazards models and multivariate additive models. The framework

reduces the implementation of new model classes to the coding of some standard

derivatives of the log likelihood.},
 author = {Wood, Simon N. and Pya, Natalya and S{\"a}fken, Benjamin},
 year = {2016},
 title = {Smoothing Parameter and Model Selection for General Smooth Models},
 url = {http://dx.doi.org/10.1080/01621459.2016.1180986},
 keywords = {phd;stat},
 pages = {1548--1563},
 volume = {111},
 number = {516},
 issn = {0162-1459},
 journal = {Journal of the American Statistical Association},
 doi = {10.1080/01621459.2016.1180986},
 howpublished = {refereed}
}

Downloads: 0