001、ROH 检测
[s20213040583@admin1 test]$ ls ## 测试文件 outcome.map outcome.ped ## plink软件检测 [s20213040583@admin1 test]$ plink --file outcome --homozyg --homozyg-density 50 --homozyg-gap 1000 --homozyg-kb 500 --homozyg-snp 100 --homozyg-window-het 1 --homozyg-window-snp 50 --homozyg-window-threshold 0.05 --homozyg-window-missing 5 --out roh &> /dev/null [s20213040583@admin1 test]$ ls outcome.map outcome.ped roh.hom roh.hom.indiv roh.hom.summary roh.log roh.nosex [s20213040583@admin1 test]$ head -n 3 roh.hom ## 查看前几行 FID IID PHE CHR SNP1 SNP2 POS1 POS2 KB NSNP DENSITY PHOM PHET SUN 3518 -9.000 2 s32258.1 s29579.1 244061250 248787856 4726.607 106 44.591 0.991 0.009 SUN 3589 -9.000 2 OAR2_13707962.1 OAR2_19604517.1 14351231 19320770 4969.540 104 47.784 0.981 0.019