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

Мне нужен вывод в виде списка матриц 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

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


r matrix mapply
person Pulsar_534511    schedule 15.02.2018    source источник
comment
Используйте Map вместо mapply., то есть Map(`-`, list1, list2)   -  person Sotos    schedule 15.02.2018
comment
Спасибо! Я не знал о функции Map.   -  person Pulsar_534511    schedule 15.02.2018

Ответы (1)


arrow_upward
1
arrow_downward

Как комментирует @Sotos, используйте Map, который является оболочкой для mapply, чтобы выводить результаты в виде списка, а не связанной матрицы. На самом деле вы все еще можете использовать mapply, но использовать аргумент SIMPLIFY = FALSE. ПРИМЕЧАНИЕ. SIMPLIFY пишется ЗАГЛАВНЫМИ БУКВАМИ, в отличие от другого аргумента семейства apply с тем же именем.

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

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

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

all.equal(sub1, sub2)
# [1] TRUE
person Parfait    schedule 15.02.2018
comment
Спасибо за ответ! Я не использовал ЗАГЛАВНЫЕ БУКВЫ в аргументе SIMPLIFY. - person Pulsar_534511; 15.02.2018