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

Эконометрика. Тихомиров

.pdf
Скачиваний:
498
Добавлен:
27.03.2016
Размер:
6 Мб
Скачать
☆

математическое ожидание М(yt|yt 0), так и вероятность P(yt 0), и, наоборот,

при отрицательном i с ростом фактора хit эти показатели уменьшаются.

Вместе с тем заметим, что эффект одновременного увеличения математического ожидания и вероятности при увеличении некоторого независимого фактора хi на практике может и не иметь место. В частности,

как показали Фин и Шмидт (Fin and Schmidt, 1984), независимая переменная

хi, увеличивающая вероятность нецензурированного наблюдения (P(yt 0)), не всегда увеличивает и математическое ожидание переменной (М(yt|yt 0)). В

качестве примера они приводят потери от пожаров в зданиях. Вероятность возникновения пожара в старом здании выше, следовательно P(yt 0)/ хit 0 (хit – возраст t-го здания), но так как старое здание стоит дешевле, то и пожар в нем приносит меньше убытков, т. е. М(yt|yt 0)/ хit 0. Таким образом, в

данной задаче предполагается, что коэффициент i при факторе “возраст здания” имеет разные знаки в функциях вероятности и математического ожидания. В рамках tobit-модели это учесть невозможно.

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

На основе probit-модели определяется вероятность нецензурированного

(или цензурированного) наблюдения при данном наборе факторов xt.

P[уt 0]= ( xt); zt =1,

 

P[уt=0]=1– ( xt); zt =0,

(10.167)

где ( xt) – интегральная функция закона нормального распределения,

определяющая вероятность нецензурированного наблюдения; – вектор

параметров модели, zt – переменная-индикатор, принимающая значение 1

для нецензурированного наблюдения и значение 0 – для цензурированного.

Далее на основе модели усеченной регрессии определяется математическое ожидание нецензурированного наблюдения. В соответствии с выражением (10.150) математическое ожидание нецензурированной переменной может быть представлено в следующем виде:

M[уt |zt =1]= xt + t.

(10.168)

Заметим, что если = / , то модель (10.167)–(10.168) сводится к tobit-

модели.

10.4.3. Модели случайно усеченных выборок (selection-model)

Предположим, что переменные у и z имеют двумерное распределение с коэффициентом корреляции . Найдем распределение у по случайной выборке (у, z) условии, что уровень переменной z превышает определенное значение (z b). Интуиция подсказывает, что если у и z положительно коррелированы, то усечение z должно подвинуть распределение у вправо.

Нахождение распределения у связано с определением, во-первых, вида функции плотности случайно усеченного распределения переменных у и z, и,

во-вторых, математического ожидания и дисперсии случайно усеченной переменной у при условии, что у и z подчинены закону двумерного нормального распределения.

Усеченная совместная плотность у и z согласно выражению (10.136) при любом распределении этих переменных определяется следующим выражением:

f ( y, z| b)

f ( y, z)

.

(10.168)

 

 

P(z b)

 

Если у и z распределены согласно двумерному нормальному закону с математическими ожиданиями y и z и стандартными отклонениями y и z, а коэффициент их парной корреляции равен , то в соответствии с выражениями (10.142)–(10.143) условные математическое ожидание и дисперсия у при усечении z определяются следующим образом:

M[y|z b]= y+ y ( z);

(10.169)

D[y|z b]= y2 [1– 2 ( z)],

(10.170)

где

 

z =(b– y)/ z;

(10.171)

( z)= ( z)/[1– ( z)];

(10.172)

( z)= ( z) [ ( z)– z].

(10.173)

Заметим, что при усечении сверху, т. е. z b, математическое ожидание и дисперсия переменной у также определяется согласно выражениям (10.169) и (10.170) при ( z)=– ( z)/ ( z).

Из выражения (10.169) следует, что при усечении “снизу” условное математическое ожидание у смещается в направлении корреляции переменных у и z, если усечение проводится “сверху”, то – в направлении противоположном корреляции. Случайное усечение уменьшает дисперсию, т.

к. ( z) и 2 принадлежат интервалу (0,1).

Рассмотрим два примера, иллюстрирующих случайное усечение.

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

если ее потенциальный доход будет превышать некоторый критический уровень (для каждой женщины свой). Допустим, переменная zt представляет собой разность между потенциальным и критическим доходом, и

зависимость между переменной zt и влияющими на нее факторами xt1 можно представить следующим образом:

zt= xt (1)+ t(1),

(10.174)

где xt(1) – вектор независимых факторов, влияющих на разность доходов

(например, возраст, образование, количество детей и т. д.); – вектор параметров модели; t(1) – ошибка модели.

Для всех женщин, у которых zt 0, требуется определить желательное количество рабочих часов yt. Предположим, что зависимость между переменной yt и влияющими на нее факторами хt(2) также можно описать линейной эконометрической моделью:

yt= xt(2)+ t(2),

