code: def __init__(self, method="mle"):
      if method == "mle":
          self.params = [
              1.062980233899,
              0.768972932892,
              -0.264824839032,
              -0.279936544064,
              0.756963578430,
              0.231557444097,
          ]
          self.aic = 686.468309958027
          self.bic = 711.118536383063
          self.arroots = [
              1.0077 + 0j,
              0.3044 - 0.9793j,
              0.3044 + 0.9793j,
              -1.2466 + 0j,
          ]
          self.maroots = [-4.3186 + 0.0j]
          self.bse = [
              2.781653916478,
              0.063404432598,
              0.091047664068,
              0.084679571389,
              0.054747989396,
              0.098952817806,
          ]
          self.cov_params = [
              [
                  7.7376,
                  0.0080220,
                  -0.0039840,
                  0.0064925,
                  0.0022936,
                  -0.0098015,
              ],
              [
                  0.0080220,
                  0.0040201,
                  -0.0054843,
                  0.0046548,
                  -0.0029922,
                  -0.0047964,
              ],
              [
                  -0.0039840,
                  -0.0054843,
                  0.0082897,
                  -0.0072913,
                  0.0043566,
                  0.0067289,
              ],
              [
                  0.0064925,
                  0.0046548,
                  -0.0072913,
                  0.0071706,
                  -0.0043610,
                  -0.0057962,
              ],
              [
                  0.0022936,
                  -0.0029922,
                  0.0043566,
                  -0.0043610,
                  0.0029973,
                  0.0036193,
              ],
              [
                  -0.0098015,
                  -0.0047964,
                  0.0067289,
                  -0.0057962,
                  0.0036193,
                  0.0097917,
              ],
          ]
          self.hqic = 696.389304738311
          self.llf = -336.234154979014
          self.resid = residsc_mle[:, 2]
          self.fittedvalues = yhatc_mle[:, 2]
          self.pvalues = [
              0.7024,
              7.50e-34,
              0.0036,
              0.0009,
              1.77e-43,
              0.0193,
          ]
          self.tvalues = [0.3821, 12.13, -2.909, -3.306, 13.83, 2.340]
          self.sigma2 = 0.915487643192**2
          self.forecast = forecast_results["fc41c"]
          self.forecasterr = forecast_results["fe41c"]
      elif method == "css":
          # NOTE: params, bse, cov_params, tvalues taken from R;
          #   commented-out versions below are from
          #   [TODO: Finish this sentence]
          # self.params = [-0.077068926631, 0.763816531155, -0.270949972390,
          #                -0.284496499726, 0.757135838677, 0.225247299659]
          self.params = [
              -2.234160612756,
              0.763815335585,
              -0.270946894536,
              -0.284497190744,
              0.757136686518,
              0.225260672575,
          ]
          self.aic = 668.907200379791
          self.bic = 693.444521131318
          self.arroots = [
              1.0141 + 0.0000j,
              0.3036 - 0.9765j,
              0.3036 + 0.9765j,
              -1.2455 + 0.0000j,
          ]
          self.maroots = [-4.4396 + 0.0000j]
          # self.bse = [0.076048453921, 0.067854052128, 0.098041415680,
          #             0.090698349822, 0.057331126067, 0.099985455449]
          self.bse = [
              2.1842857865614,
              0.0644148863289,
              0.0923502391706,
              0.0860004491012,
              0.0558014467639,
              0.1003832271008,
          ]
          # self.cov_params = [
          #   [0.0057834, 0.00052477, -0.00079965,
          #    0.00061291, -0.00013618, -0.0018963],
          #   [0.00052477, 0.0046042, -0.0062505,
          #    0.0053416, -0.0032941, -0.0047957],
          #   [-0.00079965, -0.0062505, 0.0096121,
          #    -0.0084500, 0.0047967, 0.0064755],
          #   [0.00061291, 0.0053416, -0.0084500,
          #    0.0082262, -0.0048029, -0.0057908],
          #   [-0.00013618, -0.0032941, 0.0047967,
          #    -0.0048029, 0.0032869, 0.0035716],
          #   [-0.0018963, -0.0047957, 0.0064755,
          #    -0.0057908, 0.0035716, 0.0099971]
          # ]
          self.cov_params = [
              [
                  4.77110439737413,
                  -0.00908682223670,
                  0.00330914414276,
                  -0.00684678121434,
                  -0.00232348925409,
                  0.00950558295301,
              ],
              [
                  -0.00908682223670,
                  -0.00562941039954,
                  0.00852856667488,
                  -0.00749429397372,
                  -0.00304322809665,
                  -0.00494984519949,
              ],
              [
                  0.00330914414276,
                  -0.00562941039954,
                  0.00852856667488,
                  -0.00749429397372,
                  0.00443590637587,
                  0.00693146988144,
              ],
              [
                  -0.00684678121434,
                  0.00482359594764,
                  -0.00749429397372,
                  0.00739607724561,
                  -0.00448059420947,
                  -0.00600908311031,
              ],
              [
                  -0.00232348925409,
                  -0.00304322809665,
                  0.00443590637587,
                  -0.00448059420947,
                  0.00311380146095,
                  0.00373734623817,
              ],
              [
                  0.00950558295301,
                  -0.00494984519949,
                  0.00693146988144,
                  -0.00600908311031,
                  0.00373734623817,
                  0.01007679228317,
              ],
          ]
          self.hqic = 678.787238280001
          self.llf = -327.453600189896
          self.resid = residsc_css[4:, 2]
          self.fittedvalues = yhatc_css[4:, 2]
          self.pvalues = [
              0.3109,
              2.15e-29,
              0.0057,
              0.0017,
              8.06e-40,
              0.0243,
          ]
          # self.tvalues = [-1.013, 11.26, -2.764, -3.137, 13.21, 2.253]
          self.tvalues = [
              -1.02283347101,
              11.85774561000,
              -2.93390571556,
              -3.30808959392,
              13.56840602577,
              2.24400708246,
          ]
          self.sigma2 = 0.915919923456**2
  
imports: import os
  from numpy import genfromtxt
label_desc: Use image node_id 1 to create a new Y_arma41c object  with example: obj = Y_arma41c(method)