将观星者扩展到多路湾

Question

我正在使用观星者为我的学士论文创建回归输出。由于我的数据结构，我必须使用聚类模型(下面的代码(。我正在使用 multiwaycov 包中的 vcovclust 命令，它运行良好。但是，观星者不支持它。您知道另一种创建像观星者一样漂亮的输出的方法吗？或者您是否知道由观星者提供的其他包/命令来聚类模型？

model1.1.2 <- lm(leaflet ~ partisan + as.factor(gender) + age + as.factor(education) + meaning + as.factor(polintrest), data = voxit)
summary(model1.1.2)
#clustering
vcov_clust1.1.2 <- cluster.vcov(model1.1.2, cbind(voxit$id, voxit$projetx))
coeftest(model1.1.2, vcov_clust1.1.2)

Answer 1

您可以手动向观星者提供调整后的 p 值和 se 值。

# model1 and model2 are both objects returned from coeftest()
# Capture them in an object and extract the ses (2nd column) and ps (4th column) in a list
ses <- list(model1[,2], model2[,2])
ps <- list(model1[,4], model2[,4])
# you can then run your normal stargazer command and supply
# the se- and p-values manually to the stargazer function
stargazer(model1, model2, type = "text", se = ses, p = ps, p.auto = F)

希望这有帮助！