[R] How to optimize this codes ?

Daren Tan daren76 at hotmail.com
Thu Dec 4 17:02:49 CET 2008


 
How to optimize the for-loop to be reasonably fast for sample.size=100000000 ? You may want to change sample.size=1000 to have an idea what I am achieving. 
 
 
set.seed(143)
A <- matrix(sample(0:1, sample.size, TRUE), ncol=10, dimnames=list(NULL, LETTERS[1:10]))
 
B <- list()
for(i in 1:10) { 
  B[[i]] <- apply(combn(LETTERS[1:10], i), 2, function(x) { sum(apply(data.frame(A[,x]), 1, all)) })
}
 
pdf("test.pdf")
boxplot(B, xlab = "Number of predictors", ylab = "Number of agreements across predictors")
dev.off()



More information about the R-help mailing list