Displaying 1 result from an estimated 1 matches for "residdf".
Did you mean:
residf
2009 Jul 16
0
how to get means and confidence limits after glmmPQL or lmer
...lly 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, speciesmeans, upperci)
result
I also explored lmer, but haven't found a way to extract means and
CIs on the original scale for the levels of treatment.
li...