[R] random effect question and glm

Michael Dewey info at aghmed.fsnet.co.uk
Fri Nov 24 13:13:42 CET 2006


At 21:52 23/11/2006, Aimin Yan wrote:
>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))

Well, try them and see. Then if you cannot understand the output tell us
a) what you found
b) how it differed from what you expected

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

See above


>thanks,
>
>Aimin Yan
>
>

Michael Dewey
http://www.aghmed.fsnet.co.uk



More information about the R-help mailing list