72239

Логит- и пробит-модели

Лекция

Экономическая теория и математическое моделирование

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

Русский

2014-11-19

532 KB

39 чел.

4.6 Логит- и пробит-модели

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

  1.  Метрической шкале (являться количественными).
  2.  Порядковой шкале.
  3.  Номинальной шкале.

В этой главе мы рассмотрим, в частности,  эконометрические модели, которые используют в том случае, когда ЗАВИСИМАЯ  переменная измерена в номинальной или порядковой шкале.

Предположим, что результирующий показатель у, «поведение» которого существенно зависит от количественных объясняющих переменных (матрица Х)

Х = (1, х1, ,...,хm-1) ,

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

yi  может быть приписано значение равное 0, если i-й индивидуум оказался безработным в обследуемом периоде времени, и 1 — в противном случае.

То есть:

В подобных ситуациях вектор

Y = (у1, y2,…,yn)T

исходных статистических данных зависимой переменной будет состоять только из: «0» или «1

 Можно ли построить линейную регрессионную модель, описывающую зависимость у от х в данном случае?

Ответ:  вряд ли.  Неясно, как интерпретировать в этом случае оцененные при помощи регрессии значения у,  которые будут уже измерены в метрической непрерывной (количественной) шкале и могут принимать различные значения.

Поэтому для исследования статистической связи между у и Х строят некоторую специальную регрессионную модель зависимости вероятности

Р{у = 1|Х}

от линейной функции наблюдаемых факторов.

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

Если

>,

то женщине выгоднее пойти работать, так как  получаемый дополнительный доход перевешивает уменьшение времени на детей и домашние дела.

Если

<,

то женщина не выходит на работу.

Обозначим разность

-=

и предположим, что эта величина является линейной функцией от наблюдаемых характеристик хi: величины заработной платы, возраста, наличия детей и т.д.

Предположим, что

=0+1х1+...m-1xm- 1   +

Тогда 

Р{у = 1}= Р{у*0}=

= Р{0+1х1+...m-1xm-1  + 0}=

= Р{ - (0+1х1+...m-1xm-1)}=

Р{- <0+1х1+...m-1xm-1}=

=F(0+1х1+...m-1xm-1).

Таким образом,

Р{у = 1/х}= F(0+1х1+...m-1xm-1),

где, F-это функция распределения случайной величины (- )

Во многих экономических задачах, скрытую переменную вводят сразу, считая, что:

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

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

Р{у = 1|Х}

подбирают функции, область значений которых определяется отрезком [0,1], а линейная функция

z=0+1х1+...m-1xm-1

играет роль аргумента этой функции, т.е.

Р{у =1|Х} =F(0+1х1+...m-1xm-1) =F(z)

(4.6.1)

причем функции F(z)(какую бы мы не выбрали!!!) должны удовлетворять следующим требованием

(4.6.2)

График одной из функций данного вида приведен на рисунке:

.

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

ОПРЕДЕЛЕНИЕ 4.6.1. Логит-модель. Модель вида (4.6.1) при условиях (4.6.2) для функции F(z), называется логит-моделью, если в качестве F(z) рассматривается логистическая функция, т. е.

(4.6.3)

Нетрудно проверить, что эта функция удовлетворяет условиям (4.6.2), ее часто обозначают через Λ(z).  

ОПРЕДЕЛЕНИЕ.4.1.2 Пробит-модель. Модель вида (4.6.1) называется пробит-моделью, если в качестве F(z) рассматривается  функция распределения стандартного нормального распределения, т.е.

(4.6.4)

Так же, как и функция Λ(z), нормальная  функция распределения F(z) удовлетворяет всем условиям (4.6.2) и является симметричной относительно z = 0.

Figure 4.1: Relationship between latent variable yand Pr(y = 1) for the BRM.

Рисунок заимствован из книги «REGRESSION MODELS FORCATEGORICAL DEPENDENT VARIABLES USING STATA»

Оценивание параметров в логит- и пробит- моделях.

Роль функции регрессии играет функция Λ(z) в логит-модели и функция F(z) в пробит-модели. Для того, чтобы оценить параметры j, например, в логит-модели, используют метод максимального правдоподобия.

Что такое функция правдоподобия?

Согласно определению, функцией правдоподобия называют функцию вида:

L=L(y1,y2,..,yn)=P{Y=y1}P{Y=y2}P{Y=yn} =

(например!)=P{Y=1}P{Y=0}…P{Y=0}

(4.6.5)

В нашем случае, исходя из (4.6.1), мы можем записать эту функцию в виде:

L(y1,y2,…,yn)=

(4.6.6)

Неизвестные параметры I и для логит и для пробит модели ищут из условия максимума функции (4.6.6). Полученные при этом оценки вi  называют оценками максимального правдоподобия.

