Как бы я рассчитал стандартизированные невязки из функции sarimax
модели arima?
Допустим, у нас есть некоторая базовая модель:
import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
import seaborn as sns
sns.set(style='ticks', context='poster')
from statsmodels.tsa.statespace.sarimax import SARIMAX
from statsmodels.tsa.seasonal import seasonal_decompose
import seaborn as sns
#plt.style.use("ggplot")
import pandas_datareader.data as web
import pandas as pd
import statsmodels.api as sm
import scipy
import statsmodels.stats.api as sms
import matplotlib.pyplot as plt
import datetime
model = SARIMAX(df, order = (6, 0, 0), trend = "c");
model_results = model.fit(maxiter = 500);
print(model_results.summary());
Мне нужен стандартизатор, поэтому, когда мы используем model_results.plot_diagnostics(figsize = (16, 10));
функция, а затем только основные plot
остатки функции должны выглядеть одинаково.