fork download
copy 
  1. library(magrittr)
  2. library(data.table)
  3. Sys.setlocale(locale = "cht") # windows OS
  4. d <-
  5.   fread(
  6.     "id  date  time   sth
  7.      個體A  2012/10/11  20:00   實驗資料OOXX
  8.      個體A  2012/10/11  23:00   實驗資料OOXX
  9.      個體A  2012/10/12  03:00   實驗資料OOXX
  10.      個體B  2012/12/11  05:00   實驗資料OOXX
  11.      個體B  2012/12/11  11:05   實驗資料OOXX
  12.      個體B  2012/12/11  13:00   實驗資料OOXX
  13.      個體B  2012/12/11  18:00   實驗資料OOXX
  14.      個體B  2012/12/11  20:00   實驗資料OOXX"
  15.   )
  16. # 定義日期時間欄位
  17. d[, datetime := paste(date, time) %>% as.POSIXct()] 
  18.  
  19. # 先按 id 再按 datetime 排序
  20. setorder(d, id, datetime) 
  21.  
  22. # 在各 id 之內，計算所有 datetime 的兩兩距離（秒），以 complete 法分群，
  23. # 以大於 6 小時（即 10800 秒）當不同群切點，
  24. # 在 group 欄記錄結果。
  25. # h 故意加上一個小正數是為了使正好 7 小時的情況仍在同一組
  26. d[, group :=
  27.     datetime %>%
  28.     dist %>%
  29.     hclust(method = "complete") %>% 
  30.     cutree(h = 7 * 60 * 60 + 0.0001), by = id]
  31. d
  32. #       id       date  time          sth            datetime group
  33. # 1: 個體A 2012/10/11 20:00 實驗資料OOXX 2012-10-11 20:00:00     1
  34. # 2: 個體A 2012/10/11 23:00 實驗資料OOXX 2012-10-11 23:00:00     1
  35. # 3: 個體A 2012/10/12 03:00 實驗資料OOXX 2012-10-12 03:00:00     2
  36. # 4: 個體B 2012/12/11 05:00 實驗資料OOXX 2012-12-11 05:00:00     1
  37. # 5: 個體B 2012/12/11 11:05 實驗資料OOXX 2012-12-11 11:05:00     2
  38. # 6: 個體B 2012/12/11 13:00 實驗資料OOXX 2012-12-11 13:00:00     2
  39. # 7: 個體B 2012/12/11 18:00 實驗資料OOXX 2012-12-11 18:00:00     3
  40. # 8: 個體B 2012/12/11 20:00 實驗資料OOXX 2012-12-11 20:00:00     3
  41.  
  42. # 在各 id 及各 group 內取最後一列
  43. d[, .SD[.N,], by = .(id, group)]
  44. #       id group       date  time          sth            datetime
  45. # 1: 個體A     1 2012/10/12 03:00 實驗資料OOXX 2012-10-12 03:00:00
  46. # 2: 個體B     1 2012/12/11 05:00 實驗資料OOXX 2012-12-11 05:00:00
  47. # 3: 個體B     2 2012/12/11 13:00 實驗資料OOXX 2012-12-11 13:00:00
  48. # 4: 個體B     3 2012/12/11 20:00 實驗資料OOXX 2012-12-11 20:00:00
  49.  
Success #stdin #stdout #stderr 0.34s 43380KB
comments (?)
stdin
copy 
Standard input is empty
stdout
copy 
[1] ""
      id       date  time          sth            datetime group
1: 個體A 2012/10/11 20:00 實驗資料OOXX 2012-10-11 20:00:00     1
2: 個體A 2012/10/11 23:00 實驗資料OOXX 2012-10-11 23:00:00     1
3: 個體A 2012/10/12 03:00 實驗資料OOXX 2012-10-12 03:00:00     1
4: 個體B 2012/12/11 05:00 實驗資料OOXX 2012-12-11 05:00:00     1
5: 個體B 2012/12/11 11:05 實驗資料OOXX 2012-12-11 11:05:00     2
6: 個體B 2012/12/11 13:00 實驗資料OOXX 2012-12-11 13:00:00     2
7: 個體B 2012/12/11 18:00 實驗資料OOXX 2012-12-11 18:00:00     3
8: 個體B 2012/12/11 20:00 實驗資料OOXX 2012-12-11 20:00:00     3
      id group       date  time          sth            datetime
1: 個體A     1 2012/10/12 03:00 實驗資料OOXX 2012-10-12 03:00:00
2: 個體B     1 2012/12/11 05:00 實驗資料OOXX 2012-12-11 05:00:00
3: 個體B     2 2012/12/11 13:00 實驗資料OOXX 2012-12-11 13:00:00
4: 個體B     3 2012/12/11 20:00 實驗資料OOXX 2012-12-11 20:00:00
stderr
copy 
Warning message:
In Sys.setlocale(locale = "cht") :
  OS reports request to set locale to "cht" cannot be honored
https://ideone.com/CD6sc6
language:
R (R 3.5.2)
created:
4 years ago
visibility:
public

Share or Embed source code

Discover > Sphere Engine API

The brand new service which powers Ideone!

Discover > IDE Widget

Widget for compiling and running the source code in a web browser!