This is a test version of Biostars. For the public version, visit https://www.biostars.org.
PLINK2 recode ped flag issue

Hi,

I am trying to covert .ped and .map files to a VCF file using plink2, as in the examples below:

plink2 --ped test.ped  --map .test.map --recode vcf --out testVCF
plink2 --file test --recode vcf --out testVCF

However, I am thrown the errors:

Error: Unrecognized flag ('--ped').
Error: Unrecognized flag ('--file').

Can anyone explain why this is?

I have been able to do the conversion using plink, however the REF alleles are incorrectly assigned (as is a well documented issue), and I do not have access to a reference file for the ~400K markers to define the reference alleles. I believe that plink2 is meant to have addressed this issue, and should correctly assign REF alleles.

Thanks in advance for any insight.

plink2 plink ped vcf
• 5,086 views • link
updated by tarek.mohamed tarek.mohamed • written by emyli • emyli

1 answer

Hi emily,

Try this plink 1.9-beta-4.6

plink --ped test.ped --map test.map --recode-vcf --out test_vcf
• 0 views • link
written by tarek.mohamed tarek.mohamed

Thanks for the suggestion, I believe this is the version I have been using (binary download ofStable (beta 4.6, 15 Aug) for OS X (64-bit) from https://www.cog-genomics.org/plink2). I tried again using your modified command and am getting the same result, I am using checkVCF.py to assess the validity of the output vcf, and am told most of the sites have inconsistent reference alleles.

• 0 views • link
written by emyli • emyli

What do you mean by REF alleles are incorrectly assigned. Can you attach an example.

T

• 0 views • link
written by tarek.mohamed tarek.mohamed

Log in to answer this question.

More posts like this