mvn_scatter_plot.clj

(use '(incanter core stats charts))

;; generate a sample from a bi-variate normal distribution
(def mvn-samp (sample-mvn 1000 :mean [7 5] :sigma (matrix [[2 1.5] [1.5 3]])))
(def x (sel mvn-samp :cols 0))
(def y (sel mvn-samp :cols 1))

;; add regression line to scatter plot
(def lm (linear-model y x))

;; plot scatter-plot of points
(doto (scatter-plot x y)
      view
      (add-lines x (:fitted lm)))

One small change needed to run this on current Incanter: Replace sample-multivariate-normal with sample-mvn

Thanks.