apply() vs. sweep() в R

Question

Я писал заметки для сравнения apply() и sweep() и обнаружил следующие странные отличия. Чтобы сгенерировать тот же результат, sweep() требует MARGIN = 1, тогда как apply требует MARGIN = 2. Также аргумент, определяющий матрицу, имеет верхний регистр X в apply(), но строчный в sweep().

my.matrix <- matrix(seq(1,9,1), nrow=3) 
row.sums <- rowSums(my.matrix) 
apply.matrix <- apply(X = my.matrix, MARGIN = 2, FUN = function (x) x/row.sums) 
sweep.matrix <- sweep(x = my.matrix, MARGIN = 1, STATS = rowSums(my.matrix), FUN="/") 
apply.matrix - sweep.matrix ##yup same matrix 

Не работает функция «наложения типа»? Является ли это еще одним R quirk или я потерял рассудок?

Answer 1

Обратите внимание, что для apply,

Если каждый вызов «ПОТЕХИ» возвращает вектор длины «N», то «применить» возвращает массив размерности «с (п, тусклый (X) [ MARGIN])»если 'п> 1'

в вашем примере, MARGIN может (и должен) быть установлен в 1 в обоих случаях; но возвращаемое значение от apply должно быть перенесено. Легче всего это увидеть, если исходная матрица не квадратная:

my.matrix <- matrix(seq(1,12,1), nrow=4) 
apply.matrix <- t(apply(X = my.matrix, MARGIN = 1, FUN = function(x) x/sum(x))) 
sweep.matrix <- sweep(x = my.matrix, MARGIN = 1, STATS = rowSums(my.matrix), FUN="/") 
all.equal(apply.matrix, sweep.matrix) 
# [1] TRUE 

Также см этот ответ Can you implement 'sweep' using apply in R?, который говорит, что очень то же самое.