[R-sig-ME] Trying to align lmer output with Robinson's BLUP is a good thing

[Ravi Varadhan]

Hi,

I am trying to understand the model fitting of lmer() with the aid of GK Robinson's 1991 paper "BLUP is a good thing."  I am trying to align the lmer() output with Robinson's notation.

I have a few questions in this regard.

  1.  How can I get matrices X and Z from lmer?
  2.  How do I obtain the variances G and R?
  3.  For a new subject (x*, z*), is it possible to obtain BLUP predictions of x*\beta + z*u as shown in Robinson's paper? When I try to do this using predict.merMod, it complains that new levels are not allowed for conditional prediction.

Thank you,
Ravi


	[[alternative HTML version deleted]]