r - 具有嵌套和交叉效应的 G-Study

Question

我在为以下数据设置 G-study 时遇到了一些问题”：

rater<- rep(1:4,each=12)
stu <- rep(1:8,each=6)
item<-rep(1:6,8)
score<- sample(2:4,48,replace=T)
dat<-as.data.frame(cbind(rater,stu,item,score))

学生与项目交叉并嵌套在评分者中。我假设这些影响暂时是随机的。我不知道如何编写回归公式。

我尝试了以下公式

library(gtheory)
gstudy(dat, score~(1|rater/stu:item))
gstudy(dat,score~(1|rater/item:stu))

我不断收到一条错误消息：

错误：每个分组因子的水平数必须 < 观察数。

我不确定函数中是否需要设置其他参数gstudy。

总结一下，我有两个问题：

1. 您如何为与项目交叉并嵌套在评分者中的学生编写公式？
2. 我需要在gstudy函数中包含其他参数吗？

PS：我之前没有使用过这个gtheory 包。

标签： rnestedpackagerandom-effectsgeneralization