def forecast(self, horiz=1, u=None):
        if u is None:
            u = np.zeros((horiz, self.nvars))
        return vargenerate(self.arhat, u, initvalues=self.y)
calculates forcast for horiz number of periods at end of sample
dataset/ETHPy150Open statsmodels/statsmodels/statsmodels/tsa/varma_process.py/_Var.forecast