[英]R datatable grouping based on condition and getting count based on the conditions
我有一個這樣的數據表:
timestamp type status
05-01-2020 12:07:08 A 1
05-01-2020 12:36:05 A 1
05-01-2020 13:34:25 A 1
05-01-2020 23:45:02 A 1
05-01-2020 23:55:02 B 1
05-01-2020 13:44:33 B 2
06-01-2020 01:07:08 A 1
06-01-2020 10:23:05 A 1
06-01-2020 12:11:08 A 2
06-01-2020 22:06:12 B 2
07-01-2020 00:01:05 A 2
07-01-2020 02:17:09 A 1
07-01-2020 12:36:05 B 1
07-01-2020 12:07:08 B 1
07-01-2020 12:36:05 A 1
07-01-2020 12:36:05 A 1
08-01-2020 12:36:05 B 2
08-01-2020 12:36:05 B 1
08-01-2020 12:36:05 B 1
09-01-2020 12:36:05 B 1
09-01-2020 12:07:08 B 2
09-01-2020 12:36:05 B 1
11-01-2020 12:07:08 A 1
11-01-2020 12:36:05 A 1
我正在嘗試使用rleid()
按日期和類型對其進行分組。
dt <- dt[, group_id := rleid(as.IDate(timestamp),type,status = 1)][]
現在我想得到兩個計數。
一是統計每天每個組內滿足條件的實例數。
date type count
05-01-2020 A 4
05-01-2020 B 1
06-01-2020 A 2
07-01-2020 A 3
07-01-2020 B 2
08-01-2020 B 2
09-01-2020 B 2
11-01-2020 A 2
第二個是找到每天滿足條件的組數。
date type count
05-01-2020 A 1
05-01-2020 B 1
06-01-2020 A 1
07-01-2020 A 2
07-01-2020 B 1
08-01-2020 B 1
09-01-2020 B 2
11-01-2020 A 1
1) 統計每個組內每天滿足條件的實例數。
library(data.table)
setDT(df)
df[, .(count = sum(status == 1)), .(timestamp, type)]
# timestamp type count
#1: 05-01-2020 A 4
#2: 05-01-2020 B 1
#3: 06-01-2020 A 2
#4: 06-01-2020 B 0
#5: 07-01-2020 A 3
#6: 07-01-2020 B 2
#7: 08-01-2020 B 2
#8: 09-01-2020 B 2
#9: 11-01-2020 A 2
如果不需要,您可以刪除 0 計數。
2)查找每天滿足條件的組數。
使用type
和status
的rleid
創建一個新列 ( count_N
),對於status = 1
,計算每個timestamp
和type
的唯一值。
df[, count_N := rleid(type, status), timestamp]
df[status == 1, .(count = uniqueN(count_N)), .(timestamp, type)]
# timestamp type count
#1: 05-01-2020 A 1
#2: 05-01-2020 B 1
#3: 06-01-2020 A 1
#4: 07-01-2020 A 2
#5: 07-01-2020 B 1
#6: 08-01-2020 B 1
#7: 09-01-2020 B 2
#8: 11-01-2020 A 1
我們可以先使用 as.POSIXct 將“時間戳”轉換為日期時間as.POSIXct
,然后將其轉換為Date
class
library(data.table)
setDT(dt)[, timestamp := as.POSIXct(timestamp,
format = '%m-%d-%Y %H:%M:%S')][, date := as.IDate(timestamp)]
dt[status == 1, .N, .(date, type)]
#. date type N
#1: 2020-05-01 A 4
#2: 2020-05-01 B 1
#3: 2020-06-01 A 2
#4: 2020-07-01 A 3
#5: 2020-07-01 B 2
#6: 2020-08-01 B 2
#7: 2020-09-01 B 2
#8: 2020-11-01 A 2
對於第二種情況
dt[, grp := rleid(type, status, date)]
dt[status == 1, .(count = uniqueN(grp)), .(date, type)]
# date type count
#1: 2020-05-01 A 1
#2: 2020-05-01 B 1
#3: 2020-06-01 A 1
#4: 2020-07-01 A 2
#5: 2020-07-01 B 1
#6: 2020-08-01 B 1
#7: 2020-09-01 B 2
#8: 2020-11-01 A 1
dt <- structure(list(timestamp = structure(c(1588349228, 1588350965,
1588354465, 1588391102, 1588391702, 1588355073, 1590988028, 1591021385,
1591027868, 1591063572, 1593576065, 1593584229, 1593621365, 1593619628,
1593621365, 1593621365, 1596299765, 1596299765, 1596299765, 1598978165,
1598976428, 1598978165, 1604250428, 1604252165), class = c("POSIXct",
"POSIXt"), tzone = ""), type = c("A", "A", "A", "A", "B", "B",
"A", "A", "A", "B", "A", "A", "B", "B", "A", "A", "B", "B", "B",
"B", "B", "B", "A", "A"), status = c(1L, 1L, 1L, 1L, 1L, 2L,
1L, 1L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 2L, 1L,
1L, 1L)), class = "data.frame", row.names = c(NA, -24L),
index = structure(integer(0), "`__status`" = c(1L,
2L, 3L, 4L, 5L, 7L, 8L, 12L, 13L, 14L, 15L, 16L, 18L, 19L, 20L,
22L, 23L, 24L, 6L, 9L, 10L, 11L, 17L, 21L)))
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.