全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
2872 0
2017-10-27
本人对月度税收做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))”,出现以下警告,但是还是可以得出预测结果:
预测警告.png
②对预测数做plot时间序列图,出现:'xlim'值不能是无限的
plot图.png
想请教下是什么原因所致,是否与我用疏系数模型预测导致预测数据出现问题?
二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群