Mercurial > repos > miller-lab > genome_diversity
diff population_structure.r @ 17:a3af29edcce2
Uploaded Miller Lab Devshed version a51c894f5bed
author | miller-lab |
---|---|
date | Fri, 28 Sep 2012 11:57:18 -0400 |
parents | 2c498d40ecde |
children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/population_structure.r Fri Sep 28 11:57:18 2012 -0400 @@ -0,0 +1,19 @@ +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()