next up previous
Next: MCMC in pedigree Up: Early applications of Previous: Gene dropping

Assessing the power of a genetic study

 

Another early application of simulation methods is in estimating the power of a linkage study. The idea is to simulate marker data on the pedigrees available for a genetic study under an hypothesized genetic model before genotyping pedigree members.

The method presented here was introduced by Ploughman and Boehnke [12]. Ott [10] presents a similar idea. The first step is to simulate genotypes conditional on the observed phenotypes. This is done by peeling the pedigree up to one of the ancestors, saving all the and terms computed along the way. The procedure gives the probability distribution of the genotype of the last individual conditional on the phenotype data on all the people. The genotype of the last individual is generated from that distribution. Then, the genotype of the next before last person in the peeling order is simulated given the genotype of the last individual and the phenotypes of the rest of the pedigree members. The operation is iterated on the other individuals in the pedigree. This is called reverse peeling.

Once the disease locus genotypes of all the individuals have been generated, the genotypes at a marker having a recombination fraction with the disease locus can be simulated by first sampling the founder genotypes according to the marker allele frequencies and then simulating the gene flow down the pedigree according to .



Simon Cawley
Wed Apr 22 19:50:08 PDT 1998