beta.clj

(ns various-distributions.beta
  (:require [incanter.charts :as c]
            [incanter.core :as i]
            [incanter.distributions :as d]
            [incanter.stats :as s]))

(defn show-beta-distribution []
  (let [dist (d/beta-distribution 2 5)
        n 100000
        x (repeatedly n #(d/draw dist))
        pdf-fn #(d/pdf dist %)]
    (let [nbins 50
          series-label "gamma"
          [x-min x-max] (s/quantile x :probs [0.01 0.99])]
      (-> (c/histogram x :nbins nbins :density true :series-label series-label)
          (c/set-x-range x-min x-max)
          (c/add-function pdf-fn x-min x-max :series-label series-label)
          i/view))))