def testNonMatchingPlots2HLAAgainstPolymerase(self): """ Examine the correlation between light matter scores and BLAST bit scores. """ scoreTypeX = 'Light matter score' scoreTypeY = 'Bit score' for parameterSet in testArgs.parameterSets: affinity = _AFFINITY[parameterSet] dirName = makeOutputDir( DATASET, parameterSet, '%s-%s' % (FILESYSTEM_NAME[scoreTypeX], FILESYSTEM_NAME[scoreTypeY])) for subject in SUBJECTS: lmScores = [] bitScores = [] for query in QUERIES: if query.id != subject.id: bitScores.append(BIT_SCORES[query.id][subject.id]) lmScore = getScore(affinity, query.id, subject.id) lmScores.append(lmScore) plot(lmScores, bitScores, subject.id, scoreTypeX, scoreTypeY, dirName)
def testPlotsHA(self): """ Examine the correlation between our scores and blast bit scores. """ scoreTypeX = 'Light matter score' scoreTypeY = 'Bit score' for parameterSet in testArgs.parameterSets: affinity = _AFFINITY[parameterSet] dirName = makeOutputDir( DATASET, parameterSet, '%s-%s' % (FILESYSTEM_NAME[scoreTypeX], FILESYSTEM_NAME[scoreTypeY])) allBitScores = [] allLmScores = [] for queryId in BIT_SCORES: bitScores = [] lmScores = [] for subjectId in BIT_SCORES[queryId]: if queryId != subjectId: lmScore = getScore(affinity, queryId, subjectId) lmScores.append(lmScore) bitScores.append(BIT_SCORES[queryId][subjectId]) allBitScores.extend(bitScores) allLmScores.extend(lmScores) plot(lmScores, bitScores, queryId, scoreTypeX, scoreTypeY, dirName) plot(allLmScores, allBitScores, 'all', scoreTypeX, scoreTypeY, dirName)
def testPlotsHA(self): """ Examine the correlation between BLAST bit scores and Z scores. """ scoreTypeX = 'Bit score' scoreTypeY = 'Z score' for parameterSet in testArgs.parameterSets: dirName = makeOutputDir( DATASET, parameterSet, '%s-%s' % (FILESYSTEM_NAME[scoreTypeX], FILESYSTEM_NAME[scoreTypeY])) allZScores = [] allBitScores = [] for queryId in BIT_SCORES: zScores = [] bitScores = [] for subjectId in BIT_SCORES[queryId]: if queryId != subjectId: bitScores.append(BIT_SCORES[queryId][subjectId]) zScores.append(Z_SCORES[queryId][subjectId]) allZScores.extend(zScores) allBitScores.extend(bitScores) plot(bitScores, zScores, queryId, scoreTypeX, scoreTypeY, dirName) plot(allBitScores, allZScores, 'all', scoreTypeX, scoreTypeY, dirName)
def testPlots2HLA(self): """ Examine the correlation between BLAST bit scores and Z scores. """ scoreTypeX = 'Bit score' scoreTypeY = 'Z score' for parameterSet in testArgs.parameterSets: dirName = makeOutputDir( DATASET, parameterSet, '%s-%s' % (FILESYSTEM_NAME[scoreTypeX], FILESYSTEM_NAME[scoreTypeY])) for subject in SUBJECTS: bitScores = [] zScores = [] for query in QUERIES: if query.id != subject.id: zScores.append(Z_SCORES[query.id][subject.id]) bitScores.append(BIT_SCORES[query.id][subject.id]) plot(bitScores, zScores, subject.id, scoreTypeX, scoreTypeY, dirName)
def testPlots4PH0(self): """ Examine the correlation between light matter scores and Z scores. """ scoreTypeX = 'Light matter score' scoreTypeY = 'Z score' for parameterSet in testArgs.parameterSets: affinity = _AFFINITY[parameterSet] dirName = makeOutputDir( DATASET, parameterSet, '%s-%s' % (FILESYSTEM_NAME[scoreTypeX], FILESYSTEM_NAME[scoreTypeY])) for subject in SUBJECTS: lmScores = [] zScores = [] for query in QUERIES: if query.id != subject.id: zScores.append(Z_SCORES[query.id][subject.id]) lmScore = getScore(affinity, query.id, subject.id) lmScores.append(lmScore) plot(lmScores, zScores, subject.id, scoreTypeX, scoreTypeY, dirName)