statsmodels.tsa.varma_process.VarmaPoly.hstackarma_minus1

VarmaPoly.hstackarma_minus1()

stack ar and lagpolynomial vertically in 2d array

this is the Kalman Filter representation, I think