triangular.clj

(ns more-distributions.triangular
  (:require [incanter-contrib.distributions :as icd]
            [incanter.charts :as c]
            [incanter.core :as i]
            [incanter.distributions :as d]
            [incanter.stats :as s]))

(defn show-triangular-distribution []
  (let [dist (icd/triangular-distribution 0.0 0.3 1.0)
        n 100000
        x (repeatedly n #(d/draw dist))
        pdf-fn #(d/pdf dist %)]
    (let [nbins 50
          series-label "triangular"
          [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 "triangular-pdf")
          i/view))))