对于我正在教授的统计课程（在R和Stata中），我想在R中复制Stata的robust命令。可以使用{{1}中的coeftest命令轻松完成此操作}套件。使用来自lmtest库的数据的示例，我得到以下内容
car这复制了Stata所谓的稳健回归。但是，现在我想使用library(car)
library(lmtest)
library(sandwich)
mod <- lm(prestige ~ education + income, data=Prestige)
coeftest(mod, vcov = vcovHC(mod, "HC1"))
包，并使用异方差稳健的标准误差来获得此回归的预测值（和置信区间）。根据软件包的文档，我应该能够使用effects参数将协方差矩阵更改为健壮矩阵。但是，当我尝试如下操作时
vcov.我收到错误消息
effect("education", mod, vcov.=vcovHC(mod, "HC1"))
我也尝试了Error in vcov.(mod) : could not find function "vcov."
命令，因为它与hccm包来自同一作者，但是我得到了相同的错误消息。我不知道我的错误在哪里。有什么想法吗？