Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Бардасов - Эконометрика.pdf
Скачиваний:
423
Добавлен:
06.03.2016
Размер:
2 Мб
Скачать
☆

§ 3. Дисперсии и стандартные ошибки коэффициентов

Знание дисперсий и стандартных ошибок позволяет анализировать точность оценок, строить доверительные интервалы для теоретических коэффициентов, проверять соответствующие гипотезы. Выборочные дисперсии эмпирических коэффициентов регрессии имеют вид:

где z jj

n

∑ei2

2

i=1

2

 

S b j =

 

z jj =S e z jj ,

 

n−m−1

(3.12)

- j-ый диагональный элемент матрицы

Z =(X T X )−1,

j =0,1, 2, ..., m.

 

Отметим, что здесь в матрицах X T X и (X T X )−1 первая строка

истолбец обозначены цифрой 0.

Впримере 3.1

z00 =7,310816; z11 =0,001593; z22 =0,043213 .

Как и в случае парной регрессии, Sb j =

2

называется стан-

Sb j

дартной ошибкой коэффициента регрессии.

 

 

В примере 3.1

 

 

 

 

n

 

 

 

 

 

∑ei2

 

 

 

 

 

i=1

 

=81,61831;

Sb0 =24,4; Sb1 =0,361;

Sb2 =1,88 .

 

n−m−1

 

 

 

 

 

Отметим, что рассмотренная выше парная линейная регрессия также может быть изложена в матричном виде.

§ 4. Проверка статистической значимости коэффициентов уравнения регрессии

Построение эмпирического уравнения регрессии — начальный этап эконометрического анализа. Первое же построенное по выборке уравнение регрессии очень редко является удовлетвори-

71

тельным по тем или иным характеристикам. Поэтому следующей важнейшей задачей эконометрического анализа будет проверка качества уравнения регрессии. Проверка статистического качества оцененного уравнения регрессии проводится по следующим направлениям:

•проверка статистической значимости коэффициентов уравнения регрессии;

•проверка общего качества уравнения регрессии;

•проверка свойств данных, выполнимость которых предполагалась при оценивании уравнения (проверка выполнимости предпосылок МНК).

Как и в случае парной регрессии, статистическая значимость коэффициентов множественной линейной регрессии с m объясняющими переменными проверяется на основе t-статистики:

t =

b j

,

(3.13)

 

 

Sb j

 

имеющей в данной ситуации распределение Стьюдента с числом степеней свободы ν=n−m−1 (n — объем выборки, m — количество объясняющих переменных в модели). При требуемом уровне значимости наблюдаемое значение t-статистики сравнивается с критической точкой t α2 , n−m−1 распределения Стьюдента.

Если t >t α2 , n−m−1 , то коэффициент bj считается статистически значимым. В противном случае ( t ≤t α2 , n−m−1) коэффициент bj счи-

тается статистически незначимым (статистически близким к нулю). Это означает, что фактор Xj линейно не связан с зависимой переменной Y. Его наличие среди объясняющих переменных не оправдано со статистической точки зрения. Он не оказывает скольконибудь серьезного влияния на зависимую переменную, а лишь искажает реальную картину взаимосвязи. Если коэффициент bj статистически незначим, рекомендуется исключить из уравнения регрессии переменную Xj. Это не приведет к существенной потере качества модели, но сделает ее более конкретной.

72

 

Впримере3.1 при α=0,05

(5%) t

 

α

,n−m−1 =t

0,05

; 7 =2,365 . Тогда

 

 

 

 

95,5

 

 

 

 

0,818

2

 

2

 

 

 

7,68

 

tb0

=

b0

=

=3,91;tb1

=

b1

=

=2,26

;tbj =

b2

=

=4,09.

Sb0

24,4

Sb1

0,361

Sb2

1,88

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, при уровне значимости α=0,05 (5%) коэф-

