У меня есть следующая модель:
mod.log.diab <- glm(data = data_train, INFEC ~ DIABETES, family = binomial(link = 'logit'))
Это дает следующее ИЛИ:
>exp(mod.log.diab$coefficients)
(Intercept) DIABETESno
2.826087 1.469822
Это значение больше 1, поэтому можно сказать, что переменная значима. Но когда я применяю тест Вальда, я получаю следующий ответ:
>regTermTest(mod.log.diab, "DIABETES")
Wald test for DIABETES
in glm(formula = INFEC ~ DIABETES, family = binomial(link = "logit"),
data = data_train)
F = 2.34702 on 1 and 1560 df: p= 0.12573
Значение p больше 0,05, поэтому мы отклоняем H0, тест говорит, что это несущественно.
Как я должен это интерпретировать, у кого правда ИЛИ или тест?
Вопрос еще не закончен, потому что, если я обновлю модель с более понятными переменными:
mod.log.cat <- glm(data = train.cat, INFEC ~ DIABETES + EDAD + HCTO, family = binomial(link = 'logit'))
Я получаю следующее ИЛИ:
>exp(mod.log.mult$coefficients)
(Intercept) DIABETESno EDAD HCTO
1.3275634 1.0751666 0.9817025 1.0558760
И следующий результат теста Вальда:
Wald test for EDAD
in glm(formula = INFEC ~ DIABETES + EDAD + HCTO, family = binomial(link = "logit"),
data = train.cat)
F = 18.90695 on 1 and 1558 df: p= 1.4615e-05
Wald test for DIABETES
in glm(formula = INFEC ~ DIABETES + EDAD + HCTO, family = binomial(link = "logit"),
data = train.cat)
F = 0.4122714 on 1 and 1558 df: p= 0.52091
Wald test for HCTO
in glm(formula = INFEC ~ DIABETES + EDAD + HCTO, family = binomial(link = "logit"),
data = train.cat)
F = 41.87676 on 1 and 1558 df: p= 1.2986e-10
На этот раз ДИАБЕТ не отвергает тест Уолда, как в простой модели. Но странный факт: все отношения нечетных чисел почти одинаковы, поэтому они близки к 1. Каков вывод? Как возможно с OR, близким к 1, некоторые переменные значимы, а другие нет?
Спасибо.
DIABETES
не имеет значения для $ \ alpha = 0,05 $. Обратите внимание, что в первом случае $ p = 0,12573 $. (2) То, что вы называете ИЛИ, является размером параметра и не соответствует значимости. (3) Корреляция междуDIABETES
иINFEC
, кажется, следует из двух других искажающих факторовHCTO
иEDAD
, которые сильно коррелированы сами по себе. - person tmrlvi   schedule 14.12.2019