r语言确定最优滞后阶数