Hi, all:
I am using extract.lme.cov to extract the covariance matrix of lme. But
the results are not expected.
For example,
b <- lme(travel~1,Rail,~1|Rail)
The default correlation for lme is no correlation within groups.
>extract.lme.cov(b,Rail)
The part of covariance matrix looks like:
1 2 3 4 5 6
1 631.4778 615.3111 615.3111 0.0000 0.0000 0.0000
2 615.3111 631.4778 615.3111 0.0000 0.0000 0.0000
3 615.3111 615.3111 631.4778 0.0000 0.0000 0.0000
4 0.0000 0.0000 0.0000 631.4778 615.3111 615.3111
5 0.0000 0.0000 0.0000 615.3111 631.4778 615.3111
6 0.0000 0.0000 0.0000 615.3111 615.3111 631.4778
Where does the covariance come from? Maybe I'm missing sth here?
Thanks,
Hua Li
[[alternative HTML version deleted]]
______________________________________________
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
