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

我在为以下数据进行G研究时遇到了一些问题：

我一直收到一条错误消息：

错误：每个分组因子的级别数必须小于观察数

我不确定是否需要在

gstudy

功能中设置其他参数

总而言之，我有两个问题：

你是如何为交叉项目和嵌套在评分员中的学生编写公式的

是否需要在

gstudy

函数中包含其他参数

---

注：我以前没有用过这个。

相关的，可能重复的不是我要问的。相关的，可能重复的不是我要问的。

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))