互助问答第297期：关于预测反米尔斯比率的问题
关于预测反米尔斯比率的问题
麻烦问一下各位老师同学，做了似不相关双边probit模型后，怎样预测对应的反米尔斯比率，目的是把两个方程的反米尔斯放到第二阶段的回归中，希望各位老师同学解惑，谢谢大家！
//help biprobit
//Examples
//webuse school ,giving up for firewall
cd D:\biprobit_heckprobit_mill
use school,clear
//Bivariate probit regression
biprobit private vote logptax loginc years
//Seemingly unrelated bivariate probit regression
biprobit (private = logptax loginc years) (vote = logptax years)
//Seemingly unrelated bivariate probit regression with robust standard errors
biprobit (private = logptax loginc years) (vote = logptax years), vce(robust)
****************************help heckprobit******************************************
//[R] heckprobit -- Probit model with sample selection
//(View complete PDF manual entry)
//Fit a probit model with sample selection
heckprobit private years logptax, sel(vote=years loginc logptax)
//Replay results, but display legend of coefficients rather than the statistics for the coefficients
heckprobit, coeflegend
//stata中heckprobit命令的查询可知rho即为实际的逆米尔兹比
//https://bbs.pinggu.org/thread-2660588-1-1.html
****probit 提取inverse mills ratio
probit private vote logptax loginc years,nolog // possibility of being open region predict gw, xb g lambda=normalden(gw)/normal(gw) if private==1
replace lambda=-normalden(gw)/normal(-gw) if private==0 // generate inverse mills ratio (selection model)
//怎么在stata中计算mills ratio
probit private vote logptax loginc years,nolog
predict gw, xb
replace lambda=normalden(gw,0,1)/normal(gw) if private==1
replace lambda=-normalden(gw,0,1)/normal(-gw) if private==0
*y是0-1变量，x1-xn是解释变量，lambda即IMR
//http://bbs.pinggu.org/thread-1087688-1-1.html?winzoom=1
往期回顾：
互助问答第296期：关于PSM-DID的问题
如果您在计量学习和实证研究中遇到问题，请及时发到邮箱szlw58@126.com，专业委员会有30多名编辑都会看，您的问题会得到及时关注！请您将问题描述清楚，任何有助于把问题描述清楚的细节都能使我们更方便地回答您的问题，提问细则参见：实证研究互助平台最新通知（点击文末阅读原文查看详情）
鲜活的事例更有助于提高您的研究水平，呆板的教科书让人生厌。如果您喜欢，请提出您的问题，也请转发推广！
（欢迎转发，欢迎分享；转载请注明出处，引用和合作请留言。本文作者拥有所有版权，原创文章最早发表于“学术苑”。任何侵权行为将面临追责！）
学术指导：张晓峒老师 Ben Lambert
本期解答人：谢杰老师
编辑：陈波
统筹：左川 李丹丹
技术：刘子瑗
全文完，感谢您的耐心阅读
请顺手点个“在看”吧~
阅读原文
