static int test16_core(float anno[16], anno16_t *a)
{
extern double kt_fisher_exact(int n11, int n12, int n21, int n22, double *_left, double *_right, double *two);
double left, right;
int i;
a->p[0] = a->p[1] = a->p[2] = a->p[3] = 1.;
for (i=0; i<4; i++) a->d[i] = anno[i];
a->depth = anno[0] + anno[1] + anno[2] + anno[3];
a->is_tested = (anno[0] + anno[1] > 0 && anno[2] + anno[3] > 0);
if (a->depth == 0) return -1;
a->mq = (int)(sqrt((anno[9] + anno[11]) / a->depth) + .499);
kt_fisher_exact(anno[0], anno[1], anno[2], anno[3], &left, &right, &a->p[0]);
for (i = 1; i < 4; ++i)
a->p[i] = ttest(anno[0] + anno[1], anno[2] + anno[3], anno+4*i);
return 0;
}
int main(int argc, char **argv) {
for (;;) {
ttest();
}
}
