R合计开/关时间

czfnxgou 于 2023-04-03 发布在其他

关注(0)|答案(1)|浏览(73)

我有一个数据集，其中包含时间戳和流量数据（以加仑/分钟为单位），每2分钟记录一次数据。

df <- structure(list(dt = structure(c(1519891200, 1519891320, 1519891440, 1519891560, 
                1519891680, 1519891800, 1519891920, 1519892040, 1519892160, 1519892280,
                1519892400, 1519892520, 1519892640, 1519892760, 1519892880), 
       class = c("POSIXct", "POSIXt"), tzone = ""), 
       gpm = c(0, 0, 0, 50, 50, 50, 50, 50, 0, 0, 80, 80, 80, 0, 0)), 
       .Names = c("dt", "gpm"), 
       row.names = c(NA, 15L), 
       class = "data.frame")

#                     dt gpm
# 1  2018-03-01 03:00:00   0
# 2  2018-03-01 03:02:00   0
# 3  2018-03-01 03:04:00   0
# 4  2018-03-01 03:06:00  50
# 5  2018-03-01 03:08:00  50
# 6  2018-03-01 03:10:00  50
# 7  2018-03-01 03:12:00  50
# 8  2018-03-01 03:14:00  50
# 9  2018-03-01 03:16:00   0
# 10 2018-03-01 03:18:00   0
# 11 2018-03-01 03:20:00  80
# 12 2018-03-01 03:22:00  80
# 13 2018-03-01 03:24:00  80
# 14 2018-03-01 03:26:00   0
# 15 2018-03-01 03:28:00   0

我想将数据聚合到一个dataframe中，其中包含on，off时间和该时间段的平均gpm。每次连续出现2个以上的零时，它都应该算作一个新事件。
最终数据应该如下所示：

#                    on                 off avg
# 1 2018-03-01 03:06:00 2018-03-01 03:14:00  50
# 2 2018-03-01 03:20:00 2018-03-01 03:24:00  80

来源：https://stackoverflow.com/questions/49949186/r-aggregate-on-off-time

1条答案

按热度按时间

ctehm74n1#

在dplyr中使用cumsum：

library(dplyr)
df %>% mutate(id = cumsum((gpm + c(0, lag(gpm)[-1])) == 0)) %>% 
  filter(gpm != 0) %>% group_by(id) %>% 
  mutate_at(vars(dt), .funs = list(on = min, off = max)) %>% 
  summarise_all(list(~mean(., na.rm=TRUE))) %>%
  select(-dt,-id,-gpm,gpm)

或者在data.table中使用rle和cumsum：

library(data.table)
setDT(df)[, id := {
  rid <-rle(+(df$gpm + c(0,lag(df$gpm)[-1]) != 0 ))
  rid$values <- cumsum(rid$values) * rid$values
  inverse.rle(rid) * +(df$gpm!=0)
        }][gpm !=0
          ][ , list(on = min(dt), off = max(dt), gpm = mean(gpm)) , 
               by = .(id)][,id:=NULL][]

输出将为：

#                    on                 off gpm
# 1 2018-03-01 03:06:00 2018-03-01 03:14:00  50
# 2 2018-03-01 03:20:00 2018-03-01 03:24:00  80

赞(0）回复(0）举报 2023-04-03

我来回答

R合计开/关时间

1条答案

相关问题

热门标签

最新问答