Книги по разным темам Pages:     | 1 |   ...   | 3 | 4 | 5 | 6 | 7 |   ...   | 35 |

4 2 0 -2 --4 -10 20 30 40 50 60 70 80 90 100 10 20 30 40 50 60 70 80 90 ARMA((1, 4), 1) ARMA(1, (1,4)) Зметим, что для первой модели уравнение a(z) = 0 принимает вид 1 - 2/3z + 1/z4 = 0, т.е. z4 - 32z + 48 = 0; корни этого уравнения z1 = 2, z2 = 2, z3 = - 2 + i8, z4 = - - i8 лежат вне единичного круга, что обеспечивает стационарность рассматриваемого процесса. Во второй модели a(z) = 0 принимает вид 1 - 0.4 z = 0;

корень этого уравнения z = 2.5 > 1, так что и эта модель стационарна.

Кроме рассмотренных примеров аддитивных сезонных моделей, употребляются также и мультипликативные спецификации, например, (1 - a1L)Xt = (1+ b1L)(1+ b4L4) t, (1 - a1L)(1 - a4L4) Xt = (1+ b1L) t.

Первая дает Xt = a1 XtЦ1 + t + b1 tЦ1+ b4 tЦ4 + b1b4 tЦ5, а вторая Xt = a1 XtЦ1 + a4 XtЦ4 - a1a4 XtЦ5 + t + b1 tЦ1.

В первой модели допускается взаимодействие составляющих скользящего среднего на лагах 1 и 4 (т.е. значений tЦ1 и tЦ4), а во второй - взаимодействие авторегрессионных составляющих на лагах 1 и 4 (т.е. значений XtЦ1 и XtЦ4 ). Конечно, эти две модели являются частными случаями аддитивных моделей Xt = a1 XtЦ1 + t + b1 tЦ1+ b4 tЦ4 + b5 tЦ5, Xt = a1 XtЦ1 + a4 XtЦ4 + a5 XtЦ5 + t + b1 tЦc b5 = b1b4, a5 = - a1a4. При приближенном выполнении последних соотношений (по крайней мере, если гипотезы о наличии таких соотношений не отвергаются), естественно перейти от оценивания аддитивной модели к оцениванию мультипликативной модели, опять следуя принципу УэкономностиФ модели (Уparsimony modelФ). Впрочем, каких-либо теоретических оснований, ведущих к предпочтению одной формы сезонности перед другой (мультипликативной или аддитивной), не существует.

Более подробно с сезонными ARMA моделями можно ознакомиться, например, в книге [Enders (1995)].

www.iet.ru/mipt/2/text/curs_econometrics.htm Эконометрика. Введение в регрессионный анализ временных рядов. В.П.Носко www.iet.ru www.iet.ru/mipt/2/text/curs_econometrics.htm Эконометрика. Введение в регрессионный анализ временных рядов. В.П.Носко www.iet.ru Глава 3. Подбор стационарной модели ARMA для ряда наблюдений Если мы предполагаем, что некоторый наблюдаемый временной ряд x1, x2, Е, xT порождается моделью ARMA, то при этом возникает проблема подбора конкретной модели из этого класса, решение которой предусматривает три этапа:

1. идентификация модели;

2. оценивание модели;

3. диагностика модели.

На этапе идентификации производится выбор некоторой частной модели из всего класса ARMA, т.е. выбор значений p и q. Используемые при этом процедуры являются не вполне точными, что может при последующем анализе привести к выводу о непригодности идентифицированной модели и необходимости замены ее альтернативной моделью. На этом же этапе делаются предварительные грубые оценки коэффициентов a1, a2, Е, ap, b1, b2, Е, bq идентифицированной модели.

На втором этапе производится уточнение оценок коэффициентов модели с использованием эффективных статистических методов. Для оцененных коэффициентов вычисляются приближенные стандартные ошибки, дающие возможность, при дополнительных предположениях о распределениях случайных величин X1, X2, Е, строить доверительные интервалы для этих коэффициентов и проверять гипотезы об их истинных значениях с целью уточнения спецификации модели.

На третьем этапе применяются различные диагностические процедуры проверки адекватности выбранной модели имеющимся данным (misspecification tests).

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

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

