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