如何用R对本数据进行时间序列分析？

附件xls为月度数据，具有明显的趋势性和周期性，doc为相关截图。

期待大家一起会诊，看如何才能找出最佳阶数P、D、Q


利用附件中数据做时间序列分析，用的是arima模型。初始参数只用arima(0,0,0)(0,0,0)，显示acf和pacf图如下


再用arima(0,1,0)(0,0,0)，显示acf和pacf图如下：


发现还是不行，在12，24，36等12倍数的系数明显较大，于是再调整为arima(0,1,0)(0,1,0)，显示acf和pacf图如下

从图上来看，仍然未达到稳定要求，如果系统自己选的话，模型参数为ARIMA(0,1,1) (0,1,1)，显示acf和pacf图如下，但仍有个别acf、pacf过大，超出阴影部分

各位良医谁能一诊？