Extract variants with same rsid

54 views
Skip to first unread message

Caterina

unread,
Jul 31, 2024, 9:08:38 AM7/31/24
to plink2-users
I have a few variants with same rsid but multiallelic variants:

1 10690948 rs17035482 C G
1 10690948 rs17035482 C T

How can I use the --extract flag to specify which variant I'd like to keep?

Thanks for your help :)

Christopher Chang

unread,
Jul 31, 2024, 9:56:19 AM7/31/24
to plink2-users
You must use something like plink2's --set-all-var-ids flag to first assign these variants unique IDs; otherwise --extract cannot distinguish between these variants.  (--recover-var-ids can be used later when you want to restore rsIDs.)

Caterina

unread,
Jul 31, 2024, 11:06:19 AM7/31/24
to plink2-users
Could you please give me an example of the format for set-all-var-ids?
Something like 
CHR:POS:REF:ALT?

Caterina

unread,
Jul 31, 2024, 11:08:20 AM7/31/24
to plink2-users
Would that be @:#:\$r:\$a ?

Christopher Chang

unread,
Jul 31, 2024, 11:54:40 AM7/31/24
to plink2-users
That should work.
Reply all
Reply to author
Forward
0 new messages