Dear R-list, let us say I want to report the reliability of group means for "HRS" (ICC2):
library(multilevel) data("bh1996") mult.icc(bh1996[,c("HRS","LEAD","COHES")],grpid=bh1996$GRP) Variable ICC1 ICC2 1 HRS 0.129237 0.91713 2 LEAD 0.147462 0.92804 3 COHES 0.048049 0.79008 the result is 0.91713. However, if I use GmeanRel I get a different result: tmod<-lme(HRS~1,random=~1|GRP,data=bh1996) mean((GmeanRel(tmod))$MeanRel) [1] 0.89586 What is the difference between both estimates? Thanks for helping me! Best wishes Alain [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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.