   (reduce + (map (partial log-likelihood beta-seq)
                  label-seq
                  feature-mat)))
-
-;; (defn probability-calc
-;;   "returns a vector of probabilities for each observation"
-;;   [beta-seq feature-mat]
-;;   (map (partial logistic-prob beta-seq)
-;;        feature-mat))
-
-(defn copy-row
-  [^DoubleMatrix row]
-  (.copy (to-double-matrix [[]]) row))
 
 (defn in-place-apply-fn
   [f ^DoubleMatrix row-mat]

           (dec iter)
           diff))))))

(defn estimated-probabilities
  "returns the set of probabilities, after applying the parameter
  values estimated over the training data; both `label-seq` and

    (.mmul (.muliRowVector (.transpose feature-mat) transformed-row)
           feature-mat)))

;; (info-matrix b feat)
;; #<DoubleMatrix [1.665334536937734E-14, 1.9984014443252805E-14, 2.3314683517128275E-14, 2.664535259100374E-14, 2.997602166487921E-14; 1.9984014443252805E-14, 4.539580775988847E-5, 9.079161549979293E-5, 1.3618742323969737E-4, 1.8158323097960185E-4; 2.331468351712827E-14, 9.079161549979293E-5, 1.8158323097627118E-4, 2.723748464527494E-4, 3.6316646192922767E-4; 2.664535259100374E-14, 1.3618742323969737E-4, 2.723748464527494E-4, 4.085622696658014E-4, 5.447496928788534E-4; 2.9976021664879214E-14, 1.8158323097960185E-4, 3.6316646192922767E-4, 5.447496928788534E-4, 7.263329238284793E-4]>

(defn beta-update
  "returns a vector of updates for the parameter vector; the
