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