Mercurial > repos > miller-lab > genome_diversity
comparison population_structure.r @ 0:2c498d40ecde
Uploaded
| author | miller-lab |
|---|---|
| date | Mon, 09 Apr 2012 12:03:06 -0400 |
| parents | |
| children |
comparison
equal
deleted
inserted
replaced
| -1:000000000000 | 0:2c498d40ecde |
|---|---|
| 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() |