Stata 回归表格输出之 esttab 详解

您所在的位置:网站首页 stata回归分析 Stata 回归表格输出之 esttab 详解

Stata 回归表格输出之 esttab 详解

2023-10-29 18:58| 来源: 网络整理| 查看: 265

1. 简介

esttab 命令是由瑞士波恩大学社会学研究所(University of Bern, Institute of Sociology)的 Ben Jann 教授编写的 Stata用户外部命令,主要用于生成满足用户需求的回归表格,这类命令已经成为量化实证分析中的基础性技能,兼具效率、规范与美观。本文是对该命令的详细介绍。

2. 安装ssc install estout, replace //estout 是一个通用程序,是 esttab 命令实现其功能的引擎。3. 使用 esttab 整理回归结果3.1 基本语法及使用esttab [ namelist ] [ using filename ] [, options estout_options ]

该命令的使用逻辑如下:

首先,将目标回归模型的估计结果进行存储;其次,使用esttab命令将存储好的估计值或统计量编辑在一个回归表格中;

通过下面的示例进一步掌握:

** 保存模型结果之方法一:使用 eststo 命令存储回归模型。通过在回归模型前加上 eststo 前缀,esttab 命令能够自动找到储存的信息,并自动生成对应每个回归模型的变量。 sysuse auto, clear eststo: regress price weight mpg //_est_est1 eststo: regress price weight mpg foreign //_est_est2 esttab // 为节省篇幅,提高阅读效率,除非必要,结果不再一一展示,下文不再说明 eststo clear //删除保存的模型结果,_est_est1 和 _est_est2 从变量窗口中被消失 ** 保存模型结果之方法二:使用 Stata 官方自带的命令进行结果存储,即 estimates store(或缩写形式:est store) sysuse auto, clear regress price weight mpg estimates store model1 regress price weight mpg foreign estimates store model2 esttab model1 model2 *注:此时若使用 eststo clear,保存的结果无法删除3.2 标准误、P值及概要统计量

esttab的输出结果默认展示带有 t 统计量的的点估计,并在表格注脚处进行标记,一并默认输出的还有样本量。实证研究中,我们往往需要输出估计量的标准误以及其他特定的统计量,这时需要在esttab后进行更为细致的设定。例如,我们需要报告标准误而非 t 统计量,并且需要汇报调整后的判决系数(adj. R-sq)。通过下面的例子,我们可以发现 esttab 命令输出回归结果的设定思路。

sysuse auto, clear eststo: regress price weight mpg eststo: regress price weight mpg foreign esttab, se ar2 // 跑完回归模型,可以通过 ereturn list 获得该模型的的 scalars、macros、matrices以及functions *上面的命令等价于 esttab, se scalars(r2_a) //r2_a 是存储在模型估计结果的 scalars中的一个标量参数 *也可以输出 p 值,置信区间(ci)或是任何存储在估计结果中的参数统计量,它们保存在 scalars 中。例如,我们需要输出p值、F统计量以及模型自由度。 esttab, p scalars(F df_m df_r) eststo clear //注意,如果把该命令误写为 est sto clear,它表示将上一个回归模型命名为“clear”并储存起来,使用时要理解原理,不断试错才能避免出错3.3 标准化系数为什么要进行标准化?

同一回归模型中,即便两个自变量的单位一致(例如教育年限和工作经历都以年为计数单位),其回归系数也无法直接进行比较。事实上,研究中涉及的自变量往往具有不同的测度单位,回归系数也会受到影响。

多元回归模型经常涉及各自变量对因变量的相对作用大小进行比较,进而从多个因素中找出首要和次要因素,这时便可以采用标准化的回归系数(standardized coefficients)。所谓标准化回归系数,是将自变量转为一个无量纲的变量,使得不同标准化回归系数之间具有可比性。

随机变量的标准化处理方式

标准化回归系数是通过标准化回归模型估计得到的,需要将模型中所有的自变量和因变量都进行标准化。

