G-Study с вложенными и перекрестными эффектами

У меня возникли проблемы с настройкой 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.