Question

我有一个原始的观察数据集，仅在工作日的上午6点到晚上9点之间以5分钟的间隔进行观察。这些没有用于绘图等的日期时间信息，所以我试图创建一个日期时间向量，以便将其添加到我的数据中。即：

X425 X432 X448 1 0.07994814 0.1513559 0.1293103 2 0.08102852 0.1436480 0.1259074

到这个

X425 X432 X448 2010-05-24 06:00 0.07994814 0.1513559 0.1293103 2010-05-24 06:05 0.08102852 0.1436480 0.1259074

我已经解决了以下问题：

# using lubridate and xts
library(xts)
library(lubridate)

# sequence of 5 min intervals from 06:00 to 21:00
sttime <- hms("06:00:00")
intervals <- sttime + c(0:180) * minutes(5)

# sequence of days from 2010-05-24 to 2010-11-05
dayseq <- timeBasedSeq("2010-05-24/2010-11-05/d")

# add intervals to dayseq
dayPlusTime <- function(days, times) {
  dd <- NULL
  for (i in 1:2) {
    dd <- c(dd,(days[i] + times))}
  return(dd)
}

obstime <- dayPlusTime(dayseq, intervals)`

但是obstime将作为列表出现。 days[1] + times有效，所以我猜这与POSIXct对象连接在一起制作dd的方式有关，但是我无法弄清楚我做错了什么otr下一步去哪里。

任何帮助表示赞赏

Answer 1

base替代方案：

# create some dummy dates
dates <- Sys.Date() + 0:14

# select non-weekend days
wd <- dates[as.integer(format(dates, format = "%u")) %in% 1:5]

# create times from 06:00 to 21:00 by 5 min interval
times <- format(seq(from = as.POSIXct("2015-02-18 06:00"),
                    to = as.POSIXct("2015-02-18 21:00"),
                    by = "5 min"),
                format = "%H:%M")

# create all date-time combinations, paste, convert to as.POSIXct and sort 
wd_times <- sort(as.POSIXct(do.call(paste, expand.grid(wd, times))))

Answer 2

其中一个问题是，当分钟超过hour时，您的间隔向量不会更改60。

以下是一种可以做到这一点的方法：

#create the interval vector
intervals<-c()
for(p in 6:20){
  for(j in seq(0,55,by=5)){
    intervals<-c(intervals,paste(p,j,sep=":"))
  }      
}
intervals<-c(intervals,"21:0")

#get the days
dayseq <- timeBasedSeq("2010-05-24/2010-11-05/d")


#concatenate everything and format to POSIXct at the end
obstime<-strptime(unlist(lapply(dayseq,function(x){paste(x,intervals)})),format="%Y-%m-%d %H:%M", tz="GMT")

在R中为部分时间创建非周末时间间隔的向量

2 个答案: