method
NegativeBinomial.
__call__
Negative binomial variance function
mean parameters
variance = mu + alpha*mu**2