I guess there is a mistake in the function princomp (package mva).
I believe that the scores are computed wrongly under cor=TRUE condition:
Lines 42-43:
if (scores && !missing(x))
scr <- scale(z, center = TRUE, scale = FALSE) %*% edc$vectors
should be replaced by something like
if (scores && !missing(x) && !cor)
scr <- scale(z, center = TRUE, scale = FALSE) %*% edc$vectors
if (scores && !missing(x) && cor)
scr <- scale(z) %*% edc$vectors
Sincerely
************************************************
Ya'acov Ritov
Department of Statistics
The Hebrew University of Jerusalem,
Jerusalem 91905
Israel
Tel. (972-2) 652-9929 (Okay way)
Fax (972-2) 652-9996
http://pluto.mscc.huji.ac.il/~yaacov/jr.html
AND THIS IS ANOTHER PROOF THAT IF YOU KEEP AN OPEN MIND, SOMEBODY WILL FILL IT
WITH RUBBISH (an anonymous BBC correspondent)
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !) To:
r-devel-request@stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._