| Previous changeset 0:26b8fea8a9e7 (2012-03-31) Next changeset 2:b0e04c112633 (2012-03-31) |
|
Commit message:
Uploaded |
|
added:
bed_size_stat.py |
| b |
| diff -r 26b8fea8a9e7 -r b1da797d49a5 bed_size_stat.py --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/bed_size_stat.py Sat Mar 31 14:02:27 2012 -0400 |
| [ |
| @@ -0,0 +1,18 @@ +''' +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') |