Несколько сезонных участков разложения

Question

У меня есть следующий кадр данных:

print(df)

  dt_op          Product     Quantity     
2018-12-15         3            10
2018-12-15         2            8
2018-12-16         3            5
(...)

Я хотел бы построить график сезонной декомпозиции с несколькими продуктами, представленными в виде линий разных цветов, используя:

from pylab import rcParams
import statsmodels.api as sm
import matplotlib.pyplot as plt

y  = df.set_index(["dt_op"])['Quantity'].resample('W').sum()
rcParams['figure.figsize'] = 11, 9

decomposition = sm.tsa.seasonal_decompose(y, model='additive')
fig = decomposition.plot()
plt.show()

Я пытался с:

codes = np.unique(df.cod_id)
for i in codes:

    y  = df.loc[(df.cod_id == i)].set_index(["dt_op"])['Quantity'].resample('W').sum()
    print(i)
    rcParams['figure.figsize'] = 11, 9

    decomposition = sm.tsa.seasonal_decompose(y, model='additive')
    fig = decomposition.plot()

plt.show()

Но дает не один график , а n графиков для n продуктов. Как я могу создать одну диаграмму с несколькими строками и по общей оси X? 1012 *