將Stargazer擴展爲multiwaycov

Question

我使用stargazer爲我的學士論文創建迴歸輸出。由於我的數據結構，我必須使用羣集模型（下面的代碼）。我使用multiwaycov軟件包中的vcovclust命令，該命令完美工作。但是，stargazer不支持它。你知道另一種創建輸出的方式，就像觀星者一樣好嗎？或者你知道其他的包/命令來聚集模型，這是由占星師支持？

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值提供給stargazer。

# 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, t.auto = F, p.auto = F) 

希望這會有所幫助！