- library(forecast)
- library(tseries)
- pr_arimax=matrix(ncol=1,nrow=18,0)
- for(i in 0: 17){
- y_arimax<-ts(HK_data [c((1+i):(26+i)),2],frequency=7)
- x1<-as.matrix(HK_data[c((1+i):(26+i)),-c(1,2)])
- fit_arimax<-auto.arima(y_arimax,xreg=x1)
- x2<-as.matrix( HK_data[(27+i),-c(1,2)])
- temp=forecast(fit_arimax,xreg =x2)$mean
- temp=as.data.frame(temp)
- pr_arimax[i+1,]=temp[1,]
- }
- mape_arimax=mean(abs((pr_arimax-y)/y))
- rmse_arimax=rmse(as.data.frame(y),as.data.frame(pr_arimax))


雷达卡




京公网安备 11010802022788号







