Вывод условных отклонений и невязок

Этот пример показывает, как вывести условные отклонения из подобранной модели условного отклонения. Стандартизированные невязки вычисляются с помощью предполагаемых условных отклонений, чтобы проверить подгонку модели.

Шаг 1. Загрузите данные.

Загрузите датские номинальные данные возврата запаса, включенные в тулбокс.

load Data_Danish
y = DataTable.RN;
T = length(y);

figure
plot(y)
xlim([0,T])
title('Danish Nominal Stock Returns')

Figure contains an axes. The axes with title Danish Nominal Stock Returns contains an object of type line.

Серия возврата, по-видимому, имеет ненулевое среднее смещение и волатильность кластеризации.

Шаг 2. Подбор модели EGARCH (1,1

).

Задайте, а затем подгоните модель EGARCH (1,1) к номинальной серии возвратов запаса. Включите среднее смещение и примите Гауссово инновационное распределение .

Mdl = egarch('Offset',NaN','GARCHLags',1,...
    'ARCHLags',1,'LeverageLags',1);
EstMdl = estimate(Mdl,y);
 
    EGARCH(1,1) Conditional Variance Model with Offset (Gaussian Distribution):
 
                     Value       StandardError    TStatistic     PValue  
                   __________    _____________    __________    _________

    Constant         -0.62723       0.74401        -0.84304       0.39921
    GARCH{1}          0.77419       0.23628          3.2766     0.0010507
    ARCH{1}           0.38636       0.37361          1.0341       0.30107
    Leverage{1}    -0.0024989       0.19222          -0.013       0.98963
    Offset            0.10325      0.037727          2.7368     0.0062047

Шаг 3. Вывод условных отклонений.

Вывод условных отклонений с помощью подобранной модели.

v = infer(EstMdl,y);

figure
plot(v)
xlim([0,T])
title('Inferred Conditional Variances')

Figure contains an axes. The axes with title Inferred Conditional Variances contains an object of type line.

Конечные условные отклонения показа повышенную волатильность в конце возврата ряда.

Шаг 4. Вычислите стандартизированные невязки.

Вычислите стандартизированные невязки для модели подгонки. Вычесть предполагаемое среднее смещение и разделить на квадратный корень процесса условного отклонения.

res = (y-EstMdl.Offset)./sqrt(v);

figure
subplot(2,2,1)
plot(res)
xlim([0,T])
title('Standardized Residuals')

subplot(2,2,2)
histogram(res,10)

subplot(2,2,3)
autocorr(res)

subplot(2,2,4)
parcorr(res)

Figure contains 4 axes. Axes 1 with title Standardized Residuals contains an object of type line. Axes 2 contains an object of type histogram. Axes 3 with title Sample Autocorrelation Function contains 4 objects of type stem, line. Axes 4 with title Sample Partial Autocorrelation Function contains 4 objects of type stem, line.

Стандартизированные невязки не показывают остаточной автокорреляции. Существуют несколько невязок больше, чем ожидалось для Гауссова распределения, но предположение нормальности не неразумно.

См. также

Объекты

Функции

Похожие темы