[R] random effect question and glm

Aimin Yan aiminy at iastate.edu
Thu Nov 23 22:52:16 CET 2006


consider p as random effect with 5 levels, what is difference between these
two models?

  > p5.random.p <- lmer(Y
~p+(1|p),data=p5,family=binomial,control=list(usePQL=FALSE,msV=1))
  > p5.random.p1 <- lmer(Y
~1+(1|p),data=p5,family=binomial,control=list(usePQL=FALSE,msV=1))

in addtion, I try these two models, it seems they are same.
what is the difference between these two model. Is this a cell means model?
m00 <- glm(sc ~aa-1,data = p5)
m000 <- glm(sc ~1+aa-1,data = p5)


thanks,

Aimin Yan



More information about the R-help mailing list