NA значения и R агрегатная функция

Question

Вот простой кадр данных с отсутствующим значением:

M = data.frame(Name = c('name','name'), Col1 = c(NA,1) , Col2 = c(1,1))

Когда я применяю агрегат М следующим образом:

aggregate(.~Name, M, FUN=sum, na.rm=TRUE)

результат является:

RowName Col1 Col2 
name 1 1 

Таким образом, Первая строка ire игнорируется. Но если я

aggregate(M[,2:3], by=list(M$Name), FUN=sum, na.rm=TRUE)

результат

Group.1 Col1 Col2 
name 1 2 

Так только (1,1) запись игнорируется.

Это вызвало серьезную головную боль отладки в одном из моих кодов, поскольку я думал, что эти два вызова эквивалентны. Есть ли веская причина, почему метод ввода «формула» обрабатывается по-разному?

Спасибо.

Answer 1

Хороший вопрос, но, на мой взгляд, это не должно приводить к отладке головной боли основных, поскольку она достаточно четко документирована в нескольких местах на странице руководства для aggregate.

Во-первых, в разделе использования:

## S3 method for class 'formula' 
aggregate(formula, data, FUN, ..., 
      subset, na.action = na.omit) 

Позже, в описании:

na.action: функция, которая указывает на то, что должно произойти, когда данные содержат НС значения. По умолчанию следует игнорировать отсутствующие значения в данных переменных.

---

Я не могу ответить, почему режима формулы был написан по-другим --- это то, что авторы функции должны отвечать --- но используя приведенную выше информацию, вы, вероятно, можете использовать следующий :

aggregate(.~Name, M, FUN=sum, na.rm=TRUE, na.action=NULL) 
# Name Col1 Col2 
# 1 name 1 2 

Answer 2

Если вы хотите формулу версия эквивалентной попробовать это:

M = data.frame(Name = rep('name',5), Col1 = c(NA,rep(1,4)) , Col2 = rep(1,5)) 
aggregate(. ~ Name, M, function(x) sum(x, na.rm=TRUE), na.action = na.pass)