comparison population_structure.r @ 12:4b6590dd7250

Uploaded
author miller-lab
date Wed, 12 Sep 2012 17:10:26 -0400
parents 2c498d40ecde
children
comparison
equal deleted inserted replaced
11:d4ec09e8079f 12:4b6590dd7250
1 library(RColorBrewer)
2
3 args = commandArgs(trailingOnly=TRUE)
4 q_file = args[[1]]
5 output_file = args[[2]]
6 populations = args[[3]]
7
8 tbl <- read.table(q_file)
9
10 if ( populations >= 3 && populations <= 12 ) {
11 colors = brewer.pal(populations, 'Paired')
12 } else {
13 colors = rainbow(populations)
14 }
15
16 pdf(file=output_file, onefile=TRUE, width=7, height=3)
17 barplot(t(as.matrix(tbl)), col=colors, xlab="Individual #", ylab="Ancestry", border=NA)
18
19 dev.off()