Modified Simpson O(n3) algorithm for the full sibship reconstruction problem
Journal Publication ResearchOnline@JCUAbstract
Motivation: The problem of reconstructing full sibling groups from DNA marker data remains a significant challenge for computational biology. A recently published heuristic algorithm based on Mendelian exclusion rules and the Simpson index was successfully applied to the full sibship reconstruction (FSR) problem. However, the so-called SIMPSON algorithm has an unknown complexity measure, questioning its applicability range. Results: We present a modified version of the SIMPSON (MS) algorithm that behaves as O(n3) and achieves the same or better accuracy when compared with the original algorithm. Performance of the MS algorithm was tested on a variety of simulated diploid population samples to verify its complexity measure and the significant improvement in efficiency (e.g. 100 times faster than SIMPSON in some cases). It has been shown that, in theory, the SIMPSON algorithm runs in non-polynomial time, significantly limiting its usefulness. It has been also verified via simulation experiments that SIMPSON could run in O(na), where a > 3. Availability: Computer code written in Java is available upon request from the first author.
Journal
Bioinformatics
Publication Name
N/A
Volume
21
ISBN/ISSN
1367-4811
Edition
N/A
Issue
20
Pages Count
6
Location
N/A
Publisher
Oxford University Press
Publisher Url
N/A
Publisher Location
Oxford, United Kingdom
Publish Date
N/A
Url
N/A
Date
N/A
EISSN
N/A
DOI
10.1093/bioinformatics/bti642