如何计算矩阵列表中特定行的总和

时间:2018-05-28 22:55:30

标签: r lapply purrr

我在矩阵中有数据,矩阵存储在列表中,我想要每个矩阵中特定行的总和。

一些示例数据

A1<-matrix(0:9, nrow=5, ncol=2)
A2<-matrix(10:19, nrow=5, ncol = 2)
A3<-matrix(20:29, nrow=5, ncol = 2)
Mylist<-list(A1, A2, A3)

我可以用

得到每个矩阵中所有行的总和
lapply(Mylist, function(x) apply(x, 1, sum) )

但我只希望特定行的总和,可能是第1行,可能是第4行,具体取决于我想要查看的内容。我知道我可以通过上面的代码生成的结果读取它,但我想要一个更清晰的解决方案,只给我结果。感谢

2 个答案:

答案 0 :(得分:2)

您可以使用purrr:map()

如果您知道输出类型(在这种情况下,似乎是所有整数),您可以更具体,如map_int()。使用map(),您将获得一个列表,并使用特定的map版本,例如map_int(),您将获得一个向量。

library(tidyverse)

ix <- 3 # let's say we want the sum of the third row

map_int(Mylist, ~sum(.x[ix, ]))
[1] 9 29 49

如果你关心的行索引是每个矩阵的变化,你可以改用map2(),这需要两个相同长度的输入:

ixs <- c(1, 2, 3)

map2_int(Mylist, ixs, ~sum(.x[.y, ]))
[1]  5 27 49

或者,如果您需要在基数R中工作,您可以只获取ix所需的索引(此处为sum()),而您内部不需要apply() lapply()

lapply(Mylist, function(x) sum(x[ix, ]))

[[1]]
[1] 9

[[2]]
[1] 29

[[3]]
[1] 49

答案 1 :(得分:0)

one.row.sum <- function(df, row.num) lapply(Mylist, function(df) sum(df[row.num, ]))

one.row.sum(Mylist, 1)
[[1]]
[1] 5

[[2]]
[1] 25

[[3]]
[1] 45