尝试在R中每小时创建一个时间序列。
我有一个收集每小时车辆数量的数据框,它看起来像:
> head(df)
# A tibble: 6 x 8
interval cars vans trucks total `mean speed` `% occupation` hour
<dttm> <int> <int> <int> <int> <dbl> <dbl> <int>
1 2017-10-09 00:00:00 7 0 0 7 7.37 1. 0
2 2017-10-09 01:00:00 24 0 0 24 16.1 3. 1
3 2017-10-09 02:00:00 27 0 0 27 18.1 2. 2
4 2017-10-09 03:00:00 47 3 0 50 31.5 3. 3
5 2017-10-09 04:00:00 122 1 5 128 48.0 16. 4
6 2017-10-09 05:00:00 353 6 2 361 66.3 20. 5
> tail(df,1)
# A tibble: 1 x 8
interval cars vans trucks total `mean speed` `% occupation` hour
<dttm> <int> <int> <int> <int> <dbl> <dbl> <int>
1 2018-03-15 20:00:00 48 0 2 50 31.5 5. 20
查看starting a daily time series in R的答案,清楚地解释了如何在白天创建ts 我已将此df转换为时间序列:
ts2Start <- df$interval[1]
ts2End <- df$interval[nrow(df)]
indexPerHour <- seq(ts2Start, ts2End, by = 'hour')
由于我们一年365天,每天24小时,因此创建了ts:
> df.ts <- ts(df$total, start = c(2017, as.numeric(format(indexPerHour[1], '%j'))),
+ frequency=24*365)
其中
as.numeric(format(indexPerHour[1], '%j')))
返回282
为了验证我正在做什么,我检查了从索引获得的日期是否与数据框中的第一行相同
head(date_decimal(index(df.ts)),1)
但是我的第一个日期/时间应该是:“2017-10-09 00:00:00” 我得到:“2017-01-12 16:59:59 UTC”
看起来df.ts
系列中的第一个索引已经开始于~282 / 24
我不明白我做错了什么。启动参数如何在ts()
?
我还查了帖子:How to Create a R TimeSeries for Hourly data
建议使用xts
包。
问题在于,我只是从一本使用tslm()
并且似乎不支持xts
对象的书中学习。
我可以使用ts()
创建每小时时间序列吗?
答案 0 :(得分:0)
您应该改用 xts 库。例如:
time_index <- seq(from = as.POSIXct("2016-01-01 00:00:00"),
to = as.POSIXct("2018-10-01 00:00:00"), by = "hour")
traff = xts(df, order.by = time_index)```