bkfsmooth(y, x, df, smoother = "spline", w = rep(1, length(y)))NULL the smoothing parameter is selected by cross-validationNowadays, this function interfaces the smooth.spline in stats library. It will become not dependent soon enough.
Hastie, T. J., Tibshirani, R. J.(1990) Generalized Additive Models. Chapman and Hall, London
Junger, W. L. (2004) Modelo Poisson-Gama Semi-Paramétrico: Uma Abordagem de Penalização por Rugosidade. MSc Thesis. Rio de Janeiro, PUC-Rio, Departamento de Engenharia Elétrica
pgam, predict.pgam