2017-06-14 730 views
0

我已经下载了1000G数据集的vcf格式。使用Plink 2.0我已经将它们转换成二进制格式。 现在我需要合并1-22条染色体。 我使用这个脚本:合并染色体Plink

${BIN}plink2 \ 
--bfile /mnt/jw01-aruk-home01/projects/jia_mtx_gwas_2016/common_files/data/clean/thousand_genomes/from_1000G_web/chr1_1000Gv3 \ 
--make-bed \ 
--merge-list /mnt/jw01-aruk-home01/projects/jia_mtx_gwas_2016/common_files/data/clean/thousand_genomes/from_1000G_web/chromosomes_1000Gv3.txt \ 
--out /mnt/jw01-aruk-home01/projects/jia_mtx_gwas_2016/common_files/data/clean/thousand_genomes/from_1000G_web/all_chrs_1000G_v3 \ 
--noweb 

但是,我得到这个错误

Error: --merge-list only accepts 1 parameter.

的chromosomes_1000Gv3.txt具有这种格式与染色体2-22文件:

chr2_1000Gv3.bed chr2_1000Gv3.bim chr2_1000Gv3.fam

chr3_1000Gv3.bed chr3_1000Gv3.bim chr3_1000Gv3.fam

....

任何建议可能是什么问题?

谢谢

回答

0

--merge-list不能与-bfile结合使用。您只能在一个plink命令中使用--bfile/- bmerge或--merge-list。