如何在R语言面板固定效应模型的stargazer表格中添加全局截距

在R的plm包中，固定效应（within）模型会通过组内变换消去全局截距，因此模型对象本身不会存储该截距值，而stargazer默认也不会自动添加这一项。要复刻Stata中显示全局截距的效果，可按以下步骤操作：

步骤1：计算全局截距及其标准误 #

使用within_intercept()函数并开启se=TRUE参数，同时获取截距的估计值和标准误：

library(stargazer)
library(plm)

data("Hedonic", package = "plm")
mod_fe <- plm(mv ~ age + crim, data = Hedonic, index = "townid")

# 计算全局截距和标准误
int_result <- within_intercept(mod_fe, se = TRUE)
int_est <- int_result$estimate  # 截距估计值
int_se <- int_result$se         # 截距标准误

步骤2：手动将截距项加入stargazer表格 #

通过stargazer的coef和se参数，将截距项与原模型的系数、标准误合并后传入，让表格显示全局截距：

# 提取原模型的系数和标准误
orig_coef <- coef(mod_fe)
orig_se <- sqrt(diag(vcov(mod_fe)))

# 合并截距项与原模型统计量
new_coef <- c(`(Global Intercept)` = round(int_est, 3), orig_coef)
new_se <- c(round(int_se, 3), orig_se)

# 生成带全局截距的表格
stargazer(mod_fe, type = "text", 
          coef = list(new_coef),
          se = list(new_se),
          omit.stat = "adj.rsq")  # 可选：移除调整R²，与Stata输出对齐

输出效果示例 #

运行上述代码后，表格会新增全局截距行，效果如下：

Dependent variable:    
             ---------------------------
                         mv             
----------------------------------------
(Global Intercept)     10.260***        
                       (0.102)          
                                        
age                   -0.004***         
                       (0.001)          
                                        
crim                  -0.009***         
                       (0.002)          
                                        
----------------------------------------
Observations             506            
R2                      0.147           
F Statistic    35.431*** (df = 2; 412)  

Note:        *p<0.1; **p<0.05; ***p<0.01

补充说明 #

• within_intercept()的计算逻辑是基于组内变换后的数据还原全局截距，和Stata中xtreg, fe输出的截距完全一致
• 如果不需要显示标准误，也可以用add.lines参数直接插入截距行，但这种方式无法自动添加显著性标记，因此更推荐上述合并系数的方法

内容的提问来源于stack exchange，提问作者Emily