code: def __init__(self, method="mle"):
      if method == "mle":
          self.params = [
              4.773779823083,
              0.591149657917,
              0.322267595204,
              -0.702933089342,
              0.116129490967,
              0.323009574097,
          ]
          self.aic = 720.814886758937
          self.bic = 745.465113183973
          self.arroots = [1.6916 + 0.0000j]
          # TODO: had to change order in maroots?
          self.maroots = [
              1.1071 - 0.7821j,
              1.1071 + 0.7821j,
              -1.2868 - 0.1705j,
              -1.2868 + 0.1705j,
          ]
          self.bse = [
              0.160891073193,
              0.151756542096,
              0.152996852330,
              0.140231020145,
              0.064663675882,
              0.065045468010,
          ]
          self.cov_params = [
              [
                  0.025886,
                  0.00026606,
                  -0.00020969,
                  -0.00021435,
                  4.2558e-05,
                  5.2904e-05,
              ],
              [
                  0.00026606,
                  0.023030,
                  -0.021269,
                  -0.018787,
                  0.0015423,
                  0.0011363,
              ],
              [
                  -0.00020969,
                  -0.021269,
                  0.023408,
                  0.018469,
                  -0.0035048,
                  -0.0010750,
              ],
              [
                  -0.00021435,
                  -0.018787,
                  0.018469,
                  0.019665,
                  -0.00085717,
                  -0.0033840,
              ],
              [
                  4.2558e-05,
                  0.0015423,
                  -0.0035048,
                  -0.00085717,
                  0.0041814,
                  0.0014543,
              ],
              [
                  5.2904e-05,
                  0.0011363,
                  -0.0010750,
                  -0.0033840,
                  0.0014543,
                  0.0042309,
              ],
          ]
          self.hqic = 730.735881539221
          self.llf = -353.407443379469
          self.resid = residsc_mle[:, 1]
          self.fittedvalues = yhatc_mle[:, 1]
          self.pvalues = [
              1.82e-193,
              9.80e-05,
              0.0352,
              5.37e-07,
              0.0725,
              6.84e-07,
          ]
          self.tvalues = [29.67, 3.895, 2.106, -5.013, 1.796, 4.966]
          self.sigma2 = 0.990262659233**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 = [1.502401748545, 0.683090744792, 0.197636417391,
          #                -0.763847295045, 0.137000823589, 0.304781097398]
          self.params = [
              4.740785760452,
              0.683056278882,
              0.197681128402,
              -0.763804443884,
              0.136991271488,
              0.304776424257,
          ]
          self.aic = 719.977407193363
          self.bic = 744.599577468616
          self.arroots = [1.4639 + 0.0000j]
          self.maroots = [
              1.1306 - 0.7071j,
              1.1306 + 0.7071j,
              -1.3554 - 0.0896j,
              -1.3554 + 0.0896j,
          ]
          # self.bse = [0.534723749868, 0.111273280223, 0.119840296133,
          #             0.111263606843, 0.070759105676, 0.061783181500]
          self.bse = [
              0.1750455599911,
              0.0942341854820,
              0.0999988749541,
              0.0929630759694,
              0.0628352649371,
              0.0645444272345,
          ]
          # self.cov_params = [
          #   [0.28593, -0.059175, 0.053968,
          #    0.046974, 0.00085168, 0.0028000],
          #   [-0.059175, 0.012382, -0.011333,
          #    -0.0098375, -0.00012631, -0.00058518],
          #   [0.053968, -0.011333, 0.014362,
          #    0.010298, -0.0028117, -0.00011132],
          #   [0.046974, -0.0098375, 0.010298,
          #    0.012380, 0.00031018, -0.0021617],
          #   [0.00085168, -0.00012631, -0.0028117,
          #    0.00031018, 0.0050069, 0.00079958],
          #   [.0028000, -0.00058518, -0.00011132,
          #    -0.0021617, 0.00079958, 0.0038172]]
          self.cov_params = [
              [
                  0.030640948072601,
                  -1.61599091345e-03,
                  0.001707084515950,
                  0.001163372764659,
                  -1.78587340563e-04,
                  0.000116062673743,
              ],
              [
                  -0.001615990913449,
                  8.88008171345e-03,
                  -0.007454252059003,
                  -0.006468410832237,
                  5.66645379098e-05,
                  -0.000381880917361,
              ],
              [
                  0.001707084515950,
                  -7.45425205900e-03,
                  0.009999774992092,
                  0.005860013051220,
                  -2.27726197200e-03,
                  0.000757683049669,
              ],
              [
                  0.001163372764659,
                  -6.46841083224e-03,
                  0.005860013051220,
                  0.008642133493695,
                  4.40550745987e-04,
                  -0.002170706208320,
              ],
              [
                  -0.000178587340563,
                  5.66645379098e-05,
                  -0.002277261972002,
                  0.000440550745987,
                  3.94827051971e-03,
                  0.000884171120090,
              ],
              [
                  0.000116062673743,
                  -3.81880917361e-04,
                  0.000757683049669,
                  -0.002170706208320,
                  8.84171120090e-04,
                  0.004165983087027,
              ],
          ]
          self.hqic = 729.888235701317
          self.llf = -352.988703596681
          self.resid = residsc_css[1:, 1]
          self.fittedvalues = yhatc_css[1:, 1]
          self.pvalues = [
              0.0050,
              8.31e-10,
              0.0991,
              6.64e-12,
              0.0528,
              8.09e-7,
          ]
          # self.tvalues = [2.810, 6.139, 1.649, -6.865, 1.936, 4.933]
          self.tvalues = [
              27.08315344127,
              7.24849772286,
              1.97683352430,
              -8.21621311385,
              2.18016541548,
              4.72196341831,
          ]
          self.sigma2 = 0.998687642867**2
  
imports: import os
  from numpy import genfromtxt
label_desc: Use image node_id 1 to create a new Y_arma14c object  with example: obj = Y_arma14c(method)