Displaying 1 result from an estimated 1 matches for "burnagecat1".
2009 Jul 16
0
how to get means and confidence limits after glmmPQL or lmer
...e if I've
backtransformed correctly.
I've also tried directly calculating the means and CIs, but I'm not
sure if the first line below is actually the right variance to use in
estimating the CIs.
variance<-mm$sigma #IS THIS THE RIGHT VARIANCE?
speciessize<- tapply(species, burnagecat1 ,length)
speciesmeans<- (tapply(predict(mm, type="response"), burnagecat1 ,mean))
upperci<- exp(log(speciesmeans) +
qt(0.95,residdf)*sqrt(variance/speciessize))
lowerci<- exp(log(speciesmeans)- qt(0.95,residdf)*sqrt(variance/speciessize))
result<-cbind(lowerci, speciesme...