计算群体pi

作者: 花生学生信 | 来源:发表于2021-07-12 07:30 被阅读0次

1、提取3k材料的名称

'''

bcftools query -l 3k.snp.vcf > 3k.list

'''

#2、挑选要计算pi的subset

#3、使用vcftools 计算pi #以30kb的窗口进行计算

vcftools --vcf 3k.snp.vcf --keep group1.list --window-pi 30000 --out group1.int

导入R

#读入数据

p1 <- read.table("group1.int.windowed.pi", header = T)

#按照群体命名

p1$pop <- "group1"

###合并三个文件

###df <- rbind(p1,p2,p3)

#画图

library(ggplot2)

options(scipen = 200)

ggplot(p1, aes(x=pop, y = PI, fill = pop)) +

geom_boxplot() +theme_classic()

如果是展示多个，可以使用grid或者facet_wrap进行展示

计算群体pi
1、提取3k材料的名称 ''' bcftools query -l 3k.snp.vcf > 3k.list ''...
【群体遗传学】 π （pi）的计算
杂合度 heterozygosity 某个位点的第个等位基因的样本频率为，那么该位点所有等位基因的频率和应该是1。...
vcftool 计算pi、theta w、TajimaD
（1）划窗口计算pi、TajimaD (2)awk 计算pi、theta w以及TajimaD 注：若有侵权，联系删除
pi分析
单个群体的基因组选择信号分析(pi分析) Pi值 Tajima D; Pi主要用来衡量每个site的nucleot...
canvas杂笔记-圆形
js中，Math.PI 计算角度弧度，Math.PI = 3.14 = 180° 绘制圆形： 2*Math.PI ...
Spark --- 启动、运行、关闭过程
计算PI值流程分析
【PI】函数使用技巧
Excel的PI函数主要是用来计算计算圆周率的近似值。本文介绍 Microsoft Excel中PI函数的公式语法...
python 简单脚本练习3
计算蛋白质 MW, pI, EC(消光系数)
核苷酸pi值计算方法
哈喽，今天小编给大家分享的是核苷酸pi值计算方法，主要是介绍序列比对（MEGA和MZAFFT版的两种方法）计算pi...
2018-10-25
作业1 计算半径为5的圆的面积和周长并打印出来。（pi为3.14）面积：pi*r*r;Perimeter(周长);...

网友评论

本文标题：计算群体pi

本文链接：https://www.haomeiwen.com/subject/ersjpltx.html

计算群体pi