stata分组回归命令

您所在的位置:网站首页 查询美金汇率的网站 stata分组回归命令

stata分组回归命令

#stata分组回归命令| 来源: 网络整理| 查看: 265

1. 问题

在实证分析中,我们经常需要在模型中加入行业虚拟变量、年度虚拟变量等,以便控制不可观测的行业个体效应或年度个体效应。然而,在正式报告结果时,我们无需报告这些虚拟变量的系数,否则结果表格会变得非常冗长。

简言之,在估计模型时,我们需要加入这些虚拟变量,而在最终呈现结果时,只在表格中进行标注,说明我们已经控制了这些虚拟变量,而受限于篇幅,没有呈现这些变量的估计系数。

2. 解决方法 1

事实上,Stata 里有多个命令可以帮我们处理这种情形。这里以 esttab 命令为例进行说明。

2.1 esttab 命令输出结果的原理

所有这些输出结果的命令,背后的基本原理都是相同的。在完成一个模型的估计后,Stata 会把相应的估计系数、标准误、R2 等统计量统一存储在内存中,称之为“返回值 (return values)”。而诸如 esttab、outreg2 等呈现估计结果的命令,无非是把内存中存储的这些统计量以一种比较标准的格式呈现在屏幕上,或输出到 word 或 excel 文档中。

例如,使用 regress 命令完成 OLS 估计后,进而输入 ereturn list 即可在屏幕上呈现出内存中存储的统计量:

. sysuse "nlsw88.dta", clear (NLSW, 1988 extract) . reg wage ttl_exp married Source | SS df MS Number of obs = 2,246 -------------+---------------------------------- F(2, 2243) = 85.79 Model | 5284.82149 2 2642.41074 Prob ; F = 0.0000 Residual | 69083.1459 2,243 30.7994409 R-squared = 0.0711 -------------+---------------------------------- Adj R-squared = 0.0702 Total | 74367.9674 2,245 33.1260434 Root MSE = 5.5497 ------------------------------------------------------------------------------ wage | Coef. Std. Err. t P;|t| [95% Conf. Interval] -------------+---------------------------------------------------------------- ttl_exp | 0.330 0.025 12.95 0.000 0.280 0.379 married | -0.291 0.245 -1.19 0.235 -0.771 0.189 _cons | 3.823 0.383 9.99 0.000 3.072 4.574 ------------------------------------------------------------------------------ . ereturn list //内存中的返回值 scalars: e(N) = 2246 e(df_m) = 2 e(df_r) = 2243 e(F) = 85.79411399585418 e(r2) = .0710631427762836 e(rmse) = 5.549724397561928 e(mss) = 5284.821485260705 e(rss) = 69083.14591378947 e(r2_a) = .0702348441965032 e(ll) = -7034.513553793961 e(ll_0) = -7117.294949670642 e(rank) = 3 macros: ... ... matrices: e(b) : 1 x 3 e(V) : 3 x 3 functions: ... ...

在上面的例子中,内存中的返回值包括四种类型,分别是:单值(scalars)、暂元(macros)、矩阵(matrices),以及函数(functions)。这些看似凌乱的返回值,都可以采用 esttab 等命令,以一种非常规整的方式重新呈现在屏幕上,或输出到 word 文档中。例如,当我们执行如下命令后,屏幕上会呈现如下结果:

. esttab, nogap scalar(N df_m F r2) ---------------------------- (1) wage ---------------------------- ttl_exp 0.330*** (12.95) married -0.291 (-1.19) _cons 3.823*** (9.99) ---------------------------- N 2246 df_m 2 F 85.79 r2 0.0711 ---------------------------- t statistics in parentheses * p


【本文地址】


今日新闻


推荐新闻


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