Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Tutorials
Tools
Jobs
Forum
Planet
All »
Home
Log In
Question: PLINK2 recode ped flag issue
0
gravatar for emyli
17 months ago by
emyli • 0
emyli • 0 wrote:

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.
ped plink plink2 vcf • 961 views
ADD COMMENT • link •
modified 17 months ago by tarek.mohamed • 240 • written 17 months ago by emyli • 0
0
gravatar for tarek.mohamed
17 months ago by
tarek.mohamed • 240
tarek.mohamed • 240 wrote:

Hi emily,

Try this plink 1.9-beta-4.6

plink --ped test.ped --map test.map --recode-vcf --out test_vcf
ADD COMMENT • link written 17 months ago by tarek.mohamed • 240

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 • link written 17 months ago by emyli • 0

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

T

ADD REPLY • link modified 17 months ago • written 17 months ago by tarek.mohamed • 240
Please log in to add an answer.

Similar posts • Search »

Content
Help
Access
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1802 users visited in the last hour