将函数应用于矩阵列表的每个元素

Question

我有一个矩阵列表。

（下面是一个简化的示例，我实际上有3个矩阵的列表，第一个矩阵是2D的，而第二个和第三个矩阵是3D的）

> a <- matrix(-1:2, ncol = 2)
> b <- array(c(-2:5), dim=c(2, 2, 2))
> c_list <- list(a,b)
> c_list
[[1]]
     [,1] [,2]
[1,]   -1    1
[2,]    0    2

[[2]]
, , 1

     [,1] [,2]
[1,]   -2    0
[2,]   -1    1

, , 2

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

我想将函数max（0，c_list）应用于每个元素（无循环），以便与“ c_list”具有相同类型的对象，但负值替换为零。

> output
[[1]]
     [,1] [,2]
[1,]    0    1
[2,]    0    2

[[2]]
, , 1

     [,1] [,2]
[1,]    0    0
[2,]    0    1

, , 2

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

我已经成功地针对矩阵或具有mapply或lapply的列表进行了此操作，但没有针对矩阵列表进行此操作。

答案：索托斯的答案

output <- lapply(c_list, function(i)replace(i, i < 0, 0))

或Moody_Mudskipper的答案

output <- lapply(c_list,pmax,0)

Answer 1

使用apply和lapply ：

a <- matrix(-1:2, ncol = 2)
b <- matrix(-3:0, ncol = 2)
c <- list(a,b)

d <- lapply(c, function(m) {
  apply(m, c(1, 2), function(x) max(0, x))
})

输出：

> d
[[1]]
     [,1] [,2]
[1,]    0    1
[2,]    0    2

[[2]]
     [,1] [,2]
[1,]    0    0
[2,]    0    0

Answer 2

您可以使用pmax ，它将保留源矩阵的格式并进行矢量化，比使用max循环更快。

lapply(c_list,pmax,0)