👇 连享会 · 推文导航 | www.lianxh.cn
黄熹 (汕头大学)邮箱： 1085300826@qq.com
目录
=1. 离散模型简介
2. Logit：二值选择模型
2.1 Logit 之建模过程
2.2 Logit 之参数解读
2.3 Logit 之 Stata 应用
3. mlogit 模型
3.1 mlogit 之建模过程
3.2 mlogit 之参数解读
3.3 mlogit 之 Stata 应用
4. 小结
5. 参考资料
6. 相关推文
温馨提示： 文中链接在微信中无法生效。请点击底部「阅读原文」 。
1. 离散模型简介 被解释变量为离散变量的模型被称为离散选择模型。若作进一步区分则有：
二值选择模型 ：模型的被解释变量是二值变量 (是否读博，是否出国留学等)；多值选择模型 ：模型的被解释变量是多值离散变量 (企业多元化战略的选择，资产组风险资产的数量选择等)。针对离散选择模型，考虑到被解释变量的连续性性态得不到满足，故常常不用 OLS 进行回归 (拿到一笔数据，迫不及待地想初步一览相关性的情况除外)。针对离散选择模型，实证工作中可根据被解释变量的分布函数，依情况采用 Logit, Probit, mlogit, mprobit, ologit, oprobit, negatibe binomial, zero-inflated negative binomial, poisson, zero-inflated poisson 等模型进行回归。
由于 Probit 模型的解释不如 Logit 模型直观，因此相对 Logit 模型而言，其使用范围和频数还是相当有限的。Logit 模型在实证研究中已被广泛使用，且其理论与实际的契合度也受到了学界的高度认可。
虽然本文的出发点是介绍文献中出现频数亦日渐高企的 mlogit 模型，但是鉴于不少读者对 Logit 模型的解读存在偏差，且 Logit 模型是 mlogit 模型的源起，因此在本文中，我力求以精简的语言、从 Stata 应用与模型解读的角度出发，将先行介绍 mlogit 模型的鼻祖—— Logit 模型，而后再对 mlogit 模型进行阐释。
本文将 Logit 和 mlogit 统称为 小罗 。
2. Logit：二值选择模型 2.1 Logit 之建模过程 将企业仅有的两种战略选择作为因变量，y=1 (并购) 和 y=0 (不并购)。现金流充裕程度、财务杠杆、第一大股东持股比例、行业成长性等因素都会影响企业的并购决策，将这一系列解释变量记为向量
由于
零条件均值满足时，有 β β
该函数被称为连接函数 (link function) ，其选择具有很大的灵活性，直接目的便是：在限定因变量在 [0,1] 区间的前提下，更合理地将因变量与自变量之间的非线性关系刻画出来。此时，如果选取的分布函数
, β β
对于这种非线性模型，可先求出该分布的概率密度和似然函数，进而利用最大似然估计法 (MLE) 和数值计算可得到 β
2.2 Logit 之参数解读 从 2.1 我们知道， Logit 模型的边际效应就是逻辑分布函数对
可见，OLS 下的回归系数在这里已然变味。该模型的边际效应是一个随解释变量而变的变量，β 。Stata 中默认使用平均边际效应 (在分别计算每个样本观测值上的边际效应之后进行简单算术平均) 来计算。那么 β
为了便于说明，记
通过对数转换可得
β β
即 β
2.3 Logit 之 Stata 应用 本例采用 Stata16 自带的数据集进行实操。这是一份 1988 年美国妇女工资的资料，包含 2246 笔观察值，共 17 个变量。本例将虚拟变量是否已婚（married）作为因变量，将连续变量年龄（age）、工资（wage）、工作时数（hours）、受教育年限（grade） 以及类别变量种族（race）、职业（occupation）这 6 个变量作为自变量。
Logit 回归假设数据服从伯努利分布，并运用梯度下降法来求解参数，因此 Stata 默认汇报迭代 (Iteration) 过程。如果附加选项 nolog ，则表示不要再汇报迭代过程了，结果如下表所示。
sysuse nlsw88, clear //调入1988年美国妇女工资的资料以上结果虽然汇报的是回归系数而非几率比，但是，以对数几率比的角度来进行模型解读将会更加清晰和严谨。应时刻记住的是：这些回归系数并不是边际效应 。
分析如下：工资每增加 1 美元，对数几率比 结 婚 概 率 未 婚 概 率
结 婚 概 率 未 婚 概 率
其实，对解释变量的系数可以有一个更为直观的解释：即，如果系数为负，则表示被解释变量取1的可能性越小。如：wage 越大，married 取 1 的可能性越小，即工资越高，越不可能结婚。
同时，Logit 模型也可能出现异方差问题。在执行上述命令之后附加选项 r ，则 Stata 会自动计算稳健标准误。由以下结果可知，稳健标准误与普通标准误十分接近，因此可大致判断不存在模型设定偏误的问题。
. logit married age wage hours grade race occupation, r nolog做完基本的 Logit 回归之后，我们可以输入 margins 命令并附加选项 dydx() 来得到平均边际效应。与 OLS 回归结果相比可知，Logit 回归的平均边际效应与 OLS 回归得到的系数估计值非常接近 (见下表)。
. logit married age wage hours grade race occupation, nolog以下命令附加了选项 or ，表示汇报几率比 (的变化倍数) 而非系数。
. logit married age wage hours grade race occupation, or简单来看，几率比越大，表示结婚概率越大 (不考虑显著性水平的情况下)。具体来说，在控制其他变量不变的情况下：年龄每增加一岁，妇女结婚的几率比放大 0.978 倍；工资每上涨一美元，妇女结婚的几率比放大 0.986 倍；依此类推。
3. mlogit 模型 3.1 mlogit 之建模过程 个人在多于两个选项中进行选择，以使效用最大化。
多项选择的应用包括交通工具的选择、职业的选择、参议员的选择以及竞标企业的选择，不一而足。值得注意的是，上述这些选择集合都是无序的，因此 mlogit 模型其实是无序选项集的模型 (与此相反，如果是针对基数效用、QS 世界大学排名、债券信用评级等有序量表，则应考虑有序选择的 ologit 模型)。
多元 Logit 模型可视为对被解释变量中各类选择行为两两配对后构成的多个二元 Logit 模型实施联合估计 (两两配对，形成多个二元 Logit 模型)。
mlogit 模型可以随机效用模型作为引子进行切入。假设一个消费者
如果消费者选择了 ，进而假设消费者选择行为满足显示偏好强公理，则可知 是所有选择中效用最大的。效用可拆分为个体特征 (如性别、收入、种族、受教育程度、职业经理人性格、公司文化等) 和选择两个方面，此时再假设 个扰动项独立同冈贝尔 (Gumbel medel，不具对称性) 分布，则可通过求解这个
维的方程组得到每种选择的预测概率： 需要说明的是，二项式 Logit 模型为 ，
β 一种简化处理的方式是令基准组的系数为 β ( 注意：只有完成了这一标准化处理的假定，方能推得如下模型；推导可见 3.2 mlogit 之参数解读 的式
可见，当 β
易知，当我们将第一组选定为基准组时，
3.2 mlogit 之参数解读 Logit 模型估计的关键在于选定基准组，所有系数均是相对于基准组来进行估计的 。与 Logit 模型类似，本文从胜算比的角度出发来对模型参数进行阐释。假设基准组为第
当基准组为第一组时，第
β
β β
依此类推，当基准组为第
β β β β
本例采用 Long and Feese (2006) 提供的数据集进行实操。此数据集来自 1982 年的一项调查问卷，问卷将受访者分为服务人员 (menial)、蓝领 (blue collar)、工匠 (craft)、白领 (white collar) 和专业人士 (professional) 共 5 类；解释变量包括虚拟变量是否白人 white、连续变量受教育年限 ed 以及工龄 exper ，其均只依赖于个体，不随组别而变。
. use "nomocc2.dta", clear3.3.1 mlogit 回归与模型解读 . mlogit occ ed exper white, nolog //进行多元 Logit 回归由上表可知：严格来说，受教育年限每增加一年，服务业、蓝领、工匠、白领这四个组别相对于专业人士组的胜算比将分别放大 、
、
直觉上来说，给定其他变量，受教育年限越长，越不可能选择除专业人士以外的 4 种职业，且均高度显著；在 5% 的显著性水平上，工龄越长，越不可能选择服务业或蓝领；在 5% 的显著性水平上，白人更不可能选择服务业或工匠；依此类推。
3.3.2 基准组的设置 . mlogit occ ed exper white, base(2) nolog 由上表可知，系数估计值会随着参照组的不同选择而变化。
3.3.3 进一步解释 . mlogit occ ed exper white // qui 表示让 Stata 闭嘴，不要把结果呈现在屏幕上其中，factor change 表示胜算比变化，即
β ；第一列列示了左侧选择相对于右侧选择的胜算比，第五列计算出的则是 factor change。
3.3.4 拟合优度 (信息准则) . mlogit occ ed exper whitefitstat 的结果主要是采用信息准则来评判拟合优度 (CI 越小则表明模型拟合得越好，由于 BIC 的惩罚力度较大，八成以上的文献都是选择 BIC；事实上，BIC 通常也确实优于 AIC)。由上表汇报的 CI 指标可初步判断，模型总体上来说拟合得很好。
3.3.5 检验 IIA (仅提供命令) IIA 假定 (无关选择独立性假设) 是 mlogit 模型的前提。因此从理论的角度来说，使用 mlogit 模型之前必须检验该假定是否得到满足。除进行常规的 Hausman 检验之外，还可以附加 Small - Hsiao 检验。
*-Hausman test of IIA3.3.6 顺道唠唠 mprobit mlogit 模型和 mprobit 模型实际上来说并无多大差别，只是 mprobit 需要计算多重积分，计算耗时更长，而且无法从胜算比的角度来直观地解释模型，因此在实证分析中学者多采用 mlogit 模型。且应注意：mlogit 和 mprobit 的系数估计值并不具有可比性，真正可比的是两个模型的预测概率值。
4.小结 对于因变量是离散变量的情况，之所以不使用 OLS 而使用小罗 ( Logit 或 mlogit ) ，是因为线性概率模型不考虑曲率，线性最小二乘法产生的斜率将取决于大部分数据在曲线上的位置，而小罗能更好地描述数据的分布情况，且将预测值的范围锁定在 [0,1] 之间；
进一步，小罗利用了对数转换，以几率比 (事件发生概率和不发生概率的比值) 代替概率值，站在“几率比”这一新的视角来对模型结果进行解读，促进了乘法比较，简单清晰地描述出了模型中原本复杂的非线性关系 (如 S 形关系) ，增加了可解释性。综上，小罗因其特有的优势和基本无害的计量方法，在学界和实务界 (如券商、医疗等) 备受推崇。
但需要注意的是，mlogit 模型采用 MLE 进行估计，且具有联立方程的特征，因此对大样本的要求比 Logit 模型更高 (当然，小罗回归均是采用 MLE 来进行估计的)。
5. 参考资料 温馨提示： 文中链接在微信中无法生效。请点击底部「阅读原文」 。
陈强. 高级计量经济学及 Stata 应用[M]. 高等教育出版社，2014 Whilliam H. Greene. 计量经济分析[M]. 中国人民大学出版社，2020 连享会 推文：多元 Logit 模型详解 (mlogit) 连享会 推文：为何使用 Logit 模型？Logit 模型评价 连享会专题推文：Probit - Logit 系列 J. Scott Long & Jeremy Freese：Regression Models For Categorical Dependent Variables Using Stata
6. 相关推文 温馨提示： 文中链接在微信中无法生效。请点击底部「阅读原文」 。
Note：产生如下推文列表的命令为：lianxh logit probit, mlianxh 命令：ssc install lianxh, replace
Stata新命令：面板-LogitFE-ProbitFE reg2logit：用OLS估计Logit模型参数 Stata：多元 Logit 模型详解 (mlogit) 详解 Logit/Probit 模型中的 completely determined 问题 Stata：嵌套 Logit 模型 (Nested Logit) 
