[R-sig-ME] problem with lme4
Ben Bolker
bbolker at gmail.com
Wed Jan 22 19:30:42 CET 2014
[cc'ing to r-sig-mixed-models]
You haven't given a reproducible example, but I believe that simply
replacing model.s2 at X with getME(model.s2,"X") in the second line
should work. I think you should also be able to get the same results
(predictions with and without random values) via the 'REform' argument
to predict() (possibly switching to 'ReForm' more recently): see
`?predict.merMod
On 14-01-22 12:59 PM, Linda Heeb wrote:
> Dear Mr Bolker
>
> Since the new lme4 package this script doesn’t work any more….
>
> model.s2<-glmer(y~licking.duration+CN+(1|nest)+(1|session),data=bold.t, family=poisson)
> bold.t$F0 <- inv.logit(model.s2 at X %*% fixef(model.s2))*100
> bold.t$F1 <- fitted(model.s2)*100
> p <- qplot(licking.duration, interactions*100, facets=.~CN, data=bold.t)
> p + theme_bw() + geom_line(aes(y=F0)) + geom_line(aes(y=F1, col=mumID)) + opts(legend.position = "none")
>
>
> Fehler in plogis(x) :
> kein Slot des Namens "X" für dieses Objekt der Klasse „glmerMod“
>
>
> I hope you can tell my why and show me how I can improve it.
>
> Cheers, Linda Heeb
>
>
More information about the R-sig-mixed-models
mailing list