Efficient Simulation of P Values for Linkage Analysis
- 13 January 2004
- journal article
- research article
- Published by Wiley in Genetic Epidemiology
- Vol. 26 (2) , 88-96
- https://doi.org/10.1002/gepi.10296
Abstract
In many genetic linkage analyses, the P value is obtained through simulation since the underlying distribution of the test statistic is complex and unknown. However, this can be very computationally intensive. A “bootstrap/replicate pool” approach has been suggested that generates P values more efficiently in terms of computation by resampling sums from a small set of simulated replicates for each pedigree. The replicate pool idea has been successfully applied, but, to our knowledge, has never been theoretically studied. An entirely different method for increasing the computational efficiency of P value simulation is Besag and Clifford's sequential sampling method. We propose an algorithm which combines Besag and Clifford's method with the replicate pool method to efficiently estimate P values for linkage studies. We derive variance expressions for the P value estimates from the replicate pool method and from our proposed hybrid method, and use these to show that the hybrid estimator has a substantial advantage over the other methods in most situations. Genet Epidemiol 26: 88–96, 2004.Keywords
This publication has 10 references indexed in Scilit:
- On Estimating P Values by Monte Carlo MethodsAmerican Journal of Human Genetics, 2003
- A Note on the Calculation of Empirical P Values from Monte Carlo ProceduresAmerican Journal of Human Genetics, 2003
- Simulation-Based P Values: Response to North et al.American Journal of Human Genetics, 2003
- A Note on the Calculation of Empirical P Values from Monte Carlo ProceduresAmerican Journal of Human Genetics, 2002
- Genetic Association Mapping Based on Discordant Sib Pairs: The Discordant-Alleles TestAmerican Journal of Human Genetics, 1998
- TestsPublished by Cambridge University Press (CUP) ,1997
- A bootstrap approach to estimating power for linkage heterogeneityGenetic Epidemiology, 1993
- A multisample bootstrap approach to the estimation of maximized-over-models lod score distributionsCytogenetic and Genome Research, 1992
- Sequential Monte Carlo p-valuesBiometrika, 1991
- Pedigree analysis by computer simulationZoo Biology, 1986