rb function

Creates a model using penalized smoothing splines