以 coeftest 作为参数时，Stargazer 不显示基本统计数据（例如 obs 计数、R2）答案

【问题标题】：Stargazer doesn't display basic stats (e.g. obs count, R2) when taking coeftest as argument以 coeftest 作为参数时，Stargazer 不显示基本统计数据（例如 obs 计数、R2）
【发布时间】：2020-04-01 21:30:47
【问题描述】：

我正在使用具有双向聚集标准误差的 lmtest 运行多元 OLS 回归。系数和标准误差工作得很好，但回归表底部没有显示通常的统计数据（例如观察数、R2、adj. R2、残差 SE、F 统计）。

代码如下：

model1 <- lm(Y ~ X + Z data=df)
model2 <- lm(Y ~ W + Y + Z, data=df)

stargazer(
(coeftest(model1,vcovCL, cluster = ~ A + B)),
(coeftest(model2,vcovCL, cluster = ~ A + B)),
type='text')

这是输出：

============================================
                    Dependent variable:     
                ----------------------------

                     (1)            (2)     
--------------------------------------------
X                   0.053***                 
                   (0.010)                  

W                                  0.048***   
                                  (0.014)   

Y                                  0.058***   
                                  (0.013)   

Z                   0.089***      0.089***   
                   (0.003)        (0.003)   

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

我了解到，当我不使用 coeftest 函数时，表格会在底部显示结果以及常见的统计数据。请在下面找到没有 coeftest 函数的示例，该函数在底部生成包含所需统计信息的输出：

stargazer(model1, model2, type='text')

===================================================================================
                                          Dependent variable:                      
                    ---------------------------------------------------------------                            
                                  (1)                             (2)              
-----------------------------------------------------------------------------------
X                                0.053***                                            
                                (0.002)                                            

W                                                                0.048***            
                                                                (0.002)            

Y                                                                0.058***            
                                                                (0.002)            

Z                               0.089***                        0.089***            
                               (0.0002)                        (0.0002)            

-----------------------------------------------------------------------------------
Observations                   1,119,677                       1,119,677           
R2                               0.151                           0.151             
Adjusted R2                      0.151                           0.151             
Residual Std. Error      0.079 (df = 1119673)            0.079 (df = 1119672)      
F Statistic         66,250.650*** (df = 3; 1119673) 49,690.200*** (df = 4; 1119672)
===================================================================================
Note:                                                   *p<0.1; **p<0.05; ***p<0.01

这个其他线程 (Appending statistics to coeftest output to include in stargazer tables) 解决了同样的问题，但没有一个答案解决了这个问题。

如何在通过 stargazer 传递 coeftest 函数时生成这些常见统计数据？

【问题讨论】：

标签： r linear-regression data-science stargazer

【解决方案1】：

尝试以这种方式对标准错误进行聚类：

rob_se <- list(sqrt(diag(vcovHC(model1, type = "HC1"))),
sqrt(diag(vcovHC(model2, type = "HC1"))))


stargazer(model1, model2
          digits = 3,
          header = FALSE,
          type = "latex", 
          se = rob_se)

【讨论】：

感谢您的反馈！在这个解决方案下，我在哪里指定我想用集群标准错误的两列？我必须在每个模型中以两种方式对标准错误进行聚类。