code: def test2_predict_is_intervals_mh():
      """
      Tests prediction intervals are ordered correctly
      """
      model = pf.GASX(
          formula="y ~ x1 + x2",
          data=data,
          ar=1,
          sc=1,
          family=pf.Poisson(),
      )
      x = model.fit("M-H", nsims=400)
      predictions = model.predict_is(h=10, intervals=True)
      assert np.all(
          predictions["99% Prediction Interval"].values
          > predictions["95% Prediction Interval"].values
      )
      assert np.all(
          predictions["95% Prediction Interval"].values
          > predictions["5% Prediction Interval"].values
      )
      assert np.all(
          predictions["5% Prediction Interval"].values
          > predictions["1% Prediction Interval"].values
      )
  
imports: import numpy
  import pandas
  import pyflux
label_desc: Use image node_id 44 for calling a global function with example usage: test2_predict_is_intervals_mh() without return types