library(FactoMineR) te <- read.table("http://factominer.free.fr/libra/te.csv",header=TRUE,sep=";") n.grupos <- 4 grupos <- quantile(te[,22], seq(0,1,1/n.grupos)) Xqual <- cut(te[,22],grupos, include.lowest=TRUE) summary(Xqual) hist(te$edad,col="grey",main="Histograma of the variable edad",freq=FALSE, xlab="edad", nclass=15) vari <- te[,22] res.hcpc <- HCPC(vari, iter.max=10) max.cla = unlist(by(res.hcpc$data.clust[,1],res.hcpc$data.clust[,2],max)) breaks=c(min(vari),max.cla) aa.cuali = cut(vari, breaks, include.lowest=TRUE) summary(aa.cuali) datos.cuali <- datos for (i in 1:ncol(datos.cuali)){ vari = datos.cuali[,i] res.hcpc=HCPC(vari, nb.clust=-1, graph=FALSE) maxi = unlist(by(res.hcpc$data.clust[,1], res.hcpc$data.clust[,2],max)) breaks=c(min(vari),maxi) aa.cuali = cut(vari, breaks, include.lowest=TRUE) datos.cuali[,i] = aa.cuali }