已知方程 yt=y(t-1) +x(t) +z(t)+dx(t)+x(t-1)+z(t-1)+dz(t)+ dx(t-1)+dz(t-1)
如何编程得到BIC值最小的模型

当然不编程也可以一个个算，问题是有没有更好的办法