class_docstr: Impulse response analysis class. Computes impulse responses, asymptotic
  standard errors, and produces relevant plots
  
  Parameters
  ----------
  model : VAR instance
  
  Notes
  -----
  Using Lütkepohl (2005) notation
class_docstr_tok: ["Impulse","response","analysis","class",".","Computes","impulse","responses",",","asymptotic","standard","errors",",","and","produces","relevant","plots","Parameters","--","--","--","--","--","model",":","VAR","instance","Notes","--","--","-","Using","L\u00fctkepohl","(","2005",")","notation"]
docstr: Returns
  -------
  W: array of eigenvectors
  eigva: list of eigenvalues
  k: matrix indicating column # of largest eigenvalue for each c_i,j
docstr_tok: ["Returns","--","--","--","-","W",":","array","of","eigenvectors","eigva",":","list","of","eigenvalues","k",":","matrix","indicating","column","#","of","largest","eigenvalue","for","each","c_i",",","j"]
code: def _eigval_decomp_SZ(self, irf_resim):
      """
      Returns
      -------
      W: array of eigenvectors
      eigva: list of eigenvalues
      k: matrix indicating column # of largest eigenvalue for each c_i,j
      """
      neqs = self.neqs
      periods = self.periods
  
      cov_hold = np.zeros((neqs, neqs, periods, periods))
      for i in range(neqs):
          for j in range(neqs):
              cov_hold[i, j, :, :] = np.cov(
                  irf_resim[:, 1:, i, j], rowvar=0
              )
  
      W = np.zeros((neqs, neqs, periods, periods))
      eigva = np.zeros((neqs, neqs, periods, 1))
      k = np.zeros((neqs, neqs), dtype=int)
  
      for i in range(neqs):
          for j in range(neqs):
              (
                  W[i, j, :, :],
                  eigva[i, j, :, 0],
                  k[i, j],
              ) = util.eigval_decomp(cov_hold[i, j, :, :])
      return W, eigva, k
  
imports: import numpy
  import numpy.linalg
  import scipy.linalg
  from statsmodels.tools.decorators import cache_readonly
  import statsmodels.tsa.tsatools
  import statsmodels.tsa.vector_ar.plotting
  import statsmodels.tsa.vector_ar.util
label_desc: Use image node_id 7 for calling the IRAnalysis obj's underlying member method code with example usage: obj._eigval_decomp_SZ(irf_resim) and returns: W,  eigva,  k