20 |
#include <config.h> |
#include <config.h> |
21 |
#include <gsl/gsl_cdf.h> |
#include <gsl/gsl_cdf.h> |
22 |
#include <gsl/gsl_sf_gamma.h> |
#include <gsl/gsl_sf_gamma.h> |
23 |
|
#include <gsl/gsl_math.h> |
24 |
#include <gsl/gsl_errno.h> |
#include <gsl/gsl_errno.h> |
25 |
|
|
26 |
|
#include "beta_inc.c" |
27 |
|
|
28 |
/* |
/* |
29 |
* Lower tail. |
* Lower tail. |
30 |
*/ |
*/ |
38 |
{ |
{ |
39 |
double u = x / (r + x); |
double u = x / (r + x); |
40 |
|
|
41 |
P = gsl_sf_beta_inc (nu1 / 2.0, nu2 / 2.0, u); |
P = beta_inc_AXPY (1.0, 0.0, nu1 / 2.0, nu2 / 2.0, u); |
42 |
} |
} |
43 |
else |
else |
44 |
{ |
{ |
45 |
double u = r / (r + x); |
double u = r / (r + x); |
46 |
|
|
47 |
P = 1 - gsl_sf_beta_inc (nu2 / 2.0, nu1 / 2.0, u); |
P = beta_inc_AXPY (-1.0, 1.0, nu2 / 2.0, nu1 / 2.0, u); |
48 |
} |
} |
49 |
|
|
50 |
return P; |
return P; |