The R Project SVN R

Rev

Rev 64655 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 64655 Rev 64668
Line 36... Line 36...
36 
	return p + NR + NB + n;
 36 
	return p + NR + NB + n;
37 
#endif
 37 
#endif
38 
    if(!R_FINITE(p) || !R_FINITE(NR) || !R_FINITE(NB) || !R_FINITE(n))
 38 
    if(!R_FINITE(p) || !R_FINITE(NR) || !R_FINITE(NB) || !R_FINITE(n))
39 
	ML_ERR_return_NAN;
 39 
	ML_ERR_return_NAN;
40 
 
 40 
 
41 
    NR = R_D_forceint(NR);
 41 
    NR = R_forceint(NR);
42 
    NB = R_D_forceint(NB);
 42 
    NB = R_forceint(NB);
43 
    N = NR + NB;
 43 
    N = NR + NB;
44 
    n = R_D_forceint(n);
 44 
    n = R_forceint(n);
45 
    if (NR < 0 || NB < 0 || n < 0 || n > N)
 45 
    if (NR < 0 || NB < 0 || n < 0 || n > N)
46 
	ML_ERR_return_NAN;
 46 
	ML_ERR_return_NAN;
47 
 
 47 
 
48 
    /* Goal:  Find  xr (= #{red balls in sample}) such that
 48 
    /* Goal:  Find  xr (= #{red balls in sample}) such that
49 
     *   phyper(xr,  NR,NB, n) >= p > phyper(xr - 1,  NR,NB, n)
 49 
     *   phyper(xr,  NR,NB, n) >= p > phyper(xr - 1,  NR,NB, n)