我正在尝试运行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
包将“二十九进制$ COLDATE”从因子转换为。
一旦我将xlim
的图集设置为预测延长5年,ARIMA的预测线就不会出现。