r - 从 glm 拟合的 lapply 列表中提取 p 值

Question

我正在使用 lapply 一次对一个因变量执行多个 glm 回归一个自变量。现在我对Pr(>|z|)特别感兴趣。每个自变量。但是，我不确定如何仅报告 Pr(>|z|)使用 lapply 中的列表。

如果我一次只运行一个模型:coef(summary(fit))[,"Pr(>|z|)"]或者summary(fit)$coefficients[,4]会工作(如 here 所述)，但尝试与 lapply 类似的东西似乎不起作用。我可以使用 lapply 只获得 p 值吗？和 glm使用访问器方法还是直接从模型调用？

#mtcars dataset
vars <- names(mtcars)[2:8]
fits <- lapply(vars, function(x) {glm(substitute(mpg ~ i, list(i = as.name(x))), family=binomial, data = mtcars)})
lapply(fits,summary) # this works
lapply(fits, coefficients) # this works
#lapply(fits, summary(fits)$coefficients[,4])# this for example does not work

Answer 1

你想做:

lapply(fits, function(f) summary(f)$coefficients[,4])

但是，如果每个项目只是一个 p 值，您可能更愿意使用向量而不是列表，因此您可以使用 sapply而不是 lapply :

sapply(fits, function(f) summary(f)$coefficients[,4])