У меня возникли проблемы с настройкой 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))
Я получаю сообщение об ошибке:
Ошибка: количество уровней каждого фактора группировки должно быть <количество наблюдений. </p>
Я не уверен, есть ли другие параметры, которые необходимо установить в функции gstudy
.
Подводя итог, у меня есть два вопроса:
- Как написать формулу для студентов, скрещенных с предметами и вложенных в оценщики?
- Есть ли дополнительные параметры, которые мне нужно включитьв функции
gstudy
?
PS: я не использовал этот пакет gtheory раньше.