projects
/
pspp
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
docs
[pspp]
/
lib
/
tukey
/
ptukey.c
diff --git
a/lib/tukey/ptukey.c
b/lib/tukey/ptukey.c
index 2af46cf89378fa6840502d3e986c6e38211420b2..07923d98bd7b5a363d7969847a333950243d5a6e 100644
(file)
--- a/
lib/tukey/ptukey.c
+++ b/
lib/tukey/ptukey.c
@@
-71,7
+71,7
@@
#include <assert.h>
#include <math.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 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 */
@@
-272,7
+272,6
@@
wprob (double w, double rr, double cc)
double
ptukey (double q, double rr, double cc, double df, int lower_tail, int log_p)
{
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
/* function ptukey() [was qprob() ]:
q = value of studentized range
@@
-327,7
+326,7
@@
ptukey (double q, double rr, double cc, double df, int lower_tail, int log_p)
f(x) = .5 + erf(x / sqrt(2)) / 2 for x > 0
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.
where f(x) is standard normal c. d. f.