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) Semiparametric Poisson-Gamma models: a roughness penalty approach. MSc Dissertation. Rio de Janeiro, PUC-Rio, Department of Electrical Engineering.
pgam, predict.pgam