[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