The R Project SVN R

Rev

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

Rev 59097 Rev 60268
Line 49... Line 49...
49 
 *  where rN[j] ~ Bin(n, prob[j]) ,  sum_j rN[j] == n,  sum_j prob[j] == 1,
 49 
 *  where rN[j] ~ Bin(n, prob[j]) ,  sum_j rN[j] == n,  sum_j prob[j] == 1,
50 
 */
 50 
 */
51 
{
 51 
{
52 
    int k;
 52 
    int k;
53 
    double pp;
 53 
    double pp;
54 
    long double p_tot = 0.;
 54 
    LDOUBLE p_tot = 0.;
55 
    /* This calculation is sensitive to exact values, so we try to
 55 
    /* This calculation is sensitive to exact values, so we try to
56 
       ensure that the calculations are as accurate as possible
 56 
       ensure that the calculations are as accurate as possible
57 
       so different platforms are more likely to give the same
 57 
       so different platforms are more likely to give the same
58 
       result. */
 58 
       result. */
59 
 
 59