code: def a_test_predict_intervals_mh():
      model = pf.EGARCHMReg(formula="JPM~Rf", data=data, p=1, q=1)
      x = model.fit("M-H", nsims=400)
      predictions = model.predict(h=10, intervals=True, oos_data=data)
      assert np.all(
          predictions["99% Prediction Interval"].values + 0.000001
          >= predictions["95% Prediction Interval"].values
      )
      assert np.all(
          predictions["95% Prediction Interval"].values + 0.000001
          >= predictions["5% Prediction Interval"].values
      )
      assert np.all(
          predictions["5% Prediction Interval"].values + 0.000001
          >= predictions["1% Prediction Interval"].values
      )
  
imports: import numpy
  import pyflux
  import pandas
  from pandas.io.data import DataReader
  from datetime import datetime
label_desc: Use image node_id 20 for calling a global function with example usage: a_test_predict_intervals_mh() without return types