Giả sử tôi sẽ thực hiện một hồi quy logistic đơn biến trên một số biến độc lập, như sau:
mod.a <- glm(x ~ a, data=z, family=binominal("logistic"))
mod.b <- glm(x ~ b, data=z, family=binominal("logistic"))
Tôi đã thực hiện một so sánh mô hình (kiểm tra tỷ lệ khả năng) để xem liệu mô hình có tốt hơn mô hình null bằng lệnh này không
1-pchisq(mod.a$null.deviance-mod.a$deviance, mod.a$df.null-mod.a$df.residual)
Sau đó, tôi đã xây dựng một mô hình khác với tất cả các biến trong đó
mod.c <- glm(x ~ a+b, data=z, family=binomial("logistic"))
Để xem biến có ý nghĩa thống kê trong mô hình đa biến hay không, tôi đã sử dụng lrtest
lệnh từ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
Tôi tự hỏi nếu pchisq
phương pháp và lrtest
phương pháp là tương đương để làm kiểm tra loglikabilities? Như tôi không biết làm thế nào để sử dụng lrtest
cho mô hình logistic univate.