R语言如何按145天间隔为各巢观测值分配唯一周期编号

R 实现方案 #

方案1：tidyverse 实现（逻辑直观，适配分组操作）

核心逻辑是按Nest_ID分组后，滚动跟踪每个周期的起始累计间隔值，每当当前累计间隔与周期起始值的差≥145时，就开启新周期，最终将周期起始值映射为连续编号。

# 加载依赖包
library(dplyr)
library(purrr)

# 计算Nest_period
df_result <- df %>%
  group_by(Nest_ID) %>%
  mutate(
    # 滚动记录每个周期的起始累计间隔
    period_start = accumulate(total_interval, function(prev_start, curr_total) {
      if (curr_total - prev_start >= 145) curr_total else prev_start
    }),
    # 将周期起始值转换为连续的周期编号
    Nest_period = as.integer(factor(period_start, levels = unique(period_start)))
  ) %>%
  ungroup() %>%
  select(-period_start) # 移除临时辅助列

运行后输出的df_result和你期望的结果完全一致。

方案2：base R 实现（无需额外安装包）

如果不想依赖第三方包，可以用分组循环的方式实现：

df$Nest_period <- ave(df$total_interval, df$Nest_ID, FUN = function(x) {
  if (length(x) == 0) return(numeric(0))
  period_start <- x[1]
  period_no <- 1
  res <- rep(1, length(x))
  for (i in seq_along(x)[-1]) {
    if (x[i] - period_start >= 145) {
      period_start <- x[i]
      period_no <- period_no + 1
    }
    res[i] <- period_no
  }
  res
})

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