Does someone know how to choose/fix the knots (breakpoints) of cs() or pb() smoothing functions in the gamlss R package. I am struggling with the code without any success. The knots remain definitively at the default values. Example(s) of code would be appreciated.