在R语言中用前2或3次测量均值替换重复测量值的方法

针对你的需求，我们可以通过分组计算每组的重复次数，再根据次数取对应前N次测量的均值替换整组的BP/BP2值。以下是两种实现方案：

方法1：使用dplyr（tidyverse 生态） #

先加载dplyr包，按ID和Measurement分组后处理：

library(dplyr)

# 你的原始数据
ID <- c(3,3,3,3,3,3,5,5,5,5)
Measurement <- c(2, 2, 2, 2, 2, 19, 6, 6, 7, 8)
BP <- c(14, 14, 15, 16, 14, 13, 14, 17, 18, 20)
BP2 <- c(7, 7, 8, 9, 10, 11, 14, 7, 8, 9)
DF1 <- data.frame(ID, Measurement, BP, BP2)

# 处理数据
DF_processed <- DF1 %>%
  group_by(ID, Measurement) %>%
  mutate(
    # 计算每组的重复次数
    group_size = n(),
    # 确定取前几次测量值计算均值
    take_n = case_when(
      group_size >= 3 ~ 3,
      group_size == 2 ~ 2,
      TRUE ~ 1
    ),
    # 计算前take_n次的BP均值
    BP_mean = mean(slice_head(cur_data(), n = take_n)$BP),
    # 计算前take_n次的BP2均值
    BP2_mean = mean(slice_head(cur_data(), n = take_n)$BP2)
  ) %>%
  # 替换整组的BP和BP2值
  mutate(BP = BP_mean, BP2 = BP2_mean) %>%
  # 移除辅助计算列
  select(-group_size, -take_n, -BP_mean, -BP2_mean) %>%
  ungroup()

# 查看结果
print(DF_processed)

处理结果： #

# A tibble: 10 × 4
      ID Measurement    BP   BP2
   <dbl>       <dbl> <dbl> <dbl>
 1     3           2  14.3  8   
 2     3           2  14.3  8   
 3     3           2  14.3  8   
 4     3           2  14.3  8   
 5     3           2  14.3  8   
 6     3          19  13    11   
 7     5           6  15.5 10.5 
 8     5           6  15.5 10.5 
 9     5           7  18     8   
10     5           8  20     9   

方法2：使用data.table #

如果处理大数据集，data.table的效率更高：

library(data.table)

# 转换为data.table格式
setDT(DF1)

# 处理数据
DF_processed <- DF1[, 
                    {
                      group_size = .N
                      take_n = if (group_size >= 3) 3 else if (group_size == 2) 2 else 1
                      BP_mean = mean(head(BP, take_n))
                      BP2_mean = mean(head(BP2, take_n))
                      # 返回整组替换后的数据
                      .(Measurement = Measurement, BP = BP_mean, BP2 = BP2_mean)
                    }, 
                    by = .(ID)]

# 查看结果
print(DF_processed)

核心逻辑说明： #

1. 按ID + Measurement分组，识别重复测量的组
2. 对每组判断重复次数：≥3则取前3次均值，=2则取前2次均值，=1则保留原值
3. 用计算得到的均值替换整组的BP/BP2值

内容的提问来源于stack exchange，提问作者19056530