22.11.2012 Views

Interdisciplinary Journal of Contemporary Research in ... - Webs

Interdisciplinary Journal of Contemporary Research in ... - Webs

Interdisciplinary Journal of Contemporary Research in ... - Webs

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

ijcrb.webs.com<br />

INTERDISCIPLINARY JOURNAL OF CONTEMPORARY RESEARCH IN BUSINESS<br />

In Figure. 4. Transition probabilities. The hidden states <strong>of</strong> the HMM represent different<br />

tree topologies, and state transitions correspond to recomb<strong>in</strong>ation events. The transition<br />

probability v is the probability that on mov<strong>in</strong>g from a site <strong>in</strong> the DNA sequence<br />

alignment to an adjacent site, no topology change occurs. If a topology change does<br />

occur, we assume that, a priori, all transitions are equally likely<br />

The parameters w, θ, and v need to be estimated.<br />

3. Heuristic Parameter Estimation: HMM-Heuristic<br />

McGuire, Wright, and Prentice (2000) estimated the branch lengths w for each tree<br />

topology separately with maximum likelihood. This approach is suboptimal. For a proper<br />

estimation <strong>of</strong> the branch lengths <strong>of</strong> a recomb<strong>in</strong>ant tree, one would have to restrict the<br />

parameter estimation to the recomb<strong>in</strong>ant region. The location <strong>of</strong> this region, however, is<br />

not known <strong>in</strong> advance. Estimat<strong>in</strong>g the branch lengths from the whole DNA sequence<br />

alignment leads to seriously distorted values, as demonstrated by Husmeier and Wright<br />

(2001) because the estimation <strong>in</strong>cludes regions <strong>of</strong> the alignment for which the tree<br />

topology is <strong>in</strong>correct. A heuristic way to address this problem, suggested by McGuire,<br />

Wright, and Prentice (2000), is to estimate the branch lengths from a subregion <strong>of</strong> the<br />

alignment. The length <strong>of</strong> this region should be matched to the length <strong>of</strong> the recomb<strong>in</strong>ant<br />

region, which, however, is not known <strong>in</strong> advance. Also, this approach does not <strong>of</strong>fer a<br />

way to estimate the recomb<strong>in</strong>ation parameter v.<br />

COPY RIGHT © 2011 Institute <strong>of</strong> <strong>Interdiscipl<strong>in</strong>ary</strong> Bus<strong>in</strong>ess <strong>Research</strong><br />

JANUARY 2011<br />

VOL 2, NO 9<br />

534

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!