Alternative Methods for H1 Simulations in Genome Wide Association Studies

Details Alternative Methods for H1 Simulations in Genome Wide Association Studies Alternative Methods for H1 Simulations in Genome Wide Association Studies

2012-01-24

arxiv.org/abs/1201.5046v1

Statistics

Applications

Scientific paper

Assessing the statistical power to detect susceptibility variants plays a critical role in GWA studies both from the prospective and retrospective points of view. Power is empirically estimated by simulating phenotypes under a disease model H1. For this purpose, the "gold" standard consists in simulating genotypes given the phenotypes (e.g. Hapgen). We introduce here an alternative approach for simulating phenotypes under H1 that does not require generating new genotypes for each simulation. In order to simulate phenotypes with a fixed total number of cases and under a given disease model, we suggest three algorithms: i) a simple rejection algorithm; ii) a numerical Markov Chain Monte-Carlo (MCMC) approach; iii) and an exact and efficient backward sampling algorithm. In our study, we validated the three algorithms both on a toy-dataset and by comparing them with Hapgen on a more realistic dataset. As an application, we then conducted a simulation study on a 1000 Genomes Project dataset consisting of 629 individuals (314 cases) and 8,048 SNPs from Chromosome X. We arbitrarily defined an additive disease model with two susceptibility SNPs and an epistatic effect. The three algorithms are consistent, but backward sampling is dramatically faster than the other two. Our approach also gives consistent results with Hapgen. Using our application data, we showed that our limited design requires a biological a priori to limit the investigated region. We also proved that epistatic effects can play a significant role even when simple marker statistics (e.g. trend) are used. We finally showed that the overall performance of a GWA study strongly depends on the prevalence of the disease: the larger the prevalence, the better the power.

Affiliated with

Also associated with

No associations

Rating

Alternative Methods for H1 Simulations in Genome Wide Association Studies does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.

If you have personal experience with Alternative Methods for H1 Simulations in Genome Wide Association Studies, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Alternative Methods for H1 Simulations in Genome Wide Association Studies will most certainly appreciate the feedback.

Rate now

Comments { 0 }

Profile ID: LFWR-SCP-O-107426