WWW.DISSERS.RU

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

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

Pages:     | 1 |   ...   | 12 | 13 || 15 | 16 |   ...   | 35 |

Хотя выборочные автокорреляционные и частные автокорреляционные функции определялись у нас для стационарных рядов, посмотрим на коррелограммы, построенные по представленным данным. Для ряда GNP коррелограмма имеет вид Autocorrelation Partial Correlation AC PAC Q-Stat Prob. |*******. |******* 1 0.946 0.946 56.419 0.. |*******. |. 2 0.893 -0.021 107.52 0.. |******. |. 3 0.840 -0.024 153.55 0.. |******. |. 4 0.791 0.013 195.14 0.. |******. |. 5 0.743 -0.021 232.52 0.. |*****. |. 6 0.696 -0.022 265.90 0.www.iet.ru/mipt/2/text/curs_econometrics.htm Эконометрика. Введение в регрессионный анализ временных рядов. В.П.Носко www.iet.ru. |*****. |. 7 0.648 -0.030 295.41 0.. |*****. |. 8 0.599 -0.044 321.09 0.. |****. |. 9 0.550 -0.033 343.13 0.. |****. |. 10 0.498 -0.052 361.57 0.. |*** *|. 11 0.442 -0.073 376.44 0.. |***. |. 12 0.388 -0.034 388.08 0.. |***. |. 13 0.337 -0.002 397.06 0.. |**. |. 14 0.291 0.007 403.91 0.. |**. |. 15 0.253 0.041 409.21 0.. |**. |. 16 0.218 -0.002 413.22 0.. |*. |. 17 0.182 -0.034 416.08 0.. |*. |. 18 0.143 -0.044 417.90 0.. |*. |. 19 0.106 -0.021 418.92 0.. |*. |. 20 0.069 -0.039 419.36 0.. |.. |. 21 0.032 -0.028 419.45 0.. |.. |. 22 -0.003 -0.030 419.45 0.. |.. |. 23 -0.037 -0.021 419.59 0.*|.. |. 24 -0.066 -0.005 420.04 0.А для ряда NONDURABLE ACF PACF AC PAC Q-Stat Prob. |*******. |******* 1 0.917 0.917 42.921 0.. |******. |. 2 0.843 0.014 79.976 0.. |******. |. 3 0.774 -0.004 111.92 0.. |*****. |. 4 0.704 -0.040 138.99 0.. |*****. |. 5 0.643 0.013 162.09 0.. |****. |. 6 0.581 -0.041 181.36 0.. |**** *|. 7 0.510 -0.090 196.57 0.. |***. |. 8 0.439 -0.050 208.14 0.. |*** *|. 9 0.366 -0.066 216.39 0.. |**. |. 10 0.300 -0.012 222.07 0.. |**. |. 11 0.246 0.026 225.99 0.. |*. |. 12 0.196 -0.006 228.56 0.. |*. |. 13 0.150 -0.013 230.11 0.. |*. |. 14 0.106 -0.025 230.91 0.. |*. |. 15 0.072 0.036 231.28 0.. |.. |. 16 0.041 -0.016 231.41 0.. |.. |. 17 0.019 0.028 231.43 0.. |. *|. 18 -0.007 -0.063 231.44 0.. |.. |. 19 -0.032 -0.018 231.52 0.*|. *|. 20 -0.062 -0.073 231.85 0.Если отвлечься от видимой нестационарности этих рядов, то поведение выборочных ACF и PACF предполагает идентификацию обоих рядов как рядов типа AR(1). Имея в виду наличие у рядов выраженного линейного тренда, произведем оценивание моделей Xt = + t + a1Xt –1 + ut.

(Здесь мы используем обозначение ut, а не t, поскольку ряд ut на этот раз может и не являться белым шумом.) www.iet.ru/mipt/2/text/curs_econometrics.htm Эконометрика. Введение в регрессионный анализ временных рядов. В.П.Носко www.iet.ru Это приводит к следующим результатам.Для ряда GNP:

Dependent Variable: GNP Method: Least Squares Sample(adjusted): 1947:2 1961:Included observations: 59 after adjusting endpoints Convergence achieved after 3 iterations Variable Coefficient Std. Error t-Statistic Prob.

C 216.0630 11.30237 19.11661 0.T 5.269279 0.281754 18.70170 0.AR(1) 0.846976 0.072723 11.64665 0.Inverted AR Roots.Остатки обнаруживают явную автокоррелированность : P-значение критерия Бройша – Годфри при альтернативе AR(1) равно 0.0000. Переоценивание с включением в модель запаздывания на два квартала, дает:

