statsmodels.genmod.families.family.NegativeBinomial.resid_anscombe

NegativeBinomial.resid_anscombe(endog, mu, var_weights=1.0, scale=1.0) [source]

The Anscombe residuals

Parameters:	endog (array) – The endogenous response variable mu (array) – The inverse of the link function at the linear predicted values. var_weights (array-like) – 1d array of variance (analytic) weights. The default is 1. scale (float, optional) – An optional argument to divide the residuals by sqrt(scale). The default is 1.
Returns:	resid_anscombe – The Anscombe residuals as defined below.
Return type:	array

Notes

Anscombe residuals for Negative Binomial are the same as for Binomial upon setting \(n=-\frac{1}{\alpha}\). Due to the negative value of \(-\alpha*Y\) the representation with the hypergeometric function \(H2F1(x) = hyp2f1(2/3.,1/3.,5/3.,x)\) is advantageous

\[resid\_anscombe_i = \frac{3}{2} * (Y_i^(2/3)*H2F1(-\alpha*Y_i) - \mu_i^(2/3)*H2F1(-\alpha*\mu_i)) / (\mu_i * (1+\alpha*\mu_i) * scale^3)^(1/6) * \sqrt(var\_weights)\]

Note that for the (unregularized) Beta function, one has \(Beta(z,a,b) = z^a/a * H2F1(a,1-b,a+1,z)\)