如何用原始時間序列繪制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)

通過您的數據樣本，我可以更具體一些。