2
$\begingroup$

Can I use likelihood ratio chi-squared tests to test the goodness of fit for quasi-Poisson models?

$\endgroup$
2
  • $\begingroup$ It is possible (like, I think R will do it if you ask), but I don't think it is totally appropriate because the quasi poisson technically has no likelihood because it has no probability distribution. $\endgroup$ Commented Jul 23, 2019 at 17:56
  • $\begingroup$ Thanks @DemetriPananos. Yes, when I try to compute it in R it all works fine, which I found confusing. $\endgroup$
    – vic12
    Commented Jul 23, 2019 at 18:16

1 Answer 1

3
$\begingroup$

If I understand your question correctly: no. Goodness-of-fit tests against an unspecified alternative (e.g., testing either the sum of squares of Pearson residuals or the residual deviance against the upper tail of a $\chi^2_{n-p}$ distribution) depend on the assumption of a fixed mean-variance relationship with no adjustable scale parameters. Once you estimate a dispersion parameter (as in quasi-likelihood), you've essentially soaked up all of the lack of fit into a single parameter. You can compute a pseudo-$R^2$ of some sort, but these are not associated with a particular test.

With quasi-likelihood fits, you're in a similar situation to fitting a linear/ordinary least-squares model; there are plenty of things you can do to examine the goodness of fit (e.g. plot residuals vs fitted to look for bias; scale-location plot to look for heteroscedasticity; Cook's distance/leverage plot to look for outliers), and you can if you like test particular aspects of the fit (e.g. can you reject the hypothesis that the residuals are Normally distributed? [although many people consider this a bad idea]), but there is no overall goodness-of-fit test.

You'll notice that many statistical procedures rely on tests of relative goodness of fit (likelihood ratio test, Wald $t$/$F$ test, AIC and other information criteria ...) but they usually start with a rather vague statement that the initial fit (e.g. the most complex one) should be "good enough" ...

$\endgroup$
1
  • 2
    $\begingroup$ Thanks a lot! This helped me a lot in understanding more about quasi-likelihood fits and how to make sure that the model is a good fit. $\endgroup$
    – vic12
    Commented Jul 23, 2019 at 18:20

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.