如何在R中按分组取众数替换NA值，保留无有效值的NA

API 探路者

2026-6-27

问题描述

原始数据框：

x	y	Owner	H
1	2	A	1
2	5	A	NA
2	2	A	5
2	5	A	4
2	5	B	4
2	5	C	3
3	3	C	6
3	3	C	NA
3	3	D	NA

需按(x,y,Owner)分组，将H列的NA值替换为对应组内H的众数；若分组内仅含NA值，保留NA。目标结果如下：

x	y	Owner	H
1	2	A	1
2	5	A	4
2	2	A	5
2	5	A	4
2	5	B	4
2	5	C	3
3	3	C	6
3	3	C	6
3	3	D	NA

尝试的代码存在问题，会替换所有H值而非仅NA值：

new_data <- my_data %>% group_by('Owner', 'x', 'y') %>% mutate(H=mean(H, na.rm=TRUE))

解决方案

步骤1：定义众数计算函数

R没有内置众数函数，自定义函数处理全NA分组的情况：

get_mode <- function(x) {
  x_non_na <- x[!is.na(x)]
  if (length(x_non_na) == 0) {
    return(NA)
  }
  freq <- table(x_non_na)
  mode_val <- names(freq)[freq == max(freq)]
  as.numeric(mode_val[1]) # 多众数时取第一个，可按需调整
}

步骤2：用dplyr处理数据

library(dplyr)

new_data <- my_data %>%
  group_by(Owner, x, y) %>% # 正确分组，列名无需加引号
  mutate(H = ifelse(is.na(H), get_mode(H), H)) %>% # 仅替换NA值
  ungroup()