Dependent Variable: GNP Method: Least Squares Sample(adjusted): 1947:3 1961:Included observations: 58 after adjusting endpoints Convergence achieved after 3 iterations Variable Coefficien Std. Error t-Statistic Prob.

t C 217.7399 5.054473 43.07865 0.T 5.221538 0.140436 37.18089 0.AR(1) 1.380274 0.109452 12.61078 0.AR(2) -0.630066 0.109453 -5.756490 0.Inverted AR Roots.69 -.39i.69+.39i Модули комплексных чисел, обратных авторегрессионным корням, равны 0.7926, что говорит в пользу стационарности детрендированного ряда Xt0 = Xt – µ – t.

K построению модели для ряда GNP можно подойти и иначе. Сначала произвести детрендирование ряда, оценивая модель Xt = µ + t + ut :

Variable Coefficient Std. Error t-Statistic Prob.

C 218.4825 2.640153 82.75373 0.T 5.181995 0.075274 68.84144 0.Durbin-Watson stat 0.316211 Prob(F-statistic) 0.Приводимые в таблицах оценки константы (C) и коэффициента при переменной t (T) соответствуют оценкам µ и в представлении (Xt – µ – t ) = a1(Xt – 1 – µ – ( t – 1)) + a2(Xt – 2 – µ – ( t – 2)) + ut (a2 = 0 для первой из двух таблиц). Эти оценки получаются применением нелинейного метода наименьших квадратов. При этом обозначение AR(1) указывает на оценку для a1, а AR(2) – на оценку для a2.

www.iet.ru/mipt/2/text/curs_econometrics.htm Эконометрика. Введение в регрессионный анализ временных рядов. В.П.Носко www.iet.ru Остатки, полученные при оценивании этой модели, образуют оцененный детрендированный ряд, коррелограмма которого Autocorrelation Partial Correlation AC PAC Q-Stat Prob. |******. |****** 1 0.836 0.836 44.028 0.. |**** ****|. 2 0.531 -0.554 62.115 0.. |*. **|. 3 0.183 -0.210 64.294 0..*|.. |. 4 -0.100 0.044 64.960 0.**|.. |. 5 -0.272 -0.004 69.949 0.***|..*|. 6 -0.339 -0.082 77.846 0.***|..*|. 7 -0.350 -0.169 86.446 0.***|..*|. 8 -0.332 -0.072 94.332 0.**|.. |. 9 -0.281 0.058 100.07 0.**|..*|. 10 -0.234 -0.177 104.16 0.**|. ***|. 11 -0.234 -0.321 108.32 0.**|.. |*. 12 -0.226 0.103 112.26 0.позволяет идентифицировать этот ряд как AR(2). После этого можно строить AR(2) модель для (оцененного) детрендированного ряда Xt_detrended = Xt – 218.4825 – 5.181995 t :

Variable Coeff. Std. Error t-Statistic Prob.

AR(1) 1.379966 0.107605 12.82435 0.AR(2) -0.630426 0.107605 -5.858722 0.Объединяя результаты последних двух оцениваний, получаем оцененную модель Xt – 218.4825 – 5.181995 t = = 1.379966 (Xt–1 – 218.4825 – 5.181995(t–1)) – – 0.630426 (Xt–2 – 218.4825 – 5.181995(t–2)), или Xt = [(1– 1.379966 + 0.630426)218.+ 1.3799665.181995 – 0.6304265.1819952] + (1 – 1.379966 +0.630426) 5.181995 t + 1.379966 Xt–1 – 0.630426 Xt–2 + et = 55.338375 +1.297882 t + 1.379966 Xt–1 – 0.630426 Xt–2 + et.

В то же время, по приведенным результатам оценивания модели Xt = + t + a1Xt–1 + a2Xt–2 + ut получаем Xt – 217.7399 – 5.221538 t = = 1.380274 (Xt–1 – 217.7399 – 5.221538(t–1)) – 0.630066 (Xt–2 – 217.7399 – 5.221538(t–2)), или Xt = [(1– 1.380274 + 0.630066)217.+ 1.3802745.221538 – 0.6300665.2215382] + (1 – 1.380274 +0.630066) 5.221538 t + 1.380274 Xt –1 – 0.630066 Xt –2 + et = 55.017011 + 1.304298 t + 1.380274 Xt–1 – 0.630066 Xt–2 + et, www.iet.ru/mipt/2/text/curs_econometrics.htm Эконометрика. Введение в регрессионный анализ временных рядов. В.П.Носко www.iet.ru так что результаты, полученные при использовании обоих вариантов построения модели, практически совпадают. Поэтому можно использовать смешанный вариант – использовать детрендированный ряд для идентификации порядка модели, а оценивать идентифицированную модель вместе с включенным в нее трендом, в данном случае – оценивать модель Xt = + t + a1Xt–1 + a2Xt–2 + ut.

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

