Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

konspekt_lektsy_ochnoe

.pdf
Скачиваний:
105
Добавлен:
10.02.2015
Размер:
2.84 Mб
Скачать

уровнях /2 и 1- /2 и числами степеней свободы (T1-1) и (T2-1). Если

F

1 / 2;T 1;T 1 F

F

 

/ 2;T 1;T 1

, то гипотеза

 

2

 

2

 

2

 

 

 

 

кр.

1

2

 

 

кр.

1

2

 

 

 

1

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ется на уровне

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Поскольку критические значения удовлетворяют соотношению

 

F / 2;v ;v

2

 

F 1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

/ 2;v ;v

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

на практике проверяется только соотношение

F F

/ 2;T

1;T 1

 

кр.

 

 

1

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

принима-

при условии, что

s2

s2

.

1

2

 

 

При достаточно больших объемах наблюдений временного ряда

T 40

 

вместо критерия Фишера рекомендуют использовать стандартизированное

нормальное распределение. При выборках

40 T 100

закону

 

ется случайная величина

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

s2

 

 

1

 

 

1

 

 

 

1

 

 

 

 

 

ln

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

s2

2

 

 

 

 

v

 

 

 

 

 

 

 

v

 

 

 

2

 

 

 

Ф

 

 

2

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

1

 

1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

v

 

 

 

 

 

 

 

 

 

 

 

v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

2

 

 

 

 

 

 

N 0,1

подчиня-

При больших выборках расчетное значение случайной величины определяется так:

Ф s

s

 

s2

s2

1

 

2

 

 

1

2

 

2T1

2T2

 

 

 

В любом случае, если Ф / 2 , то гипотеза о постоянстве дисперсии принимается.

Если временной ряд разбивается на большее число частей (n>2), гипотеза о постоянстве дисперсии может быть проверена критерием Кокрена, основанном на распределении Фишера. При этом обычно объемы этих частей принимаются равными между собой, то есть T1=T2=…=Tn. Рассчитывается критерий по формуле:

181

 

K

 

s

2

 

 

 

 

max

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s2

... s2

 

 

 

1

 

 

 

n

 

s

2

max s

2

,

Здесь

max

j

 

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Табличное значение

____

 

j 1,n

.

 

критерия Кокрена определяется по формуле:

 

F 1 / n; N 1; n 1 N 1

K ;n; N 1

 

 

n 1 F 1 / n; N 1; n 1 N 1

Если

K K ;n;N 1

,то гипотеза о постоянстве дисперсии временного

 

ряда принимается на уровне .

Критерий Бартлетта также используется при проверке гипотезы о постоянстве дисперсии. Он является более мощным, чем критерий Кокрена, но и более чувствительным к отклонениям значений временного ряда от нормального закона. Здесь временной ряд также разбивается на несколько частей, причем не обязательно одинаковых по величине.

Согласно критерию, следующая величина:

 

 

 

n

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 Ti

1 ln

si

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c i 1

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

распределена приблизительно по закону 2

с n 1 степенями свободы. Здесь

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s2

Ti 1 si2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ti 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i 1

-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

средняя дисперсия на n интервалах;

 

1

 

 

 

 

 

 

 

 

 

 

 

 

1

 

n

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 n 1

T 1

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i 1 i

 

 

 

Ti 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i 1

 

 

 

 

 

 

 

 

 

 

 

 

При больших Ti c 1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для

одинаковых

размеров

подвыборок

v1 v2

... vn

v , тогда

n

 

 

 

 

 

 

 

1

 

 

2

 

1

 

n

2

 

 

 

n 1

 

 

Ti 1 T n

 

c

nv lns

 

 

 

si

 

c 1

 

 

 

 

 

 

 

i 1

 

, поэтому

 

 

 

 

 

 

 

n i 1

 

, где

 

3 n v .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

182

 

 

 

 

 

 

 

 

 

 

Если

 

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

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]