Запишем формулы для оценивания вероятности события Р{У=1/х}, она имеет вид:

Для логит модели:

Для пробит модели,

где Ф-это функция Лапласа.

Реализация процедур оценивания параметров логит и пробит моделей проводится, например, в программе Stata.

Рассмотрим, как выглядит вывод итогов в Стате.

 

  

 401K.DES

prate     mrate     totpart   totelg    age       totemp    sole      ltotemp

 Obs:  1534

     

 1. prate                       participation rate, percent

 2. mrate                       401k plan match rate

 3. totpart                     total 401k participants

 4. totelg                      total eligible for 401k plan

 5. age                         age of 401k plan

 6. totemp                      total number of firm employees

 7. sole                        = 1 if 401k is firm's sole plan

 8. ltotemp                     log of totemp

4.7. Статистическая значимость логит и пробит - моделей и факторов этих моделей

Проверка значимости модели при помощи теста отношения правдоподобия(тест Вальда), начинается с выдвижения основной гипотезы:

Для проверки данной гипотезы вычисляется выборочная статистика

Z*=LR

(4.7.1)

LR=2(lnL-lnL0).

Здесь lnL величина  максимального значения логарифма функции правдоподобия, а lnL0- величина логарифма функции правдоподобия в случае справедливости основной гипотезы.

Если основная гипотеза верна, то выборочная статистика (4.7.1)  распределена по закону 2 с (m-1) степенью свободы. Границу правосторонней  критической области  К2  ищут по таблицам критических точек хи-квадрат по  уровню значимости (1-α) и (m-1)  степени свободы. Если  выполняется неравенство:

Z*≥ К2,

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

Для моделей бинарного выбора, значимость факторов проверяется при помощи тестирования для каждого фактора    хi, i=1,…, (m-1)  гипотез вида:

Н0:  i=0,

Ha:  i≠ 0.                                            

(4.7.2)

Выборочные статистики, которые используются для тестирования этих гипотез,    имеют асимптотически нормальное распределение и называются z-статистиками. Границу двусторонней критической области ищут по таблицам Лапласа по заданному уровню значимости (1-α).

Если  выполняется неравенство:

К 1<Z*<К2

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

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

Определение 4.7.1.  Псевдо - коэффициентом детерминации называют следующую величину:

 Определение  4.7.2.  Индексом отношения правдоподобия Макфаддена (McFadden)  называют характеристику:

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

 

На лекции мы рассмотрели нелинейные регрессионные модели, в частности, модели для бинарной зависимой переменной. Эти модели мы рассмотрели для двух функций регрессий: логит (использовали логистическую функцию) и пробит (использовали функцию распределения стандартного нормального закона распределения). Оценки параметров таких функций регрессии получают при помощи метода максимального правдоподобия. Модель тестируют при помощи теста Вальда, в основе которого статистика, имеющая хи-квадрат распределение.

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

Вернемся к моделям бинарного выбора. Если мы попытаемся найти производную от P{Y=1|X}, то придем к следующему выражению:

P{Y=1|X}=

где Z= 0+1х1+...m-1xm-1.

По теореме о производной сложной функции, и из свойства плотности   (производная от функции распределения это плотность распределения f(Z)), получаем:

P{Y=1|X}=

или, используя второе обозначение для оценок параметров:

P{Y=1|X}=

P{Y=1|X}=вjf(Z)

Как и раньше, через вj обозначены оценки неизвестных параметров.

Тогда, мы можем рассуждать следующим образом: плотность распределения всегда неотрицательна, поэтому знак  производной P{Y=1|X} будет зависеть только от знака оценки параметров, но будет являться функцией всех независимых переменных. Причем, если оценка параметра будет положительной, то увеличение  переменной xj будет приводить к увеличению вероятности P{Y=1|X}, а если оценка параметра будет отрицательной, то, соответственно, к уменьшению указанной вероятности.  

Замечание. Если фактор х является бинарной переменной, то для него нельзя ввести понятие предельного эффекта.

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

Еще один вопрос для обсуждения: как после оценивания параметров логит (пробит) модели прогнозировать значение у? Поступают, например, следующим образом. Подставляют найденные значения оценок параметров  и значения хj в Z  и вычисляют значение переменной. Если Z>0, то считают, что У=1, если  Z<0, то считают, что У=0. Замечание. Мы рассмотрели ситуацию, когда переменная у была измерена в номинальной шкале, но принимала всего два значения: 0 и 1. В общем случае, когда у может принимать несколько значений, например 0, 1, 2, 3, используют множественный (по у!!) логит или пробит. Кроме того, у может быть измерен в порядковой шкале, тогда в Стате используют  порядковый логит (пробит) ologit (oprobit).

