社区贡献的命令不coefplot
应该这样使用。
不过,下面是一个玩具示例,说明如何获得所需的东西:
sysuse auto, clear
estimates clear
// code for identification strategy 1
regress price mpg trunk length turn if foreign == 0
estimates store A
regress price trunk turn if foreign == 0
estimates store B
regress price weight trunk turn if foreign == 0
estimates store C
regress price trunk if foreign == 0
estimates store D
// code for identification strategy 2
regress price mpg trunk length turn if foreign == 1
estimates store E
regress price trunk turn if foreign == 1
estimates store F
regress price weight trunk turn if foreign == 1
estimates store G
regress price trunk if foreign == 1
estimates store H
然后绘制如下图:
local gap1 : display _dup(30) " "
local gap2 : display _dup(400) " "
local label1 |`gap1'|`gap1'|`gap1'|`gap2'
local gap3 : display _dup(25) " "
local gap4 : display _dup(400) " "
local label2 DV1`gap3'DV2`gap3'DV3`gap3'DV4`gap4'
coefplot (A, offset(-0.38)) (E, offset(-0.33)) ///
(B, offset(-0.15)) (F, offset(-0.1)) ///
(C, offset(0.09)) (G, offset(0.135)) ///
(D, offset(0.3)) (H, offset(0.35)), ///
vertical keep(trunk) coeflabels(trunk = `""`label1'""`label2'""') ///
xlabel(, notick labgap(0)) legend(off)
这里的代码使用 Stata 的玩具auto
数据集为每个foreign
类别运行许多简单的回归。相同的因变量price
用于说明,但您可以在其位置使用不同的变量。trunk
这里假设变量是感兴趣的变量。
上面的玩具代码片段基本上做的是保存每组回归估计,然后模拟成对的模型。标签DV
是根据我最近在以下两个问题中展示的 hack 绘制的:
结果是一个非常好的近似值,但需要您进行实验: