mlogit模型或者ologit模型结果RRR的输出

    在使用mlogit模型或者ologit模型分析问题时, 模型中的系数 \hat{ \beta } 并不是边际影响,因此其符号、大小并不能直接反映x y 的边际影响,甚至\hat{\beta } 的显著性也不一定能说明边际影响必然显著,因此在进行多值选择模型或者定序选择模型结果汇报时,系数的参考价值有限。一种办法是:不报告\hat{\beta } ,直接汇报边际影响的大小\frac{dy}{dx} ;另一种办法是不报告\hat{\beta } ,而是汇报相对风险比的变化倍数RRR,也就是e^\hat{\beta }  的值,尤其是当x 的取值只能变化一单位(比如性别、婚姻等虚拟变量)时,则应汇报e^\hat{\beta }  。第三种办法是同时汇报\hat{\beta } \frac{dy}{dx}  e^\hat{\beta }  。下面我们mlogit模型为例进行演示。

cd "E:\数据文件夹\ "

webuse sysdsn1,clear

*********基准模型:ols及结果输出***********

reg insure age male nonwhite i.site

  est store ols

  esttab ols

  outreg2 [ols] using table1,word replace tstat bdec(3) tdec(2) rdec(3) e(F)


*********模型mlogit:追加系数*********** 

mlogit insure age male nonwhite i.site

  est store  coef

  esttab  coef //此时模型m1中汇报的是系数b 

  outreg2 [coef] using table1,word append tstat bdec(3) tdec(2) rdec(3) e(r2_p)



*********模型mlogit:追加RRR***********   

mlogit, rrr //Replay, reporting relative-risk ratios

   est store RRR

    esttab  RRR

   outreg2 [RRR] using table1, stnum(replace coef=exp(coef), replace se=coef*se) ///

                      cti(odds ratio) ///

                      word append tstat bdec(3) tdec(2) rdec(3) e(r2_p)  nonote                                                           


*********输出dy/dx:此时outreg2将变量按照列输出结果,不够好用*******

********边际效果的输出,只要在后面附加post命令就可以**********

**********然后就可以用est命令把模型输出到CSV中,进行手动粘贴复制******

margins,dydx(_all) post  //边际效应

    est store dydx

    esttab dydx

    esttab dydx using table1.csv,replace compress nogap nonotes b(3) t(2) pr2(3)         


           ***********outreg2输出边际系数的效果如下

           outreg2 [dydx] using tab113, see  tdec(2) rdec(3) e(F)         

           outreg2 [dydx] using table1, word append tstat bdec(3) tdec(2) rdec(3) e(r2_p)                                                             


****est restore 激活模型,否则想要看看跑过的模型的其他指标,还得重新再跑一遍十分浪费时间

  est restore  coef  //对coef这一模型进行激活,省的再跑一遍,浪费时间。

  listcoef, help  //汇报系数b、 e^b、e^bStdX                

  listcoef, help factor pvalue(0.10) positive



**********est+ eq(1)强制把所有模型结果并列在一起*******************

est table ols coef RRR dydx, star(.1 .05 .01) eq(1)

                //eq(1)强制把所有结果并列在一起输出在一个表格里,无论模型变量是否一致 


***************************结束************************************************

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 217,277评论 6 503
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 92,689评论 3 393
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 163,624评论 0 353
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 58,356评论 1 293
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 67,402评论 6 392
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 51,292评论 1 301
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 40,135评论 3 418
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,992评论 0 275
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 45,429评论 1 314
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,636评论 3 334
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,785评论 1 348
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,492评论 5 345
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 41,092评论 3 328
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,723评论 0 22
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,858评论 1 269
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 47,891评论 2 370
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,713评论 2 354

推荐阅读更多精彩内容