使用rle和for循环的最长连续时间段高于阈值

1 回答

• 0

  简单地说，使用诸如threshold_find之类的已定义函数来概括您的流程，并将每年子集化的数据帧传递到其中，这可以使用 by 进行处理 .

  作为面向对象的包装器， by 将数据帧切片一个或多个因子（即年份），并返回定义的函数输出的任何对象的列表，这里是最大数据帧 . 最后， do.call() 行将 by 列表中的所有数据帧绑定到一个数据帧中 .

  threshold_find <- function(df) {
    rl <- with(df, rle(threshold))      
    cs <- cumsum(rl$lengths)
  
    index <- cbind(cs[rl$values == 1] - rl$length[rl$values == 1] + 1,
                   cs[rl$values == 1])
    df <- data.frame(index)
    df$X3 <- df$X2 - df$X1+1
  
    max <- df[which.max(df$X3),]
    max      
  }
  
  finaldf <- do.call(rbind, by(df, df$year, FUN=threshold_find))
  
  finaldf
  #      X1 X2 X3
  # 2009  1 10 10
  # 2010  9 31 23
  # 2011 10 18  9
  # 2012 12 30 19
  

  回复于 2024-04-16T19:13:45+08:00