使用lapply时群集SE

Question

我使用相同的IV但不同的DV运行许多不同的模型，并希望从这些输出系数和聚类SE。 虽然使用lapply时我不确定如何提取聚类SE。

基本上，我不确定如何复制以下两行：

cl.cov1 <- cluster.vcov(model, subsetviolence$ea) # cluster-robust SEs

cl.se.1 < - sqrt（diag（cl.cov1））

使用此代码：

part1<- c("eapost", "easchool", "eapolice", "eaelect", "eahealth", "eatownbuild", "gender", "age", "housen", "married", "school", "yoruba", "hausa", "igbo", "christian", "muslim", "relintensity", "agricfishmin", "indservtrader", "indservartisan", "outstudent", "outhousework","house", "land", "cattle", "radio", "cellphone", "expenditure", "attrition")
models <- lapply(1:length(part1), function(i) lm(as.formula(paste(part1[i],"~ treat")), data=surveyviolence))

将输出放入表中。