From 2bc4ef27ad896e3acb7049035fcf30f00f9e84e0 Mon Sep 17 00:00:00 2001
From: Lei Yan
Date: Wed, 5 Feb 2014 15:42:19 -0600
Subject: On branch master

---
 wqflask/maintenance/dataset/specials.py | 13 ++++++++++---
 1 file changed, 10 insertions(+), 3 deletions(-)

(limited to 'wqflask')

diff --git a/wqflask/maintenance/dataset/specials.py b/wqflask/maintenance/dataset/specials.py
index 62311a92..ba9598fe 100644
--- a/wqflask/maintenance/dataset/specials.py
+++ b/wqflask/maintenance/dataset/specials.py
@@ -1,3 +1,4 @@
+import utilities
 import datastructure
 import genotypes
 import probesets
@@ -6,7 +7,7 @@ import probesets
 For:    Rob, GeneNetwork
 Date:   2014-02-04
 Function:
-    For BXD group, fetch probesets with given locus.
+    For BXD group, fetch probesets with given locus (mapping info).
 
 locus="rs3663871"
 """
@@ -65,10 +66,11 @@ def bxd_correlations():
     t = genotypes.load_genos(genofile)
     genostrains = t[0]
     genos = t[1]
-    print "Get %d genos" % (len(genos))
+    print "From geno file, get %d strains" % (len(genostrains))
+    print "From geno file, get %d genos" % (len(genos))
     #
     probesetfreezes = datastructure.get_probesetfreezes(inbredsetid)
-    print "Get %d probesetfreezes" % (len(probesetfreezes))
+    print "From DB, get %d probesetfreezes" % (len(probesetfreezes))
     #
     for probesetfreeze in probesetfreezes:
         #
@@ -99,6 +101,11 @@ def bxd_correlations():
             probeset = probesets.get_probeset(probesetid)
             probesetname = probeset[1]
             probesetdata = probesets.get_probesetdata(probesetdataid)
+            print probesetdata
+            probesetdata = zip(*probesetdata)
+            probesetdata = utilities.to_dic([strain.lower() for strain in probesetdata[1]], probesetdata[2])
+            print probesetdata
+            return
             #
             for geno in genos:
                 genoname = geno['locus']
-- 
cgit v1.2.3