[R] meaning of "iner" etc. relating to cca(ade4)
Stephane Dray
dray at biomserv.univ-lyon1.fr
Fri Jan 14 22:04:12 CET 2005
Hi Dominique,
please, if you have some questions about a particular package, do not send
to R-Help. There is a list about ade4, see details at:
http://pbil.univ-lyon1.fr/ADE-4/adelist.html
cca is a correspondence analysis with constraint imposed by external
information. Redundancy analysis (available in ade4 by the function pcaiv)
is PCA with constraint.
Here is an example using pcaiv:
data(rhone)
pca1 <- dudi.pca(rhone$tab, scan = FALSE, nf = 3)
iv1 <- pcaiv(pca1, rhone$disch, scan = FALSE)
iv1$param
iner inercum inerC inercumC ratio R2 lambda
6.27 6.27 5.52 5.52 0.879 0.671 3.7
4.14 10.4 4.74 10.3 0.984 0.747 3.54
iner represent eigenvalues of the unconstrained analysis:
> pca1$eig[1:2]
[1] 6.274265 4.140947
inercum is equivalent to
cumsum(pca1$eig)
prediction of sites scores by environmental variables produce a new set
score (often named LC scores), inerC represent variance of this scores:
diag(t(iv1$ls)%*%diag(iv1$lw)%*%as.matrix(iv1$ls))
Axis1 Axis2
5.515241 4.736362
inercumC is the cumsum.
ratio is simply inercum/inerCcum
Finally note that CCA or RDA find a site score that maximizes the variance
explained by environmental variables. This is a product of variance by a R2:
0.671*5.52 # lambda = R2*inerC
[1] 3.70392
At 10:11 14/01/2005, Dominique Grüter wrote:
>Does someone know the meaning of the following values relating to the
>"cca"-command?
>
> iner, inercum, inerC, inercumC
>
>Thanks!
>DOMI
>
>______________________________________________
>R-help at 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
Stéphane DRAY
--------------------------------------------------------------------------------------------------
Département des Sciences Biologiques
Université de Montréal, C.P. 6128, succursale centre-ville
Montréal, Québec H3C 3J7, Canada
Tel : (514) 343-6111 poste 1233 Fax : (514) 343-2293
E-mail : stephane.dray at umontreal.ca
--------------------------------------------------------------------------------------------------
Web http://www.steph280.freesurf.fr/
More information about the R-help
mailing list