Survey Packet LR Test

[Hanibal]

Hallo liebe Community,

Ich möchte gerne einen "Likelihood Ratio Test Of Nested Models" lrtest() mit meinem "survey.design" Objekt für den Modellvergleich machen.
Gibt es eine Funktion aus dem Packet mit welchem ich das machen kann oder was würdet ihr benutzen?

Code: Alles auswählen

formula0 <- as.formula(fall~sex.f)
formula1 <- as.formula(fall~sex.f+alter_cut)

## Crude regression
fit0 <- svyglm(formula0, 
               design = wdat_07,
               family = quasibinomial(log))
fit1 <- svyglm(formula1, 
               design = wdat_12,
               family = quasibinomial(log))

sowas in der Art:

Code: Alles auswählen

lrtest(fit0, fit1)

Danke vielmals

[EDi]

Wenn das Paket (dass ich nicht kenne) für die svyglm Objekte ein Methode für logLik() zur Verfügung stellt, sollte lrtest() funktionieren.

Ansonsten bietet das Paket doch eine ?anova.svyglm methode an, damit kann man auch einen LRT rechnen...

[Hanibal]

Danke hat geklappt, man kann es mit anova() oder AIC() machen.