Извлечение AIC из нескольких моделей регрессии

Question

У меня есть несколько бинарных моделей логистической регрессии в R (более 100). Я хотел бы перечислить все отдельные модели регрессии вместе с их AIC, нулевым отклонением, остаточным отклонением и т. Д. В этом формате

  Model      AIC        Null deviance
     reg1      155.13        ..
     reg2      154.     
     reg3      

Возможно ли иметь код, который бы достиг этого для меня, избегая ручной работы?

Спасибо

Answer 1

Вы можете извлечь параметры модели из списка моделей, используя 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