我目前正在使用 R 使用以下代码运行回归:
library(lmtest)
library(sandwich)
reg <- lm(Y ~ x1 + x1_sq + x2 + x2_sq + x1x2 + d2 + d3 + d4, df)
coeftest(reg, vcov = vcovHC(reg, type="HC1")
并且已经意识到我指定 -type="HC1"- 不符合惯例我一直只使用这个。我知道这是一个异方差一致的估计器(White,1980)并利用自由度校正,但 -sandwich- 包包括 5 个不同的估计器。有 HCO,标准的“三明治”估计器,HC1,我使用和习惯使用的那个,它包括 df 校正,使得,然后是 HC2-HC4 ,它们都采用有趣的形式:
其中是帽子值,是它们的平均值, (从这组幻灯片中获得的公式。
我目前无法理解这些估计器之间的功能差异。即使用其中一个的好处和权衡是什么?观察次数是否有截止值,哪一个比下一个更受欢迎?