我有以下数据。
DOY <- c(1:30) #Day of Year
rain <- c(0,0,0,0,0,0,0,5.5,15.5,20.1,10.4,0,0,0,0,0,0,0,0,7.6,19.2,44.5,3.4,5.5,0,0,0,0,0,0)
DF <- data.frame(DOY=DOY,rain=rain)
我只对下雨天,降雨事件感兴趣。
DF <- DF[DF$rain > 0,]
df
的输出,它们是降雨事件:
DOY rain
8 5.5
9 15.5
10 20.1
11 10.4
20 7.6
21 19.2
22 44.5
23 3.4
24 5.5
在我的示例中,我有两个降雨事件。如果连续两天都没有下雨,那总是一个新事件。
现在,我想为每个降雨事件分配一个ID
。这样的东西是我的预期输出。
DOY rain ID
8 5.5 1
9 15.5 1
10 20.1 1
11 10.4 1
20 7.6 2
21 19.2 2
22 44.5 2
23 3.4 2
24 5.5 2
答案 0 :(得分:2)
您可以使用:
cumsum(+c(TRUE, diff(DF$DOY) > 1))