Fitting an interval response variable using gamlss.family distributions
