This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Plink Error: Duplicate ID 'chr:pos:A1:A2' generated by --set-missing-var-ids

Hi,

I am trying to work with the 1000Genomes data starting from the vcf file, which has missing IDs.

So, I thought of adding the flag --set-missing-var-ids @:#:\$1:\$2,

but this gives me the error:

Error: Duplicate ID generated by --set-missing-var-ids.

So I guess first we want to find the duplicate IDs, but by chr, pos, A1, and A2.

Is there a way to do this? I know there is a --list-duplicate-vars but this does not give me IDs since I have many missing IDs.

Any help would be very appreciated!

plink duplicates missing id
• 9,837 views • link
updated by Ram Ram • written by jamespower • jamespower

Have you considered a simple awk or cut followed by sort | uniq -c? It might also be worth looking into various options that bcftools offers.

• 1 views • link
written by Ram Ram

You're much better off using plink 2.0 for this operation, since (i) it lets you specify REF/ALT alleles when constructing the IDs, avoiding the duplicate indel problem, and (ii) it forces you to specify how very long allele codes should be handled instead of silently truncating them, preventing some nasty surprises later on.

It's also reasonable to use awk/cut for this, but if you do, you'll need to be careful with very long indels.

• 1 views • link
written by chrchang523 chrchang523

Thanks! Using plink2, with --set-missing-var-ids @:#\$r:\$a, and --new-id-max-allele-len 23 missing works great!

• 1 views • link
written by jamespower • jamespower

0 answers

No answers yet.

Log in to answer this question.

More posts like this