具有嵌套和交叉效果的G-研究

Question

我无法为以下数据设置G研究”：

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函数中是否需要包含其他参数？

P.S .：我以前没有使用过这个gtheory package。