Этот пример показывает, как вывести условные отклонения из подобранной модели условного отклонения. Стандартизированные невязки вычисляются с помощью предполагаемых условных отклонений, чтобы проверить подгонку модели.
Загрузите датские номинальные данные возврата запаса, включенные в тулбокс.
load Data_Danish y = DataTable.RN; T = length(y); figure plot(y) xlim([0,T]) title('Danish Nominal Stock Returns')
Серия возврата, по-видимому, имеет ненулевое среднее смещение и волатильность кластеризации.
Задайте, а затем подгоните модель 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
Вывод условных отклонений с помощью подобранной модели.
v = infer(EstMdl,y);
figure
plot(v)
xlim([0,T])
title('Inferred Conditional Variances')
Конечные условные отклонения показа повышенную волатильность в конце возврата ряда.
Вычислите стандартизированные невязки для модели подгонки. Вычесть предполагаемое среднее смещение и разделить на квадратный корень процесса условного отклонения.
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)
Стандартизированные невязки не показывают остаточной автокорреляции. Существуют несколько невязок больше, чем ожидалось для Гауссова распределения, но предположение нормальности не неразумно.