本篇内容介绍了“如何使用XP CLR检测基因组中的选择信号”的有关知识，在实际案例的操作过程中，不少人都会遇到这样的困境，接下来就让小编带领大家学习一下如何处理这些情况吧！希望大家仔细阅读，能够学有所成！
如果不会编程，上述xpclr的输入文件还是不好准备的，这里教大家利用plink准换生成.geno 文件；由于xpclr一次只能运行一条染色体，所以这里举例第一条染色体运行：
计算
安全
数据库
网络和加速
企业服务
## 生成pop1亚群 geno文件
#样品名字变成两列方便plink转换生成.geno 文件：
awk '{print $1 "\t" $1}' pop1.txt > pop1.keep.txt
plink --vcf $workdir/00.filter/clean.sorted.vcf.gz \
--keep pop1.keep.txt --chr 1 --out Chr1.pop1 \
--recode 01 transpose -output-missing-genotype 9 \
--allow-extra-chr --set-missing-var-ids @:# --keep-allele-order
cut -d " " -f 5- Chr1.pop1.tped | awk '{print $0" "}' > Chr1.pop1.geno
## 生成pop2亚群 geno文件
awk '{print $1 "\t" $1}' $datadir/pop2.txt > pop2.keep.txt
plink --vcf $workdir/00.filter/clean.sorted.vcf.gz \
--keep pop2.keep.txt --chr 1 --out Chr1.pop2 \
--recode 01 transpose -output-missing-genotype 9 \
--allow-extra-chr --set-missing-var-ids @:# --keep-allele-order
cut -d " " -f 5- Chr1.pop2.tped | awk '{print $0" "}' > Chr1.pop2.geno
## 生成snp位置信息文件
zcat clean.sorted.vcf.gz|awk '$1=="1" {print " "$1":"$2 "\tChr01\t" $2/100000000 "\t" $2 "\t" $4 "\t" $5 }' > Chr01.snp
## 运行xpclr分析
XPCLR -xpclr Chr1.pop1.geno Chr1.pop2.geno Chr01.snp Chr01.out -w1 0.005 200 2000 Chr01 -p0 0.95
## 替换染色体名称
sed 's/^0/Chr01/' Chr01.out.xpclr.txt得到的结果文件中，每一列分别代表 chr# grid# #ofSNPs_in_window physical_pos genetic_pos XPCLR_score max_s，XPCLR_score 是算出的 XP-CLR 分数。
“如何使用XP CLR检测基因组中的选择信号”的内容就介绍到这里了，感谢大家的阅读。如果想了解更多行业相关的知识可以关注亿速云网站，小编将为大家输出更多高质量的实用文章！
亿速云「云服务器」，即开即用、新一代英特尔至强铂金CPU、三副本存储NVMe SSD云盘，价格低至29元/月。点击查看>>
免责声明：本站发布的内容（图片、视频和文字）以原创、转载和分享为主，文章观点不代表本网站立场，如果涉及侵权请联系站长邮箱：is@yisu.com进行举报，并提供相关证据，一经查实，将立刻删除涉嫌侵权内容。
