高维固定效应模型 (reghdfe) 简介

概述

reghdfe 是 Stata 中用于估计包含任意数量高维固定效应（High-Dimensional Fixed Effects, HDFE）的线性回归模型的命令。它由 Sergio Correia 开发，能够高效地吸收多个维度的固定效应，解决了传统 areg 或 xtreg 在处理多维固定效应时计算速度慢或无法估计的问题。

核心功能

1. 多维吸收：允许通过 absorb() 选项同时控制多个分类变量的固定效应（例如：个体、时间、行业、地区等）。
2. 计算优化：基于 ftools 包进行加速，利用迭代去均值算法（Iterated Within Transformation），在处理大规模面板数据时显著优于传统方法。
3. 灵活的标准误：支持聚类稳健标准误（Cluster-Robust Standard Errors），即使聚类维度与固定效应维度不同也能正确处理。
4. 兼容性：支持与 estout, outreg2 等结果输出命令无缝对接。

适用场景

• 双重差分模型 (DID) 中包含个体和时间固定效应。
• 引力模型 (Gravity Model) 中控制出口国、进口国及年份固定效应。
• 任何需要控制大量虚拟变量但不希望直接估计这些虚拟变量系数的场景。

基本语法

reghdfe depvar indepvars, absorb(fe_vars) [options]

• depvar: 因变量。
• indepvars: 核心解释变量。
• absorb(): 指定需要吸收的固定效应变量列表。
• cluster(): 指定聚类变量以计算稳健标准误。

结果展示说明

本部分描述在执行 main.do 中的回归命令后，Stata 结果窗口将输出的文字内容逻辑。

1. 输入命令

用户在 Stata 命令行或 Do-file 编辑器中运行：

reghdfe price mpg weight, absorb(rep_group turn_group) cluster(turn_group)

2. 输出内容描述

执行上述命令后，Stata 结果窗口将依次显示以下文本信息：

A. 迭代收敛信息

• 首先显示固定效应吸收过程的迭代日志。
• 格式通常为：absorbing rep_group turn_group ... followed by iteration counts (e.g., iteration 1: tolerance = ...).
• 最后显示收敛状态，如：converged in X iterations。

B. 回归统计量表

• Header: 显示方差分析表头，包括 Source, SS, df, MS。
• Model Statistics:
  • Number of obs: 显示参与回归的样本量（例如 74）。
  • F(df_m, df_d): 显示 F 统计量及其分子自由度（核心变量个数）和分母自由度（扣除固定效应和聚类调整后的自由度）。
  • Prob > F: 显示模型整体显著性的 P 值。
  • R-squared: 显示总体 R 方。
  • Adj R-squared: 显示调整后的 R 方。
  • Root MSE: 显示均方根误差。

C. 系数估计表

表格包含以下列：

• Coef.: 核心自变量 (mpg, weight) 的回归系数估计值。注意：固定效应变量 (rep_group, turn_group) 的系数不会被直接列出，因为它们已被“吸收”。
• Std. Err.: 基于 cluster(turn_group) 计算出的聚类稳健标准误。
• t: t 统计量 (Coef. / Std. Err.)。
• P>|t|: 对应的双尾检验 P 值。
• [95% Conf. Interval]: 95% 置信区间的上下限。

D. 固定效应详情注脚

• 在系数表下方，会有一段文字说明被吸收的固定效应情况。
• 内容示例：note: 2 variables omitted because of collinearity (如果有共线性)。
• 明确列出吸收的变量及其层级数，例如：
  • absorbed singletons: 0 (是否吸收了单一观测组)
  • Fixed effects: rep_group (5 levels), turn_group (5 levels)。

E. 最终提示

• 如果使用了 estimates store，不会直接在回归表中显示，但系统后台已保存结果，用户随后可通过 esttab 或 estout 调用格式化表格。
• 控制台最后返回 Stata 提示符 .，表示命令执行完毕且无报错。