如何在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