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