PLINK2 recode ped flag issue
1
0
Entering edit mode
3.5 years ago
emyli ▴ 10

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 • 2.0k views
ADD COMMENT
0
Entering edit mode
3.5 years ago
tarek.mohamed ▴ 270

Hi emily,

Try this plink 1.9-beta-4.6

plink --ped test.ped --map test.map --recode-vcf --out test_vcf
ADD COMMENT
0
Entering edit mode

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.

ADD REPLY
0
Entering edit mode

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

T

ADD REPLY

Login before adding your answer.

Traffic: 2261 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6