[R-sig-ME] How to calculate predictions

[Stuart Luppescu]

On Tue, 2012-02-28 at 04:26 +0000, Ben Bolker wrote:
> 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 ... ?

Yes, I have. Unless I'm missing something (which is very possible) it
doesn't seem that the code includes the individual random effects.

> Also, the new lme4Eigen package does have a ?predict method ...

Doesn't seem to be available on CRAN -- only on R-forge? I don't think I
can convince our sysadmin to install a non-stable package.

-- 
Stuart Luppescu -=- slu .at. ccsr.uchicago.edu        
University of Chicago -=- CCSR 
才文と智奈美の父 -=-    Kernel 3.2.1-gentoo-r2                
1st mail: One of the reasons that SAM is popular,
 is that it is popular (i.e. since everyone has
 heard of it, it makes reviewers happy). So, it
 would be nice to be able to point to publications 
 in good journals so that reviewers will be
 comfortable. (I personally, am quite comfortable
 with SAM). 2nd mail: Oops, must have been a
 Freudian slip. Actually, I am not perfectly
 comfortable with SAM. But I am quite comfortable