计算具有可变等级数的R中的等级间可靠性?
Wikipedia建议,查看评估者之间可靠性的一种方法是使用随机效应模型来计算类内相关性。类内相关的例子讨论了看 σ2ασ2α+σ2ϵσα2σα2+σϵ2\frac{\sigma_\alpha^2}{\sigma_\alpha^2+\sigma_\epsilon^2} 从模型 Yij=μ+αi+ϵijYij=μ+αi+ϵijY_{ij} = \mu + \alpha_i + \epsilon_{ij} “其中Y ij是第 i 组的第 j 个观测值,μ是未观察到的总体均值,αi是组i中所有值共享的未观察到的随机效应,而εij是未观察到的噪声项。” 这是一个有吸引力的模型,尤其是因为在我的数据中,没有任何评分者对所有事物进行了评分(尽管大多数人的评分为20+),并且事物的评分次数是可变的(通常为3-4)。 问题#0:在该示例中,“组i”(“组i”)是否是一组被评级的事物? 问题#1:如果我正在寻找评估者之间的可靠性,我是否不需要一个包含两个术语的随机效应模型,一个用于评估者,一个用于评估的事物?毕竟,两者都有可能发生变化。 问题2:如何最好地用R表达此模型? 这个问题似乎有一个不错的建议: lmer(measurement ~ 1 + (1 | subject) + (1 | site), mydata) 我看了几个 问题,而lme的“ random”参数的语法对我来说是不透明的。我阅读了lme的帮助页面,但是没有示例,我对“随机”的描述是难以理解的。 这个问题有点类似于一个长 名单的问题,与此最接近的一次。但是,大多数都没有详细介绍R。