hi Riddle, You subscript is out of bounds because this line: for (j in seq(20,100,20)){
is incorrect - it is trying to index the 20th, 40th, 60th... column of a matrix that has only 5 columns. Try for (j in 1:5) instead. I'm not sure what the purpose of your script is so can't comment on what else needs to be done... Matt On 9/17/07, Riddle Chin <[EMAIL PROTECTED]> wrote: > Hi, everyone: > R is new to me. I am writing a nested loop to simulate data for t-test. > The following code is wrong. The subscript is out of bounds. Could anyone > tell me how to revise it? Thanks, Riddle Chin. > > result<-matrix(ncol=5, nrow=1000) > colnames(result)<-c('N=20','N=40','N=60','N=80','N=100') > for (i in 1:1000){ > for (j in seq(20,100,20)){ > x<-rnorm(j,2.7,1) > result[i,j]<-t.test(x,mu=4)$p.value<=0.05 > } > } > > > > --------------------------------- > Luggage? GPS? Comic books? > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Matthew C Keller Postdoctoral Fellow Virginia Institute for Psychiatric and Behavioral Genetics ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.