Mercurial > repos > xuebing > sharplabtool
view mytools/intervalSize.py @ 9:87eb5c5ddfe9
Uploaded
author | xuebing |
---|---|
date | Fri, 09 Mar 2012 20:01:43 -0500 |
parents | f0dc65e7f6c0 |
children |
line wrap: on
line source
''' plot histogram of interval size ''' import os,sys inputfile = sys.argv[1] outputfile = sys.argv[2] rf = open('tmp.r','w') rf.write("x <- read.table('"+inputfile+"')\n") rf.write("len <- x[,3]-x[,2]\n") rf.write("pdf('"+outputfile+"')\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") os.system('rm tmp.r')