Question

我有两个矩阵列表，每个列表有12个矩阵。我想从每个列表中减去两个矩阵。这是可重复的例子：

x1 <- matrix((1:25)*0.1, ncol=5)
x2<-x1+0.3
list1<-list(x1,x2)
list2<-list(x1*2,x2*0.1)
sub<-mapply("-",list1,list2)

这将输出返回为

    > sub
      [,1] [,2]
 [1,] -0.1 0.36
 [2,] -0.2 0.45
 [3,] -0.3 0.54
 [4,] -0.4 0.63
 [5,] -0.5 0.72
 [6,] -0.6 0.81
 [7,] -0.7 0.90
 [8,] -0.8 0.99
 [9,] -0.9 1.08
[10,] -1.0 1.17
[11,] -1.1 1.26
[12,] -1.2 1.35
[13,] -1.3 1.44
[14,] -1.4 1.53
[15,] -1.5 1.62
[16,] -1.6 1.71
[17,] -1.7 1.80
[18,] -1.8 1.89
[19,] -1.9 1.98
[20,] -2.0 2.07
[21,] -2.1 2.16
[22,] -2.2 2.25
[23,] -2.3 2.34
[24,] -2.4 2.43
[25,] -2.5 2.52

我想将输出作为5x5矩阵的列表：例如

>list1[[1]]-list2[[1]]
     [,1] [,2] [,3] [,4] [,5]
[1,] -0.1 -0.6 -1.1 -1.6 -2.1
[2,] -0.2 -0.7 -1.2 -1.7 -2.2
[3,] -0.3 -0.8 -1.3 -1.8 -2.3
[4,] -0.4 -0.9 -1.4 -1.9 -2.4
[5,] -0.5 -1.0 -1.5 -2.0 -2.5

向simplify添加mapply参数会引发以下错误：

Error in .Primitive("-")(dots[[1L]][[1L]], dots[[2L]][[1L]], simplify = dots[[3L]][[1L]]) : 
  operator needs one or two arguments

如何获取清单？

Answer 1

作为@Sotos评论，使用Map作为mapply的包装器将结果输出为列表而不是绑定矩阵。实际上，您仍然可以使用mapply但使用SIMPLIFY = FALSE参数。注意：SIMPLIFY在所有CAPS中，与其他同名的同名家族参数不同。

使用默认的 TRUE ，mapply将每个差异转换为矩阵。以下是同等的电话：

sub1 <- Map(`-`, list1, list2)

sub2 <- mapply(`-`, list1, list2, SIMPLIFY = FALSE)

all.equal(sub1, sub2)
# [1] TRUE

mapply从两个不同的列表中减去矩阵并返回矩阵列表

1 个答案: