如何用原始时间序列绘制Arima模型？

Question

我绘制原始数据集

ggplot(piz, aes(Date, Price.Dollars.per.Thousand.Cubic.Feet)) + geom_line()

后来我提取了一列并进行了ARIMA建模。

fit1
Series: a1 
ARIMA(4,1,1) 
plot(forecast(fit,h=48))

我的形象

如果我选择自动曝光解决方案

> a1=ts(a1)
> autoplot(a1) + forecast::autolayer(fit1)
Error in UseMethod("autolayer") : 
  no applicable method for 'autolayer' applied to an object of class "c('ARIMA', 'Arima')".

我正在添加样本数据集

Date     Price Dollars per Thousand Cubic Feet
Jan-2002    3.1
Feb-2002    2.86
Mar-2002    3.37
Apr-2002    3.8
May-2002    3.78
Jun-2002    3.61

如何将它们一起绘制？

Answer 1

尝试这样的事情：

  db<-data.frame(x=c(1:1000)+round(runif(1000,min = 5, max=100),0),
  date=seq(as.Date("2000/1/1"), by = "day", length.out = 1000))

  fit<-auto.arima(y = db$x)

  plot(fit$x~as.Date(db$date),col="red",type="l",xaxt = "n",)
  lines(fitted(fit)~as.Date(db$date),col="blue")
  axis(1, db$date, format(db$date, "%b %d"), cex.axis = .7)

根据您的数据样本，我可能会更具体。