r - 从多个回归模型中提取 AIC
问题描述
我在 R 中有一些二元逻辑回归模型(超过 100 个)。我想以这种格式列出所有单独的回归模型及其 AIC、空偏差、残余偏差等
Model AIC Null deviance
reg1 155.13 ..
reg2 154.
reg3
是否有可能有一个代码可以为我实现这一目标,避免手动工作
谢谢
解决方案
您可以使用函数从模型列表中提取模型的参数sapply
。然后将这些参数聚合到一个数据框中。
请参阅下面的代码以及逻辑回归模型的模拟列表。
set.seed(123)
# simulation
# generates list of 100 models
xs <- replicate(
n = 100,
expr = {
x <- rnorm(1000)
z <- 1 + 2 * x
pr <- 1 / (1 + exp(-z))
y <- rbinom(1000, 1, pr)
#now feed it to glm for logistic regression
model <- glm(y ~ x, family = "binomial")
model
},
simplify = FALSE
)
# list of models created
# now let's create a data frame
# extracting parameter from the models
summary(xs[[1]])
df <- data.frame(
model = seq_along(xs),
aic = sapply(xs, AIC),
null_deviance = sapply(xs, function(x) x$null.deviance),
df_residual = sapply(xs, function(x) x$df.residual)
)
head(df)
输出:
model aic null_deviance df_residual
1 1 867.0781 1259.696 998
2 2 853.3573 1311.370 998
3 3 904.3534 1276.693 998
4 4 922.6175 1296.127 998
5 5 884.0100 1271.172 998
6 6 878.9568 1289.871 998
推荐阅读
- wordpress - CloudFront 到 EC2 源返回 502 错误。我该如何调查?
- javascript - 我需要在 javascript 代码中格式化日期值
- c# - C# LINQ 无法执行 LEFT JOIN
- c# - 如何使用 PGP 加密 SFTP 流
- google-cloud-platform - 由于链接到 Dialogflow 代理,无法删除 Google Cloud Platform 项目
- javascript - 使用 ajax 和 php 保存 gif
- python - 从线程查询自动缩放组的 EC2 cloudwatch 指标
- r - 神经网络模型误差值这么大
- c - 在C中打印每个元素具有相同数量字符的矩阵
- javascript - 如何在 HTML 中以断线显示输出?