gistfile1.clj

(defn normal->beta [n,mean,sd,min,max]
  (let [range (- max min)
        m (/ (- mean min) range)
        s (/ sd range)
        a (/ (- (* m m) (* m m m) (* m (* s s))) (* s s))
        b (/ (+ m (* -2 m m) (* m m m) (- 0 (* s s) (* m s s))) (* s s))
        data (list (incanter-stats/sample-beta n :alpha a :beta b))]
    (map (fn [d] (+ min (* d range))) (flatten data))))