'IRF Monte Carlo integrated error bands of cumulative effect'

def cum_errband_mc(self, orth=False, repl=1000, signif=0.05, seed=None, burn=100):

    model = self.model
    periods = self.periods
    return model.irf_errband_mc(orth=orth, repl=repl, T=periods, signif=signif, seed=seed, burn=burn, cum=True)
