Wrapper functions to handle stability selection (with gradient boosting) of bamlss models.
Bamlss families for modelling multivariate Gaussian distributions with Cholesky decomposition of the covariance matrix.