Mercurial > repos > xuebing > bed_size_stat

comparison bed_size_stat.py @ 5:e9bc2ab48dd5 default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author xuebing
date Sat, 31 Mar 2012 18:55:56 -0400
parents cce2eabe8b8d
children
comparison
equal deleted inserted replaced
4:a65882e1636c 5:e9bc2ab48dd5
12 rf.write("x <- read.table('"+inputfile+"')\n") 12 rf.write("x <- read.table('"+inputfile+"')\n")
13 rf.write("len <- x[,3]-x[,2]\n") 13 rf.write("len <- x[,3]-x[,2]\n")
14 rf.write("pdf('"+outputfile+"')\n") 14 rf.write("pdf('"+outputfile+"')\n")
15 if log == 'log': 15 if log == 'log':
16 rf.write("len <- log10(len+1)\n") 16 rf.write("len <- log10(len+1)\n")
17 rf.write("hist(len,breaks=100,xlab='interval size (log10)',main=paste('mean=',mean(len),sep=''))\n")
18 else:
17 rf.write("hist(len,breaks=100,xlab='interval size',main=paste('mean=',mean(len),sep=''))\n") 19 rf.write("hist(len,breaks=100,xlab='interval size',main=paste('mean=',mean(len),sep=''))\n")
18 else:
19 rf.write("hist(len,breaks=100,xlab='interval size (log10)',main=paste('mean=',mean(len),sep=''))\n")
20 rf.write("dev.off()") 20 rf.write("dev.off()")
21 rf.close() 21 rf.close()
22 os.system("R --vanilla < tmp.r") 22 os.system("R --vanilla < tmp.r")
23 os.system('rm tmp.r') 23 os.system('rm tmp.r')