Mercurial > repos > davidvanzessen > mutation_analysis

comparison aa_histogram.r @ 102:e6bc976760d4 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Tue, 21 Jun 2016 03:32:50 -0400
parents 564c4f6da203
children 603a10976e9c
comparison
equal deleted inserted replaced
101:3cffb8a38bb1 102:e6bc976760d4
42 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") 42 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")
43 m = m + expand_limits(y=c(-0.1,1)) + xlab("AA position") + ylab("Frequency") + ggtitle("AA mutation frequency") 43 m = m + expand_limits(y=c(-0.1,1)) + xlab("AA position") + ylab("Frequency") + ggtitle("AA mutation frequency")
44 44
45 print("---------------- write/print ----------------") 45 print("---------------- write/print ----------------")
46 46
47 print("writing dat_dt") #need this
47 write.table(dat_dt, paste(dirname(outfile), "/aa_histogram.txt", sep=""), sep="\t",quote=F,row.names=F,col.names=T) 48 write.table(dat_dt, paste(dirname(outfile), "/aa_histogram.txt", sep=""), sep="\t",quote=F,row.names=F,col.names=T)
48 49 print("writing png") #also need this, file is haunted
49 png(filename=outfile, width=1280, height=720) 50 png(filename=outfile, width=1280, height=720)
50 print(m) 51 print(m)
51 dev.off() 52 dev.off()