R语言使用zero-crossing过零算法实现正转负过零点检测的咨询

首先修正你提供的示例数据中的语法笔误：原x向量中-2 -4之间缺少逗号，修正后的数据集构造代码如下：

time <- c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20)
x <- c(1,2,3,-2,-4,-8,-2,0,1,2,3,4,2,1,-3,-4,-7,-4,-1,1,2)
df <- data.frame(time, x)

核心判定逻辑 #

我们要识别的正转负过零事件满足：上一时刻x值为正，当前时刻x值为非正，满足条件标记为1，其余情况标记为0，第一行无前置数据默认标记为0。如果业务规则要求严格正转负（排除0值），可以把当前值的判定条件从<=0改为<0。

实现方式 #

方式1：基础R实现（无额外依赖） #

detect_pos2neg_cross <- function(df, col_name = "x") {
  cross_flag <- integer(nrow(df))
  for (i in 2:nrow(df)) {
    prev_val <- df[[col_name]][i-1]
    curr_val <- df[[col_name]][i]
    if (prev_val > 0 && curr_val <= 0) {
      cross_flag[i] <- 1
    }
  }
  df$cross_pos2neg <- cross_flag
  return(df)
}

# 调用函数获得结果
df_result <- detect_pos2neg_cross(df)

方式2：dplyr实现（代码更简洁） #

library(dplyr)
df <- df %>%
  mutate(cross_pos2neg = case_when(
    lag(x, default = 0) > 0 & x <= 0 ~ 1,
    TRUE ~ 0
  ))

结果验证 #

运行上述代码后，你会得到的标记结果中：

• 第4行（前值为3>0，当前值为-2<=0）标记为1
• 第15行（前值为1>0，当前值为-3<=0）标记为1
• 其余行均标记为0，符合正转负过零的判定规则。

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