Hi See in text. > > Dear Sir/Madam, > > I'm getting a problem with a R-code which converts a data frame to a matrix. > > It first generate a (m^(n-m) * m) matrix A and then regenerate another > matrix B having less dimension than A which satisfy some condition. Now I > wish to assign each row of B to a vector as individual. > > My problem is when I set any choice of (n,m) except m=1 it works fine but > setting m=1 I got the error : Error in B[i, ] : incorrect number of > dimensions. > > Moreover if (n,m) is large (say, (20,8)) I got the error : Error: cannot > allocate vector of size 3.0 Gb. I know this is due to large dimension of > matrix A. How to solve this problem. > > My code is given below: > > ********************************************************************** > > n=5 > m=3 > R=numeric(0) > # Generate all possible m-tuple ( variables having range 0 to n ) in a ( > m^(n-m) * m ) matrix > > r = expand.grid(rep(list(0:(n-m)), m)) > > write.table(r,file="test.txt",row.names=FALSE,col.names=FALSE) > > a= read.table(file="test.txt",sep="",header=FALSE)
Above lines do not do any sensible things. r shall be same as a. > > A= data.matrix(a) > > #......................................................................................... > > # Generate matrix whose rowsum = n-m > > meet.crit = apply(A, 1, function(.row) any((sum(.row)) == n-m)) # > criteron for being rowsum = n No error > > cbind(A, meet.crit) # > Checking rowsum = n for each row > -m > B=A[meet.crit,] # No error > Generate matrix > > #......................................................................................... > > > for(i in 1:choose(n-1,m-1)){ > R=B[i,] > } No error. However in each cycle only ith row is added to R and therefore only last row (in this case B[6,]) is added and stays in R. Either you need to use R <- c(R, B[i,]) in your construction or better as B is matrix > class(B) [1] "matrix" you can transform it to vector easily by stripping dimensions. R<-t(B) > dim(R) <-NULL > R [1] 2 0 0 1 1 0 0 2 0 1 0 1 0 1 1 0 0 2 Regards Petr > > *************************************************************************** > > Can you please help me how to get rid of these errors. Thanking you in > advance. > > Regards > > Ritwik Bhattacharya > > > Senior Research Fellow > SQC & OR UNIT, KOLKATA > INDIAN STATISTICAL INSTITUTE > > Voice : +91 9051253944 > > This mail is scanned by Ironport > > ______________________________________________ > 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. ______________________________________________ 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.