Three ways to assess the nucleotide diversity (heterozygosity).
The first is mean pairwise difference ∏, the average number of differences between pairs of sequences.
The second is number of segregating sites S. A segregating site is a site that is polymorphic in the data. the number of such sites is usually denoted by S (sometimes using K).
Last one is number of alleles Na.

Mean pairwise difference per nucleotide pi =  ∏/L, where L is the length of sequences.

Tajima's D test the neutral mutation hypothesis θ=4Neμ=θT=θW. D=[(θT–θW)/Var(θT–θW)] , where θw=S/a, a=[1+1/2+1/3+...+1/(n–1)] (Watterson 1975), θt= ∏ (Tajima 1983)

The following content comes from https://en.wikipedia.org/wiki/Tajima's_D

Tajima's D的更多相关文章

  1. .net MVC全局定时器执行作业

    首先的一个需求是在OA系统中定时跑一些定时作业,例如发放年假等事务,之前的做法是在服务器上加入一个服务,用系统定时作业去跑服务,这样有个问题就是当系统在发布的过程中,有可能忘记启动服务而导致无法定时执 ...

  2. calculate TajimaD in perl

    #!/usr/bin/perl use strict; use warnings; =pod--------------------------------------- this perl scri ...

  3. GATK--使用转载

    http://blog.sciencenet.cn/blog-1469385-819498.html 文章目录 一.准备工作 二.流程概览 三.流程 首先说说GATK可以做什么.它主要用于从seque ...

  4. 收集vcftools所有用法

    VCFtools用来处理VCF文档. 筛选特定突变 比较文件 总结突变 转化文件格式 验证并合并文件 取突变交集和差集 Get basic file statistics input可以为VCF或BC ...

  5. VCFtools

    The C++ executable module examples This page provides usage examples for the executable module. Exte ...

  6. SR4R数据库：水稻4个SNP集的筛选及其应用

    目录 前言 四个SNP集 hapmapSNPs tagSNPs fixedSNPs barcodeSNPs hapmapSNPs的指标统计 tagSNPs的群体结构验证 tagSNPs的遗传多样性 t ...

随机推荐

  1. net.js

    var http = require('http'); var https = require('https'); var url=require("url"); var pass ...

  2. Excel中如何截取字符串中指定字符后的部分字符

    1.如何给某列属性为时间整体加一个时间值:      场景一:假如我有一个excel中的某一列如下图所示,如何将该列的时间(用B代替整列)整体加一分钟呢?方法很简单,在空白单元格填写时间格式图中A所示 ...

  3. 记录Linux CentOS 7系统完整部署Docker容器环境教程

    笔者之前有在"详细介绍Ubuntu 16.04系统环境安装Docker CE容器的过程"文章中有介绍到利用Ubuntu系统安装Docker容器环境的过程.如果我们有使用CentOS ...

  4. cocoapod终端配置

  5. Python随笔--Pickle

    遇到问题:不能文件名为pickle.py

  6. MySql5.5安装详细说明

      双击MySql5.5安装文件开始: next 打勾,Next 选Custom,原因如下面所述: Typical(典型安装) Installs the most common program fea ...

  7. erlang证书加密

    -module(...). -include("ewp.hrl").-include("backend.hrl").-include_lib("pub ...

  8. 命令提示符出现-bash-4.1$如何解决

    故障原因:用户模式/etc/skel内容丢失 解决方法 复制etc/skel里内容到当前用户 [root@web01 ~]# su - test -bash-4.1$ cp /etc/skel/.ba ...

  9. servlet异步处理机制

    Servlet 3.0 之前,一个普通 Servlet 的主要工作流程大致如下:首先,Servlet 接收到请求之后,可能需要对请求携带的数据进行一些预处理:接着,调用业务接口的某些方法,以完成业务处 ...

  10. 模块化 Sea.js（CMD）规范 RequireJS(AMD)规范 的用法

    插入第三方库AMD CMD都 一样  如:JQ(再JQ源码里修改) 使用seajs的步骤 1.HTML里引入seajs <script src="./lib/sea.js"& ...