#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 */
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.