code: def test_predict_intervals():
      """
      Tests prediction intervals are ordered correctly
      """
      model = ARIMAX(formula="y ~ x1", data=data, ar=2, ma=2)
      x = model.fit()
      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
  from pyflux.arma import ARIMAX
label_desc: Use image node_id 10 for calling a global function with example usage: test_predict_intervals() without return types