code: def test_month_of_quarter(df_panel):
      """Test month_of_quarter for correctness, failure case of bug #4541."""
      y = load_PBS_dataset()
  
      FEATURES = ["month_of_quarter"]
      t = DateTimeFeatures(manual_selection=FEATURES)
  
      yt = t.fit_transform(y)[:13]
      expected = np.array([1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1])
      assert (expected == yt.values).all()
  
imports: import numpy
  import pandas
  import pytest
  from pandas.testing import assert_frame_equal
  from sktime.datasets import load_airline, load_longley, load_PBS_dataset
  from sktime.split import temporal_train_test_split
  from sktime.tests.test_switch import run_test_for_class
  from sktime.transformations.series.date import DateTimeFeatures
  from sktime.utils._testing.hierarchical import _make_hierarchical
label_desc: Use image node_id 27 for calling a global function with example usage: test_month_of_quarter(df_panel) without return types