Модели ARIMA: ошибка доли plot_diagnostics

Question

Я использую plot_diagnostics для построения модели ARIMA.

Мой код

mod = sm.tsa.statespace.SARIMAX(y,
                            order=(1, 1, 1),
                            seasonal_order=(1, 1, 0, 12),
                            enforce_stationarity=False,
                            enforce_invertibility=False)
results = mod.fit()
print(results.summary().tables[1])

Пока все хорошо.Но когда я пытаюсь использовать plot_diagonistic, я получаю ошибку.

results.plot_diagnostics(figsize=(30,15))
plt.show()

Ошибка:

ValueError: operands could not be broadcast together with shapes (9,) (3,) (9,) 

Не уверен, почему это произошло.Спасибо