http://genome.sph.umich.edu/w/index.php?action=history&feed=atom&title=830_-_MarkovModel%3A%3AWalkRight%28%29 2026-09-28T03:42:56Z Revision history for this page on the wiki MediaWiki 1.43.1 http://genome.sph.umich.edu/w/index.php?title=830_-_MarkovModel::WalkRight()&diff=8732&oldid=prev 2013-09-25T18:45:10Z

<p>Created page with &quot;&lt;source lang=&quot;cpp&quot;&gt; void MarkovModel::WalkRight(char * major, char * observed, char ** haplotypes, float ** freqs) { float * swap; float * vector = new float [states]...&quot;</p> <p><b>New page</b></p><div>&lt;source lang=&quot;cpp&quot;&gt;<br /> void MarkovModel::WalkRight(char * major, char * observed, char ** haplotypes, float ** freqs)<br /> {<br /> float * swap;<br /> float * vector = new float [states];<br /> float * extra = new float [states];<br /> <br /> // Initialize likelihoods at first position<br /> for (int i = 0; i &lt; states; i++)<br /> vector[i] = 1.;<br /> <br /> // Scan along chromosome<br /> for (int i = markers - 1; i &gt; 0; i--)<br /> {<br /> for (int j = 0; j &lt; states; j++)<br /> extra[j] = vector[j] * matrix[i][j];<br /> <br /> Impute(major, observed, extra, haplotypes, freqs, i);<br /> <br /> if (observed[i])<br /> Condition(vector, haplotypes, i, observed[i], E[i], freqs[observed[i]][i]);<br /> Transpose(vector, extra, R[i - 1]);<br /> <br /> swap = vector; vector = extra; extra = swap;<br /> }<br /> <br /> if (observed[0])<br /> Condition(vector, haplotypes, 0, observed[0], E[0], freqs[observed[0]][0]);<br /> Impute(major, observed, vector, haplotypes, freqs, 0);<br /> <br /> delete [] vector;<br /> delete [] extra;<br /> }<br /> &lt;/source&gt;</div>

Goncalo