logit模型边际效应输出-经管之家官网！

转自：http://www.sohu.com/a/144484890_697896

近期，有亲故问到：如何将logit模型的边际效应用esttab命令输出到Word表格中呢？在之前的推文中，我们介绍了用esttab输出回归结果、相关系数、描述性统计等等，今天再为大家介绍一下如何输出边际效应。

一

logit回归

首先，我们先做个logit回归：

webuse lbw

logit low age lwt i.race smoke ptl ht ui

http://img.mp.sohu.com/upload/20170529/1da816fc32f24a55b11be8f998573438_th.png

二

边际效应

然后计算边际效应，在这里我们将使用命令margins。margins命令可以估计指定协变量的边际效应并将其以表格的形式显示，语法：

margins [marginlist] [if] [in] [weight] [, response_options options]

http://img.mp.sohu.com/upload/20170529/5c517fc7fbd4449ea01de694f35a211b.png

简单栗子：

margins race

http://img.mp.sohu.com/upload/20170529/a31bc394875e40529454870a41e915fe_th.png

margins命令十分强大，但也复杂，参数非常多。

我们可以使用 dydx(varlist) 估计变量的边际效应；eyex(varlist)估计变量的弹性；dyex(varlist)估计变量的半弹性(d(y)/d(lnx))；eydx(varlist)估计变量的半弹性(d(lny)/d(x))。

此外，在求边际效应时，我们可以使用at()选项，指定在某种特定情况下，如指定协变量特定值来求解边际效应。比如

margins, dydx(smoke) at(age=20) //在age等于20时求smoke的边际效应

http://img.mp.sohu.com/upload/20170529/857c2061f4ab41e4acd22a24966f7b95_th.png

over(varlist)则计算在varlist每一值的情况下求解边际效应。如：

margins, dydx(smoke) over(race) //在race的每一值下求smoke的边际效应

http://img.mp.sohu.com/upload/20170529/f34fd88aedb94938a5084170e64ea416_th.png

margins命令还有很多选项，可以根据自己设定的模型和数据添加相应的选项。

三

输出结果

在这里我们使用命令margins,dydx(*)来求上面logit回归的所有变量的边际效应并将边际效应结果输出。

http://img.mp.sohu.com/upload/20170529/639ad687b0ac4aa1b4cc4616f5ac988c_th.png

如果我们想要将边际效应结果输出到Word中，我们需要使用estpost命令。我们可以输出边际效应结果的系数、t值。

estpost margins, dydx(*)

esttab using d:/Marginal.rtf, cell("b(star fmt(3)) t") pr2 ///

replace compress nogap star( * 0.10 ** 0.05 *** 0.01 ) title("Marginal Effect")

http://img.mp.sohu.com/upload/20170529/030aa52ef2824d62afd588ca6d4099a3_th.png

http://img.mp.sohu.com/upload/20170529/a7c6d6733a6e4a8497a72cab14c7683d_th.png

上面只输出了一个边际效应结果，但在实际中我们也许会想将多个边际效应结果放在一个表格里。此时，我们需要结合eststo命令。

clear

set more off

webuse lbw

logit low age lwt i.race smoke ptl ht ui

eststo:estpost margins, dydx(race)

logit low age lwt i.race smoke ptl ht ui

eststo:estpost margins,dydx(race age)

logit low age lwt i.race smoke ptl ht ui

eststo:estpost margins,dydx(*)

esttab est1 est2 est3 using d:/Marginal1.rtf, cell("b(star fmt(3))") pr2 ///

replace compress nogap star( * 0.10 ** 0.05 *** 0.01 ) title("Marginal Effect")

这个过程我们将race、race和age、以及所有变量的边际效应都放在了一个表中。其中est1 、est2、est3 为eststo命令自动储存的估计结果。

http://img.mp.sohu.com/upload/20170529/7c5a63588cde495889a646ac9e728ae3_th.png

以上就是将边际效应输出的全部过程，大家可以试一下并研究一下神奇的margins命令。

以上就是今天给大家分享的内容了，说得好就赏个铜板呗！有钱的捧个钱场，有人的捧个人场~。

应广大粉丝要求，爬虫俱乐部的推文公众号打赏功能可以开发票啦，累计打赏超过1000元我们即可给您开具发票，发票类别为“咨询费”。用心做事，只为做您更贴心的小爬虫。第一批发票已经寄到各位小主的手中，大家快来给小爬虫打赏呀~

「经管之家」APP：经管人学习、答疑、交友，就上经管之家！
免流量费下载资料----在经管之家app可以下载论坛上的所有资源，并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库，各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里，你可以找到任何学科方向、有共同话题的朋友。
经管之家（原人大经济论坛），跨越高校的围墙，带你走进经管知识的新世界。
扫描下方二维码下载并注册APP