如何在R中有效地总结列表中的矩阵（无循环）？

Question

假设我有一个函数返回2个矩阵。对于每个输入i。

f<-function(i){ 
    m1=matrix(i,ncol=5,nrow=5) 
    m2=matrix(i*10,ncol=5,nrow=5) 
    r=list(m1=m1,m2=m2) 
    return(r) 
} 
d=lapply(1:3,f) 

我怎样才能得到2件事情没有for循环？

d[[1]]$m1+d[[2]]$m1+d[[3]]$m1 
d[[1]]$m2+d[[2]]$m2+d[[3]]$m2 

我试图sum(lapply(1:3,function(x) (d[[x]]$m1)))为M1，但不起作用。

Answer 1

您可以使用Reduce：

library(purrr) 
map(transpose(d), reduce, `+`) 

这将返回：

Reduce("+", lapply(d, function(x) x$m1)) 

也可以使用transpose，map和reduce从purrr包同时获得m1和m2

$m1 
    [,1] [,2] [,3] [,4] [,5] 
[1,] 6 6 6 6 6 
[2,] 6 6 6 6 6 
[3,] 6 6 6 6 6 
[4,] 6 6 6 6 6 
[5,] 6 6 6 6 6 

$m2 
    [,1] [,2] [,3] [,4] [,5] 
[1,] 60 60 60 60 60 
[2,] 60 60 60 60 60 
[3,] 60 60 60 60 60 
[4,] 60 60 60 60 60 
[5,] 60 60 60 60 60 

Answer 2

通过内部嵌套list元素的名称（lapply(paste0("m", 1:2), ..）我们循环，通过list元素循环（lapply(d, ...），提取与names（d1[[nm]]）相匹配的元素，并使用Reduce与+（我没看到对方的回答同时发布）。

lapply(paste0("m", 1:2), function(nm) lapply(d, function(d1) Reduce(`+`, d1[[nm]]))) 
#[[1]] 
#  [,1] [,2] [,3] [,4] [,5] 
#[1,] 6 6 6 6 6 
#[2,] 6 6 6 6 6 
#[3,] 6 6 6 6 6 
#[4,] 6 6 6 6 6 
#[5,] 6 6 6 6 6 

#[[2]] 
#  [,1] [,2] [,3] [,4] [,5] 
#[1,] 60 60 60 60 60 
#[2,] 60 60 60 60 60 
#[3,] 60 60 60 60 60 
#[4,] 60 60 60 60 60 
#[5,] 60 60 60 60 60