фициенты b0 и b2 статистически значимо отличаются от нуля, а коэффициент b1 статистически незначим. Коэффициент b1 становится

статистически значимым при α≈0,058 (5,8%).

§ 5. Интервальные оценки коэффициентов теоретического уравнения регрессии

По аналогии с парной регрессией после определения точечных оценок bj коэффициентов βj (j =1, 2, ..., m) теоретического урав-

нения регрессии могут быть рассчитаны интервальные оценки указанных коэффициентов. Доверительный интервал, накрывающий с надежностью 1−α неизвестное значение параметра βj , определя-

ется неравенством:

b j −t α2 , n−m−1Sb j < βj <b j +t α2 , n−m−1Sb j ,

где t α2 , n−m−1 — критическая точка распределения Стьюдента с чис-

лом степеней свободы ν=n−m−1 (n — объем выборки, m — количество объясняющих переменных в модели) и уровнем доверия α.

В примере 3.1 при α=0,05 (5%) t

α

, n−m−1 =t

0,05

; 7 =2,365 . До-

 

 

2

2

 

верительные интервалы имеют вид:

 

 

 

95,5−2,365×24,4< β0

<95,5+2,365×24,4 ,

0,818−2,365×0,361< β1

<0,818+2,365×0,361,

7,68−2,365×1,88< β2

<7,68+2,365×1,88 .

По аналогии с парной регрессией может быть построена интервальная оценка для среднего значения предсказания:

73

ˆ

 

−

α

, n−m−1 e

 

Tp(

 

TX)−1

 

p <M

 

p

 

 

 

Tp

<ˆ

+

α

, n−m−1 e

 

 

Tp(

 

TX)−1

 

p .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y p

 

t

2

 

 

 

S

X

 

X

 

X

 

(Y

 

 

 

X

)

Y p

t 2

S

X

 

X

 

X

 

 

 

 

Пусть в примере 3.1

x p1 =60,

 

 

 

x p2 =4 , тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X Tp (X T X )−1 X p =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7,31081

 

−0,10049

 

−0,53537

 

 

1

 

 

 

 

(1

 

60

 

 

 

 

 

 

0,001593

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4)× −0,10049

 

−0,006644 × 60 =0,585;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−0,006644

0,043213

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−0,53537

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

ˆ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y p =95,5+0,818×60−7,680×15=113,9 ;

 

 

 

 

 

 

 

 

 

 

t

α

,n−m−1S e

X Tp(X T X )−1 X p =2,365× 81,6 × 0,585 =13,3 ;

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

113,9−13,3<M (Y p

 

 

X Tp)<113,9+13,3.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

X Tp(X T X )−1 X p → 1+X Tp(X T X )−1X p = 1,585 =1,259 .

§ 6. Проверка общего качества уравнения регрессии

После проверки значимости каждого коэффициента регрессии обычно проверяется общее качество уравнения регрессии. Для этой цели, как и в случае парной регрессии, используется коэффициент детерминации:

 

n

 

 

 

2

 

 

 

n

 

 

 

 

 

 

n

 

 

 

 

∑(ˆyi−

y

)

 

 

 

∑ei2

 

 

 

 

 

n∑ei2

 

 

 

R2 =

i=1

 

 

 

 

= 1−

 

 

i=1

 

 

=1−

 

 

 

i=1

 

 

(3.14)

n

 

 

)

2

n

(

 

)

2

 

 

 

n

 

2

 

∑

y −

y

 

∑

y −

y

 

n

 

2

 

 

 

