if (FALSE) {
data(rpjdl)
coa1 <- dudi.coa(rpjdl$fau, scann = FALSE)
pca1 <- dudi.pca(rpjdl$fau, scal = FALSE, scann = FALSE)
if(adegraphicsLoaded()) {
g1 <- s.distri(coa1$l1, rpjdl$fau, xax = 2, yax = 1, starSize = 0.3,
ellipseSize = 0, plab.cex = 0)
g2 <- s.distri(pca1$l1, rpjdl$fau, xax = 2, yax = 1, starSize = 0.3,
ellipseSize = 0, plab.cex = 0)
} else {
s.distri(coa1$l1, rpjdl$fau, 2, 1, cstar = 0.3, cell = 0)
s.distri(pca1$l1, rpjdl$fau, 2, 1, cstar = 0.3, cell = 0)
}
caiv1 <- pcaiv(coa1, rpjdl$mil, scan = FALSE)
plot(caiv1)
}
Run the code above in your browser using DataLab