(10.175)

где xt(2) – вектор независимых факторов, влияющих на желательное количество рабочих часов (например, семейный статус, количество детей и т.

д.); – вектор параметров модели; t(2) – ошибка модели.

Заметим, что вектора xt(1) и xt(2) могут как совпадать, так и отличаться друг от друга.

При формировании модели (10.175) возникает проблема усечения,

поскольку данные о часах работы имеются только для работающих женщин,

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

В разделе 10.3.1 рассматривалась модель миграции, в которой переменные,

влияющие на принятие решения о смене места жительства, были представлены эконометрическими моделями в зависимости от набора

соответствующих факторов. В целом модель содержала три уравнения:

чистая прибыль от переезда – Nt*=wt +ut;

(10.176)

доходы при переезде – ytp = xtp + tp;

(10.177)

доходы при “непереезде” – ytm = xtm + tm.

(10.178)

где wt, xtp и xtm – вектора независимых переменных, влияющих соответственно на чистую прибыль от переезда, и доходы в случае переезда и

“непереезда”; , и – вектора параметров; ut, tp и tm – ошибки модели.

Предположим, что совокупность мигрантов формируется из числа лиц,

желающих переехать, для которых чистая прибыль от переезда положительна. Чистая прибыль от переезда zt*, определяется согласно

выражению (10.176) как

 

zt*= wt +ut.

(10.179)

Для совокупности мигрантов формируется уравнение, связывающее величину их дохода на новом месте уt с некоторым набором факторов xt,

характеризующих, например, опыт работы, пол, образование и т. д:

уt= xt + t.

(10.180)

где – вектор параметров; t – вектор ошибки.

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

жительства уже осуществился, и индивидуум приступил к работе.

Поскольку доход на

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

взаимосвязаны, ошибки t

и ut моделей (10.179) и (10.180) взаимозависимы.

Предположим, что они распределены согласно двумерному нормальному закону с нулевыми математическими ожиданиями и коэффициентом корреляции . В этом случае в соответствии с выражениями (10.169)–(10.170)

получим:

M[yt | yt наблюдаемый доход]=M[yt |zt* 0]=M[yt | ut – wt]= = xt+M[ t | ut – wt]= xt + t ( u)= xt + t ( u), (10.181)

где u =– wt /u и t(u)= (wt /u)/ (wt /u).

Выражение (10.181) показывает, что условное математическое ожидание выборочной совокупности доходов мигрантов при условии zt 0 находится в

непосредственной и опосредованной зависимости от

факторов

xt.

Непосредственная зависимость выражается слагаемым

xt,

а

опосредованная, характеризующая влияние факторов xt на вероятность того,

что переменная zt* положительна, определяется слагаемым t(u).

На практике значение переменной zt* не наблюдается, она является латентной. Наблюдаемая переменная z принимает значение 1 (событие произошло) или 0 – в противном случае. В наших примерах: женщина работает или нет, индивидуум мигрирует или нет. С учетом этого представим модель (10.179)–(10.180) в виде совокупности двух следующих моделей:

1.

Модели селекции, определяющие выборку мигрантов

 

 

zt*= wt+ut;

(10.182)

 

zt=1, если zt*0;

(10.183)

 

zt=0, если zt*0;

(10.184)

 

P(zt =1)= (wt );

(10.185)

 

P(zt =0)=1–(wt ).

(10.186)

2.

Модели дохода мигранта

 

 

yt =xt + t.

(10.187)

yt

представляет собой значение дохода индивидуума,

фактически

сменившего место жительства случайную выборку мигрантов, лиц фактически сменивших место жительства, для которого zt =1.

В соответствии с введенным предположнием о зависимости между ошибками t и ut моделей (10.179) и (10.180) закон их совместного распределения характеризуется характеризуется следующими свойствами:

(ut, t) N(0,0,1, , ).

(10.188)

Согласно выражению (10.169) условное математическое ожидание yt при zt=1 определяется согласно выражению:

M[yt |zt=1]= xt+ ( wt ).

(10.189)

10.5. Методы оценки параметров моделей с дискретными и

ограниченными зависимыми переменными

Для оценивания параметров моделей с дискретными и ограниченными зависимыми переменными обычно применяется метод максимального правдоподобия. Однако, для оценивания моделей бинарного выбора используется так называемый «метод максимального счета». Рассмотрим особенности этих методов более подробно.

10.5.1. Метод максимального правдоподобия

Из-за специфических свойств моделей с дискретными и ограниченными зависимыми переменными, метод максимального правдоподобия имеет некоторые особенности. Покажем их на примере моделей бинарного выбора.

Предположим, что наблюдения y1, y2,..., yT независимы. Поскольку yt

может принимать только значения 0 и 1, то функция правдоподобия для бинарной модели имеет следующий вид:

L L( y

, y

2

,..., y

)

1

 

T

 

П [1 F(

y

0

t

 

x t)]

