-
计算重组率
-
2015-1-16 11:23
-
##file为bed文件:chr,start,end file1为10条染色体的cM数据:position,cM import os,sys import bisect for file in os.listdir(/share/xulab_10t/lxl/IBD_20141224_recobination): if file.endswith(_trop_trop.bed): file1=file.repla ...
-
个人分类: 科研笔记|6971 次阅读|没有评论
-
一个对fasta文件进行操作的python code
-
2015-1-15 16:22
-
##将inf1作为字典,在inf2中进行搜索,将匹配上的gene的fasat序列提取出来! import os, sys for i in range(1,11): inf1=open(chr+str(i)+.gene.one,r) inf2=open(chr+str(i)+.raw.cds,r) ouf=open(chr+str(i)+.chr ...
-
个人分类: 科研笔记|4047 次阅读|没有评论
-
关于R中ggplot2画图的一些流水记录
-
2015-1-15 15:53
-
##形式1 data - read.table(chr6_trop_temp2.bed.hist.sort.filt,header=F) pdf(chr6_trop_temp2.bed.hist.sort.filt.pdf) plot(data ,data ,cex=0.005,pch=21,type=p,col=black,xlab=location,ylab=count,main=coveragebed plot) dev.off() ##形式2 ...
-
个人分类: 科研笔记|5960 次阅读|没有评论
本页有 4 篇博文因作者的隐私设置或未通过审核而隐藏