The R Project SVN R

Rev

Rev 87804 | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 87804 Rev 87821
Line 64... Line 64...
64 
#endif
 64 
#endif
65 
    if(p < 0. || q < 0.) ML_WARN_return_NAN;
 65 
    if(p < 0. || q < 0.) ML_WARN_return_NAN;
66 
    // allowing p==0 and q==0  <==> treat as one- or two-point mass
 66 
    // allowing p==0 and q==0  <==> treat as one- or two-point mass
67 
 
 67 
 
68 
    double qbet[2];// = { qbeta(), 1 - qbeta() }
 68 
    double qbet[2];// = { qbeta(), 1 - qbeta() }
69 
    qbeta_raw(alpha, p, q, lower_tail, log_p,
 69 
    qbeta_raw(alpha, p, q, lower_tail != 0, log_p != 0,
70 
	      // log_q_cut ,      n_N
 70 
	      // log_q_cut ,      n_N
71 
	      USE_LOG_X_CUTOFF, n_NEWTON_FREE, qbet);
 71 
	      USE_LOG_X_CUTOFF, n_NEWTON_FREE, qbet);
72 
    return qbet[0];
 72 
    return qbet[0];
73 
}
 73 
}
74 
 
 74