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

Resolving the '‘==’ only defined for equally-sized data frames' Error in Generalized Additive Models with gratia in R