组意味着排除个人并拥有 NA

Question

在之前的链接中： Calculating a group mean while excluding each cases individual value

matt_k 对计算组意味着排除个人做了一个聪明的回答。他提出以下建议：

set.seed(123)
df <- data.frame(group = rep(letters[1:3], each = 3), 
             value = rnorm(9), stringsAsFactors = F)
df$loo_mean <- unlist(tapply(df$value, df$group, 
                 function(x) (sum(x) - x) / (length(x) - 1)))
df

但是代码没有正确处理 NA，因为如果该组有 NA，它会为该组的所有个人产生 NA。 谁能解决这个问题？

Answer 1

在调用函数之前省略 NA。

na.omit(df)

Answer 2

解决方案很简单。将na.rm = TRUE 添加到sum。

比较

> sum(c(1:3, NA), na.rm = FALSE)
[1] NA

和

> sum(c(1:3, NA), na.rm = TRUE)
[1] 6

有关详细信息，请参阅?sum。