I'm trying to fit an ARIMA model to monthly CPI. SO far, I've established stationarity using a DF test of the differenced CPI.
I have then proceeded to find the ACF and PACF of the differenced variable, which are as follows:
These are very peculiar to me (i.e. 3 irregular and significant lags after the first for both ACF and PACF). What can I make of this? Is this due to seasonality?
Regardless, I proceeded to fit the best ARIMA model using the auto.arima feature in R which picks an ARIMA model with the lowest information criterion. What I got was an ARIMA(2,1,2) as follows:
All coefficients are signficant, and the BIC would be lowest here, so I don't see any problems. Please correct me if I'm wrong in saying that there is nothing wrong here
I then proceeded to do some residual diagnostics, namely, doing a Box-Ljung test of the residuals with the following lags/dfs (not totally sure what this does) but anyways here are the results (and also the code I used):
So for df=1,4 and 5 (arbitrarily chosen, is there a better way to do this?) I get p > 0.05 which is what I am looking for. But for df= 36, 12 and 13 I get p <0.05 which is indicating that there are correlations among residuals indicating that the ARIMA(2,1,2) is not a good model.
Again, this is not what I'd like to see as we see a few significant correlations here.
So, should I accept this model regardless of the fact that there were a few indications that it wasn't ideal (ACF, PACF, Box-Ljung test, and ACF of residuals)? What would you do next? Does the ACF or PACF indicate seasonality perhaps?
Any help or guidance is appreciated! I'm still new to ARIMA modelling.