Я пытаюсь создать многоуровневую структуру данных для имитации мощности, используя пакет 'simr'.
Данные имеют 5 столбцов: ученик, класс, экстраверсия, опыт учителя, пол ученика. Ученики вложены в класс, экстраверсия (по центру большого среднего, диапазон от -4,21 до 4,78), а пол - предикторы уровня ученика, а опыт учителя - предиктор по уровню класса (среднее по центру, диапазон от -12 до 11). Я хочу, чтобы план данных содержал 600 наблюдений (учеников), по 20 вложенных в каждый из 30 классов.
Я могу создать фрейм данных, используя приведенный ниже код, но он эстетически неприемлем по многим стандартам. Еще одна проблема заключается в том, что я надеюсь, что экстраверсия ученика и опыт учителя могут быть нормально распределены, кажется, мой код выводит меньший диапазон ..
library(truncnorm)
set.seed(123)
pupil <- as.factor(seq(1:600))
class <- as.factor(rep(1:30,20))
extrav<-rtruncnorm(n=600,a=-4.21,b=4.78)
X <- cbind(pupil=pupil,class=class, extrav=extrav)
class<-seq(1:30)
texp<-rtruncnorm(30,-12,11)
class_texp<-cbind(class,texp)
X<-merge(X,class_texp)
sex<-rep(1:2,300)
sex<-sex-1
X<-cbind(X,sex)
View(X)