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

diplom25 / suslov_ibragimov_ekonometrika

.pdf
Скачиваний:
31
Добавлен:
24.03.2015
Размер:
5.55 Mб
Скачать

322

Глава10.Оценка параметров систем уравнений

x

s

 

 

d6

 

d5

 

d4

 

d3

 

d2

 

d1

 

p

 

Рис. 10.3

Действительно: k = 2, n = 1, r1 = 0, r2 = 1 и r1 < k − 1, r2 = k − 1.Более убедительно этот результат можно получить,используя необходимые и достаточные

условия идентификации(10.13).

Матрица H в этих условиях имеет следующий вид:

 

 

 

1

1

 

 

 

 

 

 

 

 

 

 

a11

0

 

H =

 

 

 

−b21

b22

.

 

 

 

 

 

 

c1

c2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Матрица R1 Ñпустая( rl = 0),иусловия(10.13)дляпервогоуравнения не выпол-

няются.Для второго уравнения R2 = [ 0

0 1 0 ],и матрица R2H равна [ −a11 0 ],

т.е.ее ранг равен единице,и условие(1

0.13)выполнено.А матрица,составлен-

ная из коэффициентов во всех прочих уравнениях,кроме второго,при переменных, которые исключены из второго уравнения,есть [−a11],т.е.она не вырождена.

Теперь рассматривается другая возможность:изучаемый товар входит в потребительскую корзину,и спрос на него зависит от доходов домашних хозяйств.В модель вводится переменная z2 доходов домашних хозяйств,т.е.в правую часть соотношений(10.15)добавляется слагаемое

z2 [ 0 a22 ] .

(10.17)

Если истинна модель(10.15, 10.17),то подвижной окажется линия спроса(разные домашние хозяйства имеют разные доходы),и регрессия x на p даст оценку эластичности предложения по цене(рис. 10.3).В такой ситуации не идентифицировано уравнение спроса.Уравнение предложения идентифицировано: k = 2, n = 1, r1 = 1, r2 = 0 и r1 = k − 1, r2 < k − 1.

Понятно,что можно говорить о модели,в которую входят обе отмеченные пере-

менные:и z1 и z2.ЭтоÑмодель(10.15, 10.16, 10.17).В правую часть(10.15)

10.2Взаимозависимые или одновременные уравнения

323

 

 

 

 

s1s2 s

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3s

 

 

 

 

 

 

 

 

 

 

4 s5

 

 

 

 

 

 

 

 

 

d2 d3

d4d5 d6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d1

 

 

 

 

 

 

 

Рис. 10.4

 

p

 

 

 

 

 

 

 

 

добавляется слагаемое

 

 

 

 

 

.

 

 

 

[ z1

z2

]

a11

0

 

 

 

 

 

 

 

0

a22

 

 

 

 

 

 

 

 

 

 

 

 

Вэтомслучаеидентифицированыобауравнения: k = 2, n = 1 , r1 = r2 = 1 = k−1. Но поскольку подвижны обе линииÑи спроса,и предложенияÑоблако наблюде-

ний не имеет вытянутостей(рис. 10.4),и регрессия x на p опять оказывается не значимой.Для оценки параметров регрессии требуется использовать специальные методы,рассматриваемые ниже.Впрочем,и в двух предыдущих случаях необходимо использование специальных методов оценки параметров взаимозависимых систем, т.к.обычный МНК дает смещенные и несостоятельные оценки.

Пусть теперь на предложение товара влияет еще один фактор z3,показывающий, например,количествоудобренийнаединицуплощади,скоторойсобираетсяпродукт, принимающий в дальнейшем форму товара.Тогда в правой части уравнения(10.15) возникает слагаемое

[ z1

z3

]

 

a11

0

 

,

 

 

 

 

a31

0

 

 

 

 

 

 

 

 

 

 

и первое уравнение по-прежнему остается не идентифицированным,а второе оказывается сверхидентифицированным.

Далее ряд утверждений будет иллюстрироваться на примере модели(10.15, 10.16).

Виллюстрациях эту модель удобнее записывать в сокращенном виде:

 

 

1

1

= zö1

 

 

 

 

 

[ xö pö] −β21

β22

[ α11 0 ] + [ ε1 ε2 ] .

(10.18)

Поскольку

 

 

 

 

 

 

 

 

 

 

 

