args <- commandArgs() fnameR <- args[6] d <- read.table(fnameR,skip=1) mm <- ncol(d)-1 dat <- d[,1:mm] cl <- d[,mm+1] #Standardization #dat <- sweep(dat,2,apply(dat,2,mean),FUN="-") #dat <- sweep(dat,2,apply(dat,2,sd),FUN="/") plt <- as.matrix(dat[,1:mm]) %*% eigen(cor(dat[,1:mm]))$vectors[1:mm,] panel.text<-function(x, y, labels, cex, font, ...){ usr <- par("usr"); on.exit(par(usr)) par(usr = c(0,1,0,1) ) txt <- c("PCM") cex<-2 font<-"Helvetica" text(0.5, 0.5, txt, cex = cex) } pairs( plt, diag.panel=NULL, text.panel=panel.text, pch=21, bg=c("red","blue","yellow","cyan","green","black")[unclass(cl)] ) eigen(cor(dat[,1:mm]))$values[1:mm]