R中的每小时时间序列如何ts(... start)有效?

时间:2018-05-14 12:49:50

标签: r time-series

尝试在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()创建每小时时间序列吗?

1 个答案:

答案 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)```