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