Displaying 1 result from an estimated 1 matches for "gsa_total_df_persub".
2010 Sep 30
1
AIC for tweedie glm
...nts. Many of my
values are between 0 and 1, so if I rounded to the nearest integer, I'd lose
a lot of the variation.
Here's my code:
out<-tweedie.profile(PUPAE_PER_SUB~1,p.vec=seq(1.1,1.9,length=9),method="interpolation",do.ci=TRUE,do.smooth=TRUE,do.plot=TRUE)
tweedie1<-glm(GSA_TOTAL_DF_PERSUB~GROUP_SIZE+PERIOD+SITE+PERIOD*GROUP_SIZE,family=tweedie(var.power=out$p.max,link.power=0))
This worked fine, and gave results I expected, but I don't know what the
best method is to evaluate the fit of the model. I am used to using AIC to
compare models. A site search turned up AICtweedie, wit...