This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to convert multiple vcf files to a ped file?

I had 421 vcf files, I want to convert the all vcf files to the ped file (format of plink),

and I want to merge these 421 ped files to a ped file, what is the best should I do?

Thanks a lot!

plink vcftools
• 8,370 views • link
updated by Medhat Medhat • written by line1438 • line1438

2 answers

using plink
in the directory that contains your vcf files run this command "assuming that plink in the path"

for i in *.vcf;  do f="${i%.*}"; plink --vcf $i --maf 0.05 --recode --out $f".ped";done

example:

plink --file fA --merge-list allfiles.txt --make-bed --out mynewdata
• 1 views • link
written by Medhat Medhat

Thanks for the syntax. But, please specify version of PLINK used for conversion and merging. I am working on exome VCF files from IonProton. Converted VCFs using PLINK1.9 to PED/MAP files. But when I try to merge the file-set (PED/MAP) using PLINK1.9 it gives error about multiallelic SNPs. I already removed multiallelic records from VCF files using bcftools

bcftools view -m2 -M2 -v snps input.vcf

Still it gives the same error. Any further help on this Many thanks

• 0 views • link
written by Mamoon • Mamoon

This isnt what worked for me, but I complied the snippets here for posterity.

this thread was very helpful, thank you all. Medhat's link has moved over the years, it is here now: https://zzz.bwh.harvard.edu/plink/dataman.shtml#mergelist

• 0 views • link
written by DNAvinci • DNAvinci

This is what worked for me:

gatk CombineGVCFs \
   -R Homo_sapiens_assembly38.fasta \
   --variant myfile1.g.vcf.gz
   --variant myfile2.g.vcf.gz \
   -O cohort.g.vcf.gz

gatk GenotypeGVCFs \
   -R Homo_sapiens_assembly38.fasta \
   -V cohort.g.vcf.gz \
   --annotations-to-exclude InbreedingCoeff \
   -O corhort.vcf.gz

Then following Kevin Blighe's post here: Produce PCA bi-plot for 1000 Genomes Phase III - Version 2

• 0 views • link
written by DNAvinci • DNAvinci

Log in to answer this question.

More posts like this