830 - Haplotyper::ScoreLeftConditional()
From Genome Analysis Wiki
Jump to navigationJump to searchvoid Haplotyper::ScoreLeftConditional()
{
ResetMemoryPool();
GetMemoryBlock(0);
SetupPrior(leftMatrices[0]);
ConditionOnData(leftMatrices[0], 0, current_genotypes[0]);
double theta = 0.0;
for (int i = 1; i < markers; i++)
{
GetMemoryBlock(i);
Transpose(leftMatrices[i-1], leftMatrices[i], thetas[i-1]);
ConditionOnData(leftMatrices[i], i, current_genotypes[i]);
}
}