Calculates a likelihood-ratio test for comparing nested models. Regression models for limited dependent variables are supported.
For example:
library(catregs)
data("Mize19AH")
m1 <- glm(alcB ~woman*parrole + age + race2 + race3 + race4 + income + ed1 + ed2 + ed3 + ed4,family="binomial",data=Mize19AH)
m0 <- glm(alcB ~woman + parrole + age + race2 + race3 + race4 + income + ed1 + ed2 + ed3 + ed4,family="binomial",data=Mize19AH)
lr.test(m1,m0)
## LL Full LL Reduced LR Statistic DF p-value
## -2692.32 -2693.539 2.437886 1 0.11844
## LL Full LL Reduced LR Statistic DF p-value
## -2693.539 -2692.32 2.437886 1 0.11844