search for: pupae_per_sub

...posts, it seems that my data would suit a tweedie approach. I can't use a zero-inflated Poisson model, because my data are not counts. 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...