| 95 | } |
| 96 | |
| 97 | inline double qlnorm_1(double p, double meanlog, int lower_tail, int log_p){ |
| 98 | #ifdef IEEE_754 |
| 99 | if (ISNAN(p) || ISNAN(meanlog)) |
| 100 | return p + meanlog + 1.0; |
| 101 | #endif |
| 102 | R_Q_P01_boundaries(p, 0, ML_POSINF); |
| 103 | |
| 104 | return ::exp(::Rf_qnorm5(p, meanlog, 1.0, lower_tail, log_p)); |
| 105 | } |
| 106 | |
| 107 | } // stats |
| 108 | } // Rcpp |