1,555
edits
Changes
From Genome Analysis Wiki
Jump to navigationJump to search830 - MarkovModel::Transpose() (view source)
Revision as of 14:39, 25 September 2013
, 14:39, 25 September 2013no edit summary
else
else
{
{
double sum = 0.0;
double sum = 0.0;
sum += from[i];
sum += from[i];
sum *= r * (1.0 - empiricalFlipRate) / states;
sum *= r / states;
double complement = 1. - r;
double complement = 1. - r;
{
{
sum *= 1e15;
sum *= 1e15;
complement *= 1e15;
complement *= 1e15;
}
}
for (int i = 0; i < states; i++)
for (int i = 0; i < states; i++)
to[i] = from[i] * complement + from[i^1] * flipRate + sum;
to[i] = from[i] * complement + sum;
}
}
}
}
</source>
</source>
