Mercurial > repos > davidvanzessen > shm_csr

changeset 22:0bea8c187a90 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Mon, 28 Nov 2016 10:27:22 -0500
parents c281993e515b
children 81453585dfc3
files aa_histogram.r
diffstat 1 files changed, 4 insertions(+), 3 deletions(-) [+]
line wrap: on
line diff
--- a/aa_histogram.r	Mon Nov 28 10:24:25 2016 -0500
+++ b/aa_histogram.r	Mon Nov 28 10:27:22 2016 -0500
@@ -50,9 +50,6 @@
 
         print("---------------- write/print ----------------")
 
-        png(filename=paste(outdir, "/aa_histogram_", gene, ".png", sep=""), width=1280, height=720)
-        print(m)
-        dev.off()
 
         dat.sums = data.frame(index=1:length(mutations.at.position), mutations.at.position=mutations.at.position, aa.at.position=aa.at.position)
 
@@ -60,4 +57,8 @@
         write.table(mutations.by.id.gene, paste(outdir, "/aa_histogram_count_", gene, ".txt", sep=""), sep="\t",quote=F,row.names=F,col.names=T)
         write.table(absent.aa.by.id.gene, paste(outdir, "/aa_histogram_absent_", gene, ".txt", sep=""), sep="\t",quote=F,row.names=F,col.names=T)
         write.table(dat_dt, paste(outdir, "/aa_histogram_", gene, ".txt", sep=""), sep="\t",quote=F,row.names=F,col.names=T)
+        
+        png(filename=paste(outdir, "/aa_histogram_", gene, ".png", sep=""), width=1280, height=720)
+        print(m)
+        dev.off()
 }