中的R拟合值的预报缺少日期/时间分量(Fitted values in R forecast mis

2019-10-23 11:50发布

我一直在做时间序列数据(以XTS格式)多种R型号和我一直运行到同一个问题,即没有日期/时间组件的拟合值/预测,因此我无法对他们绘制在同一图表作为原始数据。

使用下面的代码,从我看到的一切在网上我应该得到在同一图中,两个线图。 不过,我只是得到绘图命令和第二线从未出现原始图。

library("stats")
library("forecast")
model_a1 <- auto.arima(GED$Mfg.Shipments.Total..USA.)
plot(GED$Mfg.Shipments.Total..USA.)
lines(fitted(model_a1), col = 2)

如果我只是绘制残差:

plot(fitted(model_a1))

我得到一个图表,但不是具有日期/时间,它只是说,“索引”,范围为0 - 在x轴1,2000。

任何想法是什么引起的? 我找不到为什么会导致没有日期/时间组件残差的XTS时间序列东西。 另外,我想发表图片,但显然有一个愚蠢的规则:如果你不张贴至少10倍,你不能发布图片进一步解释您的问题.....

str(GED) 
An ‘xts’ object on 2018-01-01/2115-12-01 containing: 
Data: num [1:1176, 1:472] NA NA NA NA NA NA NA NA NA NA ... 
  - attr(*, "dimnames")=List of 2 
        ..$ : NULL
        ..$ : chr [1:472] "SRPM..USA." "WTI" "CU..Argentina." "CU..Australia."  
... Indexed by objects of class: [Date] TZ: UTC 
   xts Attributes: NULL

Answer 1:

不要使用日期在你的情节,使用数字序列为x轴。 您可以使用日期作为标签。 尝试是这样的:

y=GED$Mfg.Shipments.Total..USA.
n=length(y)
model_a1 <- auto.arima(y)
plot(x=1:n,y,xaxt="n",xlab="")
axis(1,at=seq(1,n,length.out=20),labels=index(y)[seq(1,n,length.out=20)],
     las=2,cex.axis=.5)
lines(fitted(model_a1), col = 2)

根据数据的结果将是类似的东西:



文章来源: Fitted values in R forecast missing date / time component