我使用输入 response
响应变量向量和 x
数据框创建了函数。该函数拟合具有聚类标准误差的线性模型,并创建一个包含系数和置信区间的数据框。
clusterLM = function(response,x) {
formula = paste(response, "~ vs*am")
model=lm.cluster(as.formula(formula), data = x ,cluster="carb")
model=bind_cols(data.frame( model$lm_res$coefficients),data.frame(confint(model)))
names(model)=c("estimate","conf.low","conf.high")
return(model)
}
然后此代码将模型作为数据框返回
vars_0 = names(mtcars)[c(3,5:7)]
addLog=function(x){ paste("log(",x,")")}
vars = names(mtcars)[c(3,5:7)]
vars=sapply(vars,addLog)
vars = set_names(vars,vars_0)
models = vars %>%
map_dfr( clusterLM,mtcars, .id="id")
我需要将该函数应用于数据框列表而不是一次一个
dfs=list(mtcars,mtcars)
最好的方法是什么?