Перейдем теперь к ряду NONDURABLE. Коррелограмма детрендированного ряда (ряда остатков от оцененной модели регрессии Xt на константу и t ) Autocorrelation Partial Correlation AC PAC Q-Stat Prob. |******. |****** 1 0.793 0.793 32.083 0.. |*****. |. 2 0.632 0.011 52.942 0.. |*** **|. 3 0.432 -0.195 62.887 0.. |** **|. 4 0.219 -0.193 65.515 0.. |*. |. 5 0.090 0.062 65.965 0.*|. *|. 6 -0.067 -0.152 66.218 0.**|. **|. 7 -0.242 -0.277 69.647 0.***|. *|. 8 -0.362 -0.084 77.505 0.****|. **|. 9 -0.510 -0.211 93.500 0.обнаруживает резко выделяющийся пик на лаге 1, так что можно попробовать оценить модель Xt = + t + a1Xt–1 + ut. Это дает следующие результаты:

Variable Coefficient Std. Error t-Statistic Prob.

C 47962.75 2862.678 16.75451 0.T 315.1909 76.44770 4.122961 0.AR(1) 0.884803 0.080824 10.94727 0.Наблюдаемые P-значения статистик Люнга – Бокса для ряда остатков превышают 0.05 при всех выборах M от 1 до 20. Проверка на отсутствие автокоррелированности по критерию Бройша – Годфри дает P-значения, большие 0.05, как при AR(1) альтернативе, так и при альтернативах AR(2), AR(3) и т.д. Наконец, при проверке нормальности P-значение статистики Jarque – Bera равно 0.648, так что по совокупности этих результатов мы могли бы говорить о пригодности оцененной модели Xt – 47962.75 – 315.1909 t = = 0.884803 (Xt–1 – 47962.75 – 315.1909 (t–1)) + et.

Обратим, однако, внимание на следующее обстоятельство. Оцененное значение 0.884803 коэффициента при Xt–1 достаточно близко к единице, и если ориентироваться на вычисленное значение 0.080824 стандартной ошибки для этого коэффициента, то при допущении отклонений от оцененного значения в пределах двух стандартных ошибок, в интервал допустимых значений 0.884803 ± 2·0.попадают и значения, большие или равные 1. Но последние соответствуют нестационарному процессу авторегрессии.

Таким образом, несмотря на то, что при полученной точечной оценке 0.коэффициента при Xt–1 построенная модель формально оказывается стационарной относительно детерминированного линейного тренда (т.е. детрендированный процесс следует стационарной AR(1) модели), мы не можем с достаточной степенью www.iet.ru/mipt/2/text/curs_econometrics.htm Эконометрика. Введение в регрессионный анализ временных рядов. В.П.Носко www.iet.ru уверенности гарантировать, что истинная модель порождения наблюдений также стационарна относительно линейного тренда.

Между тем вопрос о стационарности или нестационарности модели, порождающей наблюдаемый ряд, привлекает к себе постоянное внимание уже в течение нескольких десятков лет. Особенно это внимание усилилось после серии работ 80-х годов 20 века, в которых было введено понятие коинтеграции. С помощью этого понятия была обоснована методика построения “моделей коррекции ошибок”, в рамках которых удается моделировать наличие долговременных связей между переменными вместе с указанием краткосрочной динамики, обеспечивающей поддержание этих долговременных связей.

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

5.1. Нестационарные ARMA модели Начнем рассмотрение с наиболее простой модели – процесса AR(1) Xt = a1Xt–1 + t.

Мы уже знаем, что такой процесс является стационарным при выполнении условия – 1< a1 < 1. А как проявляется нестационарность ряда Xt при нарушении этого условия Приведем смоделированные реализации такого ряда при a1 = 0.5, a1 = 0.7, a1 = 0.9, a1 = 1, a1 = 1.05, a1 = 1.1.

