Journal
BIOINFORMATICS
Volume 23, Issue 5, Pages 545-554Publisher
OXFORD UNIV PRESS
DOI: 10.1093/bioinformatics/btl659
Keywords
-
Categories
Funding
- NCI NIH HHS [N01-CO-12400] Funding Source: Medline
- NHGRI NIH HHS [U01 HG003150, R01 HG002278] Funding Source: Medline
- NIGMS NIH HHS [T32-GM008802] Funding Source: Medline
Ask authors/readers for more resources
Motivation: Hidden Markov models (HMMs) and generalized HMMs been successfully applied to many problems, but the standard Viterbi algorithm for computing the most probable interpretation of an input sequence (known as decoding) requires memory proportional to the length of the sequence, which can be prohibitive. Existing approaches to reducing memory usage either sacrifice optimality or trade increased running time for reduced memory. Results: We developed two novel decoding algorithms, Treeterbi and Parallel Treeterbi, and implemented them in the TWINSCAN/N-SCAN gene-prediction system. The worst case asymptotic space and time are the same as for standard Viterbi, but in practice, Treeterbi optimally decodes arbitrarily long sequences with generalized HMMs in bounded memory without increasing running time. Parallel Treeterbi uses the same ideas to split optimal decoding across processors, dividing latency to completion by approximately the number of available processors with constant average overhead per processor. Using these algorithms, we were able to optimally decode all human chromosomes with N-SCAN, which increased its accuracy relative to heuristic solutions. We also implemented Treeterbi for Pairagon, our pair HMM based cDNA-to-genome aligner. Availability: The TWINSCAN/N-SCAN/PAIRAGON open source software package is available from http://genes.cse.wustl.edu. Contact: brent@cse.wustl.edu.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available