Question: Generating random sets of LD-independent SNPs from GWAS and calculating an empirical p-value
gravatar for AB
9 months ago by
AB0 wrote:


I am trying to calculate 'empirical p-value' by-

drawing 10000 random sets of 5000 SNPs each from a GWAS summary statistics, where SNPs in each set must be LD-independent (LD<0.8, 1000G Phase 3 EUR ref) and, then counting the number of SNPs in each set that overcome a certain association threshold, and are 'hits'.

Might someone have an idea on how this should be done and to plot the results?

I assume the empirical p-value can then be calculated as: Number of tests where SNPs overcome association threshold p-value/ total number of tests (10000).

Thank you all.

ADD COMMENTlink modified 9 months ago • written 9 months ago by AB0

Why not use real public GWAS data?

ADD REPLYlink written 9 months ago by zx87546.8k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2508 users visited in the last hour