Замечание. Очень часто в исследованиях приходится проводить исследования на усеченной выборке. Например, если изучают доходы домохозяйств, то бывают ситуацию, когда респондентов с очень большим доходом (например, больше 1 млн.рубл.) следует исключить из исследования, то есть

То в таких случаях используют Тобит-модели.

Yi

0

1

Pi

  1.  F(0+1х1+...m-1xm-1)

F(0+1х1+...m-1xm-1)

E[Y]

F(0+1х1+...m-1xm-1)

V[Y]

F(0+1х1+...m-1xm-1) - (F(0+1х1+...m-1xm-1))2 


 

А также другие работы, которые могут Вас заинтересовать

32649. Организация работ и выполнение проекта 29 KB
  Организация работ и выполнение проекта. Задачи управления проектом при его выполнении выполнение сводного плана проекта реализация плана проекта путем выполнения включенных в него работ; подтверждение предметной области процесс формальной приемки предметной области проекта. обеспечение качества процесс регулярной оценки выполнения работ проекта для подтверждения того что проект будет удовлетворять принятым стандартам качества. развитие команды освоение индивидуальных и групповых навыков и квалификации для улучшения выполнения...
32650. Контроль и регулирование хода выполнения проекта 37 KB
  Контроль и регулирование хода выполнения проекта. Контроль Основной целью контроля проекта является обеспечение выполнения плановых показателей и повышение общей эффективности функций планирования и контроля проекта. Содержание контроля проекта состоит в определении результатов деятельности на основе оценки и документирования фактических показателей выполнения работ и сравнения их с плановыми показателями. Система контроля проекта представляет собой часть общей системы управления проектом между элементами подсистемами которой имеются...
32651. Коммуникационная система управления проектом 26 KB
  Коммуникации Управление коммуникациями проекта управление взаимодействием информационными связями управленческая функция направленная на обеспечение своевременного сбора генерации распределения и сохранения необходимой проектной информации. Управление коммуникациями обеспечивает поддержку системы связи взаимодействий между участниками проекта передачу управленческой и отчетной информации направленной на обеспечение достижения целей проекта. Каждый участник проекта должен быть подготовлен к взаимодействию в рамках проекта в...
32652. Этапы закрытия проекта 30 KB
  Этапы закрытия проекта. Этапы закрытия Фаза завершения работ по проекту включает следующие этапы: ввод проекта в действие; достижение проектом заданных результатов; прекращение финансирования проекта; работы по закрытию проекта и внесению изменений не предусмотренных первоначальным замыслом; участие в эксплуатации объектов проекта. На практике завершающую фазу объединяют с фазой ликвидации проекта которая предусматривает вывод объектов проекта из эксплуатации ремонт модернизацию реконструкцию объекта проекта. Ввод проекта в действие...
32653. Анализ эффективности реализации и управления проектом 40 KB
  Анализ эффективности реализации проекта Реализация любого проекта направлена на достижение заранее заданных и четко сформулированных целей и связанных с ними конкретных результатов. В свою очередь планируемый результат является объектом целеполагания проекта. Цели и соответственно результаты проекта могут быть экономическими например получение прибыли снижение себестоимости продукции достижение определенного уровня эффективности работ повышение конкурентоспособности продукции социальноэкономическими создание новых рабочих мест рост...
32654. Место и роль управления проектами в современной теории управления. Сущность проектного и процессного подходов 35 KB
  Роль управления проектами в современном мире. Эти средства активно применяются не только для управления федеральными проектами и программами но и для осуществления управленческих функций внутри правительственного аппарата. В России же где организационная система и методы управления гораздо слабее чем на Западе эффект от внедрения УП окажется еще более значительным.
32655. Основные признаки, понятия и характеристики проекта 30 KB
  Основные признаки понятия и характеристики проекта. Но цели проекта ограничены во времени. При успешном завершении проекта целевая установка данная руководителю проекта заказчиком отпадает. Признак изменения Этот признак является наиболее важной характеристикой проекта так как осуществление проекта всегда несет изменения вещественной системы или предметной области в которой реализуется проект.
32656. Особенности управления программами и портфелями проектов 40.5 KB
  Особенности управления программами и портфелями проектов. При определенных условиях множество взаимосвязанных проектов объединяются в программу. Программа может быть сформулирована в терминах проектов и представлена как совокупность проектов объединенных общей целью выделенными ресурсами временем на ее выполнение технологией организацией и др. В дальнейшем мы будем рассматривать программы как совокупность проектов.
32657. Классификация проектов и программ 28 KB
  в которых осуществляется проект: социальные проекты экономические проекты организационные проекты технические проекты смешанные проекты. Сложные проекты включающие подпроекты и элементы различных типов проектов. Вид проекта: По характеру предметной области проекта: учебнообразовательные проекты; проекты исследования и развития; инновационные проекты.; инвестиционные проекты.