如何在时间序列预测图中显示日期而不是 x 轴上的周期？

Question

我有一个简单的时间序列数据集：df1，格式如下：

    Item | Date     | Var
      1  | 20000101 | 12
      2  | 20000102 | 13
      3  | 20000103 | 16
      4  | 20000104 | 18
      5  | 20000105 | 29
      6  | 20000106 | 36
      ...
      ...
      ...
      ...
      365| 20001231 | 78

我的代码是：

    varts <- ts(df1$var, frequency=7)
    comp <- decompose(varts)
    plot(comp)
    fcmodel <- HoltWinters(varts)
    plot(forecast(fcmodel, h=30, level=c(80,95)))

我得到了一个简单的预测图，但 x 轴的范围是 0 到 50 而不是日期。
我想在 x 轴上显示带有真实日期的图，这对解释我的数据很有意义。
但是我迷失在网络上的大量解释中，我找到的每一种方法都失败了。
如果有人可以帮助解决这个问题，我真的很感激。

Answer 1

library(forecast)
library(lubridate)

# generate dummy data
df1 <- data.frame(seq(1, 365, 1), seq(as.Date('2000-01-01'), as.Date('2000-12-30'), 1), round(100*runif(365), 0))
names(df1) <- c('item','date','var')

varts <- ts(df1$var, frequency=7)
comp <- decompose(varts)
plot(comp)
fcmodel <- HoltWinters(varts)

plot(forecast(fcmodel, h=30, level=c(80,95)), xaxt='n')
axis(1, at=seq(0, 80,10) , las=2, labels=seq(as.Date('2000-01-01'), as.Date('2000-12-31')+weeks(30), length.out=9) )

最后两行基本上复制了这里Replace X-axis with own values 的方法。

您可以通过编辑“at=”和“lables=”参数来编辑标签及其位置。