WWW.DISSERS.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА

   Добро пожаловать!

Pages:     | 1 |   ...   | 58 | 59 || 61 | 62 |   ...   | 82 |

t=t|t-1 N(0, t ), t = + 12 +... + q2. (16.1) t-1 t-q Здесь t-1 =(t-1, t-2,... ) — предыстория процесса {t}, а t —условная по предыстории дисперсия t, т.е. t = var(t|t-1) =E(2|t-1). Условную t дисперсию часто называют волатильностью процесса. Для того чтобы условная дисперсия оставалась положительной, требуется выполнение соотношений >и 1,..., q 0.

Данный процесс можно записать несколько иначе:

t NID(0, 1), t = tt, t = + 12 +... + q2.

t-1 t-q Аббревиатура NID означает, что t нормально распределены и независимы. Такая запись удобна тем, что нормированный случайный процесс t не зависит от предыстории.

Смысл модели ARCH состоит в том, что если абсолютная величина t оказывается большой, то это приводит к повышению условной дисперсии в последующие 16.1. Модель ARCH периоды. В свою очередь, при высокой условной дисперсии более вероятно появление больших (по абсолютной величине) значений t. Наоборот, если значения t в течение нескольких периодов близки к 0, то это приводит к понижению условной дисперсии в последующие периоды практически до уровня. В свою очередь, при низкой условной дисперсии более вероятно появление малых (по абсолютной величине) значений t. Таким образом, ARCH-процесс характеризуется инерционностью условной дисперсии (кластеризацией волатильности).

Несложно показать, что процесс ARCH не автокоррелирован:

E(tt-j) =E (E(tt-j|t-1)) = E (t-jE(t|t-1)) = 0.

Поскольку процесс имеет постоянное (нулевое) математическое ожидание и не автокоррелирован, то он является слабо стационарным в случае, если у него есть дисперсия.

Если обозначить разницу между величиной 2 и ее условным математическим t ожиданием, t, через t, то получится следующая эквивалентная запись процесса ARCH:

2 = + 12 +... + q2 + t. (16.2) t t-1 t-q Поскольку условное математическое ожидание t равно 0, то безусловное математическое ожидание также равно 0. Кроме того, как можно показать, {t} не автокоррелирован. Следовательно, квадраты процесса ARCH(q) следуют авторегрессионному процессу q-го порядка.

Если все корни характеристического уравнения 1 - 1z -... - qzq =лежат за пределами единичного круга, то у процесса ARCH(q) существует безусловная дисперсия, и он является слабо стационарным. Поскольку коэффициенты j q неотрицательны, то это условие эквивалентно условию j < 1.

j=Действительно, вычислим безусловную дисперсию стационарного ARCHпроцесса, которую мы обозначим через 2. Для этого возьмем математическое ожидание от обеих частей уравнения условной дисперсии (16.1):

E(t ) = + 1E(2 ) +... + qE(2 ).

t-1 t-q Заметим, что E t = E E(2|t-1) = E 2 = var(2) = 2, т.е. матеt t t матическое ожидание условной дисперсии равно безусловной дисперсии. Следовательно, 2 = + 12 +... + q2, 526 Модели с авторегрессионной условной...

0.ARCH(1),1 = 0.7, = 0. 0. 0.N(0; / (1 1)) 0. 0. 8 6 4 2 0 2 4 6 Рис. 16.2. Плотность ARCH(1) и плотность нормального распределения с той же дисперсией или 2 =.

1 - 1 -... - q Таким образом, для всех t безусловная дисперсия одинакова, т.е. имеет место гомоскедастичность. Однако условная дисперсия меняется, поэтому одновременно имеет место условная гетероскедастичность2.

Если не все корни приведенного выше характеристического уравнения лежат q за пределами единичного круга, т.е. если j 1, то безусловная дисперсия не j=существует, и поэтому ARCH-процесс не будет слабо стационарным3.

