code: def test_predict_intervals_mh():
      """
      Tests prediction intervals are ordered correctly
      """
      model = pf.GASX(
          formula="y ~ x1", data=data, ar=1, sc=1, family=pf.t()
      )
      x = model.fit("M-H", nsims=400)
      predictions = model.predict(
          h=10, oos_data=data_oos, intervals=True
      )
  
      assert np.all(
          predictions["99% Prediction Interval"].values
          > predictions["95% Prediction Interval"].values
      )
      assert np.all(
          predictions["95% Prediction Interval"].values
          > predictions[model.data_name].values
      )
      assert np.all(
          predictions[model.data_name].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 21 for calling a global function with example usage: test_predict_intervals_mh() without return types