shapeit填充

最新推荐文章于 2024-09-26 17:12:37 发布

BioXplorer

最新推荐文章于 2024-09-26 17:12:37 发布

阅读量849

点赞数 14

文章标签：笔记

本文链接：https://blog.csdn.net/2302_79242191/article/details/140787743

版权

使用shapeit软件进行填充

一，安装

下载地址，官网里面写得很详细。

https://mathgen.stats.ox.ac.uk/genetics_software/shapeit/shapeit.html

二，步骤

官网里面每一个参数都很详细

1.拆分染色体
for chr in {1..24};
do
plink --vcf /t.vcf \ --chr $chr \ --export vcf \ --out t.chr$chr \ --chr-set 24;
done

2.https://mathgen.stats.ox.ac.uk/genetics_software/shapeit/shapeit.html#model
shapeit \--input-vcf t.chr1.vcf \--rho 0.001 \--states 200 \--thread 8 \--window 0.8 \
--output-max gwas.phased.haps gwas.phased.sample \

##没有拆分染色体的情况，注意染色体的顺序，不然会报错
shapeit \ --input-bed t.2.bed t.2.bim t.2.fam \--rho 0.001 \--states 200 \--thread 8 \--window 0.8
--output-max gwas.phased.haps gwas.phased.sample \

#我怕跑多了会报错，所以是五个染色体一组

for chr in {1..5};do
shapeit \--input-vcf t.2.chr$chr.vcf \--rho 0.001 \--states 200 \--thread 8 \--window 0.8 \--output-max gwas.phased$chr.haps gwas.phased$chr.sample ;
done

3.变成vcf格式
shapeit -convert --input-haps gwas.phased --output-vcf gwas.phased.vcf

for chr in {1..5};do
shapeit -convert --input-haps gwas.phased$chr --output-vcf gwas.phased$chr.vcf;
done

4.vcf合并成一个染色体
for chr in {1..24};do plink --vcf gwas.phased$chr.vcf --make-bed --chr-set 24 --out gwas.phased$chr;done
制作一个文件gwas.txt，vcf文件的名称
gwas.phased1
gwas.phased2
gwas.phased3
gwas.phased4