П F (

y

1

t

 

x t).

(10.190)

Представим выражение (10.190) в несколько другой форме:

T

 

[1 F( x t)]

L= П F [( x t)]

y

 

t

 

t 1

 

 

где переменная yt принимает значение 0 или 1.

Логарифм выражения (10.191) имеет следующий вид:

1 y

t

,

 

 

 

(10.191)

l=lnL=

T

 

{y

t

t 1

 

lnF (

x t) (1 yt) ln[1 F ( x t)]}*.

(10.192)

Необходимыми условиями максимизации функции правдоподобия являются равенства нулю всех частных производных ее логарифма по параметрам , т. е.

T

f t

 

 

f t

 

 

l/ = [

yt

(1

yt)

] xt 0.

(10.193)

 

 

(1 Ft)

t 1 Ft

 

 

 

где ft=f(xt) и Ft=F(xt), т. е. функции ft и Ft имеют аргумент xt.

Подходы к решению системы (10.193), т. е. к получению оценок коэффициентов , зависят от формы функционалов f(x) и F(x). При этом заметим, что если F(xt) нелинейны, то уравнения в (10.193) также нелинейны. Для их решения (т. е. для получения оценок параметров )

используются итеративные методы, описанные в XI главе.

В частности, для logit-модели логарифм правдоподобия вылядит следующим образом:

T

 

' xt

 

 

' xt

 

ln

e

(1 y

) ln(1 ln

e

)}, (10.194)

ln L {y

 

 

t

 

' xt

t

 

' xt

 

t 1

 

 

 

 

 

1 e

 

 

1 e

 

а необходимыми условиями его максимизации являются:

* Если распределения являются симметричными, как, например, нормальное и логистическое, то 1– F( x)=

F(– x).

T

 

e

' xt

 

 

 

l/ = ( yt

 

1 e

' x

t 1

 

 

 

Для нормального распределения логарифм правдоподобия может быть записан в следующем

) xt 0. t

функции виде:

(10.195)

максимального

l=lnL= ln[1 ( x t)]+ ln ( x t).

y

0

y

1

t

 

t

 

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

l/ =

T

 

t

t

T

 

t

t

 

0,

 

 

 

 

 

 

 

x

 

x

 

 

yt 0 1 t

 

yt 1 t

 

 

 

где

(10.196)

функции

(10.197)

t =

1

 

 

1

( y

t

' x

)' ( y

t

' x

)

 

 

 

exp

 

 

 

t

 

 

t

 

;

(10.198)

2

2

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Фt =

' xt (u)du.

(10.199)

 

 

 

Рассмотрим особенности применения метода максимального правдоподобия для оценки двумерной probit-модели типа (10.72). Как было отмечено выше, вероятность того, что зависимые переменные Y1 и Y2

системы (10.72) для конкретного индивидуума принимают соответственно значения yt1 и yt2, рассчитывается как

P(Y1=yt1, Y2=yt2)= 2(wt1, wt2, pt*), (10.200)

где

wtj= qtjztj ; ztj= jxtj, j=1,2;

(10.201)

t*=qt1 qt2 ;

(10.202)

qt1=2yt1–1; qt2=2yt2–1.

(10.203)

Логарифм функции правдоподобия будет иметь следующий вид:

 

T

 

l ln L ln 2 (wt1, wt 2 , t*),

(10.204)

t 1

 

где

w2t

2(w1t, w2t, t*)=

w

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1t

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

2

(

u1

,

u2

t*

)d

u1

d

u2

,

(10.205)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

(1/2) (w2

w2

2

t*

w

 

 

w

2t

)/(1

t*

2)

 

 

 

 

 

 

 

 

 

 

e

1t

2t

 

 

1t

 

 

 

 

 

 

(

w1t

,

w2t

,

 

)

 

 

 

 

 

 

 

 

 

 

 

 

.

(10.206)

2

t*

 

2 (1

 

2

)

1/2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t*

 

 

 

 

 

 

 

Первые производные логарифма правдоподобия по параметрам j, j=1,2; и

определяются как

 

 

 

 

T

(2 y

 

1) g

 

xtj

 

 

 

 

 

lnL/ j=

 

 

tj

tj

 

,

j 1,2;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

t 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

(2 y

 

1) (2 y

 

1)

 

 

 

 

lnL/ =

 

 

t1

 

t 2

 

 

 

2

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

t 1

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

gt1 (wt1)

wt 2

t* wt1

= [(2 yt1

1) 1 xt1]

 

 

 

 

 

 

 

 

 

 

1 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t*

 

 

 

 

 

 

 

 

 

 

 

 

Ф

(2 yt 2 1) 2 ' xt 2 t* (2 yt1 1) 1

' xt1

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t*

 

 

 

 

 

 

 

 

(10.207)

(10.208)

(10.209)

и индексы 1 и 2 в qt1 перевернуты для получения qt2.