[R-sig-ME] predict in lmer()?
David Winsemius
dwinsemius at comcast.net
Sun Jul 31 15:48:40 CEST 2011
On Jul 31, 2011, at 9:35 AM, David Winsemius wrote:
>
> On Jul 31, 2011, at 8:47 AM, Dennis Murphy wrote:
>
>> Hi:
>>
>> See http://glmm.wikidot.com/faq
>>
>> Go about 2/3 of the way down until you see the section 'Predictions
>> and/or confidence (or prediction) intervals on predictions'.
>
> Aren't objects created by lmer of class "mer"? Attempting to follow
> your advice with the first example provided in ?lmer meets with
> frustration:
>
> require(lme4)
> (fm1 <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy))
> newdat0 <- expand.grid(Reaction=c(200,300,400), Days=c(0,4,8),
> Subject=c(5,10,15) )
> newdat0$pred <- predict(fm1, newdat0, level = 0)
>
> Error in UseMethod("predict") :
> no applicable method for 'predict' applied to an object of class
> "mer"
>
> I'm not even a journeyman use of ME methods, so this is just a
> question. Do you have a fix? (There are predict methods listed in
> nlme but none in the Index of lme4.
And then I scrolled a bit further down beyond the lme section to where
the lme4 section was . Sorry for the noise. And thanks for the
addition to my meager lme4 education.
--
David.
David Winsemius, MD
West Hartford, CT
More information about the R-sig-mixed-models
mailing list