code: def test_2_predict_is_intervals_bbvi():
      """
      Tests prediction intervals are ordered correctly
      """
      model = ARIMAX(formula="y ~ x1 + x2", data=data, ar=2, ma=2)
      x = model.fit("BBVI", iterations=100, quiet_progress=True)
      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[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 29 for calling a global function with example usage: test_2_predict_is_intervals_bbvi() without return types