我现在有一个时间序列数据，是ARMA（2,1）的，生成它的程序如下：
我现在需要利用proc arima在多组参数（p，q）下对这个数据进行分析，获取每一组（p.q）值所对应的AIC值，再生成一组表格，以p为行，q为列，表中是对应的AIC值。我写了一个很烂的宏来完成上述事情，因为这个宏提取每一组参数的AIC值并生成表格做得很麻烦。程序如下：
提取aic值这个步骤应该不用单独的data步，而可以用ods把每组的aic值输出到某个数据集，无奈我不会=。=
然后最后生成的表格也很牵强，很丑，有一行total去不掉。
请问各位应该怎么优化？