1,555
edits
Changes
From Genome Analysis Wiki
Jump to navigationJump to search830 - MarkovModel::CountErrors() (view source)
Revision as of 14:57, 25 September 2013
, 14:57, 25 September 2013Created page with "<source lang="cpp"> double MarkovModel::CountErrors(float * vector, char ** haplotypes, int position, char observed, double e, double freq) { if (observed == 0) re..."
<source lang="cpp">
double MarkovModel::CountErrors(float * vector, char ** haplotypes, int position, char observed, double e, double freq)
{
if (observed == 0)
return e;
double match = 0;
double mismatch = 0;
for (int i = 0; i < states; i++)
if (haplotypes[i][position] == observed)
match += vector[i];
else
mismatch += vector[i];
mismatch = (match + mismatch) * e * freq;
match *= 1.0 - e;
return mismatch / (mismatch + match);
}
</source>
double MarkovModel::CountErrors(float * vector, char ** haplotypes, int position, char observed, double e, double freq)
{
if (observed == 0)
return e;
double match = 0;
double mismatch = 0;
for (int i = 0; i < states; i++)
if (haplotypes[i][position] == observed)
match += vector[i];
else
mismatch += vector[i];
mismatch = (match + mismatch) * e * freq;
match *= 1.0 - e;
return mismatch / (mismatch + match);
}
</source>
