* Vol.30, pp.1-15, 1988.
*/
-
+
#include <config.h>
#include <assert.h>
#include <math.h>
-#define R_D__0 (log_p ? ML_NEGINF : 0.) /* 0 */
+#define R_D__0 (log_p ? -INFINITY : 0.) /* 0 */
#define R_D__1 (log_p ? 0. : 1.) /* 1 */
#define R_DT_0 (lower_tail ? R_D__0 : R_D__1) /* 0 */
#define R_DT_1 (lower_tail ? R_D__1 : R_D__0) /* 1 */
#define ME_PRECISION 8
-static inline double
+static inline double
pnorm(double x, double mu, double sigma, int lower_tail, int log_p)
{
assert (lower_tail == 1);
assert (log_p == 0);
assert (sigma == 1.0);
-
+
return gsl_cdf_gaussian_P (x - mu, sigma);
}
double
ptukey (double q, double rr, double cc, double df, int lower_tail, int log_p)
{
- const double ML_NEGINF = -1.0 / 0.0;
/* function ptukey() [was qprob() ]:
q = value of studentized range
f(x) = .5 + erf(x / sqrt(2)) / 2 for x > 0
- f(x) = erfc( -x / sqrt(2)) / 2 for x < 0
+ f(x) = erfc(-x / sqrt(2)) / 2 for x < 0
where f(x) is standard normal c. d. f.
double ans, f2, f21, f2lf, ff4, otsum, qsqz, rotsum, t1, twa1, ulen, wprb;
int i, j, jj;
-#ifdef IEEE_754
- abort (! (ISNAN (q) || ISNAN (rr) || ISNAN (cc) || ISNAN (df)));
-#endif
+ assert (! (isnan (q) || isnan (rr) || isnan (cc) || isnan (df)));
if (q <= 0)
return R_DT_0;