konspekt_lektsy_ochnoe
.pdfЕсли |
|
2 |
,n 1 |
, то гипотеза о равенстве дисперсий на рассматривае- |
|
||||
|
|
|
мых частях временного интервала принимается.
По рассмотренным параметрическим критериям следует отметить их ограниченность в применении вследствие достаточно жестких предположений нормальности закона распределения временного ряда. Кроме того, они требуют значительных вычислений. Однако реальные временные ряды могут быть распределены по закону, отличному от нормального. Поэтому на практике при проверке стационарности процессов часто используют непараметрические критерии, не имеющие ограничений по закону распределения и не столь сложные по вычислениям.
Тест Манна – Уитни используется вместо критерия Стьюдента для проверки идентичности распределений двух совокупностей, то есть временных последовательностей одного временного ряда, определенных на разных временных частях интервала 1,2,…,T. Он тестирует постоянство математического ожидания.
Пусть первая подвыборка образована T1 последовательными значениями yt, а вторая – T2 его последовательными значениями, и эти последовательности не пересекаются.
Обозначим элементы первой подвыборки символом y 1 , второй – симво-
лом y 2 . Затем объединим эти подвыборки в одну совокупность объемом (T1+T2), расположив все элементы в порядке возрастания их значений. При этом элементы подвыборок оказываются перемешанными между собой.
Если ряд стационарный, то элементы разных подвыборок довольно равномерно перемешаны друг с другом. В противном случае общая последовательность оказывается разделенной на массивы, состоящие в основном из единиц одной из совокупностей. Например, для возрастающего или убывающего временного ряда элементы подвыборок скапливаются на разных концах общей последовательности.
В тесте Манна – Уитни проверяется гипотеза о стационарности времен-
183
ного ряда на основе критерия |
u |
|
, равного числу случаев, когда элементы из |
|
|||
|
|
первой подвыборки предшествуют элементам из второй подвыборки. Значение
u рассчитывается по формулам: u |
|
R1 |
|
T1 T1 1 |
|
|
|
||||||||||||||||
|
|
2 |
или |
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
u T |
T |
T 1 |
R |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
T 1 |
1 |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
1 |
2 |
2 |
|
|
|
2 |
, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где R1 |
и R2 - суммы рангов элементов первой и второй подвыборок соответ- |
||||||||||||||||||||||
ственно, определяемых по их общей последовательности. |
|
|
|||||||||||||||||||||
|
|
|
|
Для достаточно больших последовательностей |
(T>50) случайная величи- |
||||||||||||||||||
на |
u |
|
распределена по нормальному закону с математическим ожиданием |
||||||||||||||||||||
|
|
||||||||||||||||||||||
|
|
|
|
||||||||||||||||||||
M u |
|
|
|
T |
T |
|
|
D[u |
* |
] |
T |
T |
T T 1 |
|
|
||||||||
1 |
2 |
2 |
|
|
|
1 |
2 |
|
1 |
2 |
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
и дисперсией |
|
|
|
|
|
|
|
|
12 |
|
. |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
u T1 |
T2 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
z |
|
2 |
2 |
|
|
|
|
|
|
Таким образом, случайная величина |
(u ) |
|
распределена |
||||||||||||||||
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
по закону N 0,1 . Поправка 1/2 вводится для обеспечения непрерывности z. Она
прибавляется, если z<0. Она прибавляется, если z<0, и вычитается, если z>0 Если обе подвыборки идентичны, их элементы перемешаны между собой,
тогда значение u будет находиться вблизи своего среднего значения, а значение z - около нуля. Поэтому гипотеза о стационарности процесса yt принимается на уровне значимости , если выполняется неравенство:
u1 / 2 z u / 2
Непараметрический тест Сигела – Тьюки используют вместо параметрического критерия Фишера для проверки гипотезы о постоянстве дисперсии временного ряда. Он также основан на сопоставлении рангов элементов двух подвыборок из данного интервала.
Сначала исходный временный ряд центрируется, то есть каждое значение заменяется отклонением от среднего согласно выражению yå yt y , где y -
184
среднее значение ряда yt.
Далее интервал (1,T) делится на две, желательно равные, части, где эле-
менты обозначаются соответственно |
y |
1 |
и |
y |
2 |
. Эти элементы в объединенной |
|
|
|||||
|
|
|
|
совокупности сортируются в порядке возрастания их значений. Затем каждому значению присваивается его ранг по следующему правилу: все нечетные номера получают отрицательные элементы в порядке возрастания их значений, а все четные номера – положительные элементы, но в порядке убывания их значений. Другими словами, ранг 1 получает наименьшее отрицательное значение, а ранг 2 – наибольшее положительное.
Если обозначить |
R |
сумму рангов элементов первой подвыборки, то слу- |
|||||||||||
1 |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
R |
|
T |
T T |
1 |
|
1 |
||||
|
|
|
|
1 |
|
1 |
2 |
|
|
||||
|
z |
1 |
|
|
|
|
|
2 |
|
|
2 |
||
|
|
|
|
|
|
|
|
|
|||||
|
|
T |
T |
T T |
1 |
|
|||||||
|
|
|
|
|
|||||||||
|
|
|
|
|
1 |
|
2 |
|
1 |
2 |
|
|
|
чайная величина |
|
|
|
|
|
|
|
|
12 |
|
|
оказывается распределен- |
|
|
|
|
|
|
|
|
|
|
|
|
ной по закону N 0,1 . Здесь также поправка 1/2 вводится для обеспечения непрерывности z.
Отдельную группу непараметрических тестов стационарности составляют тесты, основанные на так называемых сериальных критериях. Они анализируют закономерности серий измеренных значений временного ряда. Для их применения необходим достаточно большой объем данных, чтобы считать обнаруженные закономерности устойчивыми.
Серией называют последовательность значений временного ряда, отклоняющихся от значения некоторого признака в одну и ту же сторону. Например, при тестировании автокорреляции в остатках по методу рядов таким признаком было расчетное значение результативного признака. Во временных рядах в роли этого признака часто выступает медиана значений ряда. Тогда элементы, по значению превышающие медиану, образуют серии с положительным знаком, а элементы, не превосходящие по значению медиану – серии с отрицательным знаком.
185
Критерий Вальда – Вольфовица основан на подсчете общего числа серий. Среднее число серий рассчитывается по выражению:
M N |
|
2N N |
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
||||||
|
1 |
|
2 |
|
|
|
|
|
|
|
|
|
|
|
||||||||
s |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
N N |
|
|
|
|
, а дисперсия – по формуле: |
||||||||||||||||
|
|
|
|
|
|
|||||||||||||||||
|
|
1 |
|
|
|
|
2 |
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
D Ns |
2N N |
|
2N |
2 |
N |
|
N |
|
|
N |
|
|
||||||||||
1 |
|
|
2 |
|
|
|
|
|
|
2 |
|
|
|
|
1 |
|
2 |
|
||||
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
||||||
|
|
N |
1 |
N |
2 |
|
N |
1 |
N |
2 |
1 |
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Здесь N1 и N2 - количества элементов соответственно с положительным и |
||||||||||||||||||||||
с отрицательным знаком; (N1+N2)=T; |
N |
s |
|
|
- число серий. |
|||||||||||||||||
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Как видим, эти формулы в точности повторяют формулы метода рядов. При большом объеме временного ряда случайная величина
|
N |
s |
M N |
|
1 |
|
|
|
|
|
|
|
|||||
z |
|
|
s |
|
2 |
|
|
|
|
|
N |
|
|
N 0;1 |
|
||
|
|
|
распределена по закону |
. |
||||
|
|
|
|
|||||
|
|
|
s |
|
|
|||
|
|
|
|
|
|
|
|
Если реальный временной ряд не представляет стационарный процесс второго порядка, его нужно привести к стационарному процессу. Это делается с помощью соответствующих преобразований: взятия конечных разностей, логарифмирования цепных индексов, расчета темпов прироста и др.
Например, когда закон изменения yt близок к линейному, преобразование заключается во взятии первых разностей:
yt yt yt yt 1
Разности второго порядка:
yt yt yt yt 1 yt yt 1 yt 1 yt 2 yt 2yt 1 yt 2
применяются при законе изменения yt , близком к квадратической параболе.
При экспоненциальном росте yt логарифмируются цепные индексы:
yt ln |
yt |
ln yt ln yt 1 |
|
yt 1 |
|||
|
|
186
Расчет темпов прироста выполняется по формуле
y |
y |
t |
y |
t 1 |
|
y |
t |
1 |
||
|
|
|
|
|
||||||
t |
|
|
y |
|
|
|
y |
|
|
|
|
|
|
t 1 |
|
|
t 1 |
|
|||
|
|
|
|
|
|
|
|
Для трансформации исходного нестационарного ряда в стационарный можно использовать и другие преобразования. В каждом конкретном случае
надо исходить из примерной формы временного графика |
yt . Подходящее пре- |
образование должно обеспечивать приблизительное |
выполнение условия |
|
|
yt f yt const. |
|
Особенности конкретного стационарного процесса второго порядка полностью определяются характером его автокорреляционной функции, представ-
ляющей собой последовательность коэффициентов автокорреляции |
|
r0 , r1 , r2 , … |
|||
Здесь r0 1, остальные значения располагаются на отрезке 1;1 . |
|
|
|||
Аналогично формируется автокорреляционная функция как последова- |
|||||
тельность значений автокорреляций 0 , 1 , 2 ,… в зависимости |
от сдвига. |
||||
Между значениями двух функций существует взаимосвязь |
i |
r |
2 |
, |
i=0,1,2,…; |
|
i |
y |
|
|
0 y2 .
Модель авторегрессии–скользящего среднего (модель ARMA). Постро-
ение модели АР(k) сводится к решению двух задач:
-определение рационального порядка модели (величины k);
-оценивание параметров модели на основе уравнений Юла-Уокера.
yt 1 yt 1 2 yt 2 ... k yt k t
Система уравнений Юла-Уокера:
r1 a1 a2r1 ... ak rk 1; r2 a1r1 a2 ... ak rk 2 ;
rk a1rk 1 a2rk 2 ... ak ;
r1,r2,…rk – известные оценки коэффициентов автокорреляции;
187
a1,a2,…ak - неизвестные оценки коэффициентов модели. Модель авторегрессии первого порядка АР(1):
y |
t |
|
y |
t 1 |
|
, |
|
1 |
|
t |
|
||
a |
|
r |
|
|
|
|
1 |
1 |
|
|
|
|
Модель авторегрессии второго порядка АР(2):
y |
t |
|
y |
t 1 |
|
2 |
y |
t 2 |
|
, |
||||
|
|
|
1 |
|
|
|
|
|
t |
|
||||
r |
|
a |
a |
r ; |
|
|
|
|
|
|
||||
1 |
|
|
1 |
|
2 |
1 |
|
|
|
|
|
|
|
|
r |
|
a r a |
. |
|
|
|
|
|
|
|||||
2 |
|
|
1 1 |
|
2 |
|
|
|
|
|
|
|
||
a |
|
|
r (1 r ) |
; |
|
|
|
|
|
|||||
|
|
1 |
|
|
2 |
|
|
|
|
|
|
|||
1 |
|
|
1 r |
2 |
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
a |
|
|
r |
r2 |
|
|
|
|
|
|
|
|||
2 |
|
2 |
|
1 |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
1 r2 |
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
Модель скользящего среднего первого порядка СС(1):
yt t 1 t 1,y2 (1 12 ) 2 ,
1 12 1 1
Модель скользящего среднего второго порядка СС(2):
yt t 1 t 1 2 t 2 , |
|
|
|
|
||||||
2 (1 |
2 |
2 ) 2 |
, |
|
|
|
|
|
||
y |
1 |
2 |
|
|
|
|
|
|
|
|
|
1(1 1) ; |
|
|
2 |
; |
i |
0,i 3. |
|||
|
|
|||||||||
1 |
1 2 |
2 |
|
2 |
|
1 2 |
2 |
|
|
|
|
1 |
2 |
|
|
|
1 |
2 |
|
|
|
Простейшая модель авторегрессии - скользящего среднего АРСС(k,m) -
(AutoRegressive-MovingAverage (ARMA (k,m)) :
yt yt 1 t t 1
yt yt 1 t t 1, 1, 1
Значения автокорреляционной функции для ARMA (1,1) будут иметь вид:
(1) (1 2)( ) , 1 1 2
( ) ( 1) 1 (1), 1
188
Авторегрессионная модель проинтегрированного скользящего средне-
го (модель ARIMA).Для описания нестационарных однородных временных рядов применяется модель Бокса-Дженкинса (ARIMA –модель). Наиболее распространены ARIMA (k,m,q) – модели, со значениями параметров, не превышающими 2, q – порядок разности (дискретной производной).
Этапы методологии Бокса-Дженкинса:
1.Тестирование исходного ряда на стационарность. Анализ автокорреляционной функции. Переход к стационарному ряду путем взятия последовательных разностей (дискретные производные). Определение параметра q.
2.Исследование характера автокорреляционной функции и предположение
означениях параметров k (порядок авторегрессии) и m (порядок скользящего среднего).
3.Оценивание параметров ARIMA (k,m,q) – модели.
4.Проверка пробной модели на адекватность путем анализа ряда остатков. Для обнаружения «белого шума» в остатках применяют Q-статистику
Бокса-Пирса, H0 об отсутствии автокорреляции в остатках:
Q n rp2 ,
p 1
Q 2 ( ,v k m) H0 : 0
Критерии качества подгонки модели Бокса-Дженкинса:
КритерийАкайка (Akaike information criterion, AIC):
|
|
|
n |
2 |
|
|
k m |
|
et |
|
|
AIC |
ln |
t 1 |
|
|
|
n |
n |
|
|||
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
Выбор следует сделать в пользу модели с меньшим значением AIC.
Критерий Шварца (Swarzcriterion):
|
|
|
n |
2 |
|
|
|
( p q)ln n |
|
et |
|
|
|
SIK |
ln |
t 1 |
|
|
|
|
n |
n |
|
|
|||
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
189 |
|
|
|
|
|
|
Вопросы для самоконтроля
1.Какая модель временного ряда называется статической?
2.Когда модель временного ряда называется динамической?
3.Как определяются авторегрессионные модели?
4.Как определяется модель ARMA?
5.Как интерпретируют параметры моделей авторегрессии?
6.Что означает стационарность временного ряда?
7.Какой стационарный процесс называется «белым шумом»?
8.Какие типы включают модели стационарных временных рядов?
9.Какие типы включают модели нестационарных временных рядов?
10.Как определяется ARIMA-модель?
Лекция 20
Тема 18. Модели с лаговыми переменными Вопросы для изучения
1.Статические и динамические модели.
2.Модели с распределенным лагом.
3.Модель частичной корректировки и модель адаптивных ожиданий. Аннотация. Данная тема раскрывает особенности моделей с лаговыми пе-
ременными и методы их оценивания.
Ключевые слова. Модель с распределенным лагом, метод Алмон, метод Койка, модель частичной корректировки, модель адаптивных ожиданий.
Методические рекомендации по изучению темы
Изучить лекционную часть, где даются общие представления по данной теме.
Для закрепления теоретического материала ответить на вопросы для самоконтроля.
Для проверки усвоения темы выполнить практические задания и тест для самоконтроля.
190