PSM,RDD,Heckman,Panel模型的操作程序

凡是搞计量经济的，都关注这个号了

邮箱：[e-mail]

所有计量经济圈方法论丛的code程序, 宏微观数据库和各种软件都放在社群里.欢迎到计量经济圈社群交流访问.

计量经济圈经济社会等数据库合集, 社科研究者的大米(继续)

今天，我们“因果推断研究小组”将为计量经济圈的圈友引荐PSM, RDD，Heckman和Panel模型。关于这些模型的一些介绍，各位圈友可以参阅因果推断研究小组系列推文。

这篇文章里的数据与分位数回归, Oaxaca分解, Quaids模型, 非参数估计程序的数据和整个的do file都放在网盘里，提取方式参见文章后面说明。

**heckman two-step procedure----------------

summarize

reg whz h1 h2 h3 h4 c2

probit si i1 i2 c1 h1

predict xb 

generate invmills = normalden(xb)/normal(xb)

set seed 10100

bootstrap, rep(500): reg whz h1 h2 h3 h4 c2 invmills  //用手动计算二步heckman不一定准确

heckman whz h1 h2 h3 h4 c2, select(si = i1 i2 h1 c1) nolog

heckman whz h1 h2 h3 h4 c2, select(si = i1 i2 h1 c1) twostep //一次性完成

**PSM---------------------------

pscore d treat income, pscore(ps1) blockid(blockf1) comsup level(0.001) //propensity score matching

atts score d treat income, pscore(ps1) blockid(blockf1) comsup  //stratification matching

attr score d treat income, pscore(ps1) radius(0.001) comsup //radius matching

attk score d treat  income, pscore(ps1) comsup bootstrap reps(50) //kernal matching

search nnmatch //安装nnmatch

nnmatch score d treat income, tc(att) m(1)  //nearest neighbour matching

**RD------------------------

rdplot score cutoff, graph_options(title(RD Plot - Score-Meals Data)) ytitle(Score) xtitle(cutoff)  //以下四个为画图

rdplot score cutoff,binselect(es) graph_options(title(RD Plot - Score-Meals Data)) ytitle(Score) xtitle(cutoff)

rdplot score cutoff,binselect(qsnv) graph_options(title(RD Plot - Score-Meals Data)) ytitle(Score) xtitle(cutoff)

rdplot score cutoff, binselect(es) ci(95) graph_options(title("RD Plot: Scores Meals") ytitle(scores) xtitle(cutoff) graphregion(color(white)))

rdrobust score cutoff  //rdrobust估计

rdrobust score cutoff, all

rdbwselect score cutoff, all  //选择最优带宽

**panel data operation-------------------------

xtreg bmi aidmpc aidfpc fdid lnliv hhsize lncons frac_female frac_male lost_work ///

days_labor pa2_1 pa3_1 pa8_1 pa10_1 pa14_1 pa2_2 pa3_2 pa8_2 pa10_2 pa14_2 pa5_1 ///

pa5_2 if samp1==1 & male==1, fe cluster(hhid)  //聚类固定效应面板回归

test aidmpc=aidfpc   //检验两个系数是不是相等

global xlist "aidmpc aidfpc fdid lnliv lact_bre hhsize lncons frac_female frac_male lost_work days_labor pa2_1 pa3_1 pa8_1 pa10_1 pa14_1 pa2_2 pa3_2 pa8_2 pa10_2 pa14_2 pa5_1 pa5_2"

quietly regress bmi $xlist, vce(cluster hhid)  //聚类标准误

estimates store OLS_cluster

quietly xtreg bmi $xlist, be  //组间估计方法

estimates store BE

 xtreg bmi $xlist, fe  //固定效应估计方法

estimates store FE

quietly xtreg bmi $xlist, fe vce(robust)  //固定效应估计，稳健标准误

estimate store FE_rob

quietly xtreg bmi $xlist, re  //随机效应估计方法

estimates store RE

quietly xtreg bmi $xlist, re vce(robust)  //随机效应估计，稳健标准误

estimates store RE_rob 

estimates table OLS_cluster BE FE FE_rob RE RE_rob, b se stats(N r2 r2_o r2_b r2_w /// 

sigma_u sigma_e rho) b(%7.4f)  //把这些回归结果弄到一个表格

hausman FE RE, sigmamore  //hausman检验，即看看固定效应还是随机效应好

xtabond bmi aidmpc aidfpc fdid lnliv hhsize lncons frac_female frac_male lost_work ///

days_labor pa2_1 pa3_1 pa8_1 pa10_1 pa14_1 pa2_2 pa3_2 pa8_2 pa10_2 pa14_2 pa5_1  ///

pa5_2 if samp1==1 & male==1, lags(1) artests(1) nocons twostep vce(robust) //动态面板

需要完整版本数据和程序的，直接在公众号菜单栏回复“PSM,RDD,Heckman,Panel模型的操作程序”即可获得网盘链接（注：建议直接复制红色字再粘贴）

各位非社群群友可以参与咱们社群组建的一系列研讨会，只要也能汇报一篇与经典计量方法相关的顶刊论文。有意向的学者，可以就具体细节直接联系公共邮箱：[e-mail]

下面这些短链接文章属于合集，可以收藏起来阅读，不然以后都找不到了。

2年，计量经济圈公众号近1000篇文章，

Econometrics Circle

数据系列：空间矩阵 | 工企数据 | PM2.5 | 市场化指数 | CO2数据 |  夜间灯光 | 官员方言  | 微观数据 |

计量系列：匹配方法 | 内生性 | 工具变量 | DID | 面板数据 | 常用TOOL | 中介调节 | 时间序列 | RDD断点 | 合成控制 | 

数据处理：Stata | R | Python | 缺失值 | CHIP/ CHNS/CHARLS/CFPS/CGSS等 |

干货系列：能源环境 | 效率研究 | 空间计量 | 国际经贸 | 计量软件 | 商科研究 | 机器学习 | SSCI | CSSCI | SSCI查询 |

计量经济圈组织了一个计量社群，有如下特征：热情互助最多、前沿趋势最多、社科资料最多、社科数据最多、科研牛人最多、海外名校最多。因此，建议积极进取和有强烈研习激情的中青年学者到社群交流探讨，始终坚信优秀是通过感染优秀而互相成就彼此的。