# HG changeset patch # User davidvanzessen # Date 1456827367 18000 # Node ID cb66d6dd1e6627af1579565c4647b8278ec44495 # Parent a4317b006d70fd6d3fa84321721c7abe6743eef0 Uploaded diff -r a4317b006d70 -r cb66d6dd1e66 aa_histogram.r --- a/aa_histogram.r Tue Mar 01 05:06:47 2016 -0500 +++ b/aa_histogram.r Tue Mar 01 05:16:07 2016 -0500 @@ -37,6 +37,9 @@ m = m + annotate("segment", x = 55.5, y = -0.07, xend=65.5, yend=-0.07, colour="darkblue", size=1) + annotate("text", x = 60.5, y = -0.15, label="CDR2") m = m + annotate("segment", x = 65.5, y = -0.05, xend=104.5, yend=-0.05, colour="darkgreen", size=1) + annotate("text", x = 85, y = -0.1, label="FR3") m = m + expand_limits(y=c(-0.1,1)) + xlab("AA position") + ylab("Frequency") + ggtitle("AA mutation frequency") + +print("---------------- write/print ----------------") + write.table(dat_dt, paste(dirname(outfile), "/aa_histogram.txt", sep=""), sep="\t",quote=F,row.names=F,col.names=T) png(filename=outfile, width=1280, height=720)