This function plots the histogram and a fitted (GAMLSS family) distribution to a variable
Useful links