Efficient resampling in GENESIS or another Bioconductor GWAS package to obtain empirical p-values?
Entering edit mode
naglemi ▴ 10
Last seen 5.4 years ago

I wish to perform resampling by bootstrapping or a similar method so that I can obtain empirical p-values for my data. The purpose of this is to correct for non-normality while avoiding overcorrection that comes with transformation of data.

I'm currently using SKAT, which allows one to build null models that include resamples and then test the alternative models against these null models. This is much faster than permuting phenotype data and re-running the test. I'm looking for a way to do something similar in a GWAS tool for single-locus testing, such as GENESIS.

Ideally, I would like the amount of resampling to depend on how long the p-value is, since there's no need to resample 1 million times for a p-value of 0.5, or contrarily, 1000 resamples would be insufficient for a p-value of 1e-5, which would probably appear to be 0 without enough resampling.

I haven't been able to find a built-in ressampling function in GENESIS, but does this exist for GENESIS or a similar GWAS package that tests SNPs one at a time (unlike SKAT which tests for combined effect of SNPs)?

Thank you!

GENESIS GWAS resampling resample statistics • 1.0k views
Entering edit mode


Unfortunately GENESIS does not currently have any built in resampling functions, and I'm not aware of any GWAS package that does. Is it possible to use the SKAT code to test one variant at a time by making your aggregation units each variant individually? Sorry I can't be of more help than that right now.



Login before adding your answer.

Traffic: 916 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6