Optimal step length EM algorithm (OSLEM) for the estimation of haplotype frequency and its application in lipoprotein lipase genotyping.
BACKGROUND
Haplotype based linkage disequilibrium (LD) mapping has become a powerful and cost-effective method for performing genetic association studies, particularly in the search for genetic markers in linkage disequilibrium with complex disease loci. Various methods (e.g. Monte-Carlo (Gibbs sampling); EM (expectation maximization); and Clark's method) have been used to estimate haplotype frequencies from routine genotyping data.


RESULTS
These algorithms can be very slow for large number of SNPs. In order to speed them up, we have developed a new algorithm using numerical analysis technology, a so-called optimal step length EM (OSLEM) that accelerates the calculation. By optimizing approximately the step length of the EM algorithm, OSLEM can run at about twice the speed of EM. This algorithm has been used for lipoprotein lipase (LPL) genotyping analysis.


CONCLUSIONS
This new optimal step length EM (OSLEM) algorithm can accelerate the calculation for haplotype frequency estimation for genotyping data without pedigree information. An OSLEM on-line server is available, as well as a free downloadable version.
Zhang
Peisen
P
Columbia Genome Center, Columbia University, New York NY 10032, USA. pz6@columbia.edu
Sheng
Huitao
H
Morabia
Alfredo
A
Gilliam
T Conrad
TC
D000465:Algorithms; D001161:Arteriosclerosis / Q000201:enzymology / Q000235:genetics; D016022:Case-Control Studies; D019295:Computational Biology / Q000379:methods; D005787:Gene Frequency / Q000235:genetics*; D005838:Genotype; D006239:Haplotypes / Q000235:genetics*; D006801:Humans; D008071:Lipoprotein Lipase / Q000235:genetics*; D020641:Polymorphism, Single Nucleotide / Q000235:genetics
D008071:Lipoprotein Lipase