Difference between revisions of "830 - Haplotyper::ScoreLeftConditional()"
From Genome Analysis Wiki
Jump to navigationJump to searchLine 14: | Line 14: | ||
GetMemoryBlock(i); | GetMemoryBlock(i); | ||
− | Transpose(leftMatrices[i-1], leftMatrices[i], | + | Transpose(leftMatrices[i-1], leftMatrices[i], thetas[i-1]); |
ConditionOnData(leftMatrices[i], i, current_genotypes[i]); | ConditionOnData(leftMatrices[i], i, current_genotypes[i]); | ||
} | } | ||
} | } | ||
</source> | </source> |
Latest revision as of 13:30, 2 October 2013
void 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]);
}
}