本人对月度税收做ARIMA疏系数模型预测,对预测序列做plot图,返回:'xlim'值不能是无限的
模型代码如下:
tax<-read.csv("tax.csv",header=T)
#总税收预测
sumtax1<-tax[,2]
sumtax<-ts(sumtax1,frequency=12,start=c(2013,1))#转换为时间序列数据
sumtax
acf(sumtax,lag.max=24)
plot(diff(sumtax),ylab='First difference of sumtax',xlab='time')
acf(diff(sumtax),lag.max=24)
pacf(diff(sumtax),lag.max=24)
plot(diff(diff(sumtax),lag=12),ylab='First and Seasonal difference of sumtax',xlab='time')
acf(diff(diff(sumtax),lag=12),lag.max=24,ci.type='ma')#考虑疏系数模型1,6
pacf(diff(diff(sumtax),lag=12),lag.max=24,ci.type='white') #考虑疏系数模型1,5
arima11<-arima(sumtax, order=c(6,1,5),seasonal=list(order=c(1,1,1),period=12),transform.pars=F,fixed=c(NA,0,0,0,0,NA,NA,0,0,0,NA,NA,NA),method="ML")
arima11
#AIC=1869 虽然sar1,ama1参数显著性不高,但是去掉后模型性能降低,故不去掉
Box.test(sumtax,lag=5,type='Ljung')
Box.test(arima11$residuals,lag=5,type='Ljung')
forecast <- forecast(arima11,h=3,level=c(99.5))
forecast
plot(forecast)
plot(arima11$x,col="red") #原始数据曲线
lines(fitted(arima11),col="blue")
运行中出现的问题
①对模型未来3期进行预测时,代码为“forecast <- forecast(arima11,h=3,level=c(99.5))”,出现以下警告,但是还是可以得出预测结果:
②对预测数做plot时间序列图,出现:'xlim'值不能是无限的
想请教下是什么原因所致,是否与我用疏系数模型预测导致预测数据出现问题?