for (i = 0; i < n_smpl; ++i) {
const uint8_t *pi = PL[j] + i * PL_len[j];
double *p = pdg[j] + i * 3;
- p[0] = g_q2p[pi[b[j]->n_alleles]]; p[1] = g_q2p[pi[1]]; p[2] = g_q2p[pi[0]];
+ p[0] = g_q2p[pi[2]]; p[1] = g_q2p[pi[1]]; p[2] = g_q2p[pi[0]];
}
}
// iteration
memcpy(flast, f, 4 * sizeof(double));
freq_iter(n_smpl, pdg, f);
for (i = 0; i < 4; ++i) {
- double x = fabs(f[0] - flast[0]);
+ double x = fabs(f[i] - flast[i]);
if (x > eps) eps = x;
}
if (eps < LD_ITER_EPS) break;