R:如何编写一个for循环来读取矩阵中的每两行？
R:如何编写一个for循环来读取矩阵中的每两行？
提问于 2014-11-05 11:05:59
我想使用cor.test()计算相关统计信息。我有一个数据矩阵，其中两个待测试的对在连续的行上(我有超过1000对，所以我需要稍后对其进行更正)。我想我可以循环遍历矩阵中的每两行并执行测试(即首先测试row1和row2之间的相关性，然后测试row3和row4、row5和row6等)，但我不知道如何进行这种循环。
这就是我如何在单个对上进行测试：
d = read.table(file="cor-test-sample-data.txt", header=T, sep="\t", row.names = 1)
d = as.matrix(d)
cor.test(d[1,], d[2,], method = "spearman")回答 3
Stack Overflow用户
发布于 2014-11-05 11:21:04
你可以试一试
res <- lapply(split(seq_len(nrow(mat1)),(seq_len(nrow(mat1))-1)%/%2 +1),
function(i){m1 <- mat1[i,]
if(NROW(m1)==2){
cor.test(m1[1,], m1[2,], method="spearman")
}
else NA
})获取p-values
resP <- sapply(res, function(x) x$p.value)
indx <- t(`dim<-`(seq_len(nrow(mat1)), c(2, nrow(mat1)/2)))
names(resP) <- paste(indx[,1], indx[,2], sep="_")
resP
# 1_2 3_4 5_6 7_8 9_10 11_12 13_14
#0.89726818 0.45191660 0.14106085 0.82532260 0.54262680 0.25384239 0.89726815
# 15_16 17_18 19_20 21_22 23_24 25_26 27_28
#0.02270217 0.16840791 0.45563229 0.28533447 0.53088721 0.23453161 0.79235990
# 29_30 31_32
#0.01345768 0.01611903 或者使用mapply (假设行是偶数)
ind <- seq(1, nrow(mat1), by=2) #similar to the one used by @CathG in for loop
mapply(function(i,j) cor.test(mat1[i,], mat1[j,],
method='spearman')$p.value , ind, ind+1) 数据
set.seed(25)
mat1 <- matrix(sample(0:100, 20*32, replace=TRUE), ncol=20) Stack Overflow用户
发布于 2014-11-05 11:22:56
试一试
d = matrix(rep(1:9, 3), ncol=3, byrow = T)
sapply(2*(1:(nrow(d)/2)), function(pair) unname(cor.test(d[pair-1,], d[pair,], method="spearman")$estimate))Stack Overflow用户
发布于 2014-11-05 11:52:12
pvalues<-c()
for (i in seq(1,nrow(d),by=2)) {
pvalues<-c(pvalues,cor.test(d[i,],d[i+1,],method="spearman")$p.value)
}
names(pvalues)<-paste(row.names(d)[seq(1,nrow(d),by=2)],row.names(d)[seq(2,nrow(d),by=2)],sep="_")页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/26755653
复制相关文章
相似问题
编写一个循环来对R中的矩阵进行切片
循环来填充R中的矩阵。
R中的for循环来创建矩阵
如何编写一个for循环来遍历R中的列
如何在R中使用"for“循环编写矩阵？
添加站长 进交流群
领取专属 10元无门槛券
AI混元助手 在线答疑
关注 腾讯云开发者公众号
洞察 腾讯核心技术
剖析业界实践案例
社区富文本编辑器全新改版！诚邀体验～
全新交互，全新视觉，新增快捷键、悬浮工具栏、高亮块等功能并同时优化现有功能，全面提升创作效率和体验
腾讯云开发者
