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

Fixing Random Effects Issues in Multilevel Modeling with mgcv: A Simple Solution