Difference between revisions of "830 - Haplotyper::ScoreLeftConditional()"
From Genome Analysis Wiki
Jump to navigationJump to searchLine 10: | Line 10: | ||
double theta = 0.0; | double theta = 0.0; | ||
− | |||
for (int i = 1; i < markers; i++) | for (int i = 1; i < markers; i++) | ||
{ | { | ||
− | + | GetMemoryBlock(i); | |
− | |||
− | + | Transpose(leftMatrices[i-1], leftMatrices[i], theta[i-1]); | |
− | + | ConditionOnData(leftMatrices[i], i, current_genotypes[i]); | |
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
} | } | ||
} | } | ||
</source> | </source> |
Revision as of 13:27, 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], theta[i-1]);
ConditionOnData(leftMatrices[i], i, current_genotypes[i]);
}
}