The R Project SVN R

Rev

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

Rev 36820 Rev 37598
Line 66... Line 66...
66 
    NB -= xb;
 66 
    NB -= xb;
67 
 
 67 
 
68 
    if(!lower_tail || log_p) {
 68 
    if(!lower_tail || log_p) {
69 
	p = R_DT_qIv(p);
 69 
	p = R_DT_qIv(p);
70 
    }
 70 
    }
71 
    p *= 1 - 64*DBL_EPSILON;
 71 
    p *= 1 - 1000*DBL_EPSILON; /* was 64, but failed on FreeBSD sometimes */
72 
    sum = small_N ? term : exp(term);
 72 
    sum = small_N ? term : exp(term);
73 
 
 73 
 
74 
    while(sum < p && xr < xend) {
 74 
    while(sum < p && xr < xend) {
75 
	xr++;
 75 
	xr++;
76 
	NB++;
 76 
	NB++;