4 2 0 -2 --4 -5 10 15 20 25 30 35 40 45 50 5 10 15 20 25 30 35 40 45 a1= 0.5 a1= 0.4 -------6 -5 10 15 20 25 30 35 40 45 50 5 10 15 20 25 30 35 40 45 a1= 0.9 a1= www.iet.ru/mipt/2/text/curs_econometrics.htm Эконометрика. Введение в регрессионный анализ временных рядов. В.П.Носко www.iet.ru ----50 100 150 200 250 300 350 400 450 a1= ----------5 10 15 20 25 30 35 40 45 5 10 15 20 25 30 35 40 45 a1= 1.05 a1= 1.Во всех случаях в качестве начального значения X1 взято X1 = 0 и использовалась одна и та же последовательность значений 1, …, T, имитирующая гауссовский белый шум с дисперсией, равной единице:

---5 10 15 20 25 30 35 40 45 NOISE Однако поведение смоделированных рядов оказалось качественно различным.

Полезно проследить, как изменяется характер траектории ряда c возрастанием значений коэффициента a1 от a1 = 0.5 до a1 = 1.1. Заметим при этом, что в порождающих моделях математические ожидания Xt равны нулю.

www.iet.ru/mipt/2/text/curs_econometrics.htm Эконометрика. Введение в регрессионный анализ временных рядов. В.П.Носко www.iet.ru Модель Кол-во Среднее пересечений значение нулевого уровня Noise (белый 25 – 0.шум) AR(1) a1 = 14 – 0.0.AR(1) a1 = 8 – 0.0.AR(1) a1 = 8 – 0.0.AR(1) a1 = 1 – 3.1.AR(1) a1 = 1 – 13.1.AR(1) a1 = 1 – 59.1.При возрастании значения a1 от a1= 0 (белый шум) до a1= 1 количество пересечений нулевого уровня уменьшается, все более длинными становятся периоды, в течение которых значения ряда находятся по одну сторону от нулевого уровня. Расширенный график ряда при a1= 1, продленный до 500 наблюдений, иллюстрирует характерное свойство реализаций процесса Xt = Xt–1 + t, состоящее в том, что такой процесс, начавшись в момент t = 1 со значения Xt = x1 (в данном случае x1 = 0), в дальнейшем очень редко пересекает уровень x(“возвращается к уровню x1 ”) и, находясь в течение длительного времени по одну сторону от этого уровня (выше или ниже), может удаляться от этого уровня на значительные расстояния.

www.iet.ru/mipt/2/text/curs_econometrics.htm Эконометрика. Введение в регрессионный анализ временных рядов. В.П.Носко www.iet.ru “Повернутая вертикально”, траектория ряда напоминает траекторию движения сильно нетрезвого человека, пытающегося продвигаться вперед по прямой, но не имеющего возможности успешно выдерживать нужное направление. И это служит некоторым оправданием термина, используемого для AR(1) процесса с a1= 1:

Xt = Xt –1 + t – случайное блуждание (процесс случайного блуждания – random walk).

Далее мы рассмотрим этот процесс подробнее, а сейчас обратим внимание на поведение реализаций процесса Xt = a1Xt–1 + t при a1 = 1.05 и a1 = 1.1. Обе реализации иллюстрируют “взрывной” (“explosive”) характер поведения AR(1) процесса при a1 > 0: траектории процесса очень быстро удаляются от начального уровня на все возрастающие расстояния. В связи с этим “взрывные” модели непригодны для описания поведения макроэкономических рядов на сколь-нибудь протяженных интервалах времени.

Пониманию столь различного поведения реализаций AR(1) процесса помогает представление модели в виде Xt – Xt–1 = a1Xt–1 – Xt–1 + t = (a1 – 1)Xt–1 + t, или Xt = Xt–1 + t, где Xt = Xt – Xt–1, = a1 – 1.

При a1 = 1 имеем = a1 – 1= 0, и приращения Xt ряда Xt образуют процесс белого шума, так что условное математическое ожидание Xt при фиксированном (наблюдаемом) значении Xt–1 = xt–1 не зависит от xt–1 и равно 0. Соответственно, при фиксированном (наблюдаемом) значении Xt–1 = xt–1, условное математическое ожидание случайной величины Xt = Xt + Xt–1 равно xt–1. Если распределение случайной величины t симметрично относительно нуля (а именно таково и гауссовское распределение, которое использовалось нами при моделировании), то www.iet.ru/mipt/2/text/curs_econometrics.htm Эконометрика. Введение в регрессионный анализ временных рядов. В.П.Носко www.iet.ru наблюдаемое значение Xt = xt может с равным успехом оказаться как больше, так и меньше xt–1. Именно это и определяет “блуждающий” характер траектории ряда.

Pages:     | 1 |   ...   | 12 | 13 || 15 | 16 |   ...   | 35 |



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

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