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

Прикладна економетрика - Комашко О. В

.pdf
Скачиваний:
98
Добавлен:
24.05.2014
Размер:
2.11 Mб
Скачать

46

Якщо гіпотеза про стійкість моделі вірна, то

 

 

RSS RSS1

 

 

 

 

 

F =

 

n2

 

~ Fn

 

, n k .

(1.72)

 

RSS1

 

 

 

 

 

 

2

1

 

 

 

 

 

 

 

n1 k

1.3.Асимптотичні властивості МНК-оцінок.

Вбагатьох випадках скінченовимірні властивості оцінок найменших квадратів, описані вище, можуть не зберігатись. Так, якщо збурення не є нормально розподіленими, то і розподіл МНК-оцінки вже не буде нормальним. Якщо порушується умова, що всі xij некорельовані зі всіма εt , математичне

сподівання b не дорівнюватиме β . Більш того, модель лінійної регресії,

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

При послаблені деяких з класичних припущень або при переході до інших моделей скінченновимірні властивості оцінок як правило, невідомі. В таких випадках для характеризації інший підхід, який грунтується на асимптотичній теорії. Асимптотична теорія відповідає на питання, що трапиться, якщо гіпотетично розмір виборки стане нескінченно великим. Асимптотичні властивості використовуються для апроксимації скінченновимірних властивостей.

Консистентність.

Послідовність оцінок θˆn параметра θ , де n - розмір виборки, називається консистентною якщо θˆn збігаються за ймовірністю до справжнього значення параметра:

ρlimθˆn

=θ , тобто

 

 

 

δ > 0

lim P{|θˆn θ |> δ } = 0.

 

 

n→∞

 

 

 

На відміну від ситуації незалежної виборки тут нам знадобляться

додаткові умови

стосовно матриці

значень незалежних змінних X .

Припустимо, що

1

X / X збігається до

несингулярної матриці Σ XX . Точці

 

 

 

n

 

умови, що забезпечують збіжність n1 X / X , ми наводити не будемо. Зауважимо,

що вони виконуються у більшості практичних ситуацій з просторовими даними й з стаціонарними числовими рядами.

Стосовно збурень достатньо двох припущень

C1 : Mεi =0

C2 : Mxijεi =0.

Зауважимо, що умови консистеності не потребують припущень про рівність дисперсій і некорельованість збурень.

При виконанні уведених умов plim b = β .

Умова C2 є значно слабшою у порівнянні з умовою про некорельованість всіх xij зі всіма εt , а не тільки з εi , яка необхідна для забезпечення незміщенності. Так, наприклад, в моделі

yt = β0 + β1 xt + β2 yt 1 + εt

оінки параметрів будуть зміщеними, але будуть консистентними, якщо Mxt εt =0.

Кажучи, не зовсім строго, консистентність означає, що при зростанні розділу виборки ймовірність того, що оцінки будуть набагато відрізнятись від параметра прямує до нуля.

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

консистентність. Зручною властивістю є наступна. Якщо

plimb = β і g() є

неперервною функцією,то

plimg (b)= g(β).

Зауважимо

також,

що

для

консистентності s2 як оцінки

σ 2 до уведених

умов слід

додати

умову

про

рівність дисперсій збурень.

 

 

 

 

 

Асимптотична нормальність.

1 ). Оскільки

