d <- read.table("inp_iris.txt",header=TRUE) mm <- ncol(d)-1 dat <- d[,1:mm] rsl <- princomp(dat,cor=TRUE) plt <- rsl$scores pairs( plt, pch=21, bg=c("red","green","blue")[unclass(d$Species)] ) (rsl$sdev)^2 rsl$loadings