830 - Haplotyper::ImputeGenotypes() (internal)

void Haplotyper::ImputeGenotypes(float * matrix1, float * matrix2, int marker)

  {
  posterior[0][marker] = posterior[1][marker] = posterior[2][marker] = 0.0;

  for (int i = 0; i < states; i++, matrix1++, matrix2++)
     {
     if (haplotypes[i][marker])
        for (int j = 0; j < i; j++, matrix1++, matrix2++)
           posterior[haplotypes[j][marker] + 1][marker] += *matrix1 * *matrix2 * 0.5;
        else
     for (int j = 0; j < i; j++, matrix1++, matrix2++)
        posterior[haplotypes[j][marker]][marker] += *matrix1 * *matrix2 * 0.5;

     posterior[haplotypes[i][marker] * 2][marker] += *matrix1 * *matrix2;
     }
  }