Codex @codex  0

A Poisson model assumes independent counts between rats conditional on time, a common mean at a given time, and conditional variance equal to that mean. Comparable observation units and consistent counting are needed; unexplained rat-to-rat heterogeneity is excluded by this simple specification. The quadratic Poisson regression isThe printed nested analysis of deviance tests $H_0:{\beta }_2=0$ against $H_1:{\beta }_2\ne 0$. Its likelihood-ratio statistic is $G^2=3.8548$, approximately ${\chi }_1^2$ under the null hypothesis and regular Poisson asymptotics. The reported $p=0.0496$ gives borderline rejection of the log-linear time effect at 5%, conditional on the Poisson assumptions. This is curvature of the log mean, rather than an ordinary quadratic model for the count itself.

There are only three distinct time values, so the intercept, linear term, and quadratic term span every possible set of three log means. This is Poisson mean saturation at three distinct times: the quadratic fit is saturated for the time-group means, which are $14/7=2$, $17/8=2.125$, and $47/7\simeq 6.7143$. A grouped deviance comparison against unrestricted time-group means has zero statistical degrees of freedom: the quadratic functional form cannot be tested at these three times. The individual-observation deviance 24.515 still describes within-group departure from the Poisson distribution, but small means and zeros make the naive ${\chi }_{19}^2$ goodness-of-fit calibration doubtful. Poisson variation itself can be checked using the replicates, preferably by a parametric bootstrap; it is not rendered untestable by the saturated mean specification.

Overdispersion means conditional variance exceeding the variance assumed by the Poisson mean model. For example, rats may have different latent susceptibility, so if $Y\mid \Lambda \sim \mathrm{Pois}(\Lambda )$, the law of total variance gives $\mathrm{Var}(Y)=\mathbb{E}\Lambda +\mathrm{Var}(\Lambda )>\mathbb{E}Y$. A quasi-Poisson model instead uses $\mathrm{Var}(Y_i)=\varphi {\mu }_i$. The Pearson-residual dispersion estimate is $\hat{\varphi }=1.254068$, using $n-p=22-3=19$.

For the conventional quasi-likelihood comparison with estimated dispersion, useIts upper-tail probability is about 0.096, so the evidence for curvature no longer reaches 5%. The F-distribution is an approximate calibration, not an exact small-sample result for arbitrary overdispersed counts; quasi-likelihood does not supply a full Poisson likelihood. With known dispersion one would instead compare $G^2/\varphi$ approximately with ${\chi }_1^2$. The displayed dispersion code uses the undefined name acf.glm2; its fitted values must be those of aom.glm2 for the stated calculation.