1

1 −1 =

1

 

β22

−1

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−β21

β22

 

 

β21 + β22 β21

1

 

 

324

Глава10.Оценка параметров систем уравнений

приведенная форма модели имеет следующий вид:

[ xö pö] = zö1 [ d11 d12 ] + [ η1 η2 ] =

1

 

 

 

 

 

=

 

(öz1 [ α11β22

− α11 ] + [ ε1β22 + ε2

β21

ε2 − ε1 ]).

(10.19)

β21 + β22

Из этого соотношения видно,как d и η связаны с β и ε.

 

 

 

Дальнейшее изложение ведется в предположении,что строки матрицы

Rl Ñ

орты.

 

 

 

 

10.3.Оценка параметров отдельного уравнения

Вводятся дополнительные обозначения:

X l

Ñ N ×kl -матрица наблюдений заизучаемыми переменными xl ,входящи-

ми в l-е уравнение;

Xl

Ñ N -вектор-столбец наблюдений за l-й переменной xl ;

Xl Ñ N × (kl − 1)-матрица X l без столбца Xl наблюдений за xl ;

βl

Ñ kl -вектор-столбец параметров при изучаемых переменных в l-м урав-

нении;

 

βl

Ñ (kl − 1)-вектор-столбец βl с обратным знаком и без l-го элемента

βll = 1;

Z l Ñ N ×(nl +1)-матрицанаблюденийзанезависимымифакторами zl ,входящими в l-е уравнение,включая единичный столбец,соответствующий свободному члену;

αl Ñ (nl + 1)-вектор-столбец параметров при этих факторах вместе со свободным членом;

εl Ñ N -вектор-столбец остатков в l-м уравнении по наблюдениям. Тогда l-е уравнение регрессии можно записать следующим образом:

 

X l βl = Z l α + ε

l

(10.20)

 

 

l

 

или

 

 

 

 

X = X l β + Z l α + ε .

(10.21)

l

− l

l

l

 

Применение обычного МНК к этому уравнению дает в общем случае смещенные и несостоятельные оценки,прежде всего потому,что остатки εl скорее всего коррелированы с регрессорами Xl ,которые к тому же недетерминированы и наблюдаются с ошибками(гипотеза g2 нарушена).

10.3.Оценка параметров отдельного уравнения

 

 

 

 

 

325

Для иллюстрации справедливости этого утверждения используется модель(10.15,

10.16).Пусть эта модель истинна,и тогда регрессия x на p даст оценку −β22:

 

 

 

 

 

 

 

 

i i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−b22мнк =

% i2

.

 

 

 

 

 

(10.22)

Это выражение можно преобразовать,используя(10.18, 10.19) (чтобы не загро-

 

1

 

 

 

 

%

 

 

 

 

 

 

 

 

 

мождать записи,

 

обозначено через P ):

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

%pöii=

β22ii2

 

 

i=özi1d12i2

 

 

− b22мнк = P !i i

 

 

=

− β22 + P

!εi2i

 

=

 

 

 

 

 

 

 

 

 

 

= −β22 + P (d12 !i1εi2 + !ηi2εi2)

ηi2=

εi2

−εi1

 

 

 

β21

22

− β22

 

 

 

 

=

 

+

 

 

 

 

+ P #d12 !i1εi2 +

 

1

 

 

(!εi22 !εi1εi2)$.

 

 

 

 

 

 

 

 

 

 

β21 + β22

 

 

мнк

по математическому ожиданию никак не может равняться

Очевидно,что −b22

 

 

 

 

 

 

 

 

 

 

 

 

2

 

−β22,посколькувправойчастиполученноговыраженияимеется

εi2

,т.е.диспер-

%

сия(в математическом ожидании)остатка в уравнении по спросу,которая не равна нулю и к тому же не будет уменьшаться с ростом N .Эта оценка смещена и несостоятельна.

Еслиданное уравнение точно идентифицировано,то для оценки его параметров можно использовать косвенный метод (КМ)наименьших квадратов:с помощью МНК оцениваются параметры приведенной формы системы уравнений,через которые однозначно выражаются структурные параметры данного уравнения.

В качестве примера можно использовать оценку параметров второго уравнения модели(10.15, 10.16),которое точно идентифицировано.Действительно,параметры приведенной формы модели однозначно определяют оценку −β22,как это следует из(10.19):

 

