Mercurial > repos > elixir-it > vinyl_pca

view PCA.R @ 1:460883beb10c draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author elixir-it
date Wed, 15 Jul 2020 07:55:07 +0000
parents 0011da72f65a
children
line wrap: on
line source

args <- commandArgs(trailingOnly = TRUE)
fileIn=args[1]
nhealth=args[2]
ndis=args[3]
fileO=args[4]
Final_D=read.table(fileIn,header=T,row.names=1)
PCA=prcomp(t(Final_D))
M=PCA$x
png(fileO,width=1800,height=1800,res=160)
plot(M[,1],M[,2],col=rep(c("orange","purple"),c(nhealth,ndis)),main="PCA of gene scores",xlab="PC1",ylab="PC2",cex.main=2,cex.lab=2,pch=20,cex.axis=1.5,cex=1.5)
legend(min(M[,1]),min(M[,2])*0.8,fill=c("orange","purple"),legend=c("affected","unaffected"),cex=1.5)
dev.off()