為什么 stargazer 為我的 lm 回歸輸出提供了一個輸出表，但在我使用 lm_robust() 時卻沒有？

Question

當我使用

est <- lm_robust(x~y, data=df) stargazer(est, est, type="text")

R 說“ % Error: Unrecognized object type. ”但是當我不嘗試使用強大的標准錯誤時會正確打印我的輸出。

我只是想找到一種方法將我的lm回歸輸出組織到一個簡單的表格中以包含在我的論文中。 如果不使用stargazer包，還有其他方法可以做到這一點嗎？

Answer 1

您可以對強大的標准錯誤和stargazer使用解決方法，如下所示：

library("sandwich")
library("plm")
library("stargazer")

data("Produc", package = "plm")

# Regression    
model <- plm(log(gsp) ~ log(pcap) + log(pc) + log(emp) + unemp,
             data = Produc, 
             index = c("state","year"),
             method="pooling")

# Adjust standard errors
cov1         <- vcovHC(model, type = "HC1")
robust_se    <- sqrt(diag(cov1))

# Stargazer output (with and without RSE)
stargazer(model, model, type = "text",
          se = list(NULL, robust_se))