在R中计算相邻行百分比差异并生成新列的异常问题

问题原因分析 #

你的代码错误使用了group_by(WC_30cm_neg)，这会把所有WC值相同的行归为一组。分组后lag(WC_30cm_neg)取的是同一组内的上一行值，而同一组内WC值完全一致，导致差值diff始终为0，最终百分比列全部显示0。

正确解决方案 #

我们需要计算按时间排序后相邻行的WC值变化，不需要对WC值分组，同时要设置百分比精度，避免因数值过小被四舍五入为0。

步骤1：确保时间列格式正确

先确认Date_Time是时间类型，若读入后为字符型，先转换：

df$Date_Time <- as.POSIXct(df$Date_Time)

步骤2：计算相邻行差异与百分比

df_diff <- df %>% 
  arrange(Date_Time) %>% # 按时间排序保证顺序正确
  mutate(
    diff = WC_30cm_neg - lag(WC_30cm_neg), # 计算数值差值
    # 计算百分比差异，设置精度到小数点后两位
    increase_percent = scales::percent(
      (WC_30cm_neg - lag(WC_30cm_neg)) / lag(WC_30cm_neg),
      accuracy = 0.01
    )
  ) %>% 
  filter(!is.na(increase_percent)) # 移除无前置数据的第一行（结果为NA）

替代方案（不依赖scales包）

如果不想使用scales包，可直接格式化百分比字符串：

df_diff <- df %>% 
  arrange(Date_Time) %>% 
  mutate(
    diff = WC_30cm_neg - lag(WC_30cm_neg),
    # 手动计算并格式化百分比，保留两位小数
    increase_percent = sprintf("%.2f%%", (diff / lag(WC_30cm_neg)) * 100)
  ) %>% 
  filter(!is.na(increase_percent))

验证结果 #

运行后会得到正确的百分比差异：

• 第二行（0.3048）相对第一行（0.3051）的差异约为-0.10%
• 第五行（0.3051）相对第四行（0.3048）的差异约为0.10%

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