Fits non-crossing regression quantiles as a function of linear covariates and multiple smooth terms, including varying coefficients, via B-splines with L1-norm difference penalties.
Random intercepts and variable selection are allowed via the lasso penalties.
The smoothing parameters are estimated as part of the model fitting, see Muggeo and others (2021)