本文介绍了如何创建既具有普通标准误差又具有稳健标准误差的表?的处理方法,对大家解决问题具有一定的参考价值,需要的朋友们下面随着小编来一起学习吧!
问题描述
我经常在出版物中看到回归表,其中报告了普通标准误差(在括号中),以及下面的稳健标准误差(在括号中)。这些表还包括括号/括号旁边的星号,表示统计意义。
创建这样的回归报告最明智的方法是什么?
到目前为止,我一直在使用Stata中的estout
包。对于给定的模型,我可以有一列具有普通标准误差,而另一列具有稳健的标准误差。
例如,使用estout
,我可以执行以下操作:
eststo: qui reg ROE duality
eststo: qui reg ROE duality, vce(cluster firm)
esttab b(%9.3fc) ///
se(%9.3fc) ///
star (* 0.5 ** 0.25)
上述代码片段产生:
--------------------------------------------
(1) (2)
ROE ROE
--------------------------------------------
duality -8.090** -8.090*
(6.585) (7.067)
--------------------------------------------
N 647 647
--------------------------------------------
然而,该表浪费了列空间,因为两列的点估计将是相同的,唯一的区别是来自不同方差-协方差估计器的标准误差。我更想要的是下面这样的表格:
------------------------
(1)
ROE
-------------------------
duality -8.090
(6.585)**
[7.067]*
-------------------------
N 647
-------------------------
请注意,0.5和0.25的统计意义仅供说明,并不反映惯例。
推荐答案
您只需手动添加健壮的标准错误:
sysuse auto, clear
eststo clear
quietly regress price weight mpg, vce(robust)
matrix regtab = r(table)
matrix regtab = regtab[2,1...]
matrix rbse = regtab
eststo: quietly regress price weight mpg
estadd matrix rbse = rbse
esttab, cells(b se rbse)
-------------------------
(1)
price
b/se/rbse
-------------------------
weight 1.746559
.6413538
.777837
mpg -49.51222
86.15604
95.8074
_cons 1946.069
3597.05
4213.793
-------------------------
N 74
-------------------------
要按照您的规范设置格式,需要使用相关选项:
esttab , cells("b(fmt(a3) star)" "se(fmt(a2) par)" "rbse(fmt(a2) par([ ]))") ///
star(* 0.5 ** 0.25) addnote("Robust SE in brackets" "* p<0.5, ** p<0.25") ///
nonumbers
---------------------------
price
b/se/rbse
---------------------------
weight 1.747**
(0.64)
[0.78]
mpg -49.51
(86.2)
[95.8]
_cons 1946.1
(3597.0)
[4213.8]
---------------------------
N 74
---------------------------
Robust SE in brackets
* p<0.5, ** p<0.25
这篇关于如何创建既具有普通标准误差又具有稳健标准误差的表?的文章就介绍到这了,希望我们推荐的答案对大家有所帮助,也希望大家多多支持!