将随机变量 X 进行标准化的表达式如下:

z=\frac{X-E(X)}{\sigma(X)}

其中,E(X) 和 \sigma(X) 分别表示期望和标准差。z 也叫标准化随机变量或标准分,是一个无量纲的纯数。变量的标准化形式表示:以标准差为单位,测量观测值与均值之间的距离。经过标准化处理后的新变量(z),其均值为0,方差为1。此外,我们还应该知道,标准化处理其实也是一个对中(centering)和测度转换(rescaling)的过程,经过标准化转换,不同变量的位置和尺度得以一致。

估计标准化模型

应用中,我们可以直接将标准化后的变量进行回归,但理论上还是要尽量知其所以然,以便更好理解对标准化模型的估计,很多细节上的问题也自然迎刃而解。假设标准化处理之前的原始回归模型设定如下:

y_i=\beta_0+\beta_1x_{i1}+\beta_2x_{i2}+…+\beta_px_{ip}+\epsilon_i \tag{1}

对模型两侧的变量进行标准化转换,即:y_i^*=\frac{y_i-\bar{y}}{S_y} ,x_{ik}^*=\frac{x_{ik}-\bar{x_k}}{S_{xk}}

对(1)式进行对中处理并除以因变量的标准差 S_y ,模型(1)可做如下变换:

\begin{align} \frac{yi-\bar{y}}{S_y}=& \frac{\beta_1}{S_y}(x{i1}-\bar{x1})+\frac{\beta_2}{S_y}(x{i2}-\bar{x2})+…+\frac{\beta_p}{S_y}(x{ip}-\bar{x_p})+\frac{\epsilon_i}{S_y} \\ y_i^*=& \beta1\frac{S{x1}}{Sy}(\frac{x{i1}-\bar{x1}}{S{x1}})+\beta2\frac{S{x2}}{Sy}(\frac{x{i2}-\bar{x2}}{S{x2}})+…+\beta_p\frac{S{xp}}{Sy}(\frac{x{ip}-\bar{xp}}{S{xp}})+\frac{\epsilon_i-\bar\epsilon}{S_y} \\ y_i^*= & \beta_1^*x{i1}^*+\beta_2^*x{i2}^*+…+\beta_p^*x_{ip}^*+\epsilon_i^* \tag{2} \end{align} \begin{align} \end{align}此时,通过比较模型(1)和模型(2),可以得到标准化与非标准化回归系数之间的关系如下:

\beta_k^*=\beta_k\times \frac{S{xk}}{S_y} \tag{3}

显然,利用(3)式,我们也可以通过计算样本中变量 y 与 x_k 的标准差,在获得非标准化系数后求得标准化系数。

两种回归系数的比较

标准化回归系数处于 -1, 1 的区间内,并且可以进行标准化尺度下的变量间系数比较。

相比而言,非标准化回归系数虽然不可比,但它将数据本身的信息提供给我们,估计系数具有明确的现实意义。

由此可见,并非进行标准化处理就显得“高大上”,标准化处理与否需要根据所研究的具体问题进行选择。但是,不论选择哪一种,尤其要关注对两种回归系数的解释。同是边际效应,标准化回归系数表示自变量每增加1个标准差,因变量平均增加 \beta_k^*个标准差,具体的数学表达不再展示(提示:对(3)式进行变换)。

使用esttab汇报标准化回归系数sysuse auto, clear eststo: regress price weight mpg eststo: regress price weight mpg foreign esttab /*-------------------------------------------- (1) (2) price price -------------------------------------------- weight 1.747** 3.465*** (2.72) (5.49) mpg -49.51 21.85 (-0.57) (0.29) foreign 3673.1*** (5.37) _cons 1946.1 -5853.7 (0.54) (-1.73) -------------------------------------------- N 74 74 -------------------------------------------- t statistics in parentheses * p


【本文地址】


今日新闻


推荐新闻


CopyRight 2018-2019 办公设备维修网 版权所有 豫ICP备15022753号-3