Слайд 1ARCH-модель
Определение 1: Условной дисперсией называется дисперсия случайной переменной, обусловленная информацией о других случайных переменных, то есть дисперсия, найденная при условии наличия знаний о дисперсии в предыдущие моменты времени σt2= D(εt |εt-1,εt-2…). Определение 2: ARCH-модель первого порядка имеет вид: , (1) где - остатки, полученные после предварительной оценки какой-либо модели. Здесь дисперсия в момент времени t зависит от квадрата ошибок в момент (t-1), то есть условная дисперсия σt2 является AR-процессом квадрата ошибок модели. ARCH-модель (q) (порядок авторегрессии дисперсии – q) имеет вид: (2) здесь условная дисперсия представлена как линейная функция квадратов прошлых ошибок в моменты времени t-1, t-2, …, t-q AR-процесс временных рядов – процесс авторегрессии, когда текущие значения ряда линейно зависят от предыдущих значений.
Слайд 2Моделировать дисперсию можно лишь после того, как из данных убрали среднюю величину, поэтому полная модель ARCH(q) имеет вид: (3) Здесь первое уравнение является AR-процессом р-ого порядка, а во втором уравнении моделируется условная дисперсия как линейная функция от q квадратов прошлых значений ошибок (εt-q), полученных после оценки первого уравнения. Условие «положительности» дисперсии: β0>0, β1≥0, β2≥0,…, βq≥0. Смысл ARCH-модели: если абсолютная величина остатков модели εt первого уравнения (3) является большой, то это приведет к повышению условной дисперсии в последующие периоды (второе уравнение (3)), наоборот если остатки близки к нулю, то это приведет к снижению условной дисперсии. Это проявляется свойство, называемое кластеризацией волатильности, которую как раз и позволяет измерить ARCH-модель. Эффект кластеризации волатильности впервые было отмечено Манделбротом (1963 г.)
Слайд 3Алгоритм определения наличия ARCH-эффектов.
необходимо построить AR-модель ряда xt с ошибкой εt согласно первому уравнению из (3); определить остатки как оценки εt ; построить линейную регрессию квадрата ошибок в момент времени t на квадрат остатков модели после моделирования АР: ; протестировать коэффициент λ на отсутствие значимости с помощью теста Стьюдента, теста Фишера, χ2 теста, принимая за нулевую гипотезу: H0: λ1=0. Соответственно за альтернативную гипотезу H1: λ1≠0. 5. Если λ1 значимо отличается от 0, то модель можно специфицировать как ARCH-модель первого порядка (ARCH (1)).
Слайд 4Общая схема тестирования модели на ARCH-эффекты:
Проводится оценка модели (например, AR-модели, СС-модели, АРСС-модели или простой регрессии по времени); на основе знания ошибок модели ( – расчётное значение модели, построенной в пункте 1)), оценивается модель: Здесь модель тестируется на ARCH-эффекты p-го порядка. 3. для оценённой модели рассчитывается коэффициент детерминации R2, отвечающий за качество подгонки модели; 4. формируются гипотезы (нулевая и альтернативная): , ; 5. определяется значение статистики χ2расч =TR2, где Т – объём выборки ряда, R2 – коэффициент детерминации; 6. сравнивается χ2расч с χ2табл, определённым для степеней свободы p (p – число временных лагов в модели ARCH(p)) 7. если χ2расч > χ2табл , то H0 отклоняется, и считается, что ARCH-модель существенна на заданном уровне значимости и её порядок равен p.
Слайд 5GARCH-модель
Определение 3: GARCH-модель – модель с обобщённой авторегрессией условной гетероскедастичности. GARCH (p,q) имеет в отличие от ARCH-модели два порядка и в общем виде записывается: (4) где αi и βj >0 (i=1, 2,…, p; j=1, 2,…, q) иначе дисперсия была бы меньше нуля. GARCH- модель показывает, что текущее значение условной дисперсии является функцией от константы – p-го значения квадратов остатков из уравнения условной средней (или любого другого уравнения) и q-го значения предыдущей условной дисперсии (то есть АR-процесс q-го порядка от условной дисперсии). Наиболее популярной для предсказания изменчивости доходности финансовых активов является модель GARCH(1, 1): (5) .
Слайд 6Волатильность GARCH
Волатильность (изменчивость) не является постоянным процессом, и может изменяться во времени. Если известна точная модель описания процесса, изменяющегося во времени, то для нахождения годовой волатильности этого процесса нужно определить квадратный корень из условной дисперсии и умножить модель на , где N – количество наблюдений в год. Полученная мера волатильности будет изменяться во времени, т.е. текущая волатильность будет определяться как функция от прошлой волатильности. Для предсказания волатильности при помощи GARCH-модели, можно использовать следующую рекурсивную модель: (6) (7) Здесь εt2 - величина, неизвестная в будущем, которая при выполнении прогноза заменяется условной оценкой дисперсии σt. Таким образом, формула (7) позволяет предсказывать σt2 в момент времени (t+1), затем σt2 в момент времени (t+2) и т.д. При этом, например, σt+2 рассчитывается как условная дисперсия при условии известности значений y1, y2,…,yt и прогноза yt+1. Результат каждого расчёта является предсказанием условной дисперсии на j периодов вперёд.
Слайд 7Оценка моделей ARCH и GARCH
ARCH и GARCH процессы имеют, как правило, островершинное безусловное распределение. Так эксцессы (момент четвертого порядка) для модели ARCH (1), представленной уравнением (1), и GARCH (1;1), представленной уравнением (5), соответственно равны и . Коэффициенты асимметрии (моменты третьего порядка) для моделей волатильности равны нулю. Несмотря на это стандартным методом оценки моделей является метод максимального правдоподобия, основывающийся на нормальном распределении. При этом оценки моделей будут состоятельными, но асимптотически неэффективными (неэффективные в пределе при увеличении числа степеней свободы). Отметим, что наличие высоких куртозисов ARCH-процессов хорошо согласуется с поведением многих финансовых показателей, имеющих толстые хвосты в распределении.
Слайд 8Рассмотрим условное математическое ожидание в момент времени t: , Отсюда ошибка определяется как: . Предположим, что , где - условная дисперсия, а z –распределена по стандартизированному нормальному закону, то есть z ~ N(0, 1). Тогда εt ~ N(0, ), где для ARCH-модели: ; а для GARCH-модели: + В итоге для ARCH-модели мы имеем m+1+p+1 неизвестных параметров, а для GARCH-модели m+1+p+1+q параметров, подлежащих оцениванию. Оценка параметров производится методом максимального правдоподобия.
Слайд 9Проверка адекватности GARCH/ARCH моделей.
Качество подгонки GARCH/ARCH модели под исходные данные можно контролировать на основе близости к единице индекса детерминации (R2) или скорректированного на число степеней свободы индекса детерминации (R2Adjusted). или , здесь n – общее число наблюдений временного ряда, k – число степеней свободы модели (для GARCH k=p+q, для ARCH k=p), – остаточная или объясненная моделью дисперсия, – общая дисперсия. Для проверки достоверности оценок модели необходимо проанализировать стандартизованные остатки έ/σ, где σ – условное среднеквадратическое отклонение, рассчитанное по модели GARCH/ARCH, а έ - остатки в уравнении условного математического ожидания (первоначального уравнения). Если модель GARCH/ARCH достаточно хорошо описана, то стандартизованные остатки являются независимыми одинаково распределенными случайными величинами с нулевым математическим ожиданием и единичным средним квадратичным отклонением.
Слайд 10Идентификация модели GARCH
I этап: Расчет критерия Люнга-Бокса (LB). Рассчитывается LB статистика на основе предварительного расчета k коэффициентов автокорреляции по Т наблюдениям (ρk), с последующим возведением их в квадрат : (8) где m – максимальный лаг коэффициентов автокорреляции, Т – длина временного ряда. Выдвигается гипотеза о незначимости m лагов в исходной модели авторегрессии. Сравнивают расчетное значение LB с критическим значением χ2, определенным для степени свободы v=т . Если LBрасчетное > χ2 , то гипотезу о незначимости m лагов в исходной модели авторегрессии отклоняют на заданном уровне значимости α. II этап: Расчет критерия Люнга-Бокса по стандартизованным остаткам. Рассчитывают коэффициенты автокорреляции по стандартизованным остаткам с последующим возведением их в квадрат : (9) где m – максимальный лаг автокорреляции стандартизованных остатков. Выдвигается гипотеза о несущественности порядков p и q GARCH модели. LBрассчетное сравнивают с χ2табл, определенным для степени свободы v1=m-p-q, где т – общее число наблюдений, p и q – порядки GARCH модели. Если LBрассчетное < χ (m-p-q), то нулевую гипотезу отклоняют и считают, что GARCH модель хорошо определена, и ее порядки соответственно равны p и q.
Слайд 11Идентификация модели GARCH на основе анализа коррелограмм
После оценивания математического ожидания ряда данных (на основе АРПСС моделей, выделения компонент временного ряда или обычной регрессии) получают остаточную компоненту. Стандартизируют полученные остатки. Строят коррелограммы АКФ и ЧАКФ на стандартизированных остатках. Определяют количество лагов для коэффициентов АКФ и ЧАКФ, выходящих за границы белого шума. Полученное число является порядком ARCH-модели. Селекцию моделей ARCH и GARCH проводят исходв из минимума информационных критериев Акайке, Шварца и Ханена-Квина.