Если ряд порождается моделью ARMA(p, q), то мы будем в дальнейшем для краткости обозначать это как Xt ~ ARMA(p, q). Соответственно, если ряд порождается моделью AR(p), то Xt ~ AR(p), и если ряд порождается моделью MA(q), то Xt ~ MA(q).

3.1. Идентификация стационарной модели ARMA Основной отправной точкой для идентификации стационарной модели ARMA является различие поведения автокорреляционных (ACF) и частных автокорреляционных (PACF) функций (ACF - autocorrelation function, PACF - partial autocorrelation function) рядов, соответствующих различным моделям ARMA.

О поведении автокорреляционных функций для различных моделей ARMA мы уже говорили. Однако по поведению только автокорреляционной функции трудно идентифицировать даже порядок чистого (без MA составляющей) процесса авторегрессии. Решению этого вопроса помогает рассмотрение поведения частной www.iet.ru/mipt/2/text/curs_econometrics.htm Эконометрика. Введение в регрессионный анализ временных рядов. В.П.Носко www.iet.ru автокорреляционной функции (PACF) стационарного процесса Xt. Ее значение part(k) на лаге k определяется как значение коэффициента корреляции между случайными величинами Xt и Xt+k, очищенными от влияния случайных величин Xt+1, Е, Xt+kЦ1.

Это соответствует тому, что part(k) является коэффициентом при XtЦk в линейной комбинации случайных величин XtЦ1, Е, XtЦk, наилучшим образом приближающей случайную величину Xt. Исходя из последнего, можно показать (см., например, [Hamilton (1994)]), что part(k) определяется как решение относительно ak системы первых k уравнений Юла - Уокера (s) = a1 (sЦ1) + a2 (sЦ2) + Е + ak (sЦk), s = 1, 2, Е, k, которую в этом случае удобнее записать в виде (sЦ1) a1 + (sЦ2) a2 + Е + (sЦk) ak = (s), s = 1, 2, Е, k, подчеркивая, что неизвестными здесь являются a1, a2, Е, ak, а (1Цk), Е, (kЦ1) - известные коэффициенты. Исходя из этого и применяя известное из алгебры правило Крамера решения системы k линейных уравнений с k неизвестными, находим, что вычисление PACF можно производить по формулам part(0) = 1, part(1) = (1), 1 (1) (1) (2) (2) - (1) part(2) = =, 1 (1) 1- (1) (1) 1 (1) (1) 1 (1) (2) part(3) = (1) 1 (2) (1) 1 (1), (2) (1) (3) (2) (1) K 1 (1) (2) K (1) (1) 1 (1) K (2) (2) (1) 1 K (3) M M M O M part(k) = (k -1) (k - 2) (k - 3) K (k).

1 (1) (2) K (k -1) (1) 1 (1) K (k - 2) (2) (1) 1 K (k - 3) M M M O M (k -1) (k - 2) (k - 3) K Здесь определитель в числителе выражения для part(k) отличается от определителя в знаменателе этого выражения только заменой последнего столбца столбцом, состоящим из значений (1), (2),..., (k).

www.iet.ru/mipt/2/text/curs_econometrics.htm Эконометрика. Введение в регрессионный анализ временных рядов. В.П.Носко www.iet.ru Замечательным является тот факт, что если Xt - процесс типа AR(p), то тогда part(p) 0, part(k) = 0 для k > p.

Это позволяет по графику PACF определять порядок процесса авторегрессии и отличать процесс авторегрессии от процессов скользящего среднего и ARMA(p, q) с q > 0.

Напомним, что зануление ACF после лага q соответствует процессу MA(q). Теперь же мы видим, что зануление PACF после лага p соответствует процессу AR(p).

Поэтому идентификация этих моделей по ACF и PACF более определенна по сравнению с идентификацией моделей ARMA(p, q) с p 0, q 0.

