Fitting Generalized Additive Models in the Negative Binomial Family Using R's Gamlss Package