我想模拟我的设计,但效果不佳

时间:2016-04-18 14:12:52

标签: r for-loop simulation montecarlo

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次模拟。

我的代码有什么问题?

0 个答案:

没有答案