固定效应模型企业聚类HC1稳健标准误计算及默认聚类水平咨询

一、按企业聚类的HC稳健标准误计算 #

要在你的固定效应模型基础上计算按企业（个体）聚类的HC1稳健标准误，只需在vcovHC()中指定cluster = "group"参数——因为你的模型已经通过index=c(firm, year)和effect='individual'把企业设为个体分组维度。完整代码如下：

# 估计企业固定效应模型
model <- plm(dep ~ a *b + a + b + controls, index=c(firm, year), method='within', effect='individual')

# 计算按企业聚类的HC1稳健标准误
vcov_firm_cluster <- vcovHC(model, type="HC1", cluster = "group")

# 输出带稳健标准误的回归结果（需加载lmtest包）
library(lmtest)
coeftest(model, vcov = vcov_firm_cluster)

cluster = "group"会锁定模型的个体分组（企业），在企业层面聚类计算标准误，能同时修正异方差和企业内观测值的序列相关问题。

二、未指定聚类方式时的默认水平 #

当调用vcovHC()不写cluster参数时，默认是不做聚类的纯异方差修正：仅针对横截面异方差调整标准误，不考虑任何组内（比如企业内、年度内）的相关性问题。这种默认设置只适用于横截面数据，对面板数据而言，通常需要指定聚类维度才能得到更可靠的标准误。

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