mapply для вычитания матриц из двух разных списков и возврата списка матриц

У меня есть два списка матриц, каждый список имеет 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

Я хочу вывод в виде списка 5х5 матриц: например

>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

Как я могу получить список?

1 ответ

Решение

Как комментарии @Sotos, используйте Map который является оберткой для mapply выводить результаты в виде списка, а не в виде матрицы. На самом деле, вы все еще можете использовать mapply но используйте SIMPLIFY = FALSE аргумент. НОТА: SIMPLIFY находится во ВСЕХ КОПИЯХ в отличие от другого применяемого семейного аргумента с тем же именем.

Используя значение по умолчанию TRUE, mapply преобразует каждое различие в матрицу. Ниже приведены эквивалентные звонки:

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

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

all.equal(sub1, sub2)
# [1] TRUE
Другие вопросы по тегам