From bd48ca02770735cff482b1612399f236e8ac377f Mon Sep 17 00:00:00 2001
From: Zachary Sloan
Date: Fri, 3 Jan 2014 21:47:58 +0000
Subject: Wrote some code for the histogram figure

Changed some code in dataset to try and track down the issue with Amelie's file
---
 wqflask/base/data_set.py | 13 +++++++------
 1 file changed, 7 insertions(+), 6 deletions(-)

(limited to 'wqflask/base')

diff --git a/wqflask/base/data_set.py b/wqflask/base/data_set.py
index fbe78d5d..d421f21e 100755
--- a/wqflask/base/data_set.py
+++ b/wqflask/base/data_set.py
@@ -168,13 +168,14 @@ class Markers(object):
         
         for marker, p_value in itertools.izip(self.markers, p_values):
             marker['p_value'] = p_value
-            if math.isnan(marker['p_value']):
+            if math.isnan(marker['p_value']) or marker['p_value'] <= 0:
                 print("p_value is:", marker['p_value'])
-            marker['lod_score'] = -math.log10(marker['p_value'])
-            #Using -log(p) for the LRS; need to ask Rob how he wants to get LRS from p-values
-            marker['lrs_value'] = -math.log10(marker['p_value']) * 4.61
-        
-        
+                marker['lod_score'] = 0
+                marker['lrs_value'] = 0
+            else:
+                marker['lod_score'] = -math.log10(marker['p_value'])
+                #Using -log(p) for the LRS; need to ask Rob how he wants to get LRS from p-values
+                marker['lrs_value'] = -math.log10(marker['p_value']) * 4.61
 
 
 class HumanMarkers(Markers):
-- 
cgit 1.4.1