S=20

N=500

C=3

J=12

LCProb=N/3

RT=5

RespProb=0.7

p.x=0.8

varJX = matrix(,N,J)

LC= matrix(,N,1)

X=matrix(,N,1)

for(s in 1:S){

for(a in 1:J){

for(i in 1:N){

if(i <=LCProb){

LC=1

varJX[i,a]<-sample(c(1,2,3,4,5),size=1,replace=T,prob=c(0.8,0.05,0.05,0.05,0.05))

X<-sample(c(0,1),size=1,replace=T,prob= c(0.9,0.1))

}else{

if(i <= 2* LCProb){

LC=2

varJX[i,a]<-sample(c(1,2,3,4,5),size=1,replace=T,prob=c(0.05,0.8,0.05,0.05,0.05))

X<-sample(c(0,1),size=1,replace=T,prob= c(0.7,0.3))

}else{

LC=3

varJX[i,a]<-sample(c(1,2,3,4,5),size=1,replace=T,prob= c(0.05,0.05,0.8,0.05,0.05))

X<-sample(c(0,1),size=1,replace=T,prob= c(0.5,0.5))

}

}

}

}

varJX＝result

d <- data.frame(varJX,X)

sim <- paste("sim", s, ".txt",sep = "")

write.table(d, file = sim,row.names=F)

}这是我的代码，想生成20批随机数据，但是多加了一个循环语句就出错了。求大神指点！