摘要: > library(VennDiagram)载入需要的程辑包:grid> dg1<-read.table("clipboard",header=F)> dg2<-read.table("clipboard",header=F)> venn.diagram(list(dg_1_nip=dg1[,1],dg_2_nip=dg2[,1]),fill=c("red","blue"),"rice.tiff") 阅读全文
posted @ 2012-12-24 15:41 蒋蒋 阅读(302) 评论(0) 推荐(0) 编辑
摘要: #!/bin/bash#---------cutting file----------IN=/share1/projects/project_mwr/tophat_result#filename=$IN/Sample_ZJ-GFP-/countlist.txtfilename=$IN/Sample_ZJ-GFP+/countlist.txt#---------cuttint last lin -----cutnu=5cutnum=`expr $cutnu - 1`filetmp=`cat $filename|wc -l`filelen=`expr $filetmp - $cutnum`file 阅读全文
posted @ 2012-12-24 10:54 蒋蒋 阅读(198) 评论(0) 推荐(0) 编辑