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

R循环替换矩阵中的多个值,r,matrix,replace,R,Matrix,Replace,有人知道一种更有效的方法来做以下事情吗?我有两个矩阵,一个是整数值,另一个是数值。我还有一个包含整数向量的列表。我有一个循环,当列表中的整数值等于整矩阵时,它将数值矩阵的值设置为NA。然后我从得到的矩阵中得到行积和和。有没有办法避免在每一步创建数字矩阵的副本?或者是另一种方法?谢谢 mat1<-matrix(rpois(20*300000,6),20,300000) mat2<-matrix(runif(20*300000),20,300000) list1<-list(c(1

有人知道一种更有效的方法来做以下事情吗?我有两个矩阵,一个是整数值,另一个是数值。我还有一个包含整数向量的列表。我有一个循环,当列表中的整数值等于整矩阵时,它将数值矩阵的值设置为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)
#

mat1用
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)))
}
for(i在1中:长度(列表1)){

Copy尚未完全阅读您的代码,但只需快速浏览一下,它看起来像是
?outer
可能会对您有所帮助
mat2
中的数字是否与您的示例中的数字一样都是正数?如果是这样,您可以定义
colProds它们将始终是正数,谢谢。