我正在尝试使用'simr'软件包创建用于电源仿真的多级数据结构。
数据有5栏：学生，班级，外向性，教师经验，学生性别。小学生嵌套在课堂上，外向性（中位居中，在-4.21和4.78之间），性别是学生水平的预测指标，而老师经验是班级中的预测指标（中位，在-12至11之间）。我希望数据框包含600个观察值（学生），在30个类中的每一个中嵌套20个。
我可以使用以下代码创建数据框，但是从美学角度来看，这是许多标准所不能接受的。另一个问题是，我希望学生的外向性和教师经验可以正常分布，我的代码似乎输出的范围更小。 / p>
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)