Dear fellows,
I'm trying to extract the AIC statistic from a GLM model with quasipoisson link.
The formula I'm referring to is
AIC = -2(maximum loglik) + 2df * phi
with phi the overdispersion parameter, as reported in:
Peng et al., Model choice in time series studies os air pollution and
mortality. J R Stat Soc A, 2006; 162: pag 190.
Unfortunately, the function logLik doesn't work for a quasipoisson link.
Do you know a fast method to extract the AIC for these models?
Thanks in advance
Antonio Gasparrini
Public and Environmental Health Research Unit (PEHRU)
London School of Hygiene & Tropical Medicine
Keppel Street, London WC1E 7HT, UK
Office: 0044 (0)20 79272406 - Mobile: 0044 (0)79 64925523
http://www.lshtm.ac.uk/people/gasparrini.antonio (
http://www.lshtm.ac.uk/pehru/ )
[[alternative HTML version deleted]]
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
