theta parameter, using
gam() and gam.nbut(). The code is a minor modification of
negative.binomial from library(MASS).neg.binom(theta=NA, link="log")theta, or NA
if this is unknown.log, sqrt or identity.gam() to fit a
generalised additive mode with Negative Binomial errors.gam, gam.nbut