tsa
stack ar and lagpolynomial vertically in 2d array
this is the Kalman Filter representation, I think
VarmaPoly.hstackarma_minus1