mlrcs function

Fit multi-level restricted cubic spline model