code: def check_panel_theta_quantiles():
      """Test predict quantiles with theta on panel data with datetime index."""
      # make panel with hour of day panel and datetime index
      y = load_airline()
      y.index = pd.date_range(
          start="1960-01-01", periods=len(y.index), freq="H"
      )
      y.index.names = ["datetime"]
      y.name = "passengers"
      y = y.to_frame()
      y["hour_of_day"] = y.index.hour
      y = (
          y.reset_index()
          .set_index(["hour_of_day", "datetime"])
          .sort_index()
      )
  
      forecaster = ThetaForecaster(sp=1)
      forecaster.fit(y)
      forecaster.predict(fh=[1, 3])
      forecaster.predict_quantiles(fh=[1, 3], alpha=[0.1, 0.5, 0.9])
  
file_path: sktime/sktime/forecasting/tests/test_theta.py
imports: import numpy
  import pandas
  import pytest
  from sktime.datasets import load_airline
  from sktime.forecasting.tests._config import TEST_OOS_FHS
  from sktime.forecasting.theta import ThetaForecaster, ThetaModularForecaster
  from sktime.split import temporal_train_test_split
  from sktime.utils.validation._dependencies import _check_soft_dependencies
  from sktime.utils.validation.forecasting import check_fh
label_desc: Use image node_id 5 for calling a global function with example usage: check_panel_theta_quantiles() without return types