For the Poisson variance function , the code computes the Pearson chi-squared statistic
and the Pearson dispersion estimator
The first quantity measures goodness of fit; the second estimates the dispersion parameter, which equals one in a correctly specified Poisson regression. A standard rough calculation substitutes the residual deviance for the Pearson statistic and gives
If the reported upper-tail probability is inverted numerically, the actual Pearson statistic used by the code is about , giving . Either calculation reveals severe overdispersion.
The second fit is a Quasi-Poisson regression. It retains the logarithmic link function and mean model
but assumes only for an unknown dispersion parameter , rather than a complete Poisson distribution.