Deutsches R-Forum

Verfasst: **Sa Jan 02, 2021 1:15 pm**

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

Verfasst: **Sa Jan 02, 2021 2:05 pm**

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...

Verfasst: **Di Jan 05, 2021 5:49 pm**

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

Deutsches R-Forum

Survey Packet LR Test

Survey Packet LR Test

Re: Survey Packet LR Test

Re: Survey Packet LR Test