Mercurial > repos > davidvanzessen > argalaxy_tools
comparison aa_histogram.r @ 11:0510cf1f7cbc draft
Uploaded
| author | davidvanzessen |
|---|---|
| date | Tue, 04 Aug 2015 09:59:26 -0400 |
| parents | |
| children |
comparison
equal
deleted
inserted
replaced
| 10:edbf4fba5fc7 | 11:0510cf1f7cbc |
|---|---|
| 1 library(ggplot2) | |
| 2 | |
| 3 args <- commandArgs(trailingOnly = TRUE) | |
| 4 | |
| 5 input = args[1] | |
| 6 outfile = args[2] | |
| 7 | |
| 8 dat = read.table(input, header=F, sep=",") | |
| 9 dat=as.numeric(dat[1,]) | |
| 10 dat_sum = sum(dat) | |
| 11 | |
| 12 dat_freq = dat / dat_sum * 100 | |
| 13 dat_dt = data.frame(i=1:length(dat_freq), freq=dat_freq) | |
| 14 | |
| 15 m = ggplot(dat_dt, aes(x=i, y=freq)) + theme(axis.text.x = element_text(angle = 90, hjust = 1)) | |
| 16 m = m + geom_histogram(stat="identity", colour = "black", fill = "darkgrey", alpha=0.8) + scale_x_continuous(breaks=1:length(dat_freq), labels=1:length(dat_freq)) | |
| 17 m = m + annotate("segment", x = 0.5, y = -0.3, xend=26.5, yend=-0.3, colour="darkgreen", size=1) + annotate("text", x = 13, y = -0.2, label="FR1") | |
| 18 m = m + annotate("segment", x = 26.5, y = -0.4, xend=38.5, yend=-0.4, colour="darkblue", size=1) + annotate("text", x = 32.5, y = -0.3, label="CDR1") | |
| 19 m = m + annotate("segment", x = 38.5, y = -0.3, xend=55.5, yend=-0.3, colour="darkgreen", size=1) + annotate("text", x = 47, y = -0.2, label="FR2") | |
| 20 m = m + annotate("segment", x = 55.5, y = -0.4, xend=65.5, yend=-0.4, colour="darkblue", size=1) + annotate("text", x = 60.5, y = -0.3, label="CDR2") | |
| 21 m = m + annotate("segment", x = 65.5, y = -0.3, xend=104.5, yend=-0.3, colour="darkgreen", size=1) + annotate("text", x = 85, y = -0.2, label="FR3") | |
| 22 write.table(dat_dt, paste(dirname(outfile), "/aa_histogram.txt", sep=""), sep="\t",quote=F,row.names=F,col.names=T) | |
| 23 | |
| 24 png(filename=outfile, width=1280, height=720) | |
| 25 print(m) | |
| 26 dev.off() |