about summary refs log tree commit diff
path: root/web/webqtl/correlationMatrix/tissueCorrelationMatrix.py
diff options
context:
space:
mode:
Diffstat (limited to 'web/webqtl/correlationMatrix/tissueCorrelationMatrix.py')
-rwxr-xr-xweb/webqtl/correlationMatrix/tissueCorrelationMatrix.py132
1 files changed, 0 insertions, 132 deletions
diff --git a/web/webqtl/correlationMatrix/tissueCorrelationMatrix.py b/web/webqtl/correlationMatrix/tissueCorrelationMatrix.py
deleted file mode 100755
index 23dc14eb..00000000
--- a/web/webqtl/correlationMatrix/tissueCorrelationMatrix.py
+++ /dev/null
@@ -1,132 +0,0 @@
-# Copyright (C) University of Tennessee Health Science Center, Memphis, TN.
-#
-# This program is free software: you can redistribute it and/or modify it
-# under the terms of the GNU Affero General Public License
-# as published by the Free Software Foundation, either version 3 of the
-# License, or (at your option) any later version.
-#
-# This program is distributed in the hope that it will be useful,
-# but WITHOUT ANY WARRANTY; without even the implied warranty of
-# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
-# See the GNU Affero General Public License for more details.
-#
-# This program is available from Source Forge: at GeneNetwork Project
-# (sourceforge.net/projects/genenetwork/).
-#
-# Contact Drs. Robert W. Williams and Xiaodong Zhou (2010)
-# at rwilliams@uthsc.edu and xzhou15@uthsc.edu
-#
-# This module is used by GeneNetwork project (www.genenetwork.org)
-#
-# Created by GeneNetwork Core Team 2010/11/10
-#
-# Last updated by Ning Liu, 2011/01/26
-
-
-#tissueCorrelationMatrix: funciton part for TissueCorrelationPage.py
-from htmlgen import HTMLgen2 as HT
-from correlation import correlationFunction
-from dbFunction import webqtlDatabaseFunction
-import sys
-
-#########################################
-#      Tissue Correlation Page
-#########################################
-
-class  tissueCorrelationMatrix:
-	def __init__(self,tissueProbeSetFreezeId=None):
-		
-		#initialize parameters
-		self.tProbeSetFreezeId = tissueProbeSetFreezeId
-		self.cursor = webqtlDatabaseFunction.getCursor()
-
-
-
-	#retreive dataSet info from database table TissueProbeSetFreeze to get all TissueProbeSetFreezeId(List), Name(List) and FullName(List)
-	def getTissueDataSet(self):	
-		tissProbeSetFreezeIds,Names,fullNames = webqtlDatabaseFunction.getTissueDataSet(cursor=self.cursor)
-		return tissProbeSetFreezeIds,Names,fullNames
-
-
-	#retrieve DatasetName, DatasetFullName based on TissueProbeSetFreezeId, return DatasetName(string), DatasetFullName(string)
-	def getFullnameofCurrentDataset(self):
-	
-		DatasetName, DatasetFullName =webqtlDatabaseFunction.getDatasetNamesByTissueProbeSetFreezeId(cursor=self.cursor, TissueProbeSetFreezeId=self.tProbeSetFreezeId)		
-		return DatasetName, DatasetFullName
-
-				
-	#retrieve how many tissue used in the specific dataset based on TissueProbeSetFreezeId, return TissueCount(int)
-	def getTissueCountofCurrentDataset(self):
-	
-		TissueCount =webqtlDatabaseFunction.getTissueCountByTissueProbeSetFreezeId(cursor=self.cursor,TissueProbeSetFreezeId=self.tProbeSetFreezeId)
-		return TissueCount
-
-
-		
-	#retrieve corrArray(array), pvArray(array) for display by calling  calculation function:calZeroOrderCorrForTiss
-	def getTissueCorrPvArray(self,geneNameLst=None,dataIdDict=None):	
-		#retrieve SymbolValuePairDict(Dict), dictionary of Symbol and Value Pair.key is symbol, value is one list of expression values of one probeSet
-		symbolValuepairDict =correlationFunction.getGeneSymbolTissueValueDict(cursor=self.cursor,symbolList=geneNameLst,dataIdDict=dataIdDict)
-		corrArray,pvArray = correlationFunction.getCorrPvArray(cursor=self.cursor,priGeneSymbolList=geneNameLst,symbolValuepairDict=symbolValuepairDict)
-		return corrArray,pvArray
-
-
-		
-	#retrieve symbolList,geneIdList,dataIdList,ChrList,MbList,descList,pTargetDescList (all are list type) to 
-	#get multi lists for short and long label functions, and for getSymbolValuePairDict and 
-	#getGeneSymbolTissueValueDict to build dict to get CorrPvArray
-	def getTissueProbeSetXRefInfo(self,GeneNameLst=[]):
-		symbolList,geneIdDict,dataIdDict,ChrDict,MbDict,descDict,pTargetDescDict =correlationFunction.getTissueProbeSetXRefInfo(cursor=self.cursor,GeneNameLst=GeneNameLst,TissueProbeSetFreezeId=self.tProbeSetFreezeId)
-		return symbolList,geneIdDict,dataIdDict,ChrDict,MbDict,descDict,pTargetDescDict		
-
-
-
-	#retrieve corrArray(array), pvArray(array) for gene symbol pair
-	def getCorrPvArrayForGeneSymbolPair(self,geneNameLst=None):
-		corrArray = None
-		pvArray = None
-
-		if len(geneNameLst) == 2:
-			#retrieve SymbolValuePairDict(Dict), dictionary of Symbol and Value Pair.key is symbol, value is one list of expression values of one probeSet
-			symbolList,geneIdDict,dataIdDict,ChrDict,MbDict,descDict,pTargetDescDict =correlationFunction.getTissueProbeSetXRefInfo(cursor=self.cursor,GeneNameLst=geneNameLst,TissueProbeSetFreezeId=self.tProbeSetFreezeId)	
-			symbolValuepairDict =correlationFunction.getGeneSymbolTissueValueDict(cursor=self.cursor,symbolList=geneNameLst,dataIdDict=dataIdDict)
-			corrArray,pvArray = correlationFunction.getCorrPvArray(cursor=self.cursor,priGeneSymbolList=geneNameLst,symbolValuepairDict=symbolValuepairDict)
-
-		return corrArray,pvArray
-
-
-	#retrieve symbolCorrDict(dict), symbolPvalueDict(dict) to get all tissues' correlation value and P value; key is symbol	
-	def calculateCorrOfAllTissueTrait(self, primaryTraitSymbol=None, method='0'):
-		symbolCorrDict, symbolPvalueDict = correlationFunction.calculateCorrOfAllTissueTrait(cursor=self.cursor, primaryTraitSymbol=primaryTraitSymbol, TissueProbeSetFreezeId=self.tProbeSetFreezeId,method=method)
-		
-		return symbolCorrDict, symbolPvalueDict
-
-	#Translate GeneId to gene symbol and keep the original order.
-	def getGeneSymbolLst(self, geneSymbols=None):
-		geneSymbolLst=[]
-		geneIdLst=[]
-		#split the input string at every occurrence of the delimiter '\r', and return the substrings in an array. 
-		tokens=geneSymbols.strip().split('\r')
-
-		#Ning: To keep the original order of input symbols and GeneIds
-		for i in tokens:
-			i=i.strip()
-			if (len(i) >0) and (i not in geneSymbolLst):
-				geneSymbolLst.append(i)
-				# if input includes geneId(s), then put it/them into geneIdLst
-				if i.isdigit():
-					geneIdLst.append(i)
-
-		#Ning: Replace GeneId with symbol if applicable
-		if len(geneIdLst)>0:
-			# if input includes geneId(s), replace geneId by geneSymbol; 
-			geneIdSymbolPair =webqtlDatabaseFunction.getGeneIdSymbolPairByGeneId(cursor=self.cursor, geneIdLst =geneIdLst)
-			for geneId in geneIdLst:
-				if geneIdSymbolPair[geneId]:
-					index = geneSymbolLst.index(geneId)
-					geneSymbolLst[index] =geneIdSymbolPair[geneId]
-
-		return geneSymbolLst
-			
-		
-