| 65 | return R_D_exp(x); |
| 66 | } |
| 67 | inline double qweibull_1(double p, double shape /*, double scale [=1.0] */, int lower_tail, int log_p){ |
| 68 | #ifdef IEEE_754 |
| 69 | if (ISNAN(p) || ISNAN(shape)) |
| 70 | return p + shape + 1.0; |
| 71 | #endif |
| 72 | if (shape <= 0) return R_NaN; |
| 73 | |
| 74 | R_Q_P01_boundaries(p, 0, ML_POSINF); |
| 75 | |
| 76 | return ::pow(- R_DT_Clog(p), 1./shape); |
| 77 | } |
| 78 | |
| 79 | } // stats |
| 80 | } // Rcpp |