data.table 1.8.x mean（）函数自动删除NA？

Question

今天我发现我的程序中有一个错误，因为data.table自动删除NA表示mean

例如：

> a<-data.table(a=c(NA,NA,FALSE,FALSE), b=c(1,1,2,2))
> a

> a[,list(mean(a), sum(a)),by=b]
   b V1 V2
1: 1  0 NA // Why V1 = 0 here? I had expected NA
2: 2  0  0


> mean(c(NA,NA,FALSE,FALSE))
[1] NA
> mean(c(NA,NA))
[1] NA
> mean(c(FALSE,FALSE))
[1] 0

这是预期的行为吗？

Answer 1

这不是预期的。 看起来像是一个优化问题......

> a[,list(mean(a), sum(a)),by=b]
   b V1 V2
1: 1  0 NA
2: 2  0  0
> options(datatable.optimize=FALSE)
> a[,list(mean(a), sum(a)),by=b]
   b V1 V2
1: 1 NA NA
2: 2  0  0
> 

调查并修复于v1.8.9，很快将在CRAN上。 来自新闻 ：

j中的mean（）自v1.8.2起已经过优化，但不尊重na.rm = TRUE（默认值）。 非常感谢Colin Fang的报道。 测试补充。

v1.8.2中的新功能是：

mean（）现在自动优化，＃1231。 当存在大量组时，这可以加快分组20次。 请参阅wiki第3点 ，不再需要知道。 通过设置选项关闭优化（datatable.optimize = 0）。