code: def __init__(self, method="mle"):
      if method == "mle":
          self.params = [
              4.507728587708,
              0.788365037622,
              -0.358656861792,
              0.035886565643,
              -0.699600200796,
          ]
          self.aic = 813.417242529788
          self.bic = 834.546008036962
          self.arroots = [1.0991 - 1.2571j, 1.0991 + 1.2571j]
          self.maroots = [-1.1702 + 0.0000j, 1.2215 + 0.0000j]
          self.bse = [
              0.045346684035,
              0.078382496509,
              0.07004802526,
              0.069227816205,
              0.070668181454,
          ]
          self.cov_params = [
              [
                  0.0020563,
                  -2.3845e-05,
                  -6.3775e-06,
                  4.6698e-05,
                  5.8515e-05,
              ],
              [
                  -2.3845e-05,
                  0.0061438,
                  -0.0014403,
                  -0.0035405,
                  -0.0019265,
              ],
              [
                  -6.3775e-06,
                  -0.0014403,
                  0.0049067,
                  -0.00059888,
                  -0.0025716,
              ],
              [
                  4.6698e-05,
                  -0.0035405,
                  -0.00059888,
                  0.0047925,
                  0.0022931,
              ],
              [
                  5.8515e-05,
                  -0.0019265,
                  -0.0025716,
                  0.0022931,
                  0.0049940,
              ],
          ]
          self.hqic = 821.920952341460
          self.llf = -400.708621264894
          self.resid = residsc_mle[:, 3]
          self.fittedvalues = yhatc_mle[:, 3]
          self.pvalues = [0.0000, 8.48e-24, 3.05e-07, 0.6042, 4.17e-23]
          self.tvalues = [99.41, 10.06, -5.120, 0.5184, -9.900]
          self.sigma2 = 1.196309833136**2
      elif method == "css":
          # NOTE: params, bse, cov_params, tvalues taken from R;
          #   commented-out versions below are from
          #   [TODO: Finish this sentence]
          # self.params = [2.571274348147, 0.793030965872, -0.363511071688,
          #                0.033543918525, -0.702593972949]
          self.params = [
              4.507207454494,
              0.793055048760,
              -0.363521072479,
              0.033519062805,
              -0.702595834943,
          ]
          self.aic = 806.807171655455
          self.bic = 827.887744132445
          # self.bse = [0.369201481343, 0.076041378729, 0.070029488852,
          #             0.062547355221, 0.068166970089]
          self.bse = [
              0.0446913896589,
              0.0783060902603,
              0.0697866176073,
              0.0681463870772,
              0.068958002297,
          ]
          # self.cov_params = [
          #   [0.13631, -0.017255, -0.012852, 0.014091, 0.017241],
          #   [-0.017255, 0.0057823, -0.0020013, -0.0026493, -0.0014131],
          #   [-0.012852, -0.0020013, 0.0049041, -0.00042960, -0.0023845],
          #   [0.014091, -0.0026493, -0.00042960, 0.0039122, 0.0022028],
          #   [0.017241, -0.0014131, -0.0023845, 0.0022028, 0.0046467]
          # ]
          self.cov_params = [
              [
                  1.99732030964e-03,
                  -2.22972353619e-05,
                  -0.000009957435095,
                  4.64825632252e-05,
                  5.98134427402e-05,
              ],
              [
                  -2.22972353619e-05,
                  6.13184377186e-03,
                  -0.001435210779968,
                  -3.47284237940e-03,
                  -1.95077811843e-03,
              ],
              [
                  -9.95743509501e-06,
                  -1.43521077997e-03,
                  0.004870171997068,
                  -6.54767224831e-04,
                  -2.44459075151e-03,
              ],
              [
                  4.64825632252e-05,
                  -3.47284237940e-03,
                  -0.000654767224831,
                  4.64393007167e-03,
                  2.34032945541e-03,
              ],
              [
                  5.98134427402e-05,
                  -1.95077811843e-03,
                  -0.002444590751509,
                  2.34032945541e-03,
                  4.75520608091e-03,
              ],
          ]
          self.arroots = [1.0908 - 1.2494j, 1.0908 + 1.2494j]
          self.maroots = [-1.1694 + 0.0000j, 1.2171 + 0.0000j]
          self.hqic = 815.293412134796
          self.llf = -397.403585827727
          self.resid = residsc_css[2:, 3]
          self.fittedvalues = yhatc_css[2:, 3]
          self.pvalues = [
              3.30e-12,
              1.83e-25,
              2.09e-07,
              0.5918,
              6.55e-25,
          ]
          # self.tvalues = [6.964, 10.43, -5.191, 0.5363, -10.31]
          self.tvalues = [
              100.851808120009,
              10.127629231947,
              -5.209036989363,
              0.491868523669,
              -10.188749840927,
          ]
          self.sigma2 = 1.201409294941**2
  
imports: import os
  from numpy import genfromtxt
label_desc: Use image node_id 1 to create a new Y_arma22c object  with example: obj = Y_arma22c(method)