Q: mgcv::gam() with constraints on the coefficents

Question

I'm fitting a logit GAM model and need to impose an inequality constraints on beta's in the smooth terms. I've tried mgcv::pcls() but it seems I can't specify the independent inequalities for the different spline coefficients. What's the recommended way to achieve this?

Answer 1

Found an answer myself, posting it just in case someone will look for it. Mono.con function in mgcv generates constraint matrices Ain and bin. For several splines, Ain matrices should be just put together into a single block-diagonal matrix.