X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=bcftools%2Fprob1.c;h=dc4c16ccc7c24edafd39ec0de193eb618435beb8;hb=cac5f70d5bc208e7699f5b8e72ebd120d41a28ea;hp=b27bb8e676cbb5e45e0eb0309133d0aa7a8460e1;hpb=cd65d352809ab157b72524d8542e1ee66e5a09c6;p=samtools.git

diff --git a/bcftools/prob1.c b/bcftools/prob1.c
index b27bb8e..dc4c16c 100644
--- a/bcftools/prob1.c
+++ b/bcftools/prob1.c
@@ -283,7 +283,7 @@ static void contrast(bcf_p1aux_t *ma, double pc[4]) // mc_cal_y() must be called
 {
 	int k, n1 = ma->n1, n2 = ma->n - ma->n1;
 	long double sum1, sum2;
-	pc[0] = pc[1] = pc[2] = pc[3] = 0.;
+	pc[0] = pc[1] = pc[2] = pc[3] = -1.;
 	if (n1 <= 0 || n2 <= 0) return;
 	for (k = 0, sum1 = 0.; k <= 2*n1; ++k) sum1 += ma->phi1[k] * ma->z1[k];
 	for (k = 0, sum2 = 0.; k <= 2*n2; ++k) sum2 += ma->phi2[k] * ma->z2[k];