I am using S-Plus 8.0 (Windows), and the function “lme”
to fit a model with fixed and random components. The function works well. I
call the resulting object a.lme.
Then I define s.lme = summary(a.lme). Printing that gives me
what I need. However, I would like to simulate this many times, generating a
new y vector each time, and for each such iteration store the resulting
parameter estimates. I can do this for the estimates of the parameters (the betas) in the fixed part. However, neither of
the lists a.lme and s.lme contain the estimates of the standard deviations (or
variance) estimates, even though they are printed when I enter
s.lme or print(s.lme) ( which are the same).
How can I “grab” these estimates so that I can
store them?
Jim Stapleton