В то же время, вместо не известных нам истинных последовательностей автокорреляций (k) и частных автокорреляций part(k) мы можем довольствоваться только их состоятельными оценками - выборочной ACF, образованной выборочными автокорреляциями T -k - )(xt+k - ) (xt (k), T - k t=r(k) = = k = 1,...,T -1, T 2 (0) - ) (xt T t=где T = x = - оценка для = E(Xt), xt T t=T -k (k) = - )(xt+k - ) - оценка для (k), (xt T - k t=и выборочной PACF, образованной выборочными частными автокорреляциями rpart(k). Получить последние можно, заменяя входящие в выражения для part(k) автокорреляции (s) их оценками r(s). Однако проще поступить иначе, исходя из того, что part(k) является коэффициентом при XtЦk в линейной комбинации случайных величин XtЦ1, Е, XtЦk, наилучшим образом приближающей случайную величину Xt. Именно, можно просто оценить методом наименьших квадратов коэффициенты в модели Xt = a1 XtЦ1 + a2 XtЦ2 + Е + ak XtЦk + ut (в которой составляющая ut получается как разность ut = Xt - (a1 XtЦ1 + a2 XtЦ2 + Е + ak XtЦk), так что на нее не накладываются какие-либо предварительные ограничения).

Полученная в результате оценка коэффициента ak и есть rpart(k).

Если Xt является стационарным процессом типа ARMA(p,q) и E(X )<, то t указанные оценки , (k), r(k) и rpart(k) являются состоятельными оценками для , (k), (k) и part(k), соответственно. (См. [Hamilton (1994), p.199].) Поскольку r(k) и rpart(k) всего лишь оценки для (k) и part(k), то их наблюдаемые значения могут значительно отличаться от (k) и part(k). В частности, если при некоторых k =k1 и k =k2 в модели, порождающей наблюдения, (k1) = 0 и part(k2) = 0, то, как правило, r(k1) 0 и rpart(k2) 0, что вносит дополнительную неопределенность в задачу идентификации. Более того, характер изменения теоретической автокорреляционной функции вовсе не обязательно будет воспроизводиться в ее выборочном аналоге - выборочной автокорреляционной функции.

www.iet.ru/mipt/2/text/curs_econometrics.htm Эконометрика. Введение в регрессионный анализ временных рядов. В.П.Носко www.iet.ru Тем не менее, во многих случаях поведение теоретических ACF и PACF в какой-то мере отражается и на поведении их выборочных аналогов. Поэтому представление о поведении теоретических ACF и PACF может помочь в решении задачи идентификации соответствующих моделей в рамках общего класса моделей ARMA. В этой связи мы суммируем в следующей таблице свойства ACF и PACF для некоторых популярных моделей стационарных временных рядов.

Моде ACF PACF ль Белый (k) = 0 для k 0 part(k) = 0 для k шум, MA(0) AR(1) Экспоненциальное part(1) = a, убывание (k) = part(k) = 0, k a1 > 0 a1k AR(1) Осциллирующее part(1) = a, убывание (k) = part(k) = 0, k a1 < 0 a1k AR(p) Убывание к нулю с Зануление при k возможной осцилляцией p MA(1) Положительный Осциллирующее, пик убывание;

b1 > 0 при k = 1; part(1) > зануление при k > MA(1) Отрицательный Убывание по, пик абсолютной величине;

b1 < 0 при k = 1; part(k) < 0 при k зануление при k > MA(q) Зануление при k p ARM Экспоненциальное Осциллирующее A(1, 1) убывание с лага 1; убывание с лага a1 > 0 1;

знак (1) совпадает со знаком (a1+ b1) part(1) = (1) ARM Осциллирующее Экспоненциально A(1, 1) е убывание с лага 1;

a1 < 0 убывание с лага знак (1) совпадает 1;

со знаком (a1+ b1) part(1) = (1);

знак part(k) совпадает со знаком (1), k > ARM Осциллирующее Осциллирующее A(p, q) или прямое убывание, или прямое убывание, www.iet.ru/mipt/2/text/curs_econometrics.htm Эконометрика. Введение в регрессионный анализ временных рядов. В.П.Носко www.iet.ru начинающееся с начинающееся с лага q лага p SAR(1 Затухание на лагах, Пик на лаге, ) кратных периоду кратном периоду сезонности; зануление сезонности;

на остальных лагах зануление на остальных лагах SMA( Пик на лаге, Затухание на 1) кратном периоду лагах, кратных периоду сезонности; сезонности; зануление на остальных лагах зануление на остальных лагах Имея в виду возможность идентификации моделей AR(p) и MA(q) по графикам функций r(k) и rpart(k), желательно иметь статистические критерии для проверки гипотез о равенстве нулю тех или иных значений (k) и part(k) на основании наблюдаемых значений r(k) и rpart(k). Вопрос этот весьма сложный, и мы ограничимся только двумя приближенными рецептами, которые предполагают гауссовость инноваций (т.е., что t - гауссовский белый шум).

Если Xt - процесс типа MA(q), то при больших n q 1+ D(r(k)) 2 ( j) для k > q, T j = так что чем длинее ряд наблюдений, тем надежнее выявляются нулевые значения (k), k > q. При этом, limT E(r(k)) = (k).

Более того, при больших T и k > q распределение случайной величины r(k) близко к нормальному распределению. Отсюда вытекает, что естественный приближенный критерий проверки гипотезы H0 : УXt - процесс типа MA(q)Ф состоит в том, чтобы отвергать эту гипотезу, если q r(k) > 1+ ( j) r T j = для k > q. Уровень значимости такого критерия приближенно равен 0.05.

В частности, если q = 0, то Xt ~ MA(0) - белый шум, и гипотеза H0: УXt - белый шумФ отвергается указанным приближенным критерием при r(k) >, k > 0.

T (2) Если Xt - процесс типа AR(p), то при больших T и k > p распределение rpart(k) можно аппроксимировать нормальным распределением rpart(k) N (0, T Ц1) (так что D(rpart(k)) T - 1 ).

Следовательно, если гипотезу H0: Xt ~ AR(p) отвергать при rpart (k) >, k > p, T то получим критерий, уровень значимости которого приближенно равен 0.05.

www.iet.ru/mipt/2/text/curs_econometrics.htm Эконометрика. Введение в регрессионный анализ временных рядов. В.П.Носко www.iet.ru Имея в виду два указанных приближенных критерия, в процедурах анализа временных рядов обычно предусмотрена распечатка графиков выборочных ACF и PACF, на которые нанесены границы полосы 2/ T. В этих границах с вероятностью, близкой к 0.95, должно заключаться значение r(k), если Xt - белый шум, и значение rpart(k), если Xt ~ AR(p). Здесь следует сделать одно важное предупреждение. Именно, оба построенных критерия имеют уровень значимости, близкий к 0.05, только когда мы проверяем гипотезу H0 при некотором фиксированном k.

Что, однако, обычно происходит на практике Рассмотрим это на примере смоделированного белого шума, график которого уже приводился ранее. Всего там было получено T = 499 УнаблюденийФ x1, x2, Е, x499. В следующей таблице приведены значения выборочных автокорреляционной и частной автокорреляционной функций для значений (УлаговФ) k = 1, 2, Е, 36.

A P A P A P CF ACF CF ACF CF ACF - - 0 0 - 0.019 0.019 3.102.126 5 0.053 0.- - - - - 0.013 0.014 4 0.071 0.051 6 0.015 0.- - - - - 0.083 0.083 5 0.044 0.036 7 0.064 0.0 0 0 0 0.038.035 6.017.034 8.032.- - - - - 0.047 0.049 7 0.083 0.115 9 0.057 0.0 0 0 0 -.017.009 8.035.028 0 0.053 0.- - - - 0 0.024 0.019 9 0.049 0.085 1.011 0.0 0 0 0 0.062.053 0.069.032 2.034.0 0 0 0 0.061.069 1.041.022 3.029.0 0 - - - 0.074.073 2 0.014 0.057 4 0.042 0.0 0 - - 0 1.079.099 3 0.035 0.018 5.013.0 0 0 0 0 2.021.034 4.034.012 6.046.Рассматривая значения ACF, мы замечаем, что из полосы 2/T = 0.выбивается значение r(13) = 0.102. Означает ли это, что мы должны отвергнуть гипотезу H0: Xt - белый шум Рассматривая значения PACF, мы также обнаруживаем значения, выходящие за пределы этой полосы, что приводит к тому же вопросу.

Поскольку количество наблюдений у нас весьма велико (T = 499), можно воспользоваться утверждением об асимптотической независимости rpart(k), k = 1, 2, Е.

Pages:     | 1 |   ...   | 3 | 4 | 5 | 6 | 7 |   ...   | 35 |    Книги по разным темам