830 - MarkovModel::Condition()
From Genome Analysis Wiki
Jump to navigationJump to searchvoid MarkovModel::Condition(float * vector, char ** haplotypes, int position,
char observed, double e, double freq)
{
if (observed == 0) return;
double pmatch = (1. - e) + e * freq;
double prandom = e * freq;
for (int i = 0; i < states; i++)
if (haplotypes[i][position] == observed)
vector[i] *= pmatch;
else
vector[i] *= prandom;
}