Hi all,
 
Some of you may have seen my recent message to the r-help list
concerning a modification to lmer code that caused R to crash. Happily,
the same modifications to the lmer2 code do not cause this crash so as
far as I am concerned this is a major improvement to the original. Now,
however, I'm not sure how to get at the bits of the fitted object that I
want to play with to post-process the results as the methods aren't yet
implemented for this structure - I need, for example, the coefficients
of the fixed and random effects and the full variance-covariance matrix
of fixed _and_ random effects (together!). I'm sure there must be ways
around this issue, since I can see that there are slots with names that
look very suitable at least to hold the fixed and random effects - but I
can't work out how to get them. And I'm not sure if the full
variance-covariance matrix is going to be in there anyway (looks like it
might be for lmer objects). 
 
For example:
 
fm1 <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy)
fm2 <- lmer2(Reaction ~ Days + (Days|Subject), sleepstudy)
 
using:
 
 slotNames(fm1)
 slotNames(fm2)
 
shows me that there are slots with the ranef and fixef in both classes
of object - I can use fixef() on both objects to get the fixed effects,
but ranef() only works on fm1 at the moment - how else can I get these
objects directly? And there are obviously many more slots in the fm1
object than the fm2 object - can anyone point out if any of these slots
are the full variance-covariance matrix I need, containing both fixed
and random effects and if so, how I can get at it (or if not, how I
might go about creating it?). 
 
Thanks in advance,
 
Colin
 
 
Dr. Colin Beale
Spatial Ecologist
The Macaulay Institute
Craigiebuckler
Aberdeen
AB15 8QH
UK

Tel: 01224 498245 ext. 2427
Fax: 01224 311556
Email: c.beale@macaulay.ac.uk 


-- 
Please note that the views expressed in this e-mail are those of the
sender and do not necessarily represent the views of the Macaulay
Institute. This email and any attachments are confidential and are
intended solely for the use of the recipient(s) to whom they are
addressed. If you are not the intended recipient, you should not read,
copy, disclose or rely on any information contained in this e-mail, and
we would ask you to contact the sender immediately and delete the email
from your system. Thank you.
Macaulay Institute and Associated Companies, Macaulay Drive,
Craigiebuckler, Aberdeen, AB15 8QH.


	[[alternative HTML version deleted]]


