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.
