set.seed(123456)
reps<-20
n.iter<-188
n<-90
t.Z1_30=matrix(NA,nrow=600,ncol=1)
t.Z2_30=matrix(NA,nrow=600,ncol=1)
t.Z3_30=matrix(NA,nrow=600,ncol=1)
for(i in 1:n){
for(j in 1:reps){
`n<-90`
`nc<-3`
X<-rep(1:nc, each=n/nc)
b0<-5
b1<-0.3
Y<-b0+b1*X+rnorm(n,0,2)
}
group<-split(Y,X)
t.Z1_30[i]<-((group$'1')-5.3)/2
t.Z2_30[i]<-((group$'2')-5.6)/2
t.Z3_30[i]<-((group$'3')-5.9)/2
t.Z.frame<-rbind(t.Z1_30, t.Z2_30, t.Z3_30)
}
我想在r中模拟我的设计。我的IV是(样本中的平衡,不平衡)和(子样本大小:30,50,100,150,200)。我是r的初学者。
所以我想在(blanced)和(子样本大小:30,子组= 3,reps = 20)的情况下设置我的模拟。所以我申请了loop()和嵌套for循环,我预计我会得到1800个模拟 (3 X 30 X 20)但我可以得到1800次模拟我只得到90次模拟。
我的代码有什么问题?