1,555
edits
Changes
From Genome Analysis Wiki
Jump to navigationJump to search830 - Haplotyper::ImputeGenotypes() (internal) (view source)
Revision as of 14:29, 2 October 2013
, 14:29, 2 October 2013no edit summary
<source lang="cpp">
void Haplotyper::ImputeGenotypes(float * matrix1, float * matrix2, int marker)
void Haplotyper::ImputeGenotypes(float * matrix1, float * matrix2, int marker)
{
{
posterior[haplotypes[i][marker] * 2][marker] += *matrix1 * *matrix2;
posterior[haplotypes[i][marker] * 2][marker] += *matrix1 * *matrix2;
}
}
double sum = posterior[0][marker] = posterior[1][marker] = posterior[2][marker];
posterior[0][marker] /= sum;
posterior[1][marker] /= sum;
posterior[2][marker] /= sum;
}
}
</source>
