b_rcs function

Basis for restricted cubic splines