[R] How to extract R{i} from lme object?

Peng l12345p2000 at yahoo.com
Tue Mar 4 21:38:16 CET 2003

Hi, lme() users,

Can some one tell me how to do this.
I model Orthodont with the same G for random
variables, but different R{i}'s for boys and girls, so
that I can get sigma1_square_hat for boys and
sigma2_square_hat for girls.

The model is Y{i}=X{i}beta + Z{i}b + e{i}
b ~ iid N(0,G) and e{i} ~ iid N(0,R{i}) i=1,2
orth.lme <- lme(distance ~ Sex * age, data=Orthodont,
random=~age|Subject, weights=varIdent(form=~1|Sex),

I can see the numbers I need from summary(), but how
can I extract them? I tried several functions in nlme,
but I cannot find a correct one.


Peng Liu
Peng Liu                      |
Division of Statistics        |
Northern Illinois University  |
De Kalb, IL 60115, USA        |
E-mail: pliu at math.niu.edu     |

More information about the R-help mailing list