def err_band_sz2(self, orth=False, repl=1000, signif=0.05,
                     seed=None, burn=100, component=None):
        model = self.model
        periods = self.periods
        if orth:
            irfs = self.orth_irfs
        elif svar:
            irfs = self.svar_irfs
        else:
            irfs = self.irfs
        neqs = self.neqs
        irf_resim = model.irf_resim(orth=orth, repl=repl, T=periods, seed=seed,
                                   burn=100)

        W, eigva, k = self._eigval_decomp_SZ(irf_resim)

        if component != None:
            if np.shape(component) != (neqs,neqs):
                raise ValueError("Component array must be " + str(neqs) + " x " + str(neqs))
            if np.argmax(component) >= neqs*periods:
                raise ValueError("Atleast one of the components does not exist")
            else:
                k = component

        gamma = np.zeros((repl, periods+1, neqs, neqs))
        for p in range(repl):
            for i in range(neqs):
                for j in range(neqs):
                    gamma[p,1:,i,j] = W[i,j,k[i,j],:] * irf_resim[p,1:,i,j]

        gamma_sort = np.sort(gamma, axis=0) #sort to get quantiles
        indx = round(signif/2*repl)-1,round((1-signif/2)*repl)-1

        lower = np.copy(irfs)
        upper = np.copy(irfs)
        for i in range(neqs):
            for j in range(neqs):
                lower[:,i,j] = irfs[:,i,j] + gamma_sort[indx[0],:,i,j]
                upper[:,i,j] = irfs[:,i,j] + gamma_sort[indx[1],:,i,j]

        return lower, upper
Read key=value pairs from a file, into a dict.
dataset/ETHPy150Open statsmodels/statsmodels/statsmodels/tsa/vector_ar/irf.py/IRAnalysis.err_band_sz2--crossover with--dataset/ETHPy150Open coreemu/core/daemon/core/misc/utils.py/readfileintodict