view population_structure.r @ 1:0df7a9b89f13

add empty snifffers section to mitigate toolshed bug; new version of dpmix
author Richard Burhans <burhans@bx.psu.edu>
date Mon, 09 Apr 2012 12:40:57 -0400
parents 2c498d40ecde
children
line wrap: on
line source

library(RColorBrewer)

args = commandArgs(trailingOnly=TRUE)
q_file = args[[1]]
output_file = args[[2]]
populations = args[[3]]

tbl <- read.table(q_file)

if ( populations >= 3 && populations <= 12 ) {
    colors = brewer.pal(populations, 'Paired')
} else {
    colors = rainbow(populations)
}

pdf(file=output_file, onefile=TRUE, width=7, height=3)
barplot(t(as.matrix(tbl)), col=colors, xlab="Individual #", ylab="Ancestry", border=NA)

dev.off()