"Selmlog" 命令：multinomial logit

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

大家好：
     我想用“Selmlog” 命令分析multinomial logit 回归。
     因为selmlog 不是stata 自带程序，必须使用用户撰写的ado 文件。我安装这个网站的说明将下载的ado 文件放在指定位置 （http://www.parisschoolofeconomic ... mlog/selmlog13.html）。本命令源于文章Selection Bias Corrections Based on the Multinomial Logit Model: Monte-Carlo Comparisons", François Bourguignon, Martin Fournier, Marc Gurgand, Journal of Economic Surveys, vol. 21, n°1, 2007.

    我尝试了如下回归：
selmlog 收入 X1 X2, selec (选择变量=x1 x2 z1 z2； 其中（选择变量取值1， 2， 3）

    但是显示如下错误：Method option not specified variable _m* not found  r(111);


    有人能告诉我这是什么原因吗？谢谢。

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏2 回帖

关键词：Multinomial nomial logit multi Log Journal 文章 程序 网站 用户

你要加上下面这些，dmf(2)是correction的方法，可以改别的，参见程序帮助文件
dmf(2) bootstrap(100) mloptions(rrr level (95)) gen(rho_1)

xiaofengshen 发表于 2014-8-11 22:07
你要加上下面这些，dmf(2)是correction的方法，可以改别的，参见程序帮助文件
dmf(2) bootstrap(100) mlop ...

好的，谢谢。现在可以了。
最初是我忘记设置缺失变量了。
对了，我看有些文章当中报告Adjusted R-square 值或者F（n, n-k）的值，这个是怎么弄出来的。

我想请问一下multinomial logit model的命令怎么写啊

xiaofengshen 发表于 2014-8-11 22:07
你要加上下面这些，dmf(2)是correction的方法，可以改别的，参见程序帮助文件
dmf(2) bootstrap(100) mlop ...

请问下，加了之后运行出现Non-missing values of dependent variable for more than one outcome是为什么？

linzima 发表于 2014-8-20 15:46
好的，谢谢。现在可以了。
最初是我忘记设置缺失变量了。
对了，我看有些文章当中报告Adjusted R-squar ...

您好，最近在学用 selmlog，出现与您之前之问题：

1. . selmlog ca ca0, sel(err = lgdp openness pcgdpg resimp inf lgov totg_sd fd kaopen) dmf(0)
   
2. Non-missing values of dependent variable for more than one outcome
   
3. variable _m* not found
   
4. r(111);

复制代码

后来，您谈到您"最初是我忘记设置缺失变量了"，我猜这里是关键之处，请问您是怎么处理？谢谢！

同问，缺失变量怎么设置？

黃河泉 发表于 2017-8-19 17:28
您好，最近在学用 selmlog，出现与您之前之问题：后来，您谈到您"最初是我忘记设置缺失变量了"，我猜这里 ...

您好，请问这个问题你解决了吗，拜托能讲讲如何解决这个问题的吗？

zzcwan 发表于 2018-1-7 19:58
您好，请问这个问题你解决了吗，拜托能讲讲如何解决这个问题的吗？

我试过的结果是类似

1. // full sample
   
2. preserve
   
3. replace ca = . if err == 2 | err == 3
   
4. selmlog ca ca0, sel(err = lgdp kaopen fd lgov) lee boot(1000) showmlogit mloptions(vce(robust))
   
5. *est store m1
   
6. outreg2 using "selmlog", excel ctitle(R1) dec(4) replace
   
7. restore
   
8. 
   
9. preserve
   
10. replace ca = . if err == 1 | err == 3
    
11. selmlog ca ca0, sel(err = lgdp kaopen fd lgov) lee boot(1000) showmlogit mloptions(vce(robust))
    
12. *est store m2
    
13. outreg2 using "selmlog", excel ctitle(R2) dec(4) append
    
14. restore
    
15. 
    
16. preserve
    
17. replace ca = . if err == 1 | err == 2
    
18. selmlog ca ca0, sel(err = lgdp kaopen fd lgov) lee boot(1000) showmlogit mloptions(vce(robust) )
    
19. *est store m3
    
20. outreg2 using "selmlog", excel ctitle(R3) dec(4) append
    
21. restore

复制代码

我的 err 有三个状况 (err=1,2,3)，当你要估计 err=1 之结果，关键的地方是你必须令 err=2, err=3 的被解释变量 ca 为 missing values。

黃河泉 发表于 2018-1-8 07:06
我试过的结果是类似我的 err 有三个状况 (err=1,2,3)，当你要估计 err=1 之结果，关键的地方是你必须令 e ...