--indep-pairwise best practice?

584 views
Skip to first unread message

Oveis Jamialahmadi

unread,
Mar 20, 2020, 8:07:53 AM3/20/20
to plink2-users
Hi Christopher!

I am aware that this question has been asked probably several times; but I'm interested to know more details if possible:
Is there any particular way to assess the best parameters for --indep-pairwise? I'm aware that it depends on downstream purposes and 'kb' can be used (e.g. 500 kb 1 r2) for more stringnet purposes, but still the results would be different. As an example, I ran PLINK2 using these 3 different scenarios:
  1. --indep-pairwise 500 kb 1 0.2
  2. --indep-pairwise 50 5 0.2
  3. --indep-pairwise 1000 80 0.2
Each of them generates a different set of markers, with some being unique (e.g. approach 1 cannot detect markers above 500 kb distance). My question is: is there any approach to optimize these parameters? Let's say we want to detect all possible independent markers with a specific r2. Is this even possible or if not, is there any greedy approach to find them?


Thanks in advance,
Best/Oveis  

Christopher Chang

unread,
Mar 20, 2020, 1:38:24 PM3/20/20
to plink2-users

You can scan an entire chromosome at a time by setting the window size larger than any chromosome length (e.g. "--indep-pairwise 250000kb ..."), but this is rarely worth the computational cost.
Reply all
Reply to author
Forward
0 new messages