| Pré. | Proc. |
bdtrc
Complemented binomial distribution.SYNOPSIS:
int k, n;
double p, y, bdtrc();
y = bdtrc(k, n, p);
DESCRIPTION:
Returns the sum of the terms k+1 through n of the Binomial
probability density:
n
-- ( n ) j n-j
> ( ) p (1-p)
-- ( j )
j=k+1
The terms are not summed directly; instead the incomplete beta integral is employed, according to the formula:
y = bdtrc(k, n, p) = incbet(k+1, n-k, p).
The arguments must be positive, with p ranging from 0 to 1.
ACCURACY:
Tested at random points (a,b,p).
a,b Relative error:
arithmetic domain # trials peak rms
For p between 0.001 and 1:
IEEE 0,100 100000 6.7e-15 8.2e-16
For p between 0 and .001:
IEEE 0,100 100000 1.5e-13 2.7e-15
ERROR MESSAGES:
message condition value returned
domain x<0, x>1, n<k 0.0