Mercurial > repos > rico > my_400th_test
annotate population_structure.r @ 0:32f8fbdd754c default tip
Uploaded
| author | rico |
|---|---|
| date | Fri, 06 Apr 2012 13:59:29 -0400 |
| parents | |
| children |
| rev | line source |
|---|---|
| 0 | 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() |