−b22KM =

d11

 

 

(10.23)

 

 

.

 

 

 

 

d12

 

 

Поскольку

% i21

 

 

 

 

 

 

% i21

 

 

 

 

 

 

 

 

 

 

 

d11 =

i i1

, d12 =

 

i i1

,

 

 

 

 

%

 

 

 

что

 

 

 

 

 

 

 

 

 

 

то соотношение(10.23)означает, %

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i i1

 

 

 

 

 

−b22KM =

%ii1

,

 

 

в ка-

 

 

инструментальных переменных с

 

т.е.что(ср.с(10.22))используется метод

%

 

 

 

 

 

z1

 

честве инструментальной переменной.

326

Глава10.Оценка параметров систем уравнений

Можно записать уравнения для оценки косвенным методом в общем случае.

Сначала следует обратить внимание на то,что условия(10.11)эквивалентны требованиям

T B βl = B

,

T Aα = A ,

(10.24)

l

l

 

l l l

 

где TlB Ñ k × kl -матрица,полученная из Ik вычеркиванием столбцов,соответствующих тем изучаемым переменным,которые исключены из l-го уравнения;

TlA Ðаналогичная (n + 1) × (nl + 1)-матрица для Al .

Bl и Al имеют нулевые компоненты,соответствующие исключенным из l-го уравнения переменным.

Далее необходимо учесть,что параметры структурной формы,удовлетворяющие условиям(10.24),должны для своей идентификации еще удовлетворять соотношениям(10.10).Тем самым получается система уравнений для нахождения параметров структурной формы:

DTlB bl − TlAal = 0,

или по определению матрицы TlB :

Dl bl − TlAal = 0,

где Dl Ðоценки параметров приведенной формы уравнений для изучаемых переменных,вошедших в l-е уравнение,или,наконец,

D = Dl b

+ T Aa ,

(10.25)

l

− l

l l

 

где Dl Ñоценки параметров l-го уравнения в приведенной форме,

Dl Ñоценки параметров приведенной формы уравнений для изучаемых переменных,вошедших в правую часть l-го уравнения.

Эти матрицы коэффициентов приведенной формы представляются следующим образом:

Dl = (Z !Z )−1Z !X l , Dl = (Z !Z )−1Z !Xl , Dl = (Z !Z )−1Z !Xl .

Система уравнений(10.25)может быт ь также получена умножением обеих частей системы(10.21)слева на (Z !Z )−1Z !,т.к.третье слагаемое правой части отбрасывается(МНК-остатки должны быт ь ортогональны регрессорам),а во2-м

слагаемом (Z !Z )−1Z !Z l заменяется на TlA (т.к.по определению этой матрицы

Z l = Z T A ).

 

 

 

 

 

l

 

Dl

TlA имеетразмерность (n + 1)×

Вобщемслучае,матрицаэтойсистемы

× (kl + nl).Первый ее блок имеет

 

6

 

(7n + 1) × (kl − 1)

,второйÑ

 

размерность

 

(n + 1) × (nl + 1).

10.3.Оценка параметров отдельного уравнения

327

В случае точной идентификации и строгого выполнения условий(10.14)эта матрица квадратна и не вырождена.Система(10.25)дает единственное решениеÑоценку параметров структурной формы l-го уравнения косвенным методом наименьших квадратов.

В структурной форме со скрытым свободным членом модель(10.15+10.16)записывается следующим образом:

X P

 

1

1

 

= [ Z1

1N ]

 

a11

0

 

+ [ e1 e2 ] ,

6

7

−b21

b22

 

 

 

 

c1

c2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

аее второе,точно идентифицированное уравнение в форме(10.21) Ñ

X = P (−b22) + [ Z1

1N ]

0

+ [ e1

e2 ] .

(10.26)

c2

Как это было показано выше,обе части(10.26)умножаются на матрицу

 

Z1"

Z1

1N

−1 Z1" :

 

 

1N"

 

6

 

7

 

 

1N"

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d11

 

 

=

d12

 

( b22) +

0

 

,

 

 

d21

 

 

 

d22

 

 

 

c2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

b22) + T2Ac2,

где

T2A =

1

 

