Я использую много разных моделей с одинаковыми IV, но разными DV и хотел бы вывести из них коэффициенты и кластерные SE. Я не уверен, как извлечь кластерные SE, хотя при использовании lapply.
По сути, я не уверен, как повторить следующие две строки:
cl.cov1 <- cluster.vcov(model, subsetviolence$ea) # cluster-robust SEs
кл.с.1 <- sqrt (diag (cl.cov1)) </p>
Используя этот код:
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))
чтобы получить вывод в таблицу.