第9讲二元结果模型
- 格式:docx
- 大小:10.72 KB
- 文档页数:4
第9讲二元结果模型
第9讲离散选择模型之二元结果模型
参考书目:
1・ Long, J・ S・,and J・ Freese・ 2021・ Regression Models for Categorical
Dependent Variables Using Stata・ 2nd ed・ CollegeStation, TX: Stata Press
教学视频:
Logistic regression, part 1: Binary predictors Logistic regression, part 2:
Continuous predictors Logistic regression, part 3: Factor variables
一、离散被解释变量的例子
二元结果模型:考研或不考研:就业或待业:买房或不买房;买保险或不买保险;贷 款申请被批准或拒绝:出国或不出国;回国或不回国;战争或和半;医药实验中的生或死。
多元结果模型:对不同交通方式的选择(走路、骑车、坐车上班):对不同职业的选择。
这类模型被称为“离散选择模型” (discrete choice model) o考虑到离散彼解释变
量的特点,通常不宜用OLS进行回归。假设个体只有两种选择,比如y二1 (考研)或y二0
(不考研)。是否考研,取决丁•研究生毕业后的预期收入、个人兴趣、本科毕业后直接就业
的收入前景等。所有解释变量都包括在向量x中。
1
二、二元结果模型的微观基础
对于二元选择行为,可通过“潜变量” (latent variable)概括该行为的净收益(收益
减去成本)。
如果净收益大于0,则选择做;否则,选择不做。
y*=x' B + £
其中,净收益y*为潜变量,不可观测。选择规则为y二1,若y*>0 y二0,若y*W0
如果£为正态分布,则为Probit:如果£为逻辑分布,则为Logit。 logistic — Logistic regression, reporting odds ratios (Logistic 回归,扌畏告
优势比/比值比)
对 Logit 模型,记 p二 P(y 二1 x ),则 1~P= P(y 二0 x )。
P /(1-P)称为"几率比/优势比/比值比"(odds ratio) o那么什么是几率比?
2
举例说明,假设在检验药物疗效的随机实验中,“y =1 ”表示“生”,“y二0 ”表 示“死”;则几率比为2意味着存活的概率是死亡概率的两倍。
命令语法:
logisticdepvarindepvars [if] [in] Zweight] [,options]
命令logistic拟合logistic回归模型,其中depvar是一个0/1变量(更准确地说, 是一个0/非0变量)。logistic显示的估计结果是优势比:要想查看系数,运行 logistic 后,输入 logito 一、概述
logistic替代logit命令,通常是拟合最大似然logit模型的首选方法。Stata使用
logit和logistic意味着同样的事情:最大似然估计。并J1两种命令会得到相同的结果。
logistic命令通常优于logit命令,原因是logistic提供优势比而不是系数的估计结果。
对 T" logistic 回归的介绍,请参阅 Lemeshow&Hosmer (2021). Pagano&Gauvreau
(2000, 47O0C487),或Pampel (2000): 一个完整的,非数学的讨论,请参阅 Kleinbaum& Klein
(2021); 一个更深入的讨论,请参阅 Hosmer, Lemeshow, & Sturdivant (2021)。Gould
(2000)讨论了关 J' logistic 回归的解释。Dupont (2021)和 Hilbe (2021)用 Stata 实例讨论T logistic 回归。Vittinghoff et al. (2021)M点讨论 了模型设定问题。
Stata有一系列命令用来估计二分类或多分类因变量模型。Long和Freese (2021)专 门论述了如何使用Stata拟合此类模型。下面是一些常用估计命令的列表。help estimation
commands提供了 Stata全部估计命令的完整列表。
asclogit asmprobit asroprobit binreg
[R] [R] [R] [R]麦克法登离散选择模型(McFadden' s choice)多项式probit模型 排序probit模型二项式的广义线性模型3biprobit blogit bprobit clogit cloglog exlogistic glm giogit gprobit
heckoprobit heckprobit hetprobit ivprobit logit mecloglog meglm melogit meprobit
mlogit mprobit nlogit ologit oprobit probit rologit scobit slogit svy: cmd xtcloglog
xtgee xtlogit xtologit xtoprobit xtprobit
[R] [R1 [R] [R1 [R1 [R1 [R] [R] [R〕[R] [R] [R] [R] [R] [ME] [ME] [ME] [ME] [R]
[R] [R] [R] [R] [R] [R] [R] [R] [SVY] svy estimation [XT] [XT] [XT] [XT] [XT] [XT]二元probit模型分组数据的logit模型分组数据的probit模型条件(固定 效应)logistic回归互补重对数模型精确logistic回归广义线性模型分组数据的加权 最小二乘logistic回归分组数据的加权最小二乘probit回归有序probit的样本选择 模型probit的样本选择模型异方差probit模型具有连续内生变量的probit模型 Logistic回归分析,报告系数多层次混合效应互补双对数回归多层次混合效应广义线性 模型多层次混合效应logistic回归多层次混合效应probit回归多分类因变量 logistic回归多分类因变量probit回归嵌套logit回归有序logistic回归有序 probit 回归 Probit 回归有序秩 logistic 回归 Skewed logistic
regression Stereotype logistic regression命令的svy版随机效应和总体平均cloglog模型GEE 总体平均广义线性模型固定效应、随机效应和总体平均logit模型随机效应有序
logistic模型随机效应有序probit模型随机效应和总体平均probit模型例:估计决 定美国妇女就业与否的二元结果模型。数据集包括以下变量:虚拟变^work(l=就业), age(年龄),虚拟变量married(1=已婚),children(子女数),
4
educat ion (受教育年限)。考虑以下模型:
worki= P 0+ P lagei+ B 2marriedi+ P 3childreni+ P 4educationi+ £ i 作为对照,首
先使用OLS进行线性概率模型(LPM)估计:use womenwkl, clear (原数据是
womenwk. dta) reg work age married children education probit work age married
children education, nolog
mfx (计算probit模型在样本均值处的边际效应,与OLS估计的回归系数进行比较)
estat classification (计算预测准确的百分比)logit work age married children
education, nolog mfx
estat classification
hetprob work age married children education, het(age married children
education)nolog (p值为0.78,所以接受“同方差”的原假设。)generate age2=age*age generateagemari=age^married
generateagechr=age*children
quietlylogit work age married children education age2 agemariagechr test age2
agemariagechr (接受零假设)
quietlylogit work age married children education estimates store blogit
quietlyprobit work age married children education estimates store bprobit
quietlyregress work age married children education estimates store bols
quietlylogit work age married children education, vce(robust) estimates
store blogitr
quietlyprobit work age married children education, vce (robust)
5
感谢您的阅读,祝您生活愉快。