Załóżmy, że zamierzam wykonać jednoczynnikową regresję logistyczną dla kilku niezależnych zmiennych, takich jak to:
mod.a <- glm(x ~ a, data=z, family=binominal("logistic"))
mod.b <- glm(x ~ b, data=z, family=binominal("logistic"))
Zrobiłem porównanie modelu (test współczynnika prawdopodobieństwa), aby sprawdzić, czy model jest lepszy niż model zerowy za pomocą tego polecenia
1-pchisq(mod.a$null.deviance-mod.a$deviance, mod.a$df.null-mod.a$df.residual)
Następnie zbudowałem inny model ze wszystkimi zmiennymi
mod.c <- glm(x ~ a+b, data=z, family=binomial("logistic"))
Aby sprawdzić, czy zmienna jest istotna statystycznie w modelu wielowymiarowym, użyłem lrtest
polecenia zepicalc
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
Zastanawiam się, czy pchisq
metoda i lrtest
metoda są równoważne z wykonaniem testu wiarygodności? Jak nie wiem, jak korzystać lrtest
z univate modelu logistycznego.