对ARIMA疏系数模型预测序列做plot图，返回：'xlim'值不能是无限的

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

本人对月度税收做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'值不能是无限的

想请教下是什么原因所致，是否与我用疏系数模型预测导致预测数据出现问题？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：ARIMA 模型预测 plot Rim ima R语言 ARIMA ARIMA预测 SARIMA模型