vendredi 31 mai 2019

Compare Mixed Effects Model to Fixed Effects in R

I would like to test the significance of my random intercept. So I would like to compare my mixed model: fm1 <- glmer (BinaryOutcome ~ Age + (1|Center), family='binomial', data=DF1) to my random effects model: fm2 <- glm (BinaryOutcome ~ Age, family='binomial', data=DF1)

I tried the following code: lr.test (fm1, fm2)

which gives me some errors.

Anyone any idea how to solve this issue?




Aucun commentaire:

Enregistrer un commentaire