Еще одно свойство ARCH-процессов состоит в том, что безусловное распределение t имеет более высокий куртозис (т.е. более толстые хвосты и острую вершину), чем нормальное распределение (определение куртозиса и эксцесса см. в Приложении A.3.1). У ARCH(1) эксцесс равен E(4) t - 3 =, 4 1 - Она называется авторегрессионной, поскольку динамика квадратов ARCH-процесса описывается авторегрессией.

При этом у ARCH-процессов есть интересная особенность: они могут быть строго стационарны, не будучи слабо стационарны. Дело в том, что определение слабой стационарности требует существования конечных первых и вторых моментов ряда. Строгая же стационарность этого не требует, поэтому даже если условная дисперсия бесконечна (и, следовательно, ряд не является слабо стационарным), ряд всежеможет бытьстрогостационарным.

16.2. Модель GARCH причем при 31 1 четвертый момент распределения не существует (эксцесс равен бесконечности). Это свойство ARCH-процессов хорошо соответствует финансовым временным рядам, которые обычно характеризуются толстыми хвостами. На рисунке 16.2 изображен график плотности безусловного распределения ARCH(1). Для сравнения на графике приведена плотность нормального распределения с той же дисперсией.

Получить состоятельные оценки коэффициентов ARCH-процесса можно, используя вышеприведенное представление его квадратов в виде авторегрессии (16.2). Более эффективные оценки получаются при использовании метода максимального правдоподобия.

При применении ARCH-моделей к реальным данным было замечено, что модель ARCH(1) не дает достаточно длительных кластеров волатильности, а только порождает большое число выбросов (выделяющихся наблюдений). Для корректного описания данных требуется довольно большая длина лага q, что создает трудности при оценивании. В частности, зачастую нарушается условие неотрицательности оценок коэффициентов j. Поэтому Энгл предложил использовать модель со следующими ограничениями на коэффициенты лага: они задаются с помощью весов вида:

q +1- j wj =, 0.5q(q +1) сумма которых равна 1 и которые линейно убывают до нуля. Сами коэффициенты берутся равными j = wj. Получается модель с двумя параметрами, и :

t = + w12 +... + wq2.

t-1 t-q 16.2. Модель GARCH Модель GARCH (generalized ARCH — обобщенная модель ARCH), предложенная Т. Боллерслевом, является альтернативной модификацией модели ARCH (16.2), позволяющей получить более длинные кластеры при малом числе параметров. Модель ARMA зачастую позволяет получить более сжатое описание временных зависимостей для условного математического ожидания, чем модель AR.

Подобным же образом модель GARCH дает возможность обойтись меньшим количеством параметров по сравнению с моделью ARCH, если речь идет об условной дисперсии. В дальнейшем мы проведем прямую аналогию между моделями GARCH и ARMA.

528 Модели с авторегрессионной условной...

Для того чтобы вывести модель GARCH, используем в модели ARCH бесконечный геометрический лаг:

t = + j-12 = + 2.

t-j t-1 - L j=Применяя преобразование Койка, получим 2 t =(1 - ) + t-1 + 2.

t-Поменяв очевидным образом обозначения, получим модель GARCH(1, 1):

2 t = + t-1 + 2.

t-Модель GARCH(p, q) обобщает эту формулу:

2 2 t = + 1t-1 +... + pt-p + 12 +... + q2 = t-1 t-q p q = + jt-j + j2. (16.3) t-j j=1 j=При этом предполагается, что > 0, 1,..., p 0 и 1,..., q 0.

На практике, как правило, достаточно взять p =1 и q =1. Изредка используют GARCH(1, 2) или GARCH(2, 1).

Как и в модели ARCH, t служит условной дисперсией процесса:

t|t N(0, t ).

Рассчитаем безусловную дисперсию GARCH-процесса, предполагая, что он стационарен. Для этого возьмем математические ожидания от обеих частей уравнения 16.3 для условной дисперсии:

p q 2 E(t ) = jE(t-j) + jE(2 ), t-j j=1 j=откуда p q 2 = j2 + jj=1 j=и 2 =.

