costatis {ade4} | R Documentation |
Does the analysis of a series of pairs of ecological tables. This function uses Partial Triadic Analysis (pta) and coinertia to do the computations.
costatis(KTX, KTY, scannf = TRUE)
KTX |
an objet of class ktab |
KTY |
an objet of class ktab |
scannf |
a logical value indicating whether the eigenvalues bar plot should be displayed |
This function takes 2 ktabs. It does a PTA (partial triadic analysis: pta) on each ktab, and does a coinertia analysis (coinertia) on the compromises of the two PTAs.
a list of class coinertia, subclass dudi. See coinertia
IMPORTANT : KTX and KTY must have the same k-tables structure, the same number of columns, and the same column weights.
Jean Thioulouse jean.thioulouse@univ-lyon1.fr
Thioulouse J., Simier M. and Chessel D. (2004). Simultaneous analysis of a sequence of paired ecological tables. Ecology 85, 272-283..
Simier, M., Blanc L., Pellegrin F., and Nandris D. (1999). Approche simultanée de K couples de tableaux : Application a l'étude des relations pathologie végétale - environnement. Revue de Statistique Appliquée, 47, 31-46.
data(meau) wit1 <- withinpca(meau$mil, meau$plan$dat, scan = FALSE, scal = "total") pcafau <- dudi.pca(meau$fau, scale = FALSE, scan = FALSE, nf = 2) wit2 <- within(pcafau, meau$plan$dat, scan = FALSE, nf = 2) kta1 <- ktab.within(wit1, colnames = rep(c("S1","S2","S3","S4","S5","S6"), 4)) kta2 <- ktab.within(wit2, colnames = rep(c("S1","S2","S3","S4","S5","S6"), 4)) costatis1 <- costatis(kta1, kta2, scan = FALSE) plot(costatis1)