x-12-arima

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

需要用R来对中国宏观经济数据进行季度调整，想使用之前比较权威的x-12-arima来做，但是不知道在R里面怎么实现。。求指导

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：X-12-ARIMA ARIMA X-12 ima Rim 中国

用x12包里的x12()就好
  
x12(object,x12Parameter=new("x12Parameter"),x12BaseInfo=new("x12BaseInfo"),...)

举个例子吧——
##Notrun:
x12path("../x12a.exe")
xts<-x12(AirPassengers)
summary(xts)
xs<-x12(new("x12Single",ts=AirPassengers))
summary(xs)
xb<-x12(new("x12Batch",list(AirPassengers,AirPassengers,AirPassengers)))
summary(xb)
#Createnewbatchobjectwith4timeseries
xb<-new("x12Batch",list(AirPassengers,AirPassengers,AirPassengers,AirPassengers))
#changetheautomdltoFALSEinall4elements
xb<-setP(xb,list(automdl=FALSE))
#changethearima.modelandarima.smodelsettingforthefirsttsobject
xb<-setP(xb,list(arima.model=c(1,1,0),arima.smodel=c(1,1,0)),1)
#changethearima.modelandarima.smodelsettingforthesecondtsobject

xb<-setP(xb,list(arima.model=c(0,1,1),arima.smodel=c(0,1,1)),2)
#changethearima.modelandarima.smodelsettingforthethirdtsobject
xb<-setP(xb,list(arima.model=c(0,1,1),arima.smodel=c(1,1,1)),3)
#changethearima.modelandarima.smodelsettingforthefourthtsobject
xb<-setP(xb,list(arima.model=c(1,1,1),arima.smodel=c(1,1,1)),4)
#runx12onallseries
xb<-x12(xb)
summary(xb)
#Setautomdl=TRUEforthefirstts
xb<-setP(xb,list(automdl=TRUE),1)
#rerunx12onallseries(thebinarieswillonlyrunonthefirstone)
xb<-x12(xb)
#summarywitholdOutput
summary(xb,oldOutput=10)
#Changetheparameterandoutputofthefirstseriesbacktothefirstrun
xb<-prev(xb,index=1,n=1)
#summarywitholdOutput(---Novalidpreviousruns.---)
summary(xb,oldOutput=10)
##End(Notrun)


呼......

DM小菜鸟 发表于 2014-12-16 15:09
用x12包里的x12()就好
  
x12(object,x12Parameter=new("x12Parameter"),x12BaseInfo=new("x12BaseInfo") ...

抱歉，很久没上论坛了，现在才看到，非常感谢~~~