Skip to content
Prev 15538 / 20628 Next

understanding I() in lmer formula

> > gives me
 > >  Warning message:
 > >  vignette 'lmer' not found
 > 
 > That's surprising ... what's packageVersion("lme4") ?
[1] '1.1.14'
[actually I'm replacing left and right non-ascii quotes above]

 > Pinheiro and Bates (2000),
I'll look, but it seems that for something that old, even if it does 
contain these details, they are likely to be out of date.
I noticed http://lme4.r-forge.r-project.org/lMMwR/lrgprt.pdf
which is from 2010 says on p.15
 The vector u is available in fm01ML at re.  The vector Beta and the model matrix 
 X are in fm01ML at fe. 
and this seems no longer true - I get
  Error: no slot of name "fe" for this object of class "lmerMod" 

 >   Yes.    (So in particular this model are 3*(3+1)/2=6 parameters
 > (equivalent to s^2{1}, s^2(x),s^2(y), cov(1,x), cov(1,y), cov(x,y))

That might explain why changing (x+y+z+u+v+w||g) to single | seems to
run a lot longer (I interrupted it so I don't know how much longer).
Do you have an estimate of how number of terms affects run time?

If a paper describes various things being used as random effects
would you assume | or || ?
Perhaps it depends on the number of effects?  The example I sent
contained quite a few, and the original from which that was extracted
contained over a dozen.