To cbind, you don't need a loop (M <- matrix(1:50,nc=10)) c2way <- combn(ncol(M),2) MM <- M[,c2way] dim(MM) <- c(nrow(M),nrow(c2way),ncol(c2way)) MM
c3way <- combn(ncol(M),3) MMM <- M[,c3way] dim(MMM) <- c(nrow(M),nrow(c3way),ncol(c3way)) MMM etc. etc. Than you can (untested for icc) loop your function as in apply(MMM,3,somefun) cheers On Thu, Mar 29, 2012 at 7:28 AM, Dai, Hongying, <h...@cmh.edu> wrote: > Dear R users, > > I'm wondering how I can generate an arbitrary number of loops in R. > For instance, I can generate two "for" loops to get ICC among any two-way > combination among 10 variables. Here is the code > > n<-10 > for (i in 1:(n-1)) > { > for (j in (i+1):n) > { > icc(cbind(DATA[,i],DATA[,j])) > } > } > If I need three-way combination, then a code with three "for" loops will be: > n<-10 > for (i in 1:(n-2)) > { > for (j in (i+1):(n-1)) > { > for (k in (j+1):n) > { > icc(cbind(DATA[,i],DATA[,j],DATA[,k])) > } > } > } > But how can I write a code if I need all m=2, 3, 4,... loops for arbitrary > m-way combinations? > Thanks! > Daisy > > > ________________________________ > Electronic mail from Children's Mercy Hospitals and Clinics. This > communication is intended only for the use of the addressee. It may contain > information that is privileged or confidential under applicable law. If you > are not the intended recipient or the agent of the recipient, you are hereby > notified that any dissemination, copy or disclosure of this communication is > strictly prohibited. If you have received this communication in error, please > immediately forward the message to Children's Mercy Hospital's Information > Security Officer via return electronic mail at > informationsecurityoffi...@cmh.edu and expunge this communication without > making any copies. Thank you for your cooperation. > > [[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. ______________________________________________ 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.