有几处代码小bug解释一下,1.drop if year==. //去掉缺失值,要多打一个=号,去除空缺值;
2.
rename u_hat Neg // u代表了抑制效应(或称为技术效率损失)
rename w_hat Pos //W或V则代表了促进效应(或称为技术效率增益)
rename uw_diff_exp NI //净效应
replace NI=Pos-Neg //可能算的有差异可以备用
3.添加完善
*-转换成百分数,便于统计
foreach v of varlist Neg-NI{
replace `v' = `v'*100
}
tabstat Neg -NI, stat(mean sd p25 p50 p75) format(%6.3f) c(s) //算出U,W的双边效应估计拟合值
4.根据分组要求重新把bug修订为 **时间特征**和**地区特征**,
logout, save(hlh_z0) excel replace: /// //输出 Excel 表格
tabstat Pos Neg uw_diff, ////
by(zone) stat(mean sd p25 p50 p75) format(%6.3f) c(s) //分区域统计
**时间特征**
*-w_hat u_hat uw_diff (Panel 1-2)
logout, save(hlh_z1) excel replace: /// //输出 Excel 表格
tabstat w_hat u_hat uw_diff, ////
by(year) stat(mean sd p25 p50 p75) format(%6.3f) c(s) //促进效应、抑制效应时间特征
*-NI (Panel 3-4)
logout, save(hlh_z2) excel replace: /// //输出 Excel 表格
tabstat NI, ////
by(year) stat(mean sd p25 p50 p75) format(%6.3f) c(s) //净效应时间特征
|