Используйте difftime для группировки по интервалу - PullRequest
Новая
       12

Используйте difftime для группировки по интервалу

0 голосов
/ 08 февраля 2020

Как взять столбец ts класса "hms" "difftime" и использовать его для группировки остальной части кадра данных по интервалам? 

> example
# A tibble: 10 x 2
   ts        val
   <time>  <dbl>
 1 00'00" -0.7  
 2 00'01" -1.69 
 3 00'02"  0.03 
 4 00'03"  0.570
 5 00'04" -0.15 
 6 00'05" -0.34 
 7 00'06" -0.45 
 8 00'07"  0.77 
 9 00'08"  0.6  
10 00'09"  0.01 

> class(example$ts)
[1] "hms"      "difftime"

Я бы обычно использовал lubridate::floor_date, чтобы связать поле отметки времени с интервалами. Но если я попытаюсь сделать это напрямую, я получу ошибку: 

example %>% mutate(win_5s = floor_date(ts, unit = "5 seconds"))

Error in UseMethod("reclass_date", orig) : 
no applicable method for 'reclass_date' applied to an object of class "c('hms', 'difftime')"

Мой обходной путь до сих пор заключается в том, чтобы сначала преобразовать ts с as.POSIXct: 

example %>%
  mutate(ts2 = as.POSIXct(ts),
         window_5s = floor_date(ts2, "5 seconds")) %>%
  group_by(window_5s) %>%
  summarise(avg = mean(val))

# A tibble: 2 x 2
  window_5s              avg
  <dttm>               <dbl>
1 1970-01-01 00:00:00 -0.388
2 1970-01-01 00:00:05  0.118

Но это похоже на Я упускаю что-то в экосистеме lubridate - ts уже был распознан как объект time с правильными единицами, так что есть более прямой или "lubridatey" способ выполнить sh эту группировку, чем конвертировать в полную дату (с датой, которая не имеет отношения к делу или является правильной)?

dput для example: 

structure(list(ts = structure(c(0, 1, 2, 3, 4, 5, 6, 7, 8, 9), class = c("hms", 
"difftime"), units = "secs"), val = c(-0.7, -1.69, 0.03, 0.57, 
-0.15, -0.34, -0.45, 0.77, 0.6, 0.01)), row.names = c(NA, -10L
), class = c("tbl_df", "tbl", "data.frame"))

1 Ответ

1 голос
/ 12 февраля 2020

Существует функция hms::round_hms(), которая, кажется, работает: 

> test %>% mutate(hms::round_hms(ts, 5))

# A tibble: 10 x 3
   ts        val `hms::round_hms(ts, 5)`
   <time>  <dbl> <time>                 
 1 00'00" -0.7   00'00"                 
 2 00'01" -1.69  00'00"                 
 3 00'02"  0.03  00'00"                 
 4 00'03"  0.570 00'05"                 
 5 00'04" -0.15  00'05"                 
 6 00'05" -0.34  00'05"                 
 7 00'06" -0.45  00'05"                 
 8 00'07"  0.77  00'05"                 
 9 00'08"  0.6   00'10"                 
10 00'09"  0.01  00'10"

Если вы хотите получить ее, я думаю, вам нужна пользовательская функция, но источник для round_hms() дает хороший шаблон для выполнения этого: https://github.com/tidyverse/hms/blob/master/R/round.R

И вот оно: 

floor_hms <- function(x, secs) {
  vctrs::vec_restore(floor(as.numeric(x) / secs) * secs, x)
}

пример: 

> test %>% mutate(hms::round_hms(ts, 5), floor_hms(ts, 5))

# A tibble: 10 x 4
   ts        val `hms::round_hms(ts, 5)` `floor_hms(ts, 5)`
   <time>  <dbl> <time>                  <time>            
 1 00'00" -0.7   00'00"                  00'00"            
 2 00'01" -1.69  00'00"                  00'00"            
 3 00'02"  0.03  00'00"                  00'00"            
 4 00'03"  0.570 00'05"                  00'00"            
 5 00'04" -0.15  00'05"                  00'00"            
 6 00'05" -0.34  00'05"                  00'05"            
 7 00'06" -0.45  00'05"                  00'05"            
 8 00'07"  0.77  00'05"                  00'05"            
 9 00'08"  0.6   00'10"                  00'05"            
10 00'09"  0.01  00'10"                  00'05"
...