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)