(defn quantile-r8 [p d]
  (let [d (sort d)
        N (count d)
        h (+ (* (+ N (/ 1.0 3.0)) p) (/ 1.0 3.0))
        hfloor (Math/floor h)]
  (cond (< p (/ (/ 2.0 3.0) (+ N (/ 1.0 3.0)))) (first d)
        (>= p (/ (- N (/ 1.0 3.0)) (+ N (/ 1.0 3.0)))) (last d)
        :else (+ (nth d (- hfloor 1)) (* (- h hfloor) (- (nth d hfloor) (nth d (- hfloor 1)))))))
  )