Fits non-crossing regression quantiles as a function of linear covariates and multiple smooth terms via B-splines with L1-norm difference penalties. The smoothing parameters are estimated as part of the model fitting. Monotonicity and concavity constraints on the fitted curves are allowed. See Muggeo, Sciandra, Tomasello and Calvo (2013)
*
*