#CORRESPONDENCE ANALYSIS
M=read.table("c:/DATA/Multivariate/T12-10.DAT")
M
attach(M)
P=prop.table(M) # CREATES CORRESPONDENCE MATRIX P
P

library(ca)  # INSTALL LIBRARY {ca} from CRAN

?ca
CA <- ca(P)
CA
summary(CA) 

?plot.ca
plot(CA, arrows = c(FALSE, TRUE))


#PLOTTING SUBSETS OF POINTS
CA2=ca(P,suprow=c(1,3,5))
plot(CA2,what=c("passive","all"),pch=c(20,19),col=c("red","blue"))

#USING PACKAGE {vegan}:
library(vegan)

CCA=cca(M)
CCA
plot(CCA)