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