我可以使用 SARIMA 模型拟合一些数据pmdarima
。
import pmdarima as pm
from pmdarima.model_selection import train_test_split
import numpy as np
import matplotlib.pyplot as plt
# Load/split
y = pm.datasets.load_wineind()
train, test = train_test_split(y, train_size=150)
# Fit
model = pm.auto_arima(train, seasonal=True, m=12)
我可以根据这些数据进行预测,甚至可以查看样本内预测,从中可以计算残差。
N = test.shape[0] # predict N steps into the future
forecasts = model.predict(N)
in_sample_forecasts = model.predict_in_sample()
但 SARIMA 只是一个数学模型(据我所知)。所以我希望能够使用拟合的模型参数来完全预测其他一些系列。我可以这样做吗?
例如:
# Some other series entirely
some_other_series = train + np.random.randint(0, 5000, len(train))
# The following method does not exist but illustrates the desired functionality
forecasts = model.predict_for(some_other_series, N)