data = read.csv("descritiva.csv",header=TRUE) attach(data) summary(data) X = data[,1:5] apply(X,2,mean) var(X) cor(X) boxplot(X) pairs(X) par(mfrow=c(2,3)) for (i in 1:5){ boxplot(data[,i]~data[,6],axes=FALSE,xlab="Sexo",ylab=paste("X",i,sep="")) axis(2);box();axis(1,at=1:2,lab=c("F","M")) } pairs(X,col=1+data[,6],pch=16)