i=1(

i

 

 

i=1

i

 

 

n∑ y

i

 

− ∑ y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i=1

 

 

i=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Суть данного коэффициента как доли общего разброса значений зависимой переменной Y, объясненного уравнением регрессии,

74

подробно рассмотрена ранее. Как отмечалось, в общем случае 0≤R2 ≤1. Чем ближе этот коэффициент к единице, тем больше

уравнение регрессии объясняет поведение Y. Поэтому естественно желание построить регрессию с наибольшим R2.

Для множественной регрессии коэффициент детерминации является неубывающей функцией числа объясняющих переменных. Добавление новой объясняющей переменной никогда не уменьшает значение R2. Действительно, каждая следующая объясняющая переменная может лишь дополнить, но никак не сократить информацию, объясняющую поведение зависимой переменной. Это уменьшает (в худшем случае не увеличивает) область неопределенности в поведении переменной Y.

Иногда при расчете коэффициента детерминации для получения несмещенных оценок в числителе и знаменателе вычитаемой из единицы дроби делается поправка на число степеней свободы.

Вводится так называемый скорректированный (исправленный) коэффициент детерминации:

 

 

 

n

(n−m−1)

 

 

 

 

 

 

 

∑ei2

n−1

 

 

 

 

 

2 =1−

i=1

 

 

 

 

 

=1−(1−R2)

,

(3.15)

R

 

 

 

 

 

n

 

 

 

2

 

n−m−1

 

 

 

∑(yi−y)

 

(n−1)

 

 

 

 

 

 

 

 

 

 

 

 

 

i=1

 

 

 

 

 

 

 

 

 

 

где n — количество наблюдений, m — число объясняющих переменных.

Такая поправка учитывает соответствующие степени свободы RSS и TSS. Выражение (3.15) можно также представить в виде:

 

 

 

 

n

 

(n−k)

 

 

 

 

 

2 =1−

 

∑ei2

 

=1−(1−R2)

n−1

 

 

 

 

i=1

 

 

 

 

,

R

 

 

 

 

 

n

 

 

 

2

 

n−k

 

 

 

∑(yi−y)

 

(n−1)

 

 

 

 

 

 

 

 

 

 

 

 

i=1

 

 

 

 

 

 

 

 

где k =(m+1) — число параметров уравнения регрессии.

Из (3.15) очевидно, что R2 <R2 для т > 1. С ростом значения m скорректированный коэффициент детерминации растет медленнее, чем (обычный) коэффициент детерминации. Другими словами, он

75

корректируется в сторону уменьшения с ростом числа объясняю-

щих переменных. Нетрудно заметить, что

 

 

2 =R2 только при

R

R2=1.

 

2 может принимать отрицательные значения (например,

R

при R2=0). Обычно полагают

 

 

2 =max 1−(1−

R

2)

n−1

 

, 0 .

R

 

 

 

 

 

 

 

n−m−1

 

 

 

 

 

 

 

 

 

 

 

 

Доказано, что R2 увеличивается при добавлении новой объяс-

няющей переменной тогда и только тогда, когда t-статистика для этой переменной по модулю больше единицы. Поэтому добавление в модель новых объясняющих переменных осуществляется до тех пор, пока растет скорректированный коэффициент детерминации.

Обычно приводятся данные как по R2, так и по R2, являющиеся

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

После оценки индивидуальной статистической значимости каждого из коэффициентов регрессии обычно анализируется совокупная значимость коэффициентов. Такой анализ осуществляется на основе проверки гипотезы об общей значимости — гипотезы об одновременном равенстве нулю всех коэффициентов регрессии при объясняющих переменных:

H 0 : β1 = β2 =...= βm =0.

Если данная гипотеза не отклоняется, то делается вывод о том, что совокупное влияние всех т объясняющих переменных X 1,X 2, ...,X m модели на зависимую переменную Y можно считать

статистически несущественным, а общее качество уравнения регрессии — невысоким.

Проверка данной гипотезы осуществляется на основе дисперсионного анализа сравнения объясненной и остаточной дисперсий.

H 0 : (объясненная дисперсия) = (остаточная дисперсия),

76

H 1 : (объясненная дисперсия) > (остаточная дисперсия). Строится F-статистика:

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

∑(ˆyi −

 

)2 m

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

F =

 

i=1

,

(3.16)

 

 

 

 

 

 

n

 

 

 

 

 

 

 

∑(yi−ˆyi)2 (n−m−1)

 

 

 

 

 

 

 

 

 

i=1

 

 

n

 

 

 

)2 / m =ESS / m — объясненная

 

