The R Project SVN R

Rev

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

Rev 19799 Rev 20687
Line 35... Line 35...
35
    if(p < 0 || p > 1) ML_ERR_return_NAN;
35
    if(p < 0 || p > 1) ML_ERR_return_NAN;
36
 
36
 
37
    if (x < 0. || p == 0.) return R_DT_0;
37
    if (x < 0. || p == 0.) return R_DT_0;
38
    if (!R_FINITE(x)) return R_DT_1;
38
    if (!R_FINITE(x)) return R_DT_1;
39
    if(log_p && !lower_tail)
39
    if(log_p && !lower_tail)
40
	return log(1 - p) * (x + 1);
40
	return log1p(-p) * (x + 1);
41
    return R_DT_Cval(pow(1 - p, x + 1));
41
    return R_DT_Cval(pow(1 - p, x + 1));
42
}
42
}