         [clojure.contrib.math :only (sqrt)])
   (:require [incanter.core :as i]
             [incanter.stats :as s]))
-
-;; NOTE Currently, the NDVI and reliability sample time-series comes from the NDVI
-;; trends presentation.  We will have to start a test data file at
-;; some point, but I don't want to fuck with the directory structure
-;; right now, while I am learning to navigate the project.
 
 ;; The first few functions are more general, and could be pulled into,
 ;; say, the matrix operation namespace.  For now, we will leave them

(defn whizbang
  "force whizbang into an output map where the tuple for the reference
  period is separate from the tuples for all other time periods -
  those time periods that are used for est(imation) ... `for-est`!."
  [ref-pd start-pd end-pd t-series & cofactors]
  {:pre [(vector? t-series)
         (every? #{true} (vec (map vector? cofactors)))]}

  [f window xs]
  (pmap f (partition window 1 xs)))

(defn make-monotonic
  "move through a collection `coll` picking up the min or max values, depending
  on the value of `comparator` which can be either `min` or `max`.  This is very
  similar, I think, to the `reduce` function, but with side effects that print
  into a vector."
  [comparator coll]
  (reduce (fn [acc val]
            (conj acc
