performs a correspondence analysis.
dudi.coa(df, scannf = TRUE, nf = 2)
a data frame containing positive or null values
a logical value indicating whether the eigenvalues bar plot should be displayed
if scannf FALSE, an integer indicating the number of kept axes
returns a list of class coa
and dudi
(see dudi) containing
the sum of all the values of the initial table
Benz<U+00E9>cri, J.P. and Coll. (1973) L'analyse des donn<U+00E9>es. II L'analyse des correspondances, Bordas, Paris. 1--620.
Greenacre, M. J. (1984) Theory and applications of correspondence analysis, Academic Press, London.
# NOT RUN {
data(rpjdl)
chisq.test(rpjdl$fau)$statistic
rpjdl.coa <- dudi.coa(rpjdl$fau, scannf = FALSE, nf = 4)
sum(rpjdl.coa$eig)*rpjdl.coa$N # the same
if(adegraphicsLoaded()) {
g1 <- s.label(rpjdl.coa$co, plab.cex = 0.6, lab = rpjdl$frlab, plot = FALSE)
g2 <- s.label(rpjdl.coa$li, plab.cex = 0.6, plot = FALSE)
cbindADEg(g1, g2, plot = TRUE)
} else {
par(mfrow = c(1,2))
s.label(rpjdl.coa$co, clab = 0.6, lab = rpjdl$frlab)
s.label(rpjdl.coa$li, clab = 0.6)
par(mfrow = c(1,1))
}
data(bordeaux)
db <- dudi.coa(bordeaux, scan = FALSE)
db
score(db)
# }
Run the code above in your browser using DataLab