http://www.cqvip.com/qk/95364X/200605/24251970.html这里有一篇利用脉冲响应数据做ARMA的。
五个变量应该考虑多元回归吧（协整然后回归），ARMA是自回归滑动平均吧，也就是变量和自己的滞后期进行回归的。