Властивість асимптотичної нормальності дозволяє використовувати стандартні критерії перевірки гіпотез без припущення про нормальність збурень, оскількі відповідні статистики, приблизно, матимуть потрібні розподіли. Припустимо, що крім припущень, уведених у попередньому параграфі виконують припущення про рівність дисперсій і некорельованість збурень. Тоді n(b β) збігається за розподілом до N (0,σ 2 Σxx

невідома матриця Σxx консистентно оцінюється за допомогою n1 (X T X ), а σ 2

через s2 , то на практиці приблизним розподілом b є

N (β ,s2 (X T X )1 ),

що співпадає з точною формулою у випадку класичних нормально розподілених збурень.

Моделі в яких порушуються припущення про рівність дисперсії і про некорельованість збурень розглядаються, відповідно, у темах “гетероскедастичність” і “автокореляція ”. В цих ситуаціях поруч зі специфічними методами можна використовувати звичайний метод найменших квадратів, але при цьому слід використовувати інші оцінки коваріаційної матриці b. МНК-оцінки зберігають властивість асимптотичної нормальності: b N (β, S ),

де замість S слід підставити оцінки, які наведено у відповідних параграфах.

33

1.4. Модель лінійної регресії з гетероскедастичними збуреннями.

1.4.1. Вступ

В цьому розділі ми розглянемо регресійні моделі, в яких порушується припущення 2 – про рівність дисперсій збурень.

Проаналізуємо залежність особистого споживання від доходу. Згадаємо, що збурення в моделі лінійної регресїї можна вважати відхиленням рівня споживання конкретного домогосподарства від середнього рівня, який відповідає даному розміру доходу. Логiчно очiкувати, що для домогосподарств з більшими доходами спостерiгатиметься бiльший розкид рівнів споживання. Отже, оскільки дисперсія збурень є мірою цого розкиду, то припущення про рівність дисперсій збурень в такій моделі буде нереалістичним.

Наведений приклад показує необхідність вивчення нового класу моделей, які узагальнюють класичну модель лінійної регресії – моделей з гетероскедастичними збуреннями.

1.4.2. Опис моделi.

Нам знадобляться наступні позначення. Літерою υ позначимо вектор збурень у вихідній моделі, а літеру ε зарезервуємо для позначення збурень, які задовольняють припущенням 1 – 5 з параграфа 1.2.1.

Розглянемо модель лінійної регресії

yi = β0 xi0 + β1xi1 +K+βk 1xi,k1 + υi ,i =

 

,

(1.73)

1, n

або, якщо ми маємо модель з константою,

yi = β0 + β1xi1 +K+βk 1xi,k1 + υi ,i =

 

,

(1.74)

1, n

34

вякій вектор збурень υ = (υ1, υ2. . . υn)T має такi властивостi: 1. Eυ = 0

2. Гетероскедастичність збурень: Dυi = σ2i , i =1, n , σ2i ≠ σ2j при i j .

3.Незалежність збурень: υі та υj незалежні при i j .

4.Незалежність збурень та регресорів: xij та υі незалежні для всіх i та j.

5.(Додаткове) Збурення υi нормально розподілені для всіх i. Припущення 2 і 3 зручно записувати у матричному вигляді:

 

2

0

 

 

0

 

 

 

σ1

 

 

 

0

σ

22

 

 

0

 

(1.75)

Dυ =

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

0

 

 

 

0

2

 

 

 

 

σn

 

1.4.3. Наслiдки гетероскедастичності збурень на оцiнки методу найменших квадратiв

1.Оцiнки МНК будуть незмiщеними, але не будуть ефективними (не матимуть найменшої дисперсiї).

2.Стандартнi оцiнки коварiацiйної матрицi оцiнки МНК будуть змiщеними, i, як наслiдок, процедури перевiрки гiпотез та iнтервального оцiнювання, основанi на стандартних статистиках, будуть некоректними.

1.4.4. Зважений метод найменших квадратiв у випадку відомої коваріаційної матриці збурень

Припустимо, що коваріаційна матриця збурень вiдома з точністю до коефіціента пропорційності, тобто

Dυi= σ2 w2

, i =

 

,

(1.76)

1, n

i

 

 

 

 

де wi2 відомі, а σ2 – невідомий коефіціент пропорційності. В системі (1.73) (або

(1.74)) почленно розділимо i-те рівняння на wi( i =1, n ):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

35

y*

 

x*

x* +K+β

 

 

 

x*

+ ε

 

,i =

 

,

(1.77)

0

k1

i

1, n

i

 

i0

 

1 i1

 

 

 

 

 

 

 

 

 

i,k1

 

 

 

 

 

 

де

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yi* =

 

 

yi

,i =

 

 

 

 

,

 

 

 

 

 

 

 

(1.78)

 

 

 

 

 

1, n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

wi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xij*

=

xij

, j =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.79)

 

 

 

0, k 1,i =1, n ,

 

 

 

 

 

 

 

w j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

εi

=

υi

, i =

 

 

 

 

 

.

 

 

 

 

 

 

 

(1.80)

 

 

 

 

 

1, n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

wi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Якщо ми розглядаємо модель з константою (1.74),

то значення змінної x0*

обчислюються за такою формулою:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xi*0 =

1

, i =

 

 

 

.

 

 

 

 

 

 

(1.81)

 

 

 

 

 

1, n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Зауважимо, що оцінка β0 -коефіціента при змінній x0* є оцінкою постійного доданку у вихідній моделі.

Зазначимо, що вектори коефiцiентiв β в моделях (1.73) i (1.77) спiвпадають.

Знайдемо коварiацiйну матрицю збурень ε в моделi (1.775). Спочатку обчислимо

математичне сподiвання:

 

 

 

 

 

 

 

Eεi = E

υi

=

1

Eυi = 0 , i =

 

.

(1.82)

1, n

wi

 

 

 

wi

 

Отже,

36

 

 

= Eε2

 

υi

 

2

1

Eυ2

 

1

 

 

 

1

σ2 w2

= σ2 . (1.83)

Dε

 

= E

 

=

=

Dυ

 

=

 

w2

w2

 

w2

 

i

i

w

 

 

i

 

 

i

 

i

 

 

 

 

 

i

 

 

i

 

 

i

 

 

 

i

 

 

