- acnt[0] = acnt[1] = acnt[2] = acnt[3] = 0;
- for (i = 0; i < ma->n; ++i)
- for (j = 0; j < 4; ++j)
- acnt[j] += ma->bcnt[i * 4 + j];
- if (ma->rcnt[ma->ref] + ma->rcnt[ma->alt]) // never used...
- fr = (double)ma->rcnt[ref] / (ma->rcnt[ma->ref] + ma->rcnt[ma->alt]);
- if (acnt[ma->ref] + acnt[ma->alt])
- f0 = (double)acnt[ref] / (acnt[ma->ref] + acnt[ma->alt]);