Эконометрика. Тихомиров
.pdfУстановлено, что t не коррелированны и гомоскедастичны. Получены оценки a=0,3 и b=0,7.
Требуется:
1.Представить исходную модель в виде геометрической модели с распределенными лагами.
2.Определить реакцию дохода в году t0+3, если денежная масса в году t0
увеличилась на 1 единицу.
ГЛАВА VI. ЛИНЕЙНЫЕ МОДЕЛИ ВРЕМЕННЫХ РЯДОВ
6.1. Стационарные временные ряды
Широкий круг социально-экономических, технических и естественнонаучных процессов часто представляется набором последовательных значений показателя у1, у2,..., уt,..., уТ, зафиксированных в равноотстоящие друг от друга моменты времени t=1,2,... Т, так что интервал
(t, t+1) является постоянным. Этот набор значений уt, t=1,2,... обычно называется временным рядом (временной серией). Такой ряд представляет собой дискретный временной процесс.
Изменения значений уt во времени в реальной жизни обычно происходят под воздействием каких-либо причин, факторов. Однако в силу их многочисленности, сложности измерения, неразработанности теоретических предположений относительно взаимосвязей с переменной у и т. п. обосновать
ипостроить “подходящую” для описания процесса уt, t=1,2,...
многофакторную эконометрическую модель классического типа не всегда представляется возможным. В результате в отношении ряда уt часто выдвигается предположение, что совокупное влияние этих факторов формирует как бы внутренние закономерности в развитии процесса уt, что дает возможность применить для его описания эконометрическую модель из специфического класса моделей временных рядов.
Модели временных рядов активно применяются в исследованиях динамики значительного числа реальных процессов различной природы. Они часто используются в исследованиях динамики пассажиропотоков, складских запасов, спроса на различные виды продукции, миграционных процессов в человеческом и биологических сообществах, в радиотехнике, анализе химических процессов, моделировании природных событий: динамики числа солнечных пятен, природных катастроф и многих других процессов.
Самое широкое применение модели временных рядов нашли в исследованиях финансовых рынков, в анализе динамики финансовых
показателей, прогнозировании цен на различные товары, курсов акций,
соотношений курсов валют и т. п.
Пожалуй, общим для всех моделей временных рядов является предположение о том, что текущее значение процесса yt в значительной степени предопределено его предысторией, т. е. величина показателя yt
генерируется значениями yt–1, yt–2,... согласно характерным для этого временного ряда закономерностям. Математически это допущение может быть выражено следующим общим уравнением:
уt |
|
f ( yt 1, yt 2 ,... ) t , |
(6.1) |
где, как и ранее, t представляет собой ошибку модели в момент t.
Функция f как раз и выражает характер взаимосвязей, сложившихся в рассматриваемом временном ряду уt, t=1,2,... При удачном подборе этой функции правая “детерминированная” часть выражения (6.1) будет в некотором смысле “близка” к реальным значениям этого ряда. Как и ранее
“степень близости” обычно устанавливается по характеристикам и свойствам ряда ошибки t, t=1,2,... Здесь имеется в виду прежде всего минимальная дисперсия, соответствие белому шуму и т. п.
|
Для широкого круга процессов функция f имеет линейный вид. Например, |
|||||
y |
a1 y |
t 1 |
... an y |
t п |
t . |
Вопросы построения линейных моделей временных |
t |
|
|
|
|
||
рядов и анализа их основных свойств и будут рассмотрены в данной главе.
Линейные модели временных рядов применяются, как правило, для описания стационарных процессов. При этом обычно имеются в виду стационарные процессы второго порядка. Напомним, что стационарный процесс п-го порядка характеризуется постоянными значениями всех своих моментов порядка п и ниже на всех временных отрезках, входящих в интервал t=1,2,..., Т. У строго стационарных процессов постоянными являются моменты всех порядков. Таким образом, для любых двух интервалов времени (Т1, Т2) и (Т3, Т4) для стационарного процесса второго
должны выполняться условия, характеризующие равенство на рассматриваемых интервалах математических ожиданий, дисперсий и однопорядковых коэффициентов автокорреляций исследуемого процесса. На практике это означает, что для соответствующих оценок перечисленных показателей должны иметь место следующие соотношения:
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Т 2 |
|
|
у |
t |
|
|
Т |
4 |
у |
t |
|
|
|
|
у |
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
у |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
; |
|
|
|
|
|
(6.2) |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
T |
|
T 3 |
2 |
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
1 |
t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Т 1 T 2 |
|
|
|
|
Т 3 T 4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Т 2 |
|
( у |
|
|
|
у ) |
2 |
|
Т 4 |
( у |
|
у ) |
2 |
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
D1( y) |
|
|
t |
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
D2 ( y); |
|
|
(6.3) |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
T 2 |
|
T 1 |
T |
|
|
|
T 3 |
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
1 |
|
t |
4 |
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Т |
|
|
Т 3 |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
1 |
|
Т 2 i |
( у |
у) ( у |
t i |
у) |
Т 4 i |
|
( у |
у) ( у |
t i |
у) |
|
2 |
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
|
(6.4) |
||||||
|
|
|
|
|
|
|
|
ri |
|
(T |
|
|
|
|
i) D( y) |
|
(T |
|
|
|
|
|
|
|
|
|
|
ri ; |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
T |
|
t |
|
|
|
T |
|
i)) D( y) |
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
Т |
1 |
|
|
2 |
|
1 |
|
|
|
|
|
|
|
|
Т |
3 |
|
|
|
4 |
|
3 |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
i 1,2,..., |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
где |
у |
1 и |
у |
2 |
– оценки математических ожиданий; |
|
D1 |
( y) |
|
и |
D |
2 |
( y) |
– |
оценки |
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||
дисперсий; |
|
1 |
и |
2 |
оценки |
|
коэффициентов автокорреляции |
|
i-го |
порядка |
||||||||||||||||||||||||||||||||||
|
ri |
ri – |
|
|
||||||||||||||||||||||||||||||||||||||||
процесса |
уt |
|
на |
1-ом |
и на |
2-ом |
|
интервалах |
|
соответственно; у |
– |
среднее |
||||||||||||||||||||||||||||||||
значение процесса (оценка математического ожидания) на интервале (1,Т);
D(y) – оценка дисперсии процесса на интервале (1,Т).
Заметим, что на практике равенства (6.2)–(6.4) рассматриваются в
статистическом смысле. Иными словами, например, равенство |
у |
1 |
= |
у |
2 |
может |
|
|
в точности не выполняться. Однако гипотеза о постоянстве математического
ожидания процесса уt может быть принята, если значения |
у |
1 |
и |
у |
2 |
|
|
удовлетворяют соответствующему статистическому критерию.
Для проверки соответствия реального временного ряда уt, t=1,2,...
стационарному процессу, т. е. для проверки выполнимости условий (6.2)– (6.4), обычно используются соответствующие тесты. При этом в некоторых случаях для одного и того же условия возникает необходимость применять несколько тестов, если по результатам одного теста нельзя вывести суждение о безусловной истинности или ложности выдвинутой гипотезы. Все
множество таких тестов разделяется на три основные группы:
непараметрические, полупараметрические и параметрические тесты.
Непараметрические тесты не выдвигают заранее каких-либо сведений о законе распределения тестируемого временного ряда, его параметрах. Они исследуют взаимосвязи между порядками следования образующих его значений, выявляют наличие или отсутствие закономерностей в продолжительности и (или) чередовании их серий, образованных, например,
последовательностями единиц совокупности с одинаковыми знаками, сменой знаков у этих единиц и т.п.
Полупараметрические тесты обычно используют относительно слабые предположения о характере распределения значений временного ряда. Они,
например, относятся к общим свойствам функции распределения приростов значений ряда –симметричности, расположения квантилей. Оценки параметров распределения в таких тестах определяются по порядковым статистикам: среднее по медиане, среднеквадратическое отклонение – по размаху (абсолютной разнице между наибольшим и наименьшим значениями ряда) и т. п.
Параметрические тесты применяются при относительно строгих предположениях относительно законов распределения временного ряда, его параметров. Они, как правило, оценивают меру близости между эмпирическими характеристиками распределения временного ряда и их теоретическими аналогами. На основании величины этой меры делается вывод о целесообразности принятия или отвержения гипотезы о соответствии свойств рассматриваемого ряда стационарному процессу.
Рассмотрим некоторые достаточно часто используемые на практике тесты на стационарность более подробно.
6.1.1. Параметрические тесты стационарности
Из определения стационарного процесса второго порядка,
формализованного с помощью выражений (6.2)–(6.4), непосредственно
|
F F р*,1,2), |
(6.9) |
где F р*,1,2) – табличное значение критерия Фишера для уровня |
||
доверительной вероятности р* |
и числе степеней свободы |
1=п–1, |
2=Т1+Т2+...+Тn–п, то гипотеза о постоянстве математического ожидания временного ряда на всем интервале (1,Т) принимается с вероятностью р*. В
противном случае она отвергается.
Тестирование дисперсии.
Проверка гипотезы о постоянстве дисперсии временного ряда уt, t=1,2,...,Т
в случае разбиения исходного интервала на две части обычно осуществляется с использованием двухстороннего критерия Фишера. Обязательным условием при этом также является нормальный закон распределения значений уt .
Расчетное значение критерия Фишера определяется по следующей
формуле:
|
2 |
|
|
|
F |
s1 |
, |
(6.10) |
|
2 |
||||
|
|
|
||
|
s2 |
|
|
где |
2 |
и |
2 |
– оценки дисперсии |
ряда на первой и второй частях |
1 |
2 |
||||
|
s |
|
s |
|
|
соответственно с числом измерений Т1 |
и Т2. |
||||
Если для заданного уровня доверительной вероятности р* оказывается, что значение F удовлетворяет неравенству
F( |
1 p* |
, 1, 2) F F( |
1 p* |
, 1, 2), |
(6.11) |
|
|
||||
2 |
2 |
|
|
||
то гипотеза о постоянстве дисперсии временного ряда может быть принята, т.
е. предположение о том, что |
2 |
2 |
= |
2 |
является обоснованным с |
s1 |
= s2 |
|
вероятностью р*.
В выражении |
(6.11) |
|
|
значения |
|
|
|
1 p |
|
, 1, 2) |
и |
1 p |
, 1, 2) являются |
|||||||||||||||||
|
|
|
F( |
|
|
|
* |
F( |
* |
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
2 |
|
|
|
|
|
|
табличными (левосторонним |
|
и |
правосторонним) |
|
значениями |
критерия |
||||||||||||||||||||||||
Фишера, соответствующими |
|
вероятности |
|
ошибки |
второго |
|
рода |
( |
1 p |
) с |
||||||||||||||||||||
|
|
|
* |
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
числом степеней свободы |
|
|
1 |
=Т1–1 |
и |
|
2 |
=Т2–1. Заметим, что эти значения |
||||||||||||||||||||||
удовлетворяют следующему соотношению: |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
1 p |
|
|
|
|
|
|
|
|
|
1 p |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
F( |
* |
, |
1 |
, |
2 |
|
1 / F( |
|
* |
, |
2 |
, |
1 |
(6.12) |
|
|
|
|
|
|
||||||||||
|
|
|
) |
|
|
|
|
|
|
|
). |
|
|
|
|
|
|
|||||||||||||
|
2 |
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Вследствие этого на практике обычно проверяется только соотношение |
|
|||||||||||||||||||||||||||||
|
|
|
1 p |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
F F( |
|
|
|
* |
, |
1 |
, |
2 |
|
|
|
|
|
|
|
|
|
|
(6.13) |
|
|
|
|
|
|
||||
|
|
2 |
|
|
) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
2 |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
при условии, что 1 |
2 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
s |
s |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
При средних (Т 100) |
и больших (Т 100) объемах временного ряда вместо |
|||||||||||||||||||||||||||||
критерия Фишера для проверки гипотезы о постоянстве его дисперсии рекомендуется использовать стандартизованное нормальное распределение.
В первом случае, т. е. при средних выборках, принимается во внимание, что закону N(0,1) подчиняется случайная величина, определяемая как
|
|
1 |
ln |
s12 |
|
1 |
(1 / 1 1 / 2) |
|
|||
|
|
2 |
s22 |
2 |
|
||||||
|
|
|
|
|
|
. |
(6.14) |
||||
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|||||
|
|
|
1 / 2(1 / 1 1 / 2) |
|
|||||||
Во втором случае (при больших выборках) расчетное значение стандартизованной случайной величины оценивается следующим образом:
|
|
|
|
(s1 s2) |
s12 / 2Т 1 s22 / 2Т 2. |
(6.15) |
|
В обоих случаях, если оказывается справедливым соотношение