p q 1 - j - j j=1 j=16.2. Модель GARCH Таким образом, с точки зрения безусловной дисперсии GARCH-процесс гомоскедастичен.

Для того чтобы дисперсия была конечной, необходимо выполнение условия p q j + j < 1. В частности, для модели GARCH(1, 1) требуется 1 + 1 < 1.

j=1 j=Процесс GARCH можно записать в эквивалентной форме, если, как и выше, в уравнении 16.2 для модели ARCH, обозначить t = 2 - t :

t p m 2 = + (j + j)2 + t - jt-j, t t-j j=1 j=где m =max(p, q). (В этой записи подразумевается j =0 при j > p и j =при j >q.) Такая форма записи позволяет увидеть, что квадраты GARCH-процесса подчиняются модели ARMA(m, p).

Этот факт дает возможность получить автокорреляционную функцию квадратов GARCH-процесса. В частности, для GARCH(1, 1) автокорреляционная функция квадратов имеет вид 1(1 - 1 - 11) 1 =, 1 - 1 - =(1 + 1) -1-1, > 1.

Условие существования безусловного четвертого момента у отдельного наблю2 дения процесса GARCH(1, 1) состоит в том, что 31 +211 + 1 < 1. Если это условие выполняется, то эксцесс равен E(4) t - 3 = 2 4 1 - 1 - 211 - и является положительным. То есть GARCH-процесс (как и его частный случай — ARCH-процесс) имеет более высокий куртозис, чем нормальное распределение. В то же время, безусловное распределение отдельного наблюдения GARCHпроцесса является симметричным, поэтому все нечетные моменты, начиная с третьего, равны нулю.

Стандартным методом оценивания для моделей GARCH является метод максимального правдоподобия. Условно по предыстории t-1 отдельное наблюдение GARCH-процесса распределено нормально: t|t-1 N(0, t ). Функция правдоподобия для ряда 1,..., T, подчиняющегося GARCH-процессу, вычисляется как произведение плотностей этих условных нормальных распределений:

T 1 t L = exp -.

2 2t t=1 2t 530 Модели с авторегрессионной условной...

Максимизируя эту функцию правдоподобия по неизвестным параметрам, получим оценки максимального правдоподобия для GARCH-процесса.

При оценивании условную дисперсию t следует считать функцией параметров модели и вычислять по рекуррентной формуле 16.3. Для этих вычислений требуются «довыборочные» значения самого процесса и его условной дисперсии, а они неизвестны. Для решения этой проблемы можно использовать различные приемы. Самый простой, по-видимому, состоит в том, чтобы заменить условные дисперсии в начале ряда ( t =1,..., m) оценкой безусловной дисперсии, т.е. величиной T s2 = 2.

t T t=Оценки максимального правдоподобия являются состоятельными и асимптотически эффективными.

На практике модель GARCH дополняют какой-либо моделью, описывающей поведение условного или безусловного математического ожидания наблюдаемого ряда. Например, можно предположить, что наблюдается не t, а t плюс константа, т.е. что наблюдаемый ряд {xt} имеет постоянное безусловное математическое ожидание, к которому добавляется ошибка t в виде процесса GARCH:

xt = + t.

Можно моделировать безусловное математическое ожидание с помощью линейной регрессии, т.е.

xt = Zt + t.

Это позволяет учитывать линейный тренд, детерминированные сезонные переменные и т.п. При оценивании в функции правдоподобия вместо t используют xt - Zt.

С точки зрения прогнозирования перспективной является модель, сочетающая ARIMA с GARCH. Модель ARIMA в этом случае используется для моделирования поведения условного математического ожидания ряда, а GARCH —д лямод елирования условной дисперсии.

Важнейший вывод, который следует из анализа модели ARCH, состоит в том, что наблюдаемые изменения в дисперсии (волатильности) временного ряда могут иметь эндогенный характер, то есть порождаться определенной нелинейной моделью, а не какими-то внешними структурными сдвигами.

