FST分析

  • Fst分析表示群体的分化程度,值越大,差异越大,群体分化程度越高,受选择程度越高。
  • Fst值的取值范围是【0,1】,最大值为1表明两个群体完全分化,最小值为0表明群体间无分化,(0则认为两个种群间是随机交配的,基因型完全相似;为1则表示是完全隔离的,完全不相似

    在实际的研究中Fst值为0--0.05时说明群体间遗传分化很小,可以不做考虑;

    为0.05--0.15时,表明群体间存在中等程度的遗传分化;

    为0.15--0.25时群体间存在较大的遗传分化;

    为0.25以上的时候群体间就存在很大的遗传分化了。

*ls:列表清单

*unzip:解压文件

*Tab:补全名称

*pwd:打印目录

①在 ‘ls’ 中用 ‘unzip X.zip’ 解压文件

② 

vcftools --vcf Desert_high_low.vcf --fst-window-size 100000 --fst-window-step 25000 --weir-fst-pop Desert.txt  --weir-fst-pop high.txt --out Desert_high.win100k.step25k

##  其中--vcf 是输入所需要计算的群体的输入文件,注意是vcf格式的。  --weir-fst-pop 这个命令是输入第一个群体文件,注意是txt文件格式。即Desert.txt,此文件只包含一列,就是群体个体的ID。high.txt也是一样的,是第二个群体的个体的ID。 ###

mean FST:绝对值      weight FST:加权平均值(一般用这个)

③pwd出 --out的结果并下载

④下载文件选择weight扩展降序,选择前5%   ##  ?×5%=!(!选列的时候+1)

⑤复制的!+1在.txt文件命名X_X.!snp.txt

  • 3
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值