# HG changeset patch # User xuebing # Date 1333234556 14400 # Node ID e9bc2ab48dd549001e241b0f04f541701cfa197b # Parent a65882e1636cb34e6742f22d36ef911647e40177 Uploaded diff -r a65882e1636c -r e9bc2ab48dd5 bed_size_stat.py --- a/bed_size_stat.py Sat Mar 31 18:52:24 2012 -0400 +++ b/bed_size_stat.py Sat Mar 31 18:55:56 2012 -0400 @@ -14,9 +14,9 @@ rf.write("pdf('"+outputfile+"')\n") if log == 'log': rf.write("len <- log10(len+1)\n") - rf.write("hist(len,breaks=100,xlab='interval size',main=paste('mean=',mean(len),sep=''))\n") + rf.write("hist(len,breaks=100,xlab='interval size (log10)',main=paste('mean=',mean(len),sep=''))\n") else: - rf.write("hist(len,breaks=100,xlab='interval size (log10)',main=paste('mean=',mean(len),sep=''))\n") + rf.write("hist(len,breaks=100,xlab='interval size',main=paste('mean=',mean(len),sep=''))\n") rf.write("dev.off()") rf.close() os.system("R --vanilla < tmp.r")