Farz edelim ki birkaç bağımsız değişkende tek değişkenli bir lojistik regresyon yapacağım, bunun gibi:
mod.a <- glm(x ~ a, data=z, family=binominal("logistic"))
mod.b <- glm(x ~ b, data=z, family=binominal("logistic"))
Bu komutla modelin null modelden daha iyi olup olmadığını görmek için bir model karşılaştırma yaptım (olasılık oranı testi).
1-pchisq(mod.a$null.deviance-mod.a$deviance, mod.a$df.null-mod.a$df.residual)
Sonra içindeki tüm değişkenleri olan başka bir model kurdum
mod.c <- glm(x ~ a+b, data=z, family=binomial("logistic"))
Çok değişkenli modelde değişkenin istatistiksel olarak anlamlı olup olmadığını görmek için, lrtest
komutunu kullandım.epicalc
lrtest(mod.c,mod.a) ### see if variable b is statistically significant after adjustment of a
lrtest(mod.c,mod.b) ### see if variable a is statistically significant after adjustment of b
pchisq
Metot ve lrtest
metodun mantık testi yapmakla aynı olup olmadığını merak ediyorum. lrtest
Univate lojistik modeli için nasıl kullanılacağını bilmiyorum .