我想将两个循环转换为一个应用函数,其中一个用于循环,希望这会加快我的计算速度。我知道使用apply并不能保证更快的计算,但我想试试看(也适合学习经验来熟悉apply)R:尝试使用apply来转换for循环
我想要做的是;
计算两个矩阵的每一行的皮尔森相关系数,也得到p值。
矩阵双方都有约3000 X 100
现在我的代码看起来是这样的,它已经运行了天维...
cnt <- 1;
res_row1 <- c();
res_row2 <- c();
res_corr <- c();
res_pval <- c();
for (i in (1:dim(m1)[1])) {
for (j in (1:dim(m2)[1])) {
c <- cor.test(as.matrix(m1[i,]), as.matrix(m2[j,]));
res_row1[cnt] <- rownames(m1)[i];
# need both row names in the output files
res_row2[cnt] <- rownames(m2)[j];
res_corr[cnt] <- c$estimate;
res_pval[cnt] <- c$p.value;
# Storing the results for output
cnt<-cnt+1;
}
comp <- (i/dim(m1[1]) * 100;
cat(sprintf("Row number of file 1 = %f | %f percent complete \n", i, comp))
}
results <- cbind(res_row1, res_row2, res_corr, res_pval)
你们能帮助我吗?
你的问题是不是(真的)'for'环vs'apply',而是你在[hell]的第二个圈子里(http://www.burns-stat.com/pages/Tutor/R_inferno.pdf)。 – joran 2013-03-05 15:52:13
谢谢!刚刚在几周前开始使用R。一个不错的开眼界:) – 2013-03-05 18:51:09