Question: Efficient resampling in GENESIS or another Bioconductor GWAS package to obtain empirical p-values?
gravatar for naglemi
9 months ago by
naglemi0 wrote:

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!

ADD COMMENTlink written 9 months ago by naglemi0


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.


ADD REPLYlink written 9 months ago by mconomos50
Please log in to add an answer.


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