Repository 'shm_csr'
hg clone https://toolshed.g2.bx.psu.edu/repos/davidvanzessen/shm_csr

Changeset 22:0bea8c187a90 (2016-11-28)
Previous changeset 21:c281993e515b (2016-11-28) Next changeset 23:81453585dfc3 (2016-12-01)
Commit message:
Uploaded
modified:
aa_histogram.r
b
diff -r c281993e515b -r 0bea8c187a90 aa_histogram.r
--- a/aa_histogram.r Mon Nov 28 10:24:25 2016 -0500
+++ b/aa_histogram.r Mon Nov 28 10:27:22 2016 -0500
b
@@ -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()
 }