Hi,
I've noticed that I get different output from glmer and pglmm, when I think(?) the output should be more similar (example below). I imagine that I am specifying the models differently in some way, but it's not clear to me how. Could you possibly provide some guidance?
Thank you,
Neil
Hi, I've noticed that I get different output from glmer and pglmm, when I think(?) the output should be more similar (example below). I imagine that I am specifying the models differently in some way, but it's not clear to me how. Could you possibly provide some guidance? Thank you, Neil
library(phyr) library(lme4)
some example data
s <- c(11, 12, 0, 1, 4, 4, 2, 3, 24, 14, 5, 6, 7, 15, 14, 9, 0, 0, 2, 11) # successes f <- c(2, 2, 2, 2, 4, 6, 5, 3, 1, 1, 6, 3, 11, 3, 1, 0, 1, 5, 2, 3) # failures ID<-c("A", "B", "C", "D", "E", "F", "G", "H", "I", "J", "K", "L", "M", "N", "P", "Q", "R", "S", "T", "U")
dat.long<-data.frame(s,f,ID)
glmer
m1<-glmer(cbind(s,f) ~ 1 + (1|ID),data=dat.long,family="binomial")
PGLMM
m2<-pglmm(cbind(s,f) ~ 1 + (1|ID), data=dat.long,family="binomial", add.obs.re = FALSE,REML = FALSE) summary(m1) summary(m2) logLik(m1) m2$logLik