我有一个DateTime数据帧用于睡眠数据。我想根据大多数睡眠发生在截止时间的上午8点和晚上8点,将“混合”变量更改为“白天”睡眠或“晚上”睡眠。
#Current database
id<-c("m1","m1","m1","m2","m2","m2","m3","m4","m4")
x<-c("2020-01-03 10:00:00","2020-01-03 16:00:00","2020-01-03 19:20:00","2020-01-05 10:00:00","2020-01-05 15:20:00","2020-01-05 20:50:00","2020-01-06 06:30:00","2020-01-08 06:30:00","2020-01-08 07:50:00")
start<-strptime(x,"%Y-%m-%d %H:%M:%S")
y<-c("2020-01-03 16:00:00","2020-01-03 19:20:00","2020-01-03 20:50:00","2020-01-05 15:20:00","2020-01-05 20:50:00","2020-01-05 22:00:00","2020-01-06 07:40:00","2020-01-08 07:50:00","2020-01-08 08:55:00")
end<-strptime(y,"%Y-%m-%d %H:%M:%S")
mydata<-data.frame(id,start,end)
#Current ouput
df1<-mydata %>%
mutate_at(vars(start, end), ymd_hms) %>%
mutate(start_hour = hour(start),
end_hour = hour(end),
day.night = case_when(start_hour >= 8 & end_hour >= 8 & end_hour < 20 ~ "day",
start_hour >= 20 & (end_hour < 8 | end_hour <= 23) |
(start_hour < 8 & end_hour < 8)~ "night",
TRUE ~ "mixed"))
id start end start_hour end_hour day.night
1 m1 2020-01-03 10:00:00 2020-01-03 16:00:00 10 16 day
2 m1 2020-01-03 16:00:00 2020-01-03 19:20:00 16 19 day
3 m1 2020-01-03 19:20:00 2020-01-04 20:50:00 19 20 mixed
4 m2 2020-01-05 10:00:00 2020-01-05 15:20:00 10 15 day
5 m2 2020-01-05 15:20:00 2020-01-05 20:50:00 15 20 mixed
6 m2 2020-01-05 20:50:00 2020-01-05 22:00:00 20 22 night
7 m3 2020-01-06 06:30:00 2020-01-06 07:40:00 6 7 night
8 m4 2020-01-08 06:30:00 2020-01-08 07:50:00 6 7 night
9 m4 2020-01-08 07:50:00 2020-01-08 08:55:00 7 8 mixed
当前,当睡眠中断临界值时,新变量输出将设置为“混合”。
编辑:我想根据睡眠时间在哪里混合使用睡眠数据,即在第3行中:40分钟是白天,而50分钟是晚上,所以现在将是= night,在第5行:4hrs白天是40分钟,晚上是50分钟,所以现在=天。
使用lubridate
和dplyr
,您可以使用进行白天和黑夜的分类if_else
。我还根据评论调整了数据。
library(lubridate)
library(dplyr)
#data
id <- c("m1", "m1", "m1", "m2", "m2", "m2", "m3", "m4", "m4")
x <- c("2020-01-03 10:00:00", "2020-01-03 16:00:00", "2020-01-03 19:20:00", "2020-01-05 10:00:00", "2020-01-05 15:20:00", "2020-01-05 20:50:00", "2020-01-06 06:30:00", "2020-01-08 06:30:00", "2020-01-08 07:50:00")
start <- strptime(x, "%Y-%m-%d %H:%M:%S")
y <- c("2020-01-03 16:00:00", "2020-01-03 19:20:00", "2020-01-03 00:50:00", "2020-01-05 15:20:00", "2020-01-05 20:50:00", "2020-01-05 22:00:00", "2020-01-06 07:40:00", "2020-01-08 07:50:00", "2020-01-08 08:55:00")
end <- strptime(y, "%Y-%m-%d %H:%M:%S")
mydata <- data.frame(id, start, end)
#start and end times
daystart <- as.POSIXct('08:00:00', format = "%T")
nightstart <- as.POSIXct('20:00:00', format = "%T")
dayend <- as.POSIXct('19:59:00', format = "%T")
nightend <- as.POSIXct('07:59:00', format = "%T")
df1 <- mydata %>% mutate(start1 = as.POSIXct(sub("\\d+-\\d+-\\d+", Sys.Date(), start)),
end1 = as.POSIXct(sub("\\d+-\\d+-\\d+", Sys.Date(), end)),
day = ifelse(start1 > daystart & start1 < dayend & end1 < dayend & end1 > daystart, as.interval(start1, end1),
ifelse(start1 > daystart & start1 < dayend & end1 < dayend & end1 < daystart, as.interval(start1, dayend),
ifelse(start1 > daystart & start1 < dayend & end1 > dayend, as.interval(start1, dayend),
ifelse(end1 > daystart & end1 < dayend, as.interval(daystart, end1), 0)))),
night = ifelse(end1 > nightstart, as.interval(nightstart, end1),
ifelse(start1 < nightend & end1 > nightend, as.interval(start1, nightend),
ifelse(start1 < nightend & end1 < nightend, as.interval(start1, end1),
ifelse(start1 > nightstart & end1 < nightend, as.interval(start1, end1),
ifelse(start1 < nightstart & end1 < daystart, as.interval(nightstart, end1), 0))))),
day.night = ifelse(abs(day) > abs(night), "day", "night"))
df1 %>% select(names(mydata), day.night)
# id start end day.night
#1 m1 2020-01-03 10:00:00 2020-01-03 16:00:00 day
#2 m1 2020-01-03 16:00:00 2020-01-03 19:20:00 day
#3 m1 2020-01-03 19:20:00 2020-01-03 00:50:00 night
#4 m2 2020-01-05 10:00:00 2020-01-05 15:20:00 day
#5 m2 2020-01-05 15:20:00 2020-01-05 20:50:00 day
#6 m2 2020-01-05 20:50:00 2020-01-05 22:00:00 night
#7 m3 2020-01-06 06:30:00 2020-01-06 07:40:00 night
#8 m4 2020-01-08 06:30:00 2020-01-08 07:50:00 night
#9 m4 2020-01-08 07:50:00 2020-01-08 08:55:00 day
本文收集自互联网,转载请注明来源。
如有侵权,请联系 [email protected] 删除。
我来说两句