The R Project SVN R

Rev

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

Rev 13434 Rev 20687
Line 45... Line 45...
45 
 
 45 
 
46 
    x0 = floor(fmax2(c - 7. * sqrt(c), 0.));
 46 
    x0 = floor(fmax2(c - 7. * sqrt(c), 0.));
47 
    a0 = a + x0;
 47 
    a0 = a + x0;
48 
    lbeta = lgammafn(a0) + lgammafn(b) - lgammafn(a0 + b);
 48 
    lbeta = lgammafn(a0) + lgammafn(b) - lgammafn(a0 + b);
49 
    temp = pbeta_raw(x, a0, b, /* lower = */TRUE);
 49 
    temp = pbeta_raw(x, a0, b, /* lower = */TRUE);
50 
    gx = exp(a0 * log(x) + b * log(1. - x) - lbeta - log(a0));
 50 
    gx = exp(a0 * log(x) + b * log1p(-x) - lbeta - log(a0));
51 
    if (a0 > a)
 51 
    if (a0 > a)
52 
	q = exp(-c + x0 * log(c) - lgammafn(x0 + 1.));
 52 
	q = exp(-c + x0 * log(c) - lgammafn(x0 + 1.));
53 
    else
 53 
    else
54 
	q = exp(-c);
 54 
	q = exp(-c);
55 
 
 55