projects
/
pspp
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
Fix problem with category borders.
[pspp]
/
lib
/
tukey
/
qtukey.c
diff --git
a/lib/tukey/qtukey.c
b/lib/tukey/qtukey.c
index ec146001420eca7d4876e6af12a813efe3f206da..54c174be2074be0378ffaed41c6903106d47e92f 100644
(file)
--- a/
lib/tukey/qtukey.c
+++ b/
lib/tukey/qtukey.c
@@
-143,8
+143,8
@@
static double qinv(double p, double c, double v)
ps = 0.5 - 0.5 * p;
yi = sqrt (log (1.0 / (ps * ps)));
ps = 0.5 - 0.5 * p;
yi = sqrt (log (1.0 / (ps * ps)));
- t = yi + ((((
yi * p4 + p3) * yi + p2) * yi + p1) * yi + p0)
- / ((((
yi * q4 + q3) * yi + q2) * yi + q1) * yi + q0);
+ t = yi + ((((yi * p4 + p3) * yi + p2) * yi + p1) * yi + p0)
+ / ((((yi * q4 + q3) * yi + q2) * yi + q1) * yi + q0);
if (v < vmax) t += (t * t * t + t) / v / 4.0;
q = c1 - c2 * t;
if (v < vmax) q += -c3 / v + c4 * t / v;
if (v < vmax) t += (t * t * t + t) / v / 4.0;
q = c1 - c2 * t;
if (v < vmax) q += -c3 / v + c4 * t / v;
@@
-192,14
+192,14
@@
double qtukey(double p, double rr, double cc, double df,
}
/* df must be > 1 ; there must be at least two values */
}
/* df must be > 1 ; there must be at least two values */
- /* ^^
+ /* ^^
JMD: The comment says 1 but the code says 2.
Which is correct?
*/
assert (df >= 2);
assert (rr >= 1);
assert (cc >= 2);
JMD: The comment says 1 but the code says 2.
Which is correct?
*/
assert (df >= 2);
assert (rr >= 1);
assert (cc >= 2);
-
+
R_Q_P01_boundaries (p, 0, ML_POSINF);
R_Q_P01_boundaries (p, 0, ML_POSINF);