julie.clj

(defn probability-of-taking-this-transition
  [tid incoming-variations transitions payment-rule]                                                                                                                           
  (let [this-trans             
          (first               
            (filter #(= tid (:id %)) transitions)),
        this-variations (this-trans :variations),
        this-feasibility (this-trans payment-rule),
        change-points          
          (distinct            
            (for [feas (map payment-rule transitions),
                  :when #(or (nil? this-feasibility) (<= feas this-feasibility))]                                                                                              
              (or feas :infinity)))]
    (prn "CHANGE POINTS" change-points)
    (zipmap                    
      change-points            
      (for [change-point change-points] 
        (let [feasibles        
                (filter        
                  (fn [tran]   
                    (let [its-feasibility (tran payment-rule)]
                      (if (= :infinity change-point)    
                        (nil? its-feasibility)          
                        (or
                          (nil? its-feasibility)          
                          (>= its-feasibility change-point)))))                                                                                                                
                  transitions)]
          (prn "FEASIBLES" change-point feasibles)
          (if (some #{tid} (map :id feasibles))
            ; This line is the crux of the bug I think
            ; maybe we should multiply by the incoming variations??
            ; (/ this-variations (apply + (map :variations feasibles))))))))
            ; (/ 1 (apply + (map :variations feasibles))))))))
            (/ incoming-variations (apply + (map :variations feasibles)))                                                                                                      
            0))))))