Ticket #836: inexacts.patch

statistics.scm

@@ -1495 +1495 @@
                                   (- yi-rank mean-y-rank)))
                              average-x-ranks average-y-ranks)))
              (rs (if (> (* Lxx Lyy) 0) (/ Lxy (sqrt (* Lxx Lyy))) 0)) ; TODO: think about rs = 1
+             (rs (exact->inexact rs))   ; Do not return ratnums
              (ts (if (= 1 rs) 1 (/ (* rs (sqrt (- n 2))) (sqrt (- 1 (square rs))))))
+             (ts (exact->inexact ts))   ; Do not return ratnums
              (p (t-significance ts (- n 2) #:tails ':both)))
         (format #t "~%Spearman correlation coefficient ~f, p = ~f~%" rs p)
         (values rs p))))