16.3. Прогнозы и доверительные интервалы для модели GARCH 16.3. Прогнозы и доверительные интервалы для модели GARCH Одна из важнейших целей эконометрических моделей временных рядов — построение прогнозов. Какие преимущества дают модели с авторегрессионной условной гетероскедастичностью с точки зрения прогнозирования временных рядов по сравнению с моделями линейной регрессии или авторегрессии — скользящего среднего Оказывается, что прямых преимуществ нет, но есть ряд опосредованных преимуществ, которые в отдельных случаях могут иметь большое значение.

Рассмотрим модель линейной регрессии xt = Zt + t, t =1,..., T, в которой ошибка представляет собой GARCH-процесс. Поскольку ошибки не автокоррелированы и гомоскедастичны, то, как известно, оценки наименьших квадратов являются наилучшими в классе линейных по x несмещенных оценок. Однако наличие условной гетероскедастичности позволяет найти более эффективные (т.е. более точные) оценки среди нелинейных и смещенных оценок. Действительно, метод максимального правдоподобия дает асимптотически эффективные оценки, более точные, чем оценки МНК. В ошибку прогноза вносит свой вклад, во-первых, ошибка T +1, а во-вторых, разница между оценками параметров и истинными значениями параметров. Использование более точных оценок позволяет уменьшить в некоторой степени вторую составляющую ошибки прогноза.

В обычных моделях временного ряда с неизменными условными дисперсиями (например, ARMA) неопределенность ошибки прогноза — это возрастающая функция горизонта прогноза, которая (если не учитывать разницу между оценками параметров и истинными значениями параметров, отмеченную в предыдущем абзаце) не зависит от момента прогноза. Однако в присутствии ARCH-ошибок точность прогноза будет нетривиально зависеть от текущей информации и, следовательно, от момента прогноза. Поэтому для корректного построения интервальных прогнозов требуется иметь оценки будущих условных дисперсий ошибки.

Кроме того, в некоторых случаях полезно иметь прогнозы не только (условного) математического ожидания изучаемой переменной, но и ее (условной) дисперсии.

Это важно, например, при принятии решений об инвестициях в финансовые активы. В этом случае дисперсию (волатильность) доходности естественно рассматривать как меру рискованности финансового актива. Таким образом, сами по себе прогнозы условной дисперсии могут иметь практическое применение.

Покажем, что доверительный интервал прогноза зависит от предыстории T =(xT, xT -1,..., x1,... ).

532 Модели с авторегрессионной условной...

Реально прогноз делается на основе имеющегося ряда (x1,..., xT ), а не всей предыстории, однако различие это не столь существенно. При этом мы будем исходить из того, что нам известны истинные параметры процесса. Прогноз на периодов — это математическое ожидание прогнозируемой величины xT +, условное относительно имеющейся на момент t информации T. Онравен xT () =E(xT +|T ) =E(ZT + + T + |T ) =ZT +.

Здесь учитывается, что поскольку информация T содержится в информации T + -1 при 1, то по правилу повторного взятия математического ожидания выполнено E(T + |T ) =E (E (T + |T + -1) |T ) =E(0|T ) =0.

Таким образом, если известны истинные параметры, присутствие GARCH-ошибок не отражается на том, как строится точечный прогноз, — он оказывается таким же, как для обычной линейной регрессии. Ошибка предсказания в момент времени T на шагов вперед dT () =xT + - xT () =T +.

Условная дисперсия ошибки предсказания равна p = E d2 () T = E 2 T.

T T + Из этого следует, что она зависит как от горизонта прогноза, так и от предыстории T.

Заметим, что при t >T выполнено E 2|T = E t |T, поскольку t 2 E 2 - t T = E E 2 - t t-1 T = E(0|T ) =0.

Pages:     | 1 |   ...   | 58 | 59 || 61 | 62 |   ...   | 82 |



© 2011 www.dissers.ru - «Бесплатная электронная библиотека»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.