Lines Matching defs:chisqp
36 /* Given nu and chisqp, compute probability that chisq > chisqp. This uses,
39 * Q(nu,chisqp) =
40 * erfc( (3/2)*sqrt(nu) * ( cbrt(chisqp/nu) - 1 + 2/(9*nu) ) ) / 2
46 chisq_prob (mpfr_ptr q, long nu, mpfr_ptr chisqp)
54 mpfr_div_si (q, chisqp, nu, rnd); /* chisqp/nu */
55 mpfr_cbrt (q, q, rnd); /* (chisqp/nu)^(1/3) */
56 mpfr_sub_ui (q, q, 1, rnd); /* (chisqp/nu)^(1/3) - 1 */
59 mpfr_add (q, q, t, rnd); /* (chisqp/nu)^(1/3) - 1 + 2/(9*nu) */