R日付時間繁殖時間/時間

Question

私はそれが難しくないと確信していますが、私は時間を繁殖時間を合計する方法を見つけることができません。私は彼らの巣に到着し（AR）、そこから離れる（LV）鳥の数千のイベントデータ（繁殖）を持っています。 arriving-と-イベントを残す対応する時間の時間で与えられ、次のようになりされています

times = as.POSIXct(c("2007-07-11 22:47:21 UTC", "2007-07-11 22:58:39 UTC", "2007-07-11 22:58:48 UTC", "2007-07-11 23:57:45 UTC", "2007-07-12 02:29:52 UTC", "2007-07-12 03:46:23 UTC", "2007-07-12 03:46:36 UTC", "2007-07-12 04:28:54 UTC", "2007-07-12 04:29:03 UTC", "2007-07-12 05:36:38 UTC"), tz = "UTC") 
breeding = c("AR", "LV", "AR", "LV", "AR", "LV", "AR", "LV", "AR", "OFF") 

は今、私はこれらの

のような時間単位の休憩に基づいて、鳥が巣に費やす時間あたり何分数計算します

cut(times, breaks = "hour") 

右端も含めなければなりません。

difftimeと合計してみましたが、もちろん、時間と日付の区切りで切り取られませんでした。結果はどういうものかは次のようになります：

Hour  fraction 

22-23  12min 

23-00  57min 

00-01   0min 

01-02   0min 

02-03  31min 

03-04  46min 

04-05   1min 

05-06  24min 

ありがとうございました！

Answer 1

dplyr,tidyrおよびlubridateからの溶液。

library(dplyr) 
library(tidyr) 
library(lubridate) 

dt2 <- dt %>% 
    mutate(breeding = ifelse(breeding %in% "OFF", "LV", breeding)) %>% 
    mutate(ID = rep(1:(n()/2), each = 2)) %>% 
    spread(breeding, times) %>% 
    mutate(DiffTime = difftime(LV, AR, "mins")) 
dt2 
    ID     AR     LV  DiffTime 
1 1 2007-07-11 22:47:21 2007-07-11 22:58:39 11.30000 mins 
2 2 2007-07-11 22:58:48 2007-07-11 23:57:45 58.95000 mins 
3 3 2007-07-12 02:29:52 2007-07-12 03:46:23 76.51667 mins 
4 4 2007-07-12 03:46:36 2007-07-12 04:28:54 42.30000 mins 
5 5 2007-07-12 04:29:03 2007-07-12 05:36:38 67.58333 mins 

を次のように私たちは、あなたがその日の時間をお知りになりたい場合は、次のように、我々はさらに、日付と時間の情報を抽出することができ、時間差を計算することができます。

dt3 <- dt2 %>% 
    mutate(AR_Date = as.Date(AR), LV_Date = as.Date(LV), 
     AR_Hour = hour(AR), LV_Hour = hour(LV)) 
dt3 

ここから、データをさらに要約する方法を決定できます。

ID     AR     LV  DiffTime AR_Date LV_Date AR_Hour LV_Hour 
1 1 2007-07-11 22:47:21 2007-07-11 22:58:39 11.30000 mins 2007-07-11 2007-07-11  22  22 
2 2 2007-07-11 22:58:48 2007-07-11 23:57:45 58.95000 mins 2007-07-11 2007-07-11  22  23 
3 3 2007-07-12 02:29:52 2007-07-12 03:46:23 76.51667 mins 2007-07-12 2007-07-12  2  3 
4 4 2007-07-12 03:46:36 2007-07-12 04:28:54 42.30000 mins 2007-07-12 2007-07-12  3  4 
5 5 2007-07-12 04:29:03 2007-07-12 05:36:38 67.58333 mins 2007-07-12 2007-07-12  4  5 

Answer 2

面白かったので、私は本当に答えを出したいと思っていました。私はだと思います。だと思います。

あなたの意見をチェックしてください！

は次のようになります。

  Date hours Total 
1: 2007-07-11 22 01:10:15 
2: 2007-07-12  2 01:16:31 
3: 2007-07-12  3 00:42:18 
4: 2007-07-12  4 01:07:35 

ので、午後09時で、彼らは時間かけて過ごすことになります。鳥は2時まで到着せず、1時間以上も費やします。

最終出力はtimesです。

#---Load the required libraries 
require(chron)    #---For time manipulation 
require(data.table)   #---Calculate and store data.tables 
require(splitstackshape) #---Split columns by a delimiter 

#---Input Data 
times = as.POSIXct(c("2007-07-11 22:47:21 UTC", "2007-07-11 22:58:39 UTC", "2007-07-11 22:58:48 UTC", "2007-07-11 23:57:45 UTC", "2007-07-12 02:29:52 UTC", "2007-07-12 03:46:23 UTC", "2007-07-12 03:46:36 UTC", "2007-07-12 04:28:54 UTC", "2007-07-12 04:29:03 UTC", "2007-07-12 05:36:38 UTC"), tz = "UTC") 
breeding = c("AR", "LV", "AR", "LV", "AR", "LV", "AR", "LV", "AR", "OFF") 

#---The times need to be a data.table for cSplit to work 
times <- as.data.table(times) 

#---This splits your time input into "Date" and "Time" columns 
times <- cSplit(times, "x", sep = " ") 
setnames(times, c("x_1", "x_2"),c("Date", "Time")) 
times[] <- lapply(times, as.character) 

#---Bring in "AR" or "LV" into a fresh data.table 
Data <- cbind(times, breeding) 

#---Format the time column for computation 
Data$Time <- chron(times = Data$Time, format = 'h:m:s') 

#rm(times,breeding) 

#---Arrival times 
arrive <- Data[breeding == "AR",] 
setnames(arrive, c("Time"), c("ArriveTime")) 
arrive[,3] <- NULL 

#---Departure times 
leave <- Data[breeding != "AR",] 
setnames(leave, c("Time"), c("LeaveTime")) 
leave[,c(1,3)] <- NULL 

#---But them together 
times <- cbind(arrive, leave) 
#rm(arrive, leave) 

#---Calculate elapsed times 
times <- times[, Elapsed := LeaveTime - ArriveTime] 
#rm(times3,Data) 

#---Sum elapsed time by the hour of the Arrival time 
times <- times[, .(Total = sum(Elapsed)), by = .(Date,hours(ArriveTime))] 
times