2016-07-29 17 views
1

我有一組觀察對,我想標記它們的time之間的間隔。 (在真實數據集,這些觀察對代表進入和退出麥克風校準。)分組的dplyr :: mutate不正確地返回潤溼間隔的向量

# R version 3.2.3 
library(lubridate) ## Version 1.5.6 
library(dplyr) ## Version 0.5.0 

data <- data.frame(
    group = c(1,1,2,2,3,3), 
    type = rep(c("start", "end"), 3), 
    time = ymd_hms("2016-06-01 01:00:00") + c(0,1,3,6,12,18), 
    someAttribute = runif(6) 
) 

data 
## group type    time someAttribute 
## 1  1 start 2016-06-01 01:00:00  0.2540128 
## 2  1 end 2016-06-01 01:00:01  0.6845078 
## 3  2 start 2016-06-01 01:00:03  0.3576477 
## 4  2 end 2016-06-01 01:00:06  0.1223582 
## 5  3 start 2016-06-01 01:00:12  0.2715063 
## 6  3 end 2016-06-01 01:00:18  0.6392607 

我包括虛擬someAttribute在這個例子中強調的是,像tidyr::spread()一個簡單的解決方案將使屬性的混亂,屬於每行在data

我有一個功能,使間隔,並通過我組dplyr應用它:

makeTwoIntervals <- function(twoDatetimes) { 
    return(rep(interval(twoDatetimes[1], twoDatetimes[2]), 2)) 
} 

data2 <- data %>% group_by(group) %>% mutate(intervals = makeTwoIntervals(time)) 

data2$intervals 
## [1] 2016-06-01 01:00:00 UTC--2016-06-01 01:00:01 UTC 
## [2] 2016-06-01 01:00:00 UTC--2016-06-01 01:00:01 UTC 
## [3] 2016-06-01 01:00:00 UTC--2016-06-01 01:00:03 UTC 
## [4] 2016-06-01 01:00:00 UTC--2016-06-01 01:00:03 UTC 
## [5] 2016-06-01 01:00:00 UTC--2016-06-01 01:00:06 UTC 
## [6] 2016-06-01 01:00:00 UTC--2016-06-01 01:00:06 UTC 

這些值是不是我所期待得到的。正確的時間傳遞給我的函數,它創建正確的兩元素間隔向量來返回,但是當這個向量傳遞迴mutate時,會發生一些不好的事情。仔細看:

str(data2$intervals) 
## Formal class 'Interval' [package "lubridate"] with 3 slots 
## [email protected] .Data: num [1:6] 1 1 3 3 6 6 
## [email protected] start: POSIXct[1:2], format: "2016-06-01 01:00:00" "2016-06-01 01:00:00" 
## [email protected] tzone: chr "UTC" 

我不清楚這裏出了什麼問題。這些是我想看到的結果:

## Desired result of data2$intervals: 
## [1] 2016-06-01 01:00:00 UTC--2016-06-01 01:00:01 UTC 
## [2] 2016-06-01 01:00:00 UTC--2016-06-01 01:00:01 UTC 
## [3] 2016-06-01 01:00:03 UTC--2016-06-01 01:00:06 UTC 
## [4] 2016-06-01 01:00:03 UTC--2016-06-01 01:00:06 UTC 
## [5] 2016-06-01 01:00:12 UTC--2016-06-01 01:00:18 UTC 
## [6] 2016-06-01 01:00:12 UTC--2016-06-01 01:00:18 UTC 

任何人都可以提供一些洞察什麼出錯了,或者我怎麼能達到預期的結果?我濫用mutate,還是它只是沒有設計處理像lubridate::Interval這樣的對象?

+0

這是最近在data.table項目頁面上提交的問題,[#1777](https://github.com/Rdatatable/data.table/issues/1777),我建議在那裏使用一種解決方法。 – Arun

回答

1

這是一個基於@ Arun的data.table解決方法(#1777)一個解決辦法,但在dplyr語言:

data2 <- data %>% group_by(group) %>% mutate(ranges = list(range(time))) 
data3 <- data2 %>% mutate(intervals = list(interval(ranges[[1]][1], ranges[[1]][2]))) 
data3$intervals2 <- do.call("c", data3$intervals) 

data3$intervals2 
## [1] 2016-06-01 01:00:00 UTC--2016-06-01 01:00:01 UTC 
## [2] 2016-06-01 01:00:00 UTC--2016-06-01 01:00:01 UTC 
## [3] 2016-06-01 01:00:03 UTC--2016-06-01 01:00:06 UTC 
## [4] 2016-06-01 01:00:03 UTC--2016-06-01 01:00:06 UTC 
## [5] 2016-06-01 01:00:12 UTC--2016-06-01 01:00:18 UTC 
## [6] 2016-06-01 01:00:12 UTC--2016-06-01 01:00:18 UTC 

不能完全令人滿意,但它的作品。感謝您的提示,@阿倫。