Крім того,

cov(ε

,ε

 

) = Eε

ε

 

= E

υiυj

=

1

Eυυ

 

=

1

cov(υ

,υ

 

) =0

(1.84)

j

j

 

 

j

 

j

i

 

i

 

 

wi wj

 

wi wj

i

 

wi wj

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

З (1.83) і (1.84) випливає, що модель (1.77) є моделлю класичної лiнiйної регресiї.

Отже, оцiнки вектора параметрів регресії β, знайдені в моделi (1.77) методом найменших квадратів, мають бажанi статистичнi властивостi, тобто задовольняють теоремі Гауса–Маркова і, що є для нас головним, цими оцінками можна користуватись для статистичних висновків.

Означення.

Оцiнкою зваженого МНК коефiцiентiв моделi (3.1) називається оцінка звичайного МНК, знайдена за моделлю (3.5).

На практицi у бiльшостi випадкiв ваги wi невiдомі. Якщо не робити нiяких додаткових припущень, то їх оцiнити неможливо, тому що іх кількість дорівнює кількості спостережень.

1.4.5. Виявлення гетероскедастичності

Критерії виявлення гетероскедастичності розділяються на дві групи: загальні та регресійні.

Загальні критерії виявлення гетероскедастичності

Загальні критерії відрізняються тим, що при їх формуванні не використовуються припущення про характер гетероскедастичності. В цьому полягає іх перевага. Недоліком є те , що такі критерії лише виявляють наявність гетероскедастичності, але не дають інформації для розв’язання проблеми. В цьому початковому курсі ми розглянемо лише один з цієї групи критеріїв, а саме критерій Голфельда-Квондта.

37

Критерій Голфельда-Квондта.

Його використовують тоді, коли всі наявні спостереження можна розбити за деякою ознакою на дві групи. У випадку однієї незалежноої змінної споcтереження з найменшими значеннями можуть складати одну групу, а другу – спостереження з найбільшими значеннями незалежної змінної. Розбиття можна робити також за значеннями залежної змінної.

Нехай сукупність n спостережень розбита на дві групи обємами n1 і n2. Частину спостережень з середніми значеннями можна виключити. В цьому випадку n1 + n2 < n. Для того щоб застосувати критерій Голдфельда-Квондта, необхідно оцінити модель за методом найменших квадратів окремо на кожній підвиборці і знайти

σ$12 – оцінку дисперсії збурень за першою групою спостережень, та

σ$22 – оцінку дисперсії збурень за другою групою спостережень (див. (1.46).

У припущенні,що гетероскедастичність відсутня, статистика

 

σ$

2

 

F =

1

(1.85)

σ$

2

 

 

2

 

має розподіл Фішера з n1 k, n2 k ступенями свободи. Перевірка гіпотези виконується таким чином.

Якщо σ$12 > σ$22 , то обчислюють статистику (1.85) і порівнюють її з критичним значенням Fкр(α,n1 k,n2 k), знайденим за вибраним рівнем значущості α в

таблиці розподілу Фішера з n1 k, n2 k ступенями свободи. Якщо

σ$2

<

σ$

2

, то

 

 

1

 

 

2

 

обчислюють статистику

 

 

 

 

 

 

F =

σ$22

 

 

(3.86)

σ$2

 

 

 

1

 

 

 

 

 

38

і порівнюють її з критичним значенням Fкр(α,n2 k,n1 k), знайденим за вибраним рівнем значущості α в таблиці розподілу Фішера з n2 k, n1 k ступенями свободи.

Далі значення F-статистики (1.85) або (1.86) порівнюють з табличним. Якщо F < Fкр, то вважають, що гетероскедастичність відсутня. Якщо F Fкр, то вважають, що гетероскедастичність має місце.

Регресійні критерії виявлення гетероскедастичності

Регресійні критерії гетероскедастичності будуються на основі припущення, що дисперсія пропорційна функції від деякої відомої змінної:

wi2 = f (zi ), i = 1, n .

Критерій Глейзера.

Застосування цього критерія розглянемо на прикладі моделі

yi = α + βx i + υi, i =

 

.

(1.87)

1, n

На 1-му етапі модель (1.87) оцінюємо за методом найменших квадратів і знаходимо залишки uі, i = 1, n . На 2-му етапі будуємо регресію модуля залишків відносно однієї з таких функцій:

|uі| = γ + δxi + εі,

(1.88a)

|uі| = γ + δ

xi + εі,

(1.88b)

|uі| = γ + δ

1

 

+ εі,

(1.88c)

xi

 

 

 

|uі| = γ + δ

1

 

+ εі.

(1.88d)

 

xi

 

 

Зауважимо, що замість змінної x можна використати іншу змінну, яка вибирається, як правило, з економічних міркувань. Будується послідовно декілька

Соседние файлы в предмете Экономика