I have daily visitors data for the last 10 years. I want to do some basic tests like which is the busiest day, which is the busiest month, busiest week etc. I used auto.arima
function with argument xreg
to find out the coefficients of all the days of the week, week of the month. This is the output I got:
> summary(arima1)
Series: dailysea
ARIMA(1,1,2)
Coefficients:
ar1 ma1 ma2 Sun Mon Tue Wed Thu
-0.1250 -0.4506 -0.3712 -1466.6853 -3623.175 -3895.0555 -3722.146 -3327.4288
s.e. 0.1207 0.1117 0.0891 325.7253 386.738 379.8793 379.883 386.7512
Fri
-2146.910
s.e. 325.736
sigma^2 estimated as 7776468: log likelihood=-6808.5
AIC=13637 AICc=13637.31 BIC=13682.92
Training set error measures:
ME RMSE MAE MPE MAPE MASE ACF1
Training set 59.63838 2784.809 1952.625 -Inf Inf 0.8353728 -0.001839015
Can I use these coefficients to conclude that Saturday is the busiest followed by Sunday, Friday etc.? Also I have infinite MAPE which is not making sense to me.
auto.arima
? That means that the ARIMA(0,1,2) had a higher AIC (or whichever criterion was used) value; recall thatauto.arima
does search the "neighbouring" models, so it must have evaluated ARIMA(0,1,2) and found it inferior to ARIMA(1,1,2). This makes ARIMA(0,1,2) an inferior choice. Also recall that variable significance is not the right tool for variable selection (Rob J. Hyndman had a post on that somewhere, if I remember correctly). $\endgroup$