This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Pruning a dense SNP dataset without losing too much LD

Hi,

I have a set of ~1.7 million SNPs and I want to bring it down to ~300,000, but without losing too much LD between the remaining SNPs (because the methods I'm using after pruning rely heavily on IBD tracts).

What are the best parameters for the --indep-pairwise flag in PLINK to create a dataset with the desirable features? Any ideas?

Thanks!

plink
• 4,475 views • link
updated by Ram Ram • written by yorgos.athanasiadis • yorgos.athanasiadis

2 answers

If you’re explicitly trying to keep LD, and you’ve already performed basic quality control, —maf and —bp-space are probably the best tools to use.

• 1 views • link
written by chrchang523 chrchang523

You may need to experiment with increasing the --indep-pairwise flag a bit to get down to 300k. There is really only one parameter to play with, so just keep going until you reach your desired density.

• 0 views • link
updated by Ram Ram • written by Sean Davis Sean Davis

Log in to answer this question.

More posts like this