http://genome.sph.umich.edu/w/index.php?action=history&feed=atom&title=830_-_MarkovModel%3A%3ACountErrors%28%29 2026-09-28T08:27:11Z Revision history for this page on the wiki MediaWiki 1.43.1 http://genome.sph.umich.edu/w/index.php?title=830_-_MarkovModel::CountErrors()&diff=8737&oldid=prev 2013-09-25T18:57:26Z

<p>Created page with &quot;&lt;source lang=&quot;cpp&quot;&gt; double MarkovModel::CountErrors(float * vector, char ** haplotypes, int position, char observed, double e, double freq) { if (observed == 0) re...&quot;</p> <p><b>New page</b></p><div>&lt;source lang=&quot;cpp&quot;&gt;<br /> double MarkovModel::CountErrors(float * vector, char ** haplotypes, int position, char observed, double e, double freq)<br /> {<br /> if (observed == 0)<br /> return e;<br /> <br /> double match = 0;<br /> double mismatch = 0;<br /> <br /> for (int i = 0; i &lt; states; i++)<br /> if (haplotypes[i][position] == observed)<br /> match += vector[i];<br /> else<br /> mismatch += vector[i];<br /> <br /> mismatch = (match + mismatch) * e * freq;<br /> match *= 1.0 - e;<br /> <br /> return mismatch / (mismatch + match);<br /> }<br /> &lt;/source&gt;</div>

Goncalo