D1 = D2(

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Непосредственно в форме(10.25)при учете условий нормализации эта система записалась бы в виде:

D2b22 = −D1 + T2Ac2.

Из решения этой системы −bK22M получается таким же,как в(10.23),кроме того, получается оценка свободного члена:

cKM = d21 − d22 d11 . 2 d12

328

Глава10.Оценка параметров систем уравнений

Если уравнение не идентифицировано,переменных в системе(10.21)оказывается больше,чем уравнений,и эта систем а представляет бесконечное множество значений параметров структурной формы.Чтобы выбрать из этого множество какое-то решение,часть параметров структурной формы надо зафиксировать, т.е.сделать уравнение идентифицированным.

Длясверхидентифицированногоуравнениясистема(10.21)являетсяпереопределенной,и ее уравнения не могут выполняться как равенства.Различные методы оценки такого уравнения реализуют различные подходы к минимизации невязок по уравнениям этой системы.

Одним из таких методов является двухшаговый метод (2М)наименьших квадратов.

Напервом шагеспомощью МНК оцениваются параметры приведенной формы для переменных Xl :

Xl = Z Dl + V l ,

где V l Ñ N ×(kl −1)-матрицаостатковпоуравнениям;иопределяютсярасчетные значения этих переменных уже без ошибок:

Xlc = Z Dl .

На втором шаге с помощью МНК оцениваются искомые параметры структурной формы из уравнения:

X = X lc b

+ Z l a + e .

(10.27)

l

− l

l l

 

Для этого уравнения гипотеза g2 выполняется,т.к.регрессоры не имеют ошибок, и поэтому применим обычный МНК.

Можно определить единый оператор2M-оценивания.Поскольку

Xlc = F Xl ,

где F = Z (Z !Z )−1Z !,уравнение(10.24)записывается как:

 

 

Xl = F X l

Z l

bl + el ,

 

 

 

 

6

7

a

 

 

 

 

 

 

 

 

 

 

l

 

 

а оператор,входящий в него,как:

 

 

 

 

 

 

 

Xl! Z l

 

 

Xl! F Xl .

bl = Xl! F Xl

 

−1

 

 

 

 

 

 

 

 

 

 

al

 

 

Z l!X l

 

Z l!Z l

 

 

Z l!Xl

 

 

 

 

 

 

 

 

 

 

 

(10.28)

(10.29)

10.3.Оценка параметров отдельного уравнения

329

Оператор в такой форме получается как результат применения МНК к уравнению (10.25),т.е.результат умножения обеих частей этого уравнения слева на транспонированную матрицу регрессоров и отбрасывания компоненты остатков:

Xl! F

Xl =

Xl! F

 

F X l

Z l

bl

.

(10.30)

 

Z l!

 

 

Z l!

W

 

X al

l

 

 

 

 

 

 

 

 

 

 

 

 

Откуда следует оператор2М-оценивания в указанной форме,т.к. F Ñсимметричная идемпотентная матрица и

F Z l = F Z TlA = Z TlA = Z l.

Такой оператор оценивания сверхидентифицированного уравнения можно получить,если МНК применить к системе(10.2 1) (в этом случае она переопределена и в ее уравнениях возникают невязки),умножив предварительно обе ее части слева на Z .

Система нормальных уравнений для оценки(10.21),умноженной на Z ,записывается следующим образом:

Dl! Z "Z Dl =

Dl!

Z "Z Dl

TlA

bl

 

,

TlA!

 

 

TlA!

 

W

 

X al

 

 

 

 

 

 

 

 

 

 

 

 

и,учитывая,что

 

 

 

 

 

 

 

 

 

Dl! Z "Z Dl = X l!

F Xl, T A! Z "Z Dl = Z l! Xl и т.д.,

 

 

l

 

 

 

 

 

она преобразуется к виду(10.29).

Отсюда,в частности,следует,что для точно идентифированного уравнения 2М-оценка совпадает с КМ-оценкой,т.к.параметры структурной формы уравнения,однозначно определяемые соотношениями(10.21),удовлетворяют в этом случае и условиям(10.25).

Соотношения(10.29) Ñпервая форма записи оператора2М-оценивания.Ес- ли в(10.24)учесть,что Xlc = Xl − V l ,этот оператор можно записать в более прозрачной второй форме:

bl

 

=

Xl!

Xl − V l!V l

al

 

 

 

Z l!X l

 

 

 

 

Xl! Z l

 

−1

(Xl! − V l!)Xl

.

 

 

(10.31)

Z l!Z l

 

 

 

Z l!Xl

 

 

 

 

 

 

 

 

 

330

Глава10.Оценка параметров систем уравнений

 

Это доказывается аналогично с учетом того,что остатки V l ортогональны регрес-

 

сорам Z и,соответственно,

 

 

Z "V l = 0, X l!

V l = V l! V l, X l!cV l = 0.

 

Попытка применить оператор2М-оценивания для не идентифицированного уравнения не имеет смысла,т.к.обращаемая матрица в данном операторе вырождена.

В этом легко убедиться,т.к.

W X W X

F Xl Z l = Z Dl TlA ,

т.е.матрица наблюдений за регрессора ми в(10.25)получается умножением на Z слева матрицы системы(10.21).В последней,если уравнение не идентифицировано, Ñстолбцов больше,чем строк.Следовательно,регрессоры в(10.25)линейно связаны между собой,а матрица системы нормальных уравнений(матрица оператора оценивания)вырождена.

Для сверхидентифицированного уравнения можно использовать также метод наименьшего дисперсионного отношения (МНДО).Строгое обоснование его применимости вытекает из метода максимального правдоподобия.

Пусть bl в уравнении(10.20)оценено,и X l bl рассматривается как единая

эндогенная переменная.В результа те применения МНК определяются:

 

al = (Z l!Z l )−1Z l!X l bl ,

 

 

 

 

e = (I

N

F l )X l bl ,

где F l = Z l (Z l!Z l )−1Z l!,

(10.32)

l

 

 

 

 

 

e!e = bl!W lbl ,

где W l = X l!(I

N

F l )X l .

 

l l

 

 

 

 

 

Теперь находится остаточная сумма квадратов при условии,что все экзогенные переменные входят в l-е уравнение.Она равна bl!W bl ,где W = X l!(IN − F )X l . Тогда bl должны были бы быть оценены так,чтобы

 

bl!W l bl

λ =

bl!W bl min!

Иначе было бы трудно понять,почему в этом уравнении присутствуют не все экзогенные переменные.

Решение этой задачи приводит к следующим условиям:

(W l − λW )bl = 0.

(10.33)

10.4.Оценка параметров системы идентифицированных уравнений

331

Действительно,из условия равенства нулю первой производной:

∂λ

=

2W lbl(bl! W bl) − 2W bl(bl! W lbl )

=

2

(W l bl

λW bl ) = 0,

∂bl

bl! W bl

 

(bl! W bl)2

 

 

 

сразу следует(10.33).

Следовательно, λ находится как минимальный корень характеристического уравнения(см.ПриложениеA.1.2)

2W l − λW

2

= 0,

2

2

 

2

2

 

а bl определяетсяиз10.33сточностью допостоянного множителя,т.е.сточностью до нормировки bll = 1.

В общем случае λmin > 1,но при правильной спецификации модели

λmin −→ 1.

N →∞

Оператор

bl X l! X l − kV l!V l

= − −

a

l

Z l!X l

 

 

 

Xl! Z l

−1

(Xl!

 

 

 

 

 

 

Z l!Z l

− kV l!)Xl

Z l!Xl

позволяет получить так называемые оценки k-класса (не путать с k Ñколичеством эндогенных переменных в системе).

При k = 0,они являются обычными МНК-оценками для

l-го уравнения,что

легко проверяется;при k = 1,этоÑ 2М-оценки;при k

= λmin ÑМНДО-

оценки(принимается без доказательст ва). 2М-оценки занимают промежуточное

положение между МНК-и МНДО-оценками(т.к.

λmin > 1).Исследования пока-

зывают,что эффективные оценки получаются при

k < 1.

10.4.Оценка параметров системы идентифицированных уравнений

Из приведенной формы системы уравнений следует,что

x!ε = (B−1)!A!z!ε + (B−1)!ε!ε.

Как и прежде,в любом наблюдении E(ε) = 0, E(ε!ε) = σ2Ω,и ошибки не коррелированы по наблюдениям.Тогда

E(x!ε) = (B−1)!E(ε!ε) = σ2(B−1)!Ω,

Соседние файлы в папке diplom25