def fill_coll(sim): ret = VAR(sim, exog=self.model.exog[-nobs:]).fit(maxlags=k_ar, ) ret = self.exog_irf(self.results, self.max_n) return ret.cumsum(axis=0) if cum else ret