,求大神解答,很着急!!拜谢拜谢。。附代码与数据
- data1<-read.csv("data1.csv",header=T,sep=",",fileEncoding="UTF-8")->ad
- p<-as.numeric(c(data1[1:123,2]))
- dapts<-ts(p,frequency=360,start=c(2014,3,1))
- plot.ts(dapts,xlab="time")
- adf.test(dapts,alternative="stationary",k=5)
- dapdiff<-diff(dapts, differences=7)#去周期差分
- plot.ts(dapdiff)
- adf.test(dapdiff,alternative="stationary",k=5)
- #p=0.01<0.05 通过检验
- #acf与pacf的检验,初期定阶
- acf(dapdiff, lag.max=20)
- acf(dapdiff, lag.max=20,plot=FALSE)
- pacf(dapdiff, lag.max=20)
- pacf(dapdiff, lag.max=20,plot=FALSE)
- dap<-arima(dapdiff,order=c(4,0,3),seasonal=list(order=c(0,1,1),period=7))
- dapforecast <- forecast.Arima(dap,h=30,level=c(99.5))
。


雷达卡




京公网安备 11010802022788号







