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