From 5ec99c0ddaa605df282876613eff20b9c09656ac Mon Sep 17 00:00:00 2001 From: zsloan Date: Fri, 4 Dec 2015 22:48:27 +0000 Subject: Commented out code related to additive/dominance effect in gn1 marker regression code; gn1 code now runs without errors but still need to give output to template --- wqflask/htmlgen.py | 14 ------- .../marker_regression/marker_regression_gn1.py | 47 +++++++++++++--------- wqflask/wqflask/views.py | 17 +++++--- 3 files changed, 39 insertions(+), 39 deletions(-) delete mode 100755 wqflask/htmlgen.py (limited to 'wqflask') diff --git a/wqflask/htmlgen.py b/wqflask/htmlgen.py deleted file mode 100755 index 8406d925..00000000 --- a/wqflask/htmlgen.py +++ /dev/null @@ -1,14 +0,0 @@ -from __future__ import absolute_import, division, print_function - -class HTMLgen(object): - """A redefined HT until we manage to completely eliminate it""" - def __getattr__(self, name): - return "" - - def Item(self, *args, **kw): - print("This way of generating html is obsolete!") - return "foo" - - Href = Span = TD = Blockquote = Image = Item - -HTMLgen2 = HTMLgen() diff --git a/wqflask/wqflask/marker_regression/marker_regression_gn1.py b/wqflask/wqflask/marker_regression/marker_regression_gn1.py index 5deeb305..c838b3ec 100644 --- a/wqflask/wqflask/marker_regression/marker_regression_gn1.py +++ b/wqflask/wqflask/marker_regression/marker_regression_gn1.py @@ -626,6 +626,9 @@ class MarkerRegression(object): self.dict['body'] = TD_LR self.dict['title'] = "Mapping" + + print("AT END OF GN1 MAPPING") + def writeQTL2Text(self, filename): if self.multipleInterval: return "" @@ -756,7 +759,8 @@ class MarkerRegression(object): drawAreaHeight -= 60 #Image map - gifmap = HT.Map(name='WebQTLImageMap') + gifmap = HT.Map(name = "WebQTLImageMap") + #gifmap = None newoffset = (xLeftOffset, xRightOffset, yTopOffset, yBottomOffset) # Draw the alternating-color background first and get plotXScale @@ -1689,7 +1693,7 @@ class MarkerRegression(object): distScale = 10 else: distScale = 5 - for tickdists in range(distScale, ceil(distLen), distScale): + for tickdists in range(distScale, int(ceil(distLen)), distScale): canvas.drawLine(startPosX + tickdists*plotXScale, yZero, startPosX + tickdists*plotXScale, yZero + 7, color=pid.black, width=1*zoom) canvas.drawString(str(tickdists), startPosX+tickdists*plotXScale, yZero + 10*zoom, color=pid.black, font=MBLabelFont, angle=270) startPosX += (self.ChrLengthDistList[i]+self.GraphInterval)*plotXScale @@ -1809,7 +1813,12 @@ class MarkerRegression(object): lodm = 1.0 if self.lrsMax <= 0: #sliding scale - LRSMax = max(map(max, self.qtlresults)).lrs + if "lrs_value" in self.qtlresults[0]: + LRSMax = max([result['lrs_value'] for result in self.qtlresults]) + #LRSMax = max(map(max, self.qtlresults)).lrs_value + else: + LRSMax = max([result['lod_score'] for result in self.qtlresults]) + #LRSMax = max(map(max, self.qtlresults)).lod_score #genotype trait will give infinite LRS LRSMax = min(LRSMax, webqtlConfig.MAXLRS) if self.permChecked and not self.multipleInterval: @@ -1882,11 +1891,11 @@ class MarkerRegression(object): if self.multipleInterval: lrsEdgeWidth = 1 else: - additiveMax = max(map(lambda X : abs(X.additive), self.qtlresults[0])) - if INTERCROSS: - dominanceMax = max(map(lambda X : abs(X.dominance), self.qtlresults[0])) - else: - dominanceMax = -1 + #additiveMax = max(map(lambda X : abs(X.additive), self.qtlresults[0])) + #if INTERCROSS: + # dominanceMax = max(map(lambda X : abs(X.dominance), self.qtlresults[0])) + #else: + # dominanceMax = -1 lrsEdgeWidth = 2 if zoom == 2: @@ -1894,7 +1903,7 @@ class MarkerRegression(object): for i, qtlresult in enumerate(self.qtlresults): m = 0 startPosX = xLeftOffset - thisLRSColor = self.colorCollection[i] + thisLRSColor = self.colorCollection[0] for j, _chr in enumerate(self.genotype): LRSCoordXY = [] AdditiveCoordXY = [] @@ -1907,20 +1916,20 @@ class MarkerRegression(object): # updated by NL 06-18-2011: # fix the over limit LRS graph issue since genotype trait may give infinite LRS; # for any lrs is over than 460(LRS max in this system), it will be reset to 460 - if qtlresult[m].lrs> 460 or qtlresult[m].lrs=='inf': + if qtlresult['lrs_value'] > 460 or qtlresult['lrs_value']=='inf': Yc = yZero - webqtlConfig.MAXLRS*LRSHeightThresh/LRSMax else: - Yc = yZero - qtlresult[m].lrs*LRSHeightThresh/LRSMax + Yc = yZero - qtlresult['lrs_value']*LRSHeightThresh/LRSMax LRSCoordXY.append((Xc, Yc)) - if not self.multipleInterval and self.additiveChecked: - if additiveMax == 0.0: - additiveMax = 0.000001 - Yc = yZero - qtlresult[m].additive*AdditiveHeightThresh/additiveMax - AdditiveCoordXY.append((Xc, Yc)) - if not self.multipleInterval and INTERCROSS and self.additiveChecked: - Yc = yZero - qtlresult[m].dominance*DominanceHeightThresh/dominanceMax - DominanceCoordXY.append((Xc, Yc)) + #if not self.multipleInterval and self.additiveChecked: + # if additiveMax == 0.0: + # additiveMax = 0.000001 + # Yc = yZero - qtlresult[m].additive*AdditiveHeightThresh/additiveMax + # AdditiveCoordXY.append((Xc, Yc)) + #if not self.multipleInterval and INTERCROSS and self.additiveChecked: + # Yc = yZero - qtlresult[m].dominance*DominanceHeightThresh/dominanceMax + # DominanceCoordXY.append((Xc, Yc)) m += 1 canvas.drawPolygon(LRSCoordXY,edgeColor=thisLRSColor,closed=0, edgeWidth=lrsEdgeWidth, clipX=(xLeftOffset, xLeftOffset + plotWidth)) diff --git a/wqflask/wqflask/views.py b/wqflask/wqflask/views.py index 3b100fa6..23c74817 100755 --- a/wqflask/wqflask/views.py +++ b/wqflask/wqflask/views.py @@ -354,11 +354,12 @@ def marker_regression_page(): if key in wanted or key.startswith(('value:')): start_vars[key] = value - version = "v1" + version = "v3" key = "marker_regression:{}:".format(version) + json.dumps(start_vars, sort_keys=True) print("key is:", pf(key)) with Bench("Loading cache"): - result = Redis.get(key) + result = None # Just for testing + #result = Redis.get(key) #print("************************ Starting result *****************") #print("result is [{}]: {}".format(type(result), result)) @@ -377,14 +378,18 @@ def marker_regression_page(): indent=" ") result = template_vars.__dict__ - print("initial result:", result['qtl_results']) + #print("initial result:", result['qtl_results']) - for item in template_vars.__dict__.keys(): - print(" ---**--- {}: {}".format(type(template_vars.__dict__[item]), item)) + #for item in template_vars.__dict__.keys(): + # print(" ---**--- {}: {}".format(type(template_vars.__dict__[item]), item)) #causeerror - gn1_template_vars = marker_regression_gn1.MarkerRegression(result) + print("TESTING GN1!!!") + gn1_template_vars = marker_regression_gn1.MarkerRegression(result).__dict__ + print("gn1_template_vars:", gn1_template_vars) + causeerror + #qtl_length = len(result['js_data']['qtl_results']) #print("qtl_length:", qtl_length) -- cgit v1.2.3