On Thu, 2010-01-21 at 21:56 +0800, echo_july wrote:> hello,
> I am doing PCA in R using some habitat factors, and I used the
> function result1=rda() and result2=princomp(),then pick up scores of
> the result1 and result2 using scores(),but the scores are
> significantly different,i do not know the meaning of it.
> Best wishes!
Because scores.rda allows for scaling of the scores to emphasise aspects
of the data whereas the default scores method (which you are using to
extract the loadings from princomp) doesn't.
Are they the same if you do:
result1 <- rda(....)
scores(result1, scaling = 0)
HTH
G
--
%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%
Dr. Gavin Simpson [t] +44 (0)20 7679 0522
ECRC, UCL Geography, [f] +44 (0)20 7679 0565
Pearson Building, [e] gavin.simpsonATNOSPAMucl.ac.uk
Gower Street, London [w] http://www.ucl.ac.uk/~ucfagls/
UK. WC1E 6BT. [w] http://www.freshwaters.org.uk
%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%