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

How to show to date instead of periods at x-axis in time series forecast plot?

我有一个简单的时间序列数据集: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 轴上显示带有真实日期的图,这对解释我的数据很有意义。
但是我迷失在网络上的大量解释中,我发现的每一种方法都失败了。
如果有人可以帮助解决这个问题,我将不胜感激。

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=' 参数来编辑标签及其位置。