Lines Matching refs:pdg
44 static double est_freq(int n, const double *pdg) in est_freq() argument
50 const double *p = pdg + i * 3; in est_freq()
67 const double *pdg; member
79 const double *pdg = a->pdg + i * 3; in prob1() local
80 p *= pdg[0] * f3[0] + pdg[1] * f3[1] + pdg[2] * f3[2]; in prob1()
94 const double *pdg = _pdg + i * 3; in freq_iter() local
95 f0 += (pdg[1] * f3[1] + 2. * pdg[2] * f3[2]) in freq_iter()
96 / (pdg[0] * f3[0] + pdg[1] * f3[1] + pdg[2] * f3[2]); in freq_iter()
109 static double freqml(double f0, int beg, int end, const double *pdg) in freqml() argument
114 if (freq_iter(&f, pdg, beg, end) < EPS) break; in freqml()
117 a.beg = beg; a.end = end; a.pdg = pdg; in freqml()
132 const double *pdg = _pdg + i * 3; in g3_iter() local
133 tmp[0] = pdg[0] * g[0]; tmp[1] = pdg[1] * g[1]; tmp[2] = pdg[2] * g[2]; in g3_iter()
144 static double lk_ratio_test(int n, int n1, const double *pdg, double f3[3][3]) in lk_ratio_test() argument
149 const double *p = pdg + i * 3; in lk_ratio_test()
154 const double *p = pdg + i * 3; in lk_ratio_test()
169 double *pdg; in bcf_em1() local
177 pdg = call->pdg; in bcf_em1()
178 if (pdg == 0) return -1; in bcf_em1()
181 if ((x[0] = est_freq(n, pdg)) < 0.) return -1; // no data in bcf_em1()
182 x[0] = freqml(x[0], 0, n, pdg); in bcf_em1()
190 if (g3_iter(g, pdg, 0, n) < EPS) break; in bcf_em1()
193 double *p = pdg + i * 3; in bcf_em1()
199 x[5] = freqml(x[0], 0, n1, pdg); in bcf_em1()
200 x[6] = freqml(x[0], n1, n, pdg); in bcf_em1()
207 tmp = log(lk_ratio_test(n, n1, pdg, f3)); in bcf_em1()
215 if (g3_iter(g[1], pdg, 0, n1) < EPS) break; in bcf_em1()
217 if (g3_iter(g[2], pdg, n1, n) < EPS) break; in bcf_em1()
218 tmp = log(lk_ratio_test(n, n1, pdg, g)); in bcf_em1()
234 static int pair_freq_iter(int n, double *pdg[2], double f[4])
242 p[0] = pdg[0] + i * 3; p[1] = pdg[1] + i * 3;