Thanks, all! I've tried changing the start and end of the loop and I tried cbind initial column with the other columns to get the matrix I want. It did work! Thanks also for the note about the parentheses. It real useful. I'll try Oarray if that would fasten the iteration. Thanks again!
Deana PhD Candidate On Oct 24, 2011, at 19:02, "David Winsemius" <dwinsem...@comcast.net> wrote: > > On Oct 24, 2011, at 7:23 PM, Md Desa, Zairul Nor Deana Binti wrote: > >> Hello, >> Does anyone knows how to deal with zero subscript in R. I have this code: >> for (i in 1:nitems){ >> + for (j in 1:ncat-1) { >> + draw<-matrix(rnorm(nitems*(ncat-1),seed1,seed2),nitems,(ncat-1)) >> + d<-( sigma_d*draw ) + mu_d >> + draw<-matrix(rtnorm((nitems*(ncat-1)),mean = seed1, sd = seed2, >> lower = .1, upper = 1.5),nitems,(ncat-1)) >> + d<-(sigma_d*draw) + mu_d >> + write.matrix(cbind(b.i0,d), file = "F:/KU/MIRT >> group/MIMIC-DIF/R/cpprcode/b0d.dat", sep = " ") >> + b[i,j]<-b[i,j-1]+d[i,j] >> + } >> + } >> >> The error as following: >> >> Error in b[i, j ] <- b[i, j - 1] + d[i, j] : >> replacement has length zero >> >> I would like to to have a matrix where the first column takes from initial >> pre-assigned value (b[i,0]+d1), the second column is additive from the first >> column and a constant d2 that is (b(i,1)=b(i,0)+d1+d2, and so forth. Is >> there any way that R can read subscript of zero? > > I don't. (I seem to remember the Matrix package allows such shenanigans.) > > Why don't you instead augment b[,] with a top column of initial values, using > rbind() and then iterate from 2:(ncat-2). (Do note the parentheses and if you > do not know why, then consult ?Syntax for operatr precendence.) > > -- > David Winsemius, MD > West Hartford, CT > > ______________________________________________ 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.