normal distribution of deltas

normal_deltas.clj

; data is a list of numbers
(let [deltas (loop [res '() rem data]
               (if (empty? rem)
                 res
                 (recur (concat res
                                (map #(Math/abs (- % (first rem))) (rest rem)))
                        (rest rem))))
      total (count deltas)
      mean (/ (reduce + 0 deltas) (count deltas))
      std-dev (-> (reduce #(+ % (Math/pow (- %2 mean) 2))
                         0
                         deltas)
                 (/ (dec total))
                 Math/sqrt)]
  (println "total:" total)
  (println "mean:" mean)
  (println "std-dev:" std-dev))