Changing/swapping per-sample calls

10 views
Skip to first unread message

Matthew Maher

unread,
Aug 7, 2025, 5:05:01 PMAug 7
to plink2-users
I have a dataset I was supplied it which it's been determined that somewhere upstream there was some allele confusion and thus for a small set of variants (these are all bi-allelic), I want to change the actual per-sample calls - e.g. if it's a A/C variant, I want any AAs to become CCs and vice versa. 

Is there a PLINK* option to do this? 

I do not want --flip as that's for strand problems; i.e. the current two alleles are correct
I do not want to change what's REF and whats ALT
i.e. I think my pvar is just fine, and this should be a change within the pgen.

Thanks for any info and thanks for PLINK*

Chris Chang

unread,
Aug 7, 2025, 6:01:08 PMAug 7
to Matthew Maher, plink2-users
One way to do this with PLINK 1.9 is to create small fileset(s) with just the genotypes you want to patch in, and then use --merge-mode 3 (putting the newly-created fileset(s) last in the merge order).

--
You received this message because you are subscribed to the Google Groups "plink2-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to plink2-users...@googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/plink2-users/d2635d3d-903c-43e6-a94b-898b27024294n%40googlegroups.com.
Reply all
Reply to author
Forward
0 new messages