где ∑(ˆyi−

y

сумма квадратов в

i=1

 

 

 

 

 

 

 

 

 

 

 

 

расчете на одну независимую переменную;

 

 

n

y

−ˆy

i)

2 / (n−m−1)=RSS / (n−m−1)

— остаточная сумма

∑

i=1 (

i

 

 

 

 

 

 

 

 

 

 

квадратов в расчете на одну степень свободы.

 

 

При выполнении предпосылок МНК построенная F-статистика

имеет

распределение

Фишера с числами

 

степеней свободы

ν1 =m, ν2 =n−m−1. Поэтому, если при требуемом уровне значимости α, F набл >F α; m; n−m−1, где F α; m; n−m−1 — критическая точ-

ка распределения Фишера, то H0 отклоняется в пользу H1. Это означает, что объясненная дисперсия существенно больше остаточной дисперсии, следовательно, уравнение регрессии достаточно качественно отражает динамику изменения зависимой переменной Y. Если F набл <F α; m; n−m−1 , то нет оснований для отклоне-

ния H0. Значит объясненная дисперсия соизмерима с дисперсией, вызванной случайными факторами. Это дает основание считать, что совокупное влияние объясняющих переменных модели несущественно, следовательно, общее качество модели невысоко.

Однако на практике чаще вместо указанной гипотезы проверяют тесно связанную с ней гипотезу о статистической значимости коэффициента детерминации R2:

H 0 : R2 =0,

H 1 : R2 >0.

Для проверки данной гипотезы используется следующая F- статистика:

77

F =

R2

×

n−m−1

.

(3.17)

1−R2

m

 

 

 

 

 

Выражение (3.17) следует из (3.16), если числитель и знаменатель разделить на TSS.

Величина F при выполнении предпосылок МНК и при справедливости H0 имеет распределение Фишера, аналогичное распределению F-статистики (3.16). Показатели F и R2 равны или не равны нулю одновременно. Для проверки нулевой гипотезы

H0:F =0, R2 =0 при заданном уровне значимости α по таблицам

критических точек распределения Фишера находится критическое значение F крит =F α; m; n−m−1 . Нулевая гипотеза отклоняется, если

F набл >F крит . Это равносильно тому, что R2 >0 , т. е. R2 статисти-

чески значим.

Анализ статистики F позволяет сделать вывод о том, что для принятия гипотезы об одновременном равенстве нулю всех коэффициентов линейной регрессии коэффициент детерминации R2 не должен существенно отличаться от нуля. Его критическое значение уменьшается при росте числа наблюдений и может стать сколь угодно малым.

n

n ∑ei2

В примере 3.1 R2 =1−

 

 

 

 

i=1

 

 

 

 

=0,936 .

n

 

 

 

n

 

 

2

 

 

 

 

 

2

 

 

 

 

 

 

 

n ∑ y

−

∑ y

 

 

 

 

 

 

 

 

 

i

 

 

 

i

 

 

 

 

 

 

i=1

 

 

i=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда F =

R2

 

n−m−1

=

0,936×7

=51,2 . По таблицам кри-

1−R2

 

0,064×2

 

 

m

 

найдем F0,05; 2;7 =4,74 ;

тических точек распределения

 

 

Фишера

F0,01; 2;27 =9,55 . Поскольку Fнабл =51,2 > Fкрит как при 5%-м, так и

при 1%-м уровне значимости, то нулевая гипотеза H0 в обоих случаях отклоняется в пользу H1. Это означает, что объясненная дисперсия существенно больше остаточной, следовательно, уравнение регрессии достаточно качественно отражает динамику изменения зависимой переменной Y.

78

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