03 2017 档案

scripts of perl
该文被密码保护。

posted @ 2017-03-28 03:22 Vermont 阅读(14) 评论(0) 推荐(0)

parallel::ForkManager
摘要：use Parallel::ForkManager; my $MAX_PROCESSES=10; #申明最大进程数（一次创建的进程越多，越耗内存）； my $pm = new Parallel::ForkManager($MAX_PROCESSES); #创建一个实例； #系统中始终保持处理的进程数 阅读全文

posted @ 2017-03-22 00:33 Vermont 阅读(1277) 评论(0) 推荐(0)

NeedlenamWunsch alignment in perl
该文被密码保护。

posted @ 2017-03-21 11:09 Vermont 阅读(6) 评论(0) 推荐(0)

sort_gff.py
摘要：import sysinfile = sys.argv[1]outfile = sys.argv[2]gff_list = []fh = open(infile)for x in fh: if x[0] == '#': continue y = x.strip().split('\t') gff_l 阅读全文

posted @ 2017-03-18 03:23 Vermont 阅读(172) 评论(0) 推荐(0)

Tajima's D
摘要：Three ways to assess the nucleotide diversity (heterozygosity).The first is mean pairwise difference ∏, the average number of differences between pair 阅读全文

posted @ 2017-03-09 05:20 Vermont 阅读(1379) 评论(0) 推荐(0)

calculate TajimaD in perl
摘要：#!/usr/bin/perl use strict; use warnings; =pod this perl script is used to compute tajima's D Former of tajimaD.tmp.txt is : chr position #sampled #de 阅读全文

posted @ 2017-03-09 01:17 Vermont 阅读(845) 评论(0) 推荐(0)

ggplot2 multiply graphs on one figure
摘要：library(ggplot2)library(grid)grid.newpage()pushViewport(viewport(layout = grid.layout(1,2)))vplayout = function(x, y) viewport(layout.pos.row = x, lay 阅读全文

posted @ 2017-03-06 04:57 Vermont 阅读(182) 评论(0) 推荐(0)

manhattan plots in qqplot2
摘要：###manhattan plots in qqplot2library(ggplot2)setwd("~/ncbi/zm/XPCLR/")read.table("LW.gene.xpclr.position.txt",header=F,sep=',')->ILhead(IL)#V1 V2 V3 V 阅读全文

posted @ 2017-03-06 01:00 Vermont 阅读(309) 评论(0) 推荐(0)

导航