Searched refs:pSame (Results 1 – 1 of 1) sorted by relevance
1463 double *PDiffVector(double *pSame, rates_t *rates);5146 double *pSame = mymalloc(sizeof(double) * rates->nRateCategories); in PSameVector() local5149 pSame[iRate] = 0.25 + 0.75 * exp((-4.0/3.0) * fabs(length*rates->rates[iRate])); in PSameVector()5150 return(pSame); in PSameVector()5153 double *PDiffVector(double *pSame, rates_t *rates) { in PDiffVector() argument5157 pDiff[iRate] = (1.0 - pSame[iRate])/3.0; in PDiffVector()5190 double *pSame = PSameVector(length, rates); in PairLogLk() local5191 double *pDiff = PDiffVector(pSame, rates); in PairLogLk()5214 lkAB = pSame[iRate] * wA*wB + 0.25 * (1-wA*wB); in PairLogLk()5236 lkAB += fB[j] * (fA[j] * pSame[iRate] + (1-fA[j])* pDiff[iRate]); /* P(A|B) */ in PairLogLk()[all …]