R循环替换矩阵中的多个值答案

【问题标题】：R loop replace multiple values in matrixR循环替换矩阵中的多个值
【发布时间】：2013-10-27 15:39:21
【问题描述】：

有谁知道执行以下操作的更有效方法？我有两个矩阵，一个带有整数值，另一个带有数字。我还有一个包含整数向量的列表。当列表中的整数值等于整数矩阵时，我有一个循环将数值矩阵的值设置为 NA。然后我从结果矩阵中得到行积和总和。有没有办法避免每一步都创建数字矩阵的副本？或者也许是另一种方法？谢谢

mat1<-matrix(rpois(20*300000,6),20,300000)
mat2<-matrix(runif(20*300000),20,300000)
list1<-list(c(1,2,3),c(4,6),c(8,9,10,11))
results<-vector('numeric',length(list1))

start.time=Sys.time()
for(i in 1:length(list1)){
    copy<-mat2
    copy[mat1 %in% list1[[i]]]=NA
    results[i]=sum(apply(copy,2,prod,na.rm=T))
}
print(Sys.time()-start.time)
#

【问题讨论】：

还没有完整地浏览过你的代码，但只要快速浏览一下，?outer 可能对你有帮助
mat2 中的数字是否与您的示例一样都是正数？如果是这种情况，您可以定义 colProds <- function(x, na.rm = TRUE) exp(colSums(log(x), na.rm = na.rm)) 以利用 colSums 比 apply 快得多的优势。

标签： r matrix replace

【解决方案1】：

用colSums 替换你的apply（假设你的数字是正数，否则你需要更多的摆弄）让我的速度提高了近2倍：

for(i in 1:length(list1)){
    copy<-mat2
    copy[mat1 %in% list1[[i]]]=NA
    results[i]=sum(exp(colSums(log(copy), na.rm = T)))
}

【讨论】：

他们永远是积极的，谢谢。