Как цикл GLM модели для выбора модели (model.sel) - PullRequest
Новая
       8

Как цикл GLM модели для выбора модели (model.sel)

0 голосов
/ 05 февраля 2019

Я пытаюсь выбрать лучшую модель через AIC.Мне дали эти коды, которые мне сказали, что они автоматически создают комбинацию моделей glm, но я действительно не понимаю, как это работает.

Может кто-нибудь объяснить мне, как это работает, простыми словами и как заставить его работатьпотому что я не смог. 

library(MuMIn)

as.list(rep(NA, 44)) ->models

for (i in 1:44) {
  glm(log(Tbra+1)~bats[,i],data=bats)-> models[[i]]
}

names(models)<-names(bats)[7:50]

model.sel(models)

РЕДАКТИРОВАТЬ 1

У меня есть еще вопросы по этому коду.

  1. as.list(rep(NA, 44)) ->models

Что значит (NA, 44)?Как составлен список?Я задавался вопросом, что список должен быть сделан из базы данных, но этот код создал пустой список, который я не понимаю почему.

2. 

for (i in 1:44) {
  glm(log(Tbra+1)~bats[,i],data=bats)-> models[[i]]
}

Я понимаю, чтоglm сделан из данных "bats", но он требует, чтобы Bats был списком, верно?так как это сделано.Сначала я подумал, что, возможно, это относится к списку, созданному as.list, вроде: 

as.list(rep(NA, 44)) ->bats

  for (i in 1:44) {
      glm(log(Tbra+1)~bats[,i],data=bats)-> models[[i]]
    }

Но это ничего не значит.Есть мысли?

1 Ответ

0 голосов
/ 06 февраля 2019

Цикл составляет список моделей.Каждая модель пытается предсказать ответ log(Tbra + 1), используя одну переменную из данных.Затем функция model.sel() берет список, оценивает AIC каждого из них, а затем отображает их в таблице вместе с коэффициентами для каждой переменной, логарифмической вероятностью и т. Д.

...