'Negative binomial variance function
Parameters
mu : array-like
mean parameters
Returns
variance : array
variance = mu + alpha*mu**2'

def __call__(self, mu):

    p = self._clean(mu)
    return (p + (self.alpha * (p ** 2)))
