Adds the fits after fitting a natural cubic smoothing spline to subsets of a response to a data.frame
data.frame