[R-sig-ME] How to calculate predictions
Ben Bolker
bbolker at gmail.com
Tue Feb 28 05:26:16 CET 2012
Stuart Luppescu <slu at ...> writes:
>
> On Mon, 2012-02-27 at 16:07 -0700, Jake Westfall wrote:
>
> I believe that fitted() is giving fitted values for each observation.
>
> length(fitted(math.lme3))
> [1] 520573
> length(unique(math.lme3 <at> frame$sid))
> [1] 236994
>
> I have about 520,000 observations nested within 236,994 students. I want
> a fitted value for each student at age 15 (which may or may not be
> actually observed).
Have you looked at the code on http://glmm.wikidot.com/faq ... ?
Also, the new lme4Eigen package does have a ?predict method ...
Ben Bolker
More information about the R-sig-mixed-models
mailing list