使用基因型数据计算近交系数

2022-07-30 00:01:52 阅读：171 来源： 互联网

标签：基因型 -- PC1 近交系数计算 test home outcome root

001、plink --het （纯合度近交系数）

root@PC1:/home/test# ls
outcome.map  outcome.ped
root@PC1:/home/test# plink --file outcome --het 1> /dev/null
root@PC1:/home/test# ls
outcome.map  outcome.ped  plink.het  plink.log
root@PC1:/home/test# head plink.het      ## 最后一列F为近交系数
    FID     IID       O(HOM)       E(HOM)        N(NM)            F
   DOR1    DOR1        28943    2.669e+04        42629       0.1412
   DOR2    DOR2        29737    2.669e+04        42629       0.1911
   DOR3    DOR3        28395    2.669e+04        42629       0.1068
   DOR4    DOR4        28430    2.669e+04        42629        0.109
   DOR5    DOR5        29235    2.669e+04        42629       0.1596
   DOR6    DOR6        28082    2.669e+04        42629       0.0872
   DOR7    DOR7        27779    2.669e+04        42629      0.06819
   DOR8    DOR8        29014    2.669e+04        42629       0.1457
   DOR9    DOR9        28135    2.669e+04        42629      0.09053

002、ROH 近交系数

root@PC1:/home/test# ls
outcome.map  outcome.ped                                                        ## plink检测ROH
root@PC1:/home/test# plink --file outcome --homozyg-snp 30 --homozyg-kb 1000 --homozyg-density 1000 --homozyg-gap 1000 --homozyg-window-snp 50 --homozyg-window-het 1 --homozyg-window-missing 1 --out roh 1> /dev/null
root@PC1:/home/test# ls
outcome.map  outcome.ped  roh.hom  roh.hom.indiv  roh.hom.summary  roh.log                                 ## 统计染色体的长度
root@PC1:/home/test# cut -f 1 outcome.map | sort -k 1n,1n| uniq | while read i; do grep -w "^$i" outcome.map | sort -k 4rn,4rn | head -n 1 >>chr_len.txt; done
root@PC1:/home/test# ls
chr_len.txt  outcome.map  outcome.ped  roh.hom  roh.hom.indiv  roh.hom.summary  roh.log
root@PC1:/home/test# geno_len=$(awk '{sum += $4} END {print sum / 1000}' chr_len.txt )  ## 计算基因组的长度                                           
root@PC1:/home/test# awk -v a=$geno_len '{print $2, $(NF - 1)/a}' roh.hom.indiv | head       ## 计算基因组近交系数
IID 0
DOR1 0.141861
DOR2 0.193722
DOR3 0.115116
DOR4 0.114285
DOR5 0.162598
DOR6 0.0992099
DOR7 0.0779776
DOR8 0.14785
DOR9 0.0968709

标签：基因型,--,PC1,近交系数,计算,test,home,outcome,root
来源： https://www.cnblogs.com/liujiaxin2018/p/16534047.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

ICode9

使用基因型数据计算近交系数