Прогнозування

Моделі ARIMA в Python

James Fulton

Climate informatics researcher

Прогноз наступного значення

Візьмімо модель AR(1)

$y_t = a_1 y_{t-1} + \epsilon_t$

Передбачте наступне значення

$y_t = 0.6\; \text{x} \; 10 + \epsilon_t$

$y_t = 6.0 + \epsilon_t$

Невизначеність прогнозу

$ 5.0 < y_t < 7.0$

Моделі ARIMA в Python

Однокрокові прогнози

Моделі ARIMA в Python

Створення однокрокових прогнозів

# Make predictions for last 25 values
results = model.fit()

# Make in-sample prediction
forecast = results.get_prediction(start=-25)
Моделі ARIMA в Python

Створення однокрокових прогнозів

# Make predictions for last 25 values
results = model.fit()

# Make in-sample prediction
forecast = results.get_prediction(start=-25)

# forecast mean
mean_forecast = forecast.predicted_mean

Середнє прогнозу — це серія pandas

2013-10-28    1.519368
2013-10-29    1.351082
2013-10-30    1.218016
Моделі ARIMA в Python

Довірчі інтервали

# Get confidence intervals of forecasts
confidence_intervals = forecast.conf_int()

Метод довірчого інтервалу повертає DataFrame pandas

                                 lower y                       upper y
2013-09-28                     -4.720471                     -0.815384
2013-09-29                     -5.069875                      0.112505
2013-09-30                     -5.232837                      0.766300
2013-10-01                     -5.305814                      1.282935
2013-10-02                     -5.326956                      1.703974
Моделі ARIMA в Python

Візуалізація прогнозів

plt.figure()

# Plot prediction
plt.plot(dates, 
         mean_forecast.values, 
         color='red', 
         label='forecast')

# Shade uncertainty area
plt.fill_between(dates, lower_limits, upper_limits, color='pink')

plt.show()
Моделі ARIMA в Python

Візуалізація прогнозів

Моделі ARIMA в Python

Динамічні прогнози

Моделі ARIMA в Python

Створення динамічних прогнозів

results = model.fit()
forecast = results.get_prediction(start=-25, dynamic=True)

# forecast mean
mean_forecast = forecast.predicted_mean

# Get confidence intervals of forecasts
confidence_intervals = forecast.conf_int()
Моделі ARIMA в Python

Прогноз поза вибіркою

forecast = results.get_forecast(steps=20)

# forecast mean
mean_forecast = forecast.predicted_mean

# Get confidence intervals of forecasts
confidence_intervals = forecast.conf_int()
Моделі ARIMA в Python

Прогноз поза вибіркою

forecast = results.get_forecast(steps=20)

Моделі ARIMA в Python

Давайте потренуємось!

Моделі ARIMA в Python

Preparing Video For Download...