fixest 固定效应回归简介

什么是 fixest？

fixest 是 R 语言中一个高性能的计量经济学包，专门用于估计固定效应模型。它在处理面板数据（Panel Data）时表现卓越，能够同时吸收多个维度的固定效应（如个体效应、时间效应），且计算速度远超传统的 plm 或 lm 包。

核心优势

1. 极速计算：采用优化的 C++ 后端，即使在大样本和高维固定效应下也能秒级完成回归。
2. 语法简洁：使用 ~ 分隔解释变量与固定效应变量（例如：y ~ x | id + time），无需手动生成虚拟变量。
3. 稳健标准误：内置多种聚类标准误（Clustered Standard Errors）修正方法，一键解决序列相关和异方差问题。
4. 多功能输出：支持直接导出 LaTeX、Markdown 或 HTML 格式的回归表格。

适用场景

• 面板数据分析（个体/时间双固定效应）
• 高维固定效应模型（如控制成千上万个行业或地区效应）
• 需要快速迭代模型设定的实证研究

结果展示说明

本部分描述运行 main.R 脚本后，在 R 控制台或 RStudio 中预期的文字输出内容。

1. 输入操作

用户在 R 环境中 sourced main.R 文件，或直接逐行运行代码块。

2. 预期输出内容

A. 模型估计信息

系统首先会输出一行简短的估计信息，格式如下：

OLS estimation, Dep. Var.: y
Observations: 1,000 
Fixed-effects: id: 100, time: 10
Standard-errors: IID 

• 描述：显示使用了 OLS 估计，因变量为 y，样本量 1000，并成功吸收了 100 个个体固定效应和 10 个时间固定效应。默认标准误类型为 IID（同方差）。

B. 系数估计表

紧接着输出详细的系数表格，包含以下列：

• Estimate: 估计系数值。预期 x1 接近 2.0，x2 接近 3.0。
• Std. Error: 标准误。
• t value: t 统计量。
• Pr(>|t|): P 值，通常标记显著性星级（*** 表示 1% 水平显著）。

示例文本片段：

            Estimate Std. Error   t value   Pr(>|t|)    
x1          2.0150     0.0320    62.9700   < 2.2e-16 ***
x2          2.9840     0.0315    94.7300   < 2.2e-16 ***

- --
   Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

• 描述：展示了 x1 和 x2 的回归系数及其统计显著性。由于是模拟数据，系数应非常接近真实生成过程设定的 2 和 3，且 P 值极小，显示高度显著。

C. 拟合优度统计量

表格下方会列出模型的整体拟合指标：

RMSE: 0.985   Adj. R2: 0.8543   Within R2: 0.7621

• 描述：
• RMSE: 均方根误差，反映残差大小。
• Adj. R2: 调整后的 R 平方，考虑了固定效应后的整体解释力。
• Within R2: 组内 R 平方，专门衡量解释变量对去均值后数据的解释能力（固定效应模型的核心指标）。

3. 结论解读

输出结果表明，在控制了个体和时间固定效应后，变量 x1 每增加 1 个单位，y 平均增加约 2.01 个单位；变量 x2 每增加 1 个单位，y 平均增加约 2.98 个单位，且两者均在 1% 的水平上统计显著。模型具有较高的组内拟合优度。