Question

在计算聚合方式时，我无法处理NAs。请参阅以下代码：

messages.push_front(*new Message(messageID));

如果向量具有所有NA，我希望NA而不是NaN，即我希望输出为

tab=data.frame(a=c(1:3,1:3), b=c(1,2,NA,3,NA,NA))
tab
  a  b
1 1  1
2 2  2
3 3 NA
4 1  3
5 2 NA
6 3 NA

attach(tab)
aggregate(b, by=list(a), data=tab, FUN=mean, na.rm=TRUE)
  Group.1   x
1       1   2
2       2   2
3       3 NaN

我尝试使用自定义功能：

  Group.1   x
1       1   2
2       2   2
3       3  NA

但是，我收到以下错误：

adjmean=function(x) {if(all(is.na(x))) NA else mean(x,na.rm=TRUE)}

简而言之，如果列中包含所有NA，我希望NA作为输出而不是NaN。如果它具有很少的NA，那么它应该计算忽略NA的均值。

任何帮助都将不胜感激。

由于

Answer 1

这与你所拥有的非常接近，但用自定义函数替换mean(x, na.rm=TRUE)，该函数可以计算非NA值的平均值，或者提供NA本身：

R> with(tab, 
        aggregate(b, by=list(a), FUN=function(x) 
             if (any(is.finite(z<-na.omit(x)))) mean(z) else NA))
  Group.1  x
1       1  2
2       2  2
3       3 NA
R>

这真的是一条线，但是我将它拆开以使其适合SO显示。

你已经有了类似的想法，但我更改了函数以在所有情况下都返回合适的值。

Answer 2

您的功能没有任何问题。错误的原因是您在aggregate的默认方法中使用了一个不存在的参数：

adjmean = function(x) {if(all(is.na(x))) NA else mean(x,na.rm=TRUE)}
attach(tab)  ## Just because you did it. I don't recommend this.

## Your error
aggregate(b, by=list(a), data=tab, FUN=adjmean)
# Error in FUN(X[[i]], ...) : 
#   unused argument (data = list(a = c(1, 2, 3, 1, 2, 3), b = c(1, 2, NA, 3, NA, NA)))

## Dropping the "data" argument
aggregate(b, list(a), FUN = adjmean)
#   Group.1  x
# 1       1  2
# 2       2  2
# 3       3 NA

如果您想使用data参数，则应对formula使用aggregate方法。但是，此方法会以不同的方式处理NA，因此您需要一个额外的参数na.action。

示例：

detach(tab) ## I don't like having things attached
aggregate(b ~ a, data = tab, adjmean)
#   a b
# 1 1 2
# 2 2 2
aggregate(b ~ a, data = tab, adjmean, na.action = na.pass)
#   a  b
# 1 1  2
# 2 2  2
# 3 3 NA

R中的聚合NA

2 个答案: