kuangzhu1990 发表于 2011-8-11 12:05
连老师,您好~
看了您写的那篇“股权激励是否真的有效”的文章,受益匪浅,结合stata journal基本把PSM方 ...
那篇文章已经发表:
Lian, Y., Z. Su, Y. Gu, 2011, Evaluating the effects of equity incentives using PSM: Evidence from China, Frontiers of Business Research in China, 5 (2): 266-290.
Lian, Y., Z. Su, Y. Gu, 2011, Evaluating the effects of equity incentives using PSM.pdf
(516.14 KB)
处理过程如下:
*-----------------------------------------------------------
*
* PSM 论文数据处理程序(Given by Lian Yu-jun 连玉君)
*
* Lian, Y., Z. Su, Y. D. Gu, 2011,
* Evaluating the Effects of Equity Incentives Using PSM:
* Evidence from China,
* Frontiers of Business Research in China, 5(2):266-290.
*-----------------------------------------------------------
*=================
* PSM 分析
*=================
* -psmatch2- 命令
global x "size lev prof tobin gpay hhi5"
global xv "_Isicmen* state"
global if56 "if year>=2005"
dropvars jili_gp jili_qq
gen jili_gp = (jili_type==2)
label var jili_gp "1=股票激励,0=其它"
gen jili_qq = (jili_type==1)
label var jili_qq "1=期权激励,0=其它"
*- 比较各变量匹配前后的情况
*- 最近邻匹配
cap drop u
gen u = uniform()
sort u
psmatch2 jili $x $xv $if56, neighbor(3) logit
order id year jili _pscore- _pdif
*- 匹配前后解释变量的差异对比
pstest $x
*----------
*-Figure 01
*----------
/* 需要输出图形时,直接执行里面的命令即可
*- 比较密度函数图
*-匹配前的密度函数图
*-E-version(large)
twoway (kdensity _ps if _treat==1,lp(solid) lw(*2.5)) ///
(kdensity _ps if _treat==0,lp(dash) lw(*2.5)), ///
ytitle("Density", size(*1.1)) ///
ylabel(,angle(0) labsize(*1.1)) ///
xtitle("Propensity Score", size(*1.1)) ///
xscale(titlegap(2)) ///
xlabel(0(0.2)0.8, format(%2.1f) labsize(*1.1)) ///
legend(label(1 "Incentive") label(2 "Control") row(2) ///
position(3) ring(0) size(*1.1)) ///
scheme(s1mono)
graph export "Figs\kn01_large.wmf", ///
replace fontface("Times New Roman")
*-匹配后的密度函数图
*-E-version(large)
twoway (kdensity _ps if _treat==1,lp(solid) lw(*2.5)) ///
(kdensity _ps if _wei!=1 & _wei!=.,lp(dash) lw(*2.5)), ///
ytitle("Density", size(*1.1)) ///
ylabel(,angle(0) labsize(*1.1)) ///
xtitle("Propensity Score", size(*1.1)) ///
xscale(titlegap(2)) ///
xlabel(0(0.2)0.8, format(%2.1f) labsize(*1.1)) ///
legend(label(1 "Incentive") label(2 "Control") row(2) ///
position(3) ring(0) size(*1.1)) ///
scheme(s1mono)
graph export "Figs\kn02_large.wmf", ///
replace fontface("Times New Roman")