ai （ i=1，2，3，4），是模型AR（4）各滞后项的系数，P值显示a2不显著，此时要对系数重新估计，即估计a1，a3，a4.

请问：在Splus中，这种对AR模型特定阶数估计的程序该怎么编呢？函数arima.mle（）可以实现吗？