ARIMA прогнозирование с датами вместо последовательного времени

Question

Я пытаюсь запустить прогнозирование ARIMA, но я могу заставить модель работать только с последовательными временными рядами, а не с датами или временем. Как узнать фактическое время (годы) по оси X?

Я попробовал следующее:

fish<- arima(log10(twentynine$E.coli), 
         order=c(6,0,0),
         include.mean=F) 

plot(log10(twentynine$E.coli), xlim=c(0,140), ylim=c(1,4), ylab="log(E.coli)", xlab="Time", main="Amazon Creek at 29th Avenue",  type="o")

fish.pred<- predict(fish, n.ahead=60)
lines(fish.pred$pred, col="blue")
lines(fish.pred$pred+2*fish.pred$se,col="red")
lines(fish.pred$pred-2*fish.pred$se,col="red")
abline(h=log10(406), col="red", lwd=2)

Этот код работает, но результаты в последовательности по времени. То, что я хочу с годами, как:

plot(log10(twentynine$E.coli)~twentynine$COLDATE)

Я использую lubridate пакет для преобразования "twentynine$COLDATE" из фактора в дату.

Как только я получу сюжет с xlim для прогнозирования, продленного на 5 лет, строки прогноза от ARIMA не будут отображаться.