def __init__(self, analysisLine, genome, trackName1, trackName2, reversed=False):
#print 'IN ANALYSIS: ',analysisLine
AnalysisDefHandler.__init__(self, analysisLine, reversed)
self._genome = genome
self._setTracks(trackName1, trackName2)
self._useConvertersFromId()
self._validStatClass = None
def getOptionsBoxMcfdrDepth(cls, prevChoices):
if not prevChoices.isBasic:
if prevChoices.analysisName == cls.Q3:
return AnalysisDefHandler(REPLACE_TEMPLATES['$MCFDRv5$']
).getOptionsAsText().values()[0]
elif prevChoices.analysisName == cls.Q4:
return AnalysisDefHandler(REPLACE_TEMPLATES['$MCFDRv4$']
).getOptionsAsText().values()[0]
def getAnalysisDict():
analyses = {}
for category in ANALYSIS_SPECS.keys():
analyses[category] = dict(
zip(ANALYSIS_SPECS[category],
[AnalysisDefHandler(q) for q in ANALYSIS_SPECS[category]]))
return analyses
def doAnalysis(analysisSpec, analysisBins, trackStructure):
'''Performs an analysis,
as specified by analysisSpec object,
in each bin specified by analysisBins,
on data sets specified in tracks.
Typical usage:
analysisSpec = AnalysisSpec(AvgSegLenStat)
analysisSpec.addParameter("withOverlaps","no")
analysisBins = GlobalBinSource('hg18')
tracks = [ Track(['Genes and gene subsets','Genes','Refseq']) ]
results = doAnalysis(analysisSpec, analysisBins, tracks)
'''
# TODO: handle multiple tracks analysis
# assert len(tracks) in [1,2] #for now..
# in an API setting, exceptions should not generally be hidden.
# Maybe this should be optional.
# setupDebugModeAndLogging()
silenceRWarnings()
silenceNumpyWarnings()
if isinstance(trackStructure, TrackStructureV2):
analysisDef = AnalysisDefHandler(analysisSpec.getDefAfterChoices())
statClass = analysisDef._statClassList[0]
validStatClass = wrapClass(
statClass,
keywords=analysisDef.getChoices(filterByActivation=True))
job = StatJob(analysisBins, trackStructure, validStatClass)
else:
tracks = trackStructure
if len(tracks) > 2:
from gold.util.CommonConstants import MULTIPLE_EXTRA_TRACKS_SEPARATOR
analysisSpec.addParameter(
'extraTracks',
MULTIPLE_EXTRA_TRACKS_SEPARATOR.join([
'^'.join([quote(part) for part in x.trackName])
for x in tracks[2:]
]))
job = AnalysisDefJob(analysisSpec.getDefAfterChoices(),
tracks[0].trackName,
tracks[1].trackName if len(tracks) > 1 else None,
analysisBins,
galaxyFn=None)
res = job.run(printProgress=False) # printProgress should be optional?
return res
def _tryAnalysisDefForValidity(cls,
analysisDef,
genome,
trackName1,
trackName2,
tryReversed=True):
analysis = AnalysisDefHandler(analysisDef)
return cls._tryAnalysisForValidity(analysis,
genome,
trackName1,
trackName2,
tryReversed=tryReversed)
def _getHeader(self):
#analysisDef = self._results._analysisDef
#if analysisDef is not None:
# return str(analysisDef)
#else:
if self._results._analysisText is not None:
try:
return AnalysisDefHandler.splitAnalysisText(self._results._analysisText)[1].strip()
except AssertionError:
pass
trackName1, trackName2 = self._results.getTrackNames()
return self._results.getStatClassName() + \
(' between ' if not trackName2 in [None,[]] else ' on ') + ':'.join(trackName1) + \
((' and ' + ':'.join(trackName2)) if not trackName2 in [None,[]] else '')
def _getHeader(self):
#analysisDef = self._results._analysisDef
#if analysisDef is not None:
# return str(analysisDef)
#else:
if self._results._analysisText is not None:
try:
return AnalysisDefHandler.splitAnalysisText(
self._results._analysisText)[1].strip()
except AssertionError:
pass
trackName1, trackName2 = self._results.getTrackNames()
return self._results.getStatClassName() + \
(' between ' if not trackName2 in [None,[]] else ' on ') + ':'.join(trackName1) + \
((' and ' + ':'.join(trackName2)) if not trackName2 in [None,[]] else '')
def _filterAndParseAnalyses(cls, analysisList):
from config.Config import IS_EXPERIMENTAL_INSTALLATION
analysisDict = dict(
zip(analysisList, [AnalysisDefHandler(q) for q in analysisList]))
if not IS_EXPERIMENTAL_INSTALLATION:
experimentalAnalyses = [
analysisDef
for analysisDef, analysis in analysisDict.iteritems()
if analysis.getChoice('isExperimental') == 'True'
]
for analysisDef in experimentalAnalyses:
del analysisDict[analysisDef]
return analysisDict
def execute(cls, choices, galaxyFn=None, username=''):
cls._setDebugModeIfSelected(choices)
genome = choices.genome
genomicRegions = choices.genomicRegions
genomicRegionsTracks = choices.genomicRegionsTracks
sourceTfs = choices.sourceTfs
sourceTfsDetails = choices.sourceTfsDetails
tfTracks = choices.tfTracks
# Get Genomic Region track name:
if genomicRegions == cls.REGIONS_FROM_HISTORY:
galaxyTN = genomicRegionsTracks.split(':')
genElementTrackName = ExternalTrackManager.getPreProcessedTrackFromGalaxyTN(
genome, galaxyTN)
#queryGSuite = getGSuiteFromGalaxyTN(genomicRegionsTracks)
#queryTrackList = [Track(x.trackName, x.title) for x in queryGSuite.allTracks()]
elif genomicRegions == 'Hyperbrowser repository':
selectedGenRegTrack = TfbsTrackNameMappings.getTfbsTrackNameMappings(
genome)[genomicRegionsTracks]
if isinstance(selectedGenRegTrack, dict):
genElementTrackName = selectedGenRegTrack.values()
else:
genElementTrackName = selectedGenRegTrack
elif genomicRegions == 'Hyperbrowser repository (cell-type-specific)':
genElementTrackName = ['Private', 'Antonio'
] + genomicRegionsTracks.split(':')
else:
return
# Get TF track names:
if isinstance(tfTracks, dict):
selectedTfTracks = [
key for key, val in tfTracks.iteritems() if val == 'True'
]
else:
selectedTfTracks = [tfTracks]
queryTrackTitle = '--'.join(genElementTrackName)
trackTitles = [queryTrackTitle]
tracks = [Track(genElementTrackName, trackTitle=queryTrackTitle)]
for i in selectedTfTracks:
if sourceTfs == 'Hyperbrowser repository':
tfTrackName = TfTrackNameMappings.getTfTrackNameMappings(
genome)[sourceTfsDetails] + [i]
tracks.append(
Track(tfTrackName,
trackTitle=tfTrackName[len(tfTrackName) - 1]))
trackTitles.append(tfTrackName[len(tfTrackName) - 1])
else:
tfTrackName = i.split(':')
queryGSuite = getGSuiteFromGalaxyTN(sourceTfsDetails)
for x in queryGSuite.allTracks():
selectedTrackNames = (':'.join(x.trackName))
if i == selectedTrackNames:
tracks.append(Track(x.trackName, x.title))
trackTitles.append(x.trackName[-1])
# queryGSuite = getGSuiteFromGalaxyTN(sourceTfsDetails)
# tfTrackName = [x.trackName for x in queryGSuite.allTracks()] + [i]
# tracks += [Track(x.trackName, x.title) for x in queryGSuite.allTracks()]
# trackTitles += tfTrackName
# print tfTrackName
# print tracks
# print trackTitles
trackTitlesForStat = trackTitles
trackTitles = CommonConstants.TRACK_TITLES_SEPARATOR.join(trackTitles)
##first statistic for Q2
resultsForStatistics = OrderedDict()
similarityFunc = [ #GSuiteStatUtils.T7_RATIO_OF_OBSERVED_TO_EXPECTED_OVERLAP,
GSuiteStatUtils.T5_RATIO_OF_OBSERVED_TO_EXPECTED_OVERLAP
]
for similarityStatClassName in similarityFunc:
regSpec, binSpec = UserBinMixin.getRegsAndBinsSpec(choices)
analysisBins = GalaxyInterface._getUserBinSource(regSpec,
binSpec,
genome=genome)
mcfdrDepth = AnalysisDefHandler(
REPLACE_TEMPLATES['$MCFDR$']).getOptionsAsText().values()[0][0]
analysisDefString = REPLACE_TEMPLATES[
'$MCFDR$'] + ' -> GSuiteSimilarityToQueryTrackRankingsAndPValuesWrapperStat'
analysisSpec = AnalysisDefHandler(analysisDefString)
analysisSpec.setChoice('MCFDR sampling depth', mcfdrDepth)
analysisSpec.addParameter('assumptions',
'PermutedSegsAndIntersegsTrack_')
analysisSpec.addParameter(
'rawStatistic', GSuiteStatUtils.
PAIRWISE_STAT_LABEL_TO_CLASS_MAPPING[similarityStatClassName])
analysisSpec.addParameter(
'pairwiseStatistic', GSuiteStatUtils.
PAIRWISE_STAT_LABEL_TO_CLASS_MAPPING[similarityStatClassName]
) #needed for call of non randomized stat for assertion
analysisSpec.addParameter('tail', 'more')
analysisSpec.addParameter('trackTitles',
trackTitles) #that need to be string
analysisSpec.addParameter('queryTracksNum', str(len(tracks)))
results = doAnalysis(analysisSpec, analysisBins,
tracks).getGlobalResult()
if not similarityStatClassName in resultsForStatistics:
resultsForStatistics[similarityStatClassName] = {}
resultsForStatistics[similarityStatClassName] = results
keyTitle = [
#'Normalized ratio of observed to expected overlap (normalized Forbes similarity measure)',
'Ratio of observed to expected overlap (Forbes similarity measure)'
]
# 'Normalized Forbes coefficient: ratio of observed to expected overlap normalized in relation to the reference GSuite',
# 'Forbes coefficient: ratio of observed to expected overlap'
keyTitle = [
#GSuiteStatUtils.T7_RATIO_OF_OBSERVED_TO_EXPECTED_OVERLAP,
GSuiteStatUtils.T5_RATIO_OF_OBSERVED_TO_EXPECTED_OVERLAP
]
resultDict = AllTfsOfRegions.countStatistics(similarityFunc, choices,
genome, tracks,
trackTitlesForStat)
resultDictShow = AllTfsOfRegions.countStatisticResults(
resultDict, keyTitle, trackTitlesForStat)
# print resultsForStatistics
'''selectedTrackNames = []
if sourceTfs == 'History (user-defined)':
if selectedTfTracks.split(":")[1] == "gsuite":
gSuite = getGSuiteFromGalaxyTN(selectedTfTracks)
for track in gSuite.allTracks():
selectedTrackNames.append(track.trackName)
else:
galaxyTN = selectedTfTracks.split(':')
gRegTrackName = ExternalTrackManager.getPreProcessedTrackFromGalaxyTN(genome, galaxyTN)
selectedTrackNames.append(gRegTrackName)
else:'''
tfNameList = []
#Intersection between TF Tracks and selected region (Table 1):
n = 0
allTargetBins = []
alltfNames = []
table1 = []
for i in selectedTfTracks:
n = n + 1
#newGalaxyFn = galaxyFn.split(".")[0] + str(n) + "." + "dat"
if sourceTfs == 'Hyperbrowser repository':
tfTrackName = TfTrackNameMappings.getTfTrackNameMappings(
genome)[sourceTfsDetails] + [i]
else:
tfTrackName = i.split(':')
tfTrackName.pop(0)
#tfIntersection.expandReferenceTrack(upFlankSize, downFlankSize)
tfIntersection = TrackIntersection(genome, genElementTrackName,
tfTrackName, galaxyFn, str(n))
regFileNamer = tfIntersection.getIntersectedRegionsStaticFileWithContent(
)
targetBins = tfIntersection.getIntersectedReferenceBins()
#regSpec, targetBins = UserBinSelector.getRegsAndBinsSpec(choices)
tfHits = [i] * len(targetBins)
fixedTargetBins = [str(a).split(" ")[0] for a in targetBins]
extendedTargetBins = [
list(a) for a in zip(fixedTargetBins, tfHits)
]
allTargetBins = allTargetBins + extendedTargetBins
tfName = i
alltfNames = alltfNames + [tfName]
# Save output table:
tfNameList.append(tfName)
line = [tfName] + [len(targetBins)] + [
regFileNamer.getLink('Download bed-file')
] + [
regFileNamer.getLoadToHistoryLink('Send bed-file to History')
]
table1 = table1 + [line]
# Computing totals:
fullCase = ','.join(alltfNames)
firstColumn = [item[0] for item in allTargetBins]
uniqueAllTargetBins = list(set(firstColumn))
# Group TFs by bound region:
d1 = defaultdict(list)
for k, v in allTargetBins:
d1[k].append(v)
allTFTargetBins = dict((k, ','.join(v)) for k, v in d1.iteritems())
allTFTargetList = []
fullCaseTFTargetList = []
for key, value in allTFTargetBins.iteritems():
allTFTargetList = allTFTargetList + [[key, value]]
if value == fullCase:
fullCaseTFTargetList = fullCaseTFTargetList + [[key, value]]
analysis3 = TrackIntersection.getFileFromTargetBins(
allTFTargetList, galaxyFn, str(3))
analysis4 = TrackIntersection.getFileFromTargetBins(
fullCaseTFTargetList, galaxyFn, str(4))
# Print output to table:
title = 'TF targets and co-occupancy of ' + genElementTrackName[
-1] + ' genomic regions'
htmlCore = HtmlCore()
pf = plotFunction(tableId='resultsTable')
htmlCore.begin()
htmlCore.header(title)
htmlCore.divBegin('resultsDiv')
htmlCore.line(pf.createButton(bText='Show/Hide more results'))
# htmlCore.tableHeader(['Transcription Factor', 'Normalized ratio of observed to expected overlap (normalized Forbes similarity measure) -- Similarity to genomic regions track', 'Normalized ratio of observed to expected overlap (normalized Forbes similarity measure) -- p-value','Ratio of observed to expected overlap (Forbes similarity measure) -- Similarity to genomic regions track', 'Ratio of observed to expected overlap (Forbes similarity measure) -- p-value', 'Number of TF-Target Track Regions', 'File of TF Target Regions', 'File of TF Target Regions', 'Number of TF-co-occupied Regions', 'File of TF co-occupied Regions', 'File of TF co-occupied Regions', 'Rank of TF co-occupancy motifs', 'Rank of TF co-occupancy motifs'], sortable=True, tableId='resultsTable')
#previous ordering
# htmlCore.tableHeader(['Transcription Factor', 'Normalized Forbes index --overlap score',
# 'Normalized Forbes index --p-value',
# 'Forbes index --overlap score', 'Forbes index --p-value',
# 'Number of TF-Target Track Regions', 'File of TF Target Regions',
# 'File of TF Target Regions', 'Number of target track regions occupied by this TF',
# 'File of TF co-occupied Regions', 'File of TF co-occupied Regions',
# 'Rank of TF co-occupancy motifs', 'Rank of TF co-occupancy motifs'],
# sortable=True, tableId='resultsTable')
htmlCore.tableHeader(
[
'Transcription Factor',
'Number of TF-Target Track Regions',
'File of TF Track Regions',
'Number of target track regions occupied by this TF',
'File of TF Target Regions',
'Forbes index --overlap score',
'Forbes index --p-value',
#'Normalized Forbes index --overlap score', 'Normalized Forbes index --p-value',
'File of TF co-occupied Regions',
'Rank of TF co-occupancy motifs'
],
sortable=True,
tableId='resultsTable')
# Adding co-occupancy results to table:
n = 1000
genRegionNumElements = [
int(x) for x in getTrackRelevantInfo.getNumberElements(
genome, genElementTrackName)
]
for key0, it0 in resultsForStatistics.iteritems():
for el in tfNameList:
if el not in it0:
resultsForStatistics[key0][el] = [None, None]
resultsPlotDict = {}
resultPlotCat = []
resultsPlot = []
resultsForStatisticsProper = {}
for key0, it0 in resultsForStatistics.iteritems():
if not key0 in resultsPlotDict:
resultsPlotDict[key0] = {}
resultsPlotPart = []
for key1, it1 in it0.iteritems():
resultsPlotPart.append(it1[0])
if not key1 in resultsForStatisticsProper:
resultsForStatisticsProper[key1] = []
if not key1 in resultsPlotDict[key0]:
resultsPlotDict[key0][key1] = None
for el in it1:
resultsForStatisticsProper[key1].append(el)
resultsPlotDict[key0][key1] = it1[0]
resultPlotCat.append(tfNameList)
resultPlotCat.append(tfNameList)
#resultPlotCatPart = tfNameList
# print resultPlotCatPart
for key0, it0 in resultsPlotDict.iteritems():
resultsPlotPart = []
for el in tfNameList:
if el in it0:
resultsPlotPart.append(it0[el])
else:
resultsPlotPart.append(None)
resultsPlot.append(resultsPlotPart)
for i in table1:
thisCaseTFTargetList = []
for key, value in allTFTargetList:
if i[0] in value and ',' in value:
thisCaseTFTargetList = thisCaseTFTargetList + [[
key, value
]]
n = n + 1
thisAnalysis = TrackIntersection.getFileFromTargetBins(
thisCaseTFTargetList, galaxyFn, str(n))
thisCaseCoCountsList = []
thing = [x[1] for x in thisCaseTFTargetList]
for k in list(set(thing)):
thisCount = thing.count(k)
thisCaseCoCountsList = thisCaseCoCountsList + \
[[k, thisCount, 100*float(thisCount)/float(sum(genRegionNumElements)), 100*float(thisCount)/float(len(thisCaseTFTargetList))]]
thisCaseCoCountsList.sort(key=lambda x: x[2], reverse=True)
n = n + 1
thisCoCountsAnalysis = TrackIntersection.getOccupancySummaryFile(
thisCaseCoCountsList, galaxyFn, str(n))
thisLine = [len(thisCaseTFTargetList)] + \
[thisAnalysis.getLink('Download file')] + [thisAnalysis.getLoadToHistoryLink('Send file to History')] + \
[thisCoCountsAnalysis.getLink('Download file')] + [thisCoCountsAnalysis.getLoadToHistoryLink('Send file to History')]
newLineI = []
tfName = i[0]
newLineI.append(tfName)
for el in resultsForStatisticsProper[tfName]:
newLineI.append(el)
for elN in range(1, len(i)):
newLineI.append(i[elN])
# htmlCore.tableLine(i + thisLine)
# htmlCore.tableHeader(['Transcription Factor', 'Normalized Forbes index --overlap score',
# 'Normalized Forbes index --p-value',
# 'Forbes index --overlap score', 'Forbes index --p-value',
# 'Number of TF-Target Track Regions', 'File of TF Target Regions',
# 'File of TF Target Regions', 'Number of target track regions occupied by this TF',
# 'File of TF co-occupied Regions', 'File of TF co-occupied Regions',
# 'Rank of TF co-occupancy motifs', 'Rank of TF co-occupancy motifs'],
# sortable=True, tableId='resultsTable')
# htmlCore.tableHeader(['Transcription Factor', 'Number of TF-Target Track Regions', 'File of TF Track Regions',
# 'Number of target track regions occupied by this TF', 'File of TF Target Regions',
# 'Forbes index --overlap score', 'Forbes index --p-value',
# 'Normalized Forbes index --overlap score', 'Normalized Forbes index --p-value',
# 'File of TF co-occupied Regions', 'Rank of TF co-occupancy motifs'],
# sortable=True, tableId='resultsTable')
tl = newLineI + thisLine
# previous ordering tl - 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12
# actual ordering - 0, 5, 7, 8, 7, 3, 4, 1, 2, 9, 11
#ordering = [0, 5, 7, 8, 10, 3, 4, 1, 2, 10, 12]
ordering = [0, 3, 5, 6, 8, 1, 2, 8, 10]
#1, 2, => delete
eoList = []
for eo in ordering:
eoList.append(tl[eo])
htmlCore.tableLine(eoList)
totalCoOccupancyTargetList = []
n = 2000
for key, value in allTFTargetList:
n = n + 1
if ',' in value:
totalCoOccupancyTargetList = totalCoOccupancyTargetList + [[
key, value
]]
#newGalaxyFn = galaxyFn.split(".")[0] + str(n) + "." + "dat"
totalCoOccupancyAnalysis = TrackIntersection.getFileFromTargetBins(
totalCoOccupancyTargetList, galaxyFn, str(n))
#line = ['Total reported regions'] + [len(allTargetBins)] + [''] + [''] + [''] + [''] + ['']
#line = ['Full co-occupancy of ' + fullCase] + ['-'] + ['-'] + ['-'] + ['-'] + ['-'] + ['-'] + ['-'] + [len(fullCaseTFTargetList)] + [analysis4.getLink('Download file')] + [analysis4.getLoadToHistoryLink('Send file to History')] + ['-'] + ['-']
line = ['Full co-occupancy of ' + fullCase] + \
['-'] + \
['-'] + \
[len(fullCaseTFTargetList)] + \
['-'] + \
['-'] + \
['-'] + \
[analysis4.getLoadToHistoryLink('Send file to History')] + \
['-']
htmlCore.tableLine(line)
#line = ['Total unique regions'] + ['-'] + ['-'] + ['-'] + ['-'] + [len(allTFTargetList)] + [analysis3.getLink('Download bed-file')] + [analysis3.getLoadToHistoryLink('Send bed-file to History')] + [len(totalCoOccupancyTargetList)] + [totalCoOccupancyAnalysis.getLink('Download file')] + [totalCoOccupancyAnalysis.getLoadToHistoryLink('Send file to History')] + ['-'] + ['-']
line = ['Total unique regions'] + \
[len(allTFTargetList)] + \
['-'] + \
[len(totalCoOccupancyTargetList)] + \
[analysis3.getLoadToHistoryLink('Send bed-file to History')] + \
['-'] +\
['-'] + \
[totalCoOccupancyAnalysis.getLoadToHistoryLink('Send file to History')] + \
['-']
htmlCore.tableLine(line)
htmlCore.tableFooter()
htmlCore.divEnd()
# htmlCore.line(pf.hideColumns(indexList=[2, 4]))
#
sumRes = 0
for r in resultsPlot[0]:
if r != None:
sumRes += r
if sumRes != 0:
vg = visualizationGraphs()
result = vg.drawColumnCharts(
[resultsPlot[0]],
height=300,
categories=resultPlotCat,
legend=False,
addOptions='width: 90%; float:left; margin: 0 4%;',
#titleText=['Overlap between TFs and genomic region using normalized Forbes', 'Overlap between TFs and genomic region using Forbes'],
titleText=[
'Overlap between TFs and genomic region using Forbes'
],
xAxisRotation=90,
xAxisTitle='TF',
yAxisTitle='value')
htmlCore.line(result)
for key0, it0 in resultDictShow.iteritems():
htmlCore.divBegin('resultsDiv' + str(key0))
htmlCore.header(key0)
htmlCore.tableHeader(it0[0],
sortable=True,
tableId='resultsTable' + str(key0))
for elN in range(1, len(it0)):
htmlCore.tableLine(it0[elN])
htmlCore.tableFooter()
htmlCore.divEnd()
htmlCore.hideToggle(styleClass='debug')
htmlCore.end()
print htmlCore
def __init__(self, choices, galaxyFn=None, username=''):
'''
Is called when execute-button is pushed by web-user. Should print
output as HTML to standard out, which will be directed to a results page
in Galaxy history. If getOutputFormat is anything else than HTML, the
output should be written to the file with path galaxyFn. If needed,
StaticFile can be used to get a path where additional files can be put
(e.g. generated image files). choices is a list of selections made by
web-user in each options box.
'''
analysisDefString = REPLACE_TEMPLATES[
'$MCFDRv3$'] + ' -> CollectionBinnedHypothesisWrapperStat'
analysisSpec = AnalysisDefHandler(analysisDefString)
analysisSpec.setChoice('MCFDR sampling depth', choices.mcfdrDepth)
analysisSpec.addParameter('assumptions',
'PermutedSegsAndIntersegsTrack_')
if choices.question == "question 8":
analysisSpec.addParameter('rawStatistic',
'MultitrackRawBinnedOverlapV2Stat')
else:
analysisSpec.addParameter('rawStatistic',
'MultitrackRawSingleBinV2Stat')
# analysisSpec.addParameter('pairwiseStatistic', choices.stat)
analysisSpec.addParameter('summaryFunc', choices.summaryFunc)
analysisSpec.addParameter('tail', 'right-tail')
analysisSpec.addParameter('localBinSize', choices.binSize)
analysisSpec.addParameter('question', choices.question)
analysisBins = GlobalBinSource(choices.genome)
gsuite = getGSuiteFromGalaxyTN(choices.tracks)
tracks = [Track(x.trackName) for x in gsuite.allTracks()]
tracks = tracks[0:2]
results = doAnalysis(analysisSpec, analysisBins, tracks)
print results
def countStatistics(similarityFunc, choices, genome, tracks, trackTitles):
trackList = tracks[1:]
resultsForStatistics = OrderedDict()
llDict = OrderedDict()
trackT = trackTitles[1:]
i = 0
for tt1 in trackT:
if not tt1 in llDict:
llDict[tt1] = []
resultsForStatistics[tt1] = {}
for tt2 in range(i, len(trackT)):
llDict[tt1].append(trackT[tt2])
i += 1
#
# print 'llDict=' + str(llDict)
# print 'trackT=' + str(trackT)
# print 'trackList=' + str(trackList)
for key0, it0 in llDict.iteritems():
if len(it0) > 1:
trackCollection = []
for it1 in it0:
trackNumber = trackT.index(it1)
trackCollection.append(trackList[trackNumber])
trackTitles = CommonConstants.TRACK_TITLES_SEPARATOR.join(it0)
# print str(key0) + '- trackCollection: ' + str(trackCollection) + ' trackTitles: ' + str(trackTitles)
for similarityStatClassName in similarityFunc:
regSpec, binSpec = UserBinMixin.getRegsAndBinsSpec(choices)
analysisBins = GalaxyInterface._getUserBinSource(
regSpec, binSpec, genome=genome)
mcfdrDepth = AnalysisDefHandler(
REPLACE_TEMPLATES['$MCFDR$']).getOptionsAsText(
).values()[0][0]
analysisDefString = REPLACE_TEMPLATES[
'$MCFDR$'] + ' -> GSuiteSimilarityToQueryTrackRankingsAndPValuesWrapperStat'
analysisSpec = AnalysisDefHandler(analysisDefString)
analysisSpec.setChoice('MCFDR sampling depth', mcfdrDepth)
analysisSpec.addParameter(
'assumptions', 'PermutedSegsAndIntersegsTrack_')
analysisSpec.addParameter(
'rawStatistic',
GSuiteStatUtils.PAIRWISE_STAT_LABEL_TO_CLASS_MAPPING[
similarityStatClassName])
analysisSpec.addParameter(
'pairwiseStatistic',
GSuiteStatUtils.PAIRWISE_STAT_LABEL_TO_CLASS_MAPPING[
similarityStatClassName]
) #needed for call of non randomized stat for assertion
analysisSpec.addParameter('tail', 'more')
analysisSpec.addParameter(
'trackTitles', trackTitles) #that need to be string
#i added that one later
analysisSpec.addParameter('queryTracksNum',
str(len(trackCollection)))
results = doAnalysis(analysisSpec, analysisBins,
trackCollection).getGlobalResult()
if not similarityStatClassName in resultsForStatistics[
key0]:
resultsForStatistics[key0][
similarityStatClassName] = {}
resultsForStatistics[key0][
similarityStatClassName] = results
return resultsForStatistics
def getOptionsBoxMcfdrDepth(cls, prevChoices):
if not prevChoices.isBasic and prevChoices.analysisName == cls.Q2:
analysisSpec = AnalysisDefHandler(REPLACE_TEMPLATES['$MCFDR$'])
return analysisSpec.getOptionsAsText().values()[0]
def execute(cls, choices, galaxyFn=None, username=''):
'''
Is called when execute-button is pushed by web-user. Should print
output as HTML to standard out, which will be directed to a results
page in Galaxy history. If getOutputFormat is anything else than HTML,
the output should be written to the file with path galaxyFn. If needed,
StaticFile can be used to get a path where additional files can be put
(e.g. generated image files). choices is a list of selections made by
web-user in each options box.
'''
import warnings
#warnings.simplefilter('error')
cls._setDebugModeIfSelected(choices)
similarityStatClassName = choices.similarityFunc if choices.similarityFunc else GSuiteStatUtils.T5_RATIO_OF_OBSERVED_TO_EXPECTED_OVERLAP
summaryFunc = choices.summaryFunc if choices.summaryFunc else cls.SUMMARY_FUNC_DEFAULT
pairwiseStatName = GSuiteStatUtils.PAIRWISE_STAT_LABEL_TO_CLASS_MAPPING[similarityStatClassName]
gsuite = getGSuiteFromGalaxyTN(choices.gsuite)
tracks = [Track(x.trackName) for x in gsuite.allTracks()]
statTxt = "Average"
if(summaryFunc == "max"): statTxt = "Maximum"
if choices.analysisName == cls.Q2:
mcfdrDepth = choices.mcfdrDepth if choices.mcfdrDepth else AnalysisDefHandler(REPLACE_TEMPLATES['$MCFDR$']).getOptionsAsText().values()[0][0]
# First compute pvalue by running the statistic through a wrapper stat that computes the max per bin
#from quick.statistic.CollectionBinnedHypothesisWrapperStat import CollectionBinnedHypothesisWrapperStat
#analysisSpec = AnalysisSpec(CollectionBinnedHypothesisWrapperStat)
analysisDefString = REPLACE_TEMPLATES['$MCFDRv3$'] + ' -> CollectionBinnedHypothesisWrapperStat'
analysisSpec = AnalysisDefHandler(analysisDefString)
analysisSpec.setChoice('MCFDR sampling depth', mcfdrDepth)
analysisSpec.addParameter("rawStatistic", "GenericMaxBinValueStat")
# analysisSpec.addParameter('perBinStatistic', 'SummarizedStat')
analysisSpec.addParameter('perBinStatistic', 'MultitrackSummarizedInteractionV2Stat')
# analysisSpec.addParameter('mcSamplerClass', 'NaiveMCSamplingV2Stat')
analysisSpec.addParameter('pairwiseStatistic', 'ObservedVsExpectedStat')
analysisSpec.addParameter('summaryFunc', summaryFunc)
# analysisSpec.addParameter('evaluatorFunc','evaluatePvalueAndNullDistribution')
analysisSpec.addParameter('tail', 'right-tail')
analysisSpec.addParameter('assumptions', 'RandomGenomeLocationTrack')
#analysisSpec.addParameter('maxSamples', 10)
analysisSpec.addParameter('multitrackSummaryFunc', summaryFunc)
regSpec, binSpec = cls.getRegsAndBinsSpec(choices)
analysisBins = GalaxyInterface._getUserBinSource(regSpec,
binSpec,
choices.genome)
results = doAnalysis(analysisSpec, analysisBins, tracks)
results = results.getGlobalResult()
resultsTxt = "The highest ranking bin based on the " + statTxt.lower() + " of the Forbes similarity measure for pairs of tracks within each bin had a score of <b>%.3f</b> with p-value <b>%.6f</b>" % (results["TSMC_GenericMaxBinValueStat"], results['P-value'])
# Stat question 7
core = HtmlCore()
core.begin()
analysisSpec = AnalysisSpec(MultitrackSummarizedInteractionWrapperStat)
#analysisSpec.addParameter('pairwiseStatistic', 'ObservedVsExpectedStat')
analysisSpec.addParameter('pairwiseStatistic', GSuiteStatUtils.PAIRWISE_STAT_LABEL_TO_CLASS_MAPPING[similarityStatClassName])
analysisSpec.addParameter('summaryFunc', summaryFunc)
analysisSpec.addParameter('multitrackSummaryFunc', summaryFunc)
gsuite = getGSuiteFromGalaxyTN(choices.gsuite)
tracks = [Track(x.trackName) for x in gsuite.allTracks()]
regSpec, binSpec = cls.getRegsAndBinsSpec(choices)
analysisBins = GalaxyInterface._getUserBinSource(regSpec,
binSpec,
choices.genome)
results = doAnalysis(analysisSpec, analysisBins, tracks)
#print '<br>results: ', results, '<br><br>'
prettyResults = OrderedDict()
for key, val in results.iteritems():
if "Result" in val.keys():
prettyResults[key] = val["Result"]
else:
prettyResults[key] = "No result"
core.header(statTxt + " co-occurence between pairs of tracks within each bin")
if choices.analysisName == cls.Q2:
core.paragraph(resultsTxt)
core.divBegin(divClass='resultsExplanation')
core.paragraph('The following is a list of all bins and the <b>' + statTxt.lower() + '</b> co-occurrence of tracks within each bin.')
core.divEnd()
"""
core.paragraph('''
Suite data is coinciding the most in bin %s
''' % ('test'))
"""
visibleRows = 20
makeTableExpandable = len(prettyResults) > visibleRows
columnNames = ['Bin', 'Co-occurrence within the bin']
if choices.analysisName == cls.Q1:
shortQuestion = cls.Q1_SHORT
else:
shortQuestion = cls.Q2_SHORT
addTableWithTabularAndGsuiteImportButtons(
core, choices, galaxyFn, shortQuestion, tableDict=prettyResults,
columnNames=columnNames, sortable=True, presorted=0,
expandable=makeTableExpandable, visibleRows=visibleRows)
core.divEnd()
core.end()
print str(core)
def getOptionsBoxMcfdrDepth(cls, prevChoices):
return AnalysisDefHandler(REPLACE_TEMPLATES['$MCFDRv5$']).getOptionsAsText().values()[0]
def __init__(self, analysisLine, genome, trackName1, trackName2, reversed=False):
#print 'IN ANALYSIS: ',analysisLine
AnalysisDefHandler.__init__(self, analysisLine, reversed)
self._genome = genome
self._setTracks(trackName1, trackName2)
self._useConvertersFromId()
def _prepareAnalysisWithHypothesisTests(cls, choices, localAnalysis):
mcfdrDepth = choices.mcfdrDepth if choices.mcfdrDepth else \
AnalysisDefHandler(REPLACE_TEMPLATES['$MCFDRv5$']).getOptionsAsText().values()[0][0]
analysisDefString = REPLACE_TEMPLATES['$MCFDRv5$'] + ' -> ' + ' -> MultipleRandomizationManagerStat'
analysisSpec = AnalysisDefHandler(analysisDefString)
analysisSpec.setChoice('MCFDR sampling depth', mcfdrDepth)
analysisSpec.addParameter('rawStatistic', PairedTSStat.__name__)
analysisSpec.addParameter('pairedTsRawStatistic', ObservedVsExpectedStat.__name__)
analysisSpec.addParameter('tail', 'right-tail')
analysisSpec.addParameter('evaluatorFunc', 'evaluatePvalueAndNullDistribution')
analysisSpec.addParameter('runLocalAnalysis', 'Yes' if localAnalysis else 'No')
return analysisSpec
def execute(cls, choices, galaxyFn=None, username=''):
'''
Is called when execute-button is pushed by web-user. Should print
output as HTML to standard out, which will be directed to a results page
in Galaxy history. If getOutputFormat is anything else than HTML, the
output should be written to the file with path galaxyFn. If needed,
StaticFile can be used to get a path where additional files can be put
(e.g. generated image files). choices is a list of selections made by
web-user in each options box.
'''
import numpy
numpy.seterr(all='raise')
cls._setDebugModeIfSelected(choices)
# DebugUtil.insertBreakPoint(username=username, currentUser='******')
genome = choices.genome
analysisQuestion = choices.analysisName
similaryStatClassName = choices.similarityFunc if choices.similarityFunc else GSuiteStatUtils.T5_RATIO_OF_OBSERVED_TO_EXPECTED_OVERLAP
summaryFunc = choices.summaryFunc if choices.summaryFunc else 'average'
reverse = 'Yes' if choices.reversed else 'No'
gsuite = getGSuiteFromGalaxyTN(choices.gsuite)
regSpec, binSpec = UserBinMixin.getRegsAndBinsSpec(choices)
analysisBins = GalaxyInterface._getUserBinSource(regSpec,
binSpec,
genome=genome)
tracks = [
Track(x.trackName, trackTitle=x.title) for x in gsuite.allTracks()
]
trackTitles = CommonConstants.TRACK_TITLES_SEPARATOR.join(
[quote(x.title, safe='') for x in gsuite.allTracks()])
additionalResultsDict = OrderedDict()
additionalAttributesDict = OrderedDict()
if analysisQuestion in [cls.Q1, cls.Q2, cls.Q3]:
additionalAttributesDict = cls.getSelectedAttributesForEachTrackDict(
choices.additionalAttributes, gsuite)
#additional analysis
stats = [CountStat, CountElementStat]
additionalResultsDict = runMultipleSingleValStatsOnTracks(
gsuite, stats, analysisBins, queryTrack=None)
if analysisQuestion == cls.Q1:
analysisSpec = AnalysisSpec(
GSuiteRepresentativenessOfTracksRankingsWrapperStat)
analysisSpec.addParameter(
'pairwiseStatistic', GSuiteStatUtils.
PAIRWISE_STAT_LABEL_TO_CLASS_MAPPING[similaryStatClassName])
analysisSpec.addParameter(
'summaryFunc',
GSuiteStatUtils.SUMMARY_FUNCTIONS_MAPPER[summaryFunc])
analysisSpec.addParameter('reverse', reverse)
analysisSpec.addParameter('ascending', 'No')
analysisSpec.addParameter('trackTitles', trackTitles)
analysisSpec.addParameter('queryTracksNum', len(tracks))
results = doAnalysis(analysisSpec, analysisBins,
tracks).getGlobalResult()
gsPerTrackResultsModel = GSuitePerTrackResultModel(
results,
['Similarity to rest of tracks in suite (%s)' % summaryFunc],
additionalResultsDict=additionalResultsDict,
additionalAttributesDict=additionalAttributesDict)
if choices.leadAttribute and choices.leadAttribute != GSuiteConstants.TITLE_COL:
columnTitles, decoratedResultsDict = \
gsPerTrackResultsModel.generateColumnTitlesAndResultsDict(choices.leadAttribute)
else:
columnTitles, decoratedResultsDict = \
gsPerTrackResultsModel.generateColumnTitlesAndResultsDict()
core = HtmlCore()
core.begin()
core.divBegin(divId='results-page')
core.divBegin(divClass='results-section')
core.header(analysisQuestion)
topTrackTitle = results.keys()[0]
core.paragraph('''
The track "%s" is the most representative track of the GSuite with %s %s similarity to the rest of the tracks
as measured by "%s" track similarity measure.
''' % (topTrackTitle, results[topTrackTitle], summaryFunc,
similaryStatClassName))
addTableWithTabularAndGsuiteImportButtons(
core,
choices,
galaxyFn,
cls.Q1_SHORT,
decoratedResultsDict,
columnTitles,
gsuite=gsuite,
results=results,
gsuiteAppendAttrs=['similarity_score'],
sortable=True)
# plot
columnInd = 0
if choices.leadAttribute and choices.leadAttribute != GSuiteConstants.TITLE_COL:
columnInd = 1
res = GSuiteTracksCoincidingWithQueryTrackTool.drawPlot(
results,
additionalResultsDict,
'Similarity to rest of tracks in suite (%s)' % summaryFunc,
columnInd=columnInd)
core.line(res)
core.divEnd()
core.divEnd()
core.end()
# elif analysisQuestion == cls.Q2:
# analysisSpec = AnalysisSpec(GSuiteRepresentativenessOfTracksRankingsWrapperStat)
# analysisSpec.addParameter('pairwiseStatistic', GSuiteStatUtils.PAIRWISE_STAT_LABEL_TO_CLASS_MAPPING[similaryStatClassName])
# analysisSpec.addParameter('summaryFunc', GSuiteStatUtils.SUMMARY_FUNCTIONS_MAPPER[summaryFunc])
# analysisSpec.addParameter('reverse', reverse)
# analysisSpec.addParameter('ascending', 'Yes')
# analysisSpec.addParameter('trackTitles', trackTitles)
# results = doAnalysis(analysisSpec, analysisBins, tracks).getGlobalResult()
#
# gsPerTrackResultsModel = GSuitePerTrackResultModel(
# results, ['Similarity to rest of tracks in suite (%s)' % summaryFunc],
# additionalResultsDict=additionalResultsDict,
# additionalAttributesDict=additionalAttributesDict)
# if choices.leadAttribute and choices.leadAttribute != GSuiteConstants.TITLE_COL:
# columnTitles, decoratedResultsDict = \
# gsPerTrackResultsModel.generateColumnTitlesAndResultsDict(choices.leadAttribute)
# else:
# columnTitles, decoratedResultsDict = \
# gsPerTrackResultsModel.generateColumnTitlesAndResultsDict()
#
# core = HtmlCore()
# core.begin()
# core.divBegin(divId='results-page')
# core.divBegin(divClass='results-section')
# core.header(analysisQuestion)
# topTrackTitle = results.keys()[0]
# core.paragraph('''
# The track "%s" is the most atypical track of the GSuite with %s %s similarity to the rest of the tracks
# as measured by the "%s" track similarity measure.
# ''' % (topTrackTitle, strWithNatLangFormatting(results[topTrackTitle]), summaryFunc, similaryStatClassName))
# # core.tableFromDictionary(results, columnNames=['Track title', 'Similarity to rest of tracks in suite (' + summaryFunc+')'], sortable=False)
#
# from quick.util import CommonFunctions
# rawDataURIList = CommonFunctions.getHyperlinksForRawTableData(
# dataDict=decoratedResultsDict, colNames=columnTitles,
# tableId="resultsTable", galaxyFn=galaxyFn)
# core.tableFromDictionary(decoratedResultsDict, columnNames=columnTitles, sortable=True,
# tableId='resultsTable', addInstruction=True,
# addRawDataSelectBox=True, rawDataURIList=rawDataURIList)
# # core.tableFromDictionary(decoratedResultsDict, columnNames=columnTitles, sortable=True, tableId='resultsTable')
#
# columnInd = 0
# if choices.leadAttribute and choices.leadAttribute != GSuiteConstants.TITLE_COL:
# columnInd = 1
# res = GSuiteTracksCoincidingWithQueryTrackTool.drawPlot(
# results, additionalResultsDict,
# 'Similarity to rest of tracks in suite (%s)' % summaryFunc,
# columnInd=columnInd)
# core.line(res)
# core.divEnd()
# core.divEnd()
# core.end()
#
# if choices.addResults == 'Yes':
# GSuiteStatUtils.addResultsToInputGSuite(
# gsuite, results, ['Similarity_score'],
# cls.extraGalaxyFn[GSUITE_EXPANDED_WITH_RESULT_COLUMNS_FILENAME])
elif analysisQuestion == cls.Q3:
mcfdrDepth = choices.mcfdrDepth if choices.mcfdrDepth else \
AnalysisDefHandler(REPLACE_TEMPLATES['$MCFDR$']).getOptionsAsText().values()[0][0]
analysisDefString = REPLACE_TEMPLATES[
'$MCFDRv3$'] + ' -> GSuiteRepresentativenessOfTracksRankingsAndPValuesWrapperStat'
analysisSpec = AnalysisDefHandler(analysisDefString)
analysisSpec.setChoice('MCFDR sampling depth', mcfdrDepth)
analysisSpec.addParameter('assumptions',
'PermutedSegsAndIntersegsTrack')
analysisSpec.addParameter(
'rawStatistic',
SummarizedInteractionWithOtherTracksV2Stat.__name__)
analysisSpec.addParameter(
'pairwiseStatistic', GSuiteStatUtils.
PAIRWISE_STAT_LABEL_TO_CLASS_MAPPING[similaryStatClassName])
analysisSpec.addParameter(
'summaryFunc',
GSuiteStatUtils.SUMMARY_FUNCTIONS_MAPPER[summaryFunc])
analysisSpec.addParameter('tail', 'right-tail')
analysisSpec.addParameter('trackTitles', trackTitles)
results = doAnalysis(analysisSpec, analysisBins,
tracks).getGlobalResult()
core = HtmlCore()
gsPerTrackResultsModel = GSuitePerTrackResultModel(
results, [
'Similarity to rest of tracks in suite (%s)' % summaryFunc,
'P-value'
],
additionalResultsDict=additionalResultsDict,
additionalAttributesDict=additionalAttributesDict)
if choices.leadAttribute and choices.leadAttribute != GSuiteConstants.TITLE_COL:
columnTitles, decoratedResultsDict = \
gsPerTrackResultsModel.generateColumnTitlesAndResultsDict(choices.leadAttribute)
else:
columnTitles, decoratedResultsDict = \
gsPerTrackResultsModel.generateColumnTitlesAndResultsDict()
core.begin()
core.divBegin(divId='results-page')
core.divBegin(divClass='results-section')
core.header(analysisQuestion)
topTrackTitle = results.keys()[0]
core.paragraph('''
The track "%s" has the lowest P-value of %s corresponding to %s %s similarity to the rest of the tracks
as measured by "%s" track similarity measure.
''' % (topTrackTitle,
strWithNatLangFormatting(results[topTrackTitle][1]),
strWithNatLangFormatting(results[topTrackTitle][0]),
summaryFunc, similaryStatClassName))
# core.tableFromDictionary(results, columnNames=['Track title', 'Similarity to rest of tracks in suite (' + summaryFunc+')', 'P-value'], sortable=False)
addTableWithTabularAndGsuiteImportButtons(
core,
choices,
galaxyFn,
cls.Q3_SHORT,
decoratedResultsDict,
columnTitles,
gsuite=gsuite,
results=results,
gsuiteAppendAttrs=['similarity_score', 'p_value'],
sortable=True)
core.divEnd()
core.divEnd()
core.end()
else: # Q4
mcfdrDepth = choices.mcfdrDepth if choices.mcfdrDepth else \
AnalysisDefHandler(REPLACE_TEMPLATES['$MCFDR$']).getOptionsAsText().values()[0][0]
analysisDefString = REPLACE_TEMPLATES[
'$MCFDRv3$'] + ' -> CollectionSimilarityHypothesisWrapperStat'
analysisSpec = AnalysisDefHandler(analysisDefString)
analysisSpec.setChoice('MCFDR sampling depth', mcfdrDepth)
analysisSpec.addParameter('assumptions',
'PermutedSegsAndIntersegsTrack')
analysisSpec.addParameter('rawStatistic',
'MultitrackSummarizedInteractionV2Stat')
analysisSpec.addParameter(
'pairwiseStatistic', GSuiteStatUtils.
PAIRWISE_STAT_LABEL_TO_CLASS_MAPPING[similaryStatClassName])
analysisSpec.addParameter(
'summaryFunc',
GSuiteStatUtils.SUMMARY_FUNCTIONS_MAPPER[summaryFunc])
analysisSpec.addParameter('multitrackSummaryFunc',
'avg') # should it be a choice?
analysisSpec.addParameter('tail', 'right-tail')
results = doAnalysis(analysisSpec, analysisBins,
tracks).getGlobalResult()
pval = results['P-value']
observed = results['TSMC_MultitrackSummarizedInteractionV2Stat']
significanceLevel = 'strong' if pval < 0.01 else (
'weak' if pval < 0.05 else 'no')
core = HtmlCore()
core.begin()
core.divBegin(divId='results-page')
core.divBegin(divClass='results-section')
core.header(analysisQuestion)
core.paragraph('''
The tracks in the suite show %s significance in their collective similarity
(average similarity of a track to the rest) of %s
and corresponding p-value of %s,
as measured by "%s" track similarity measure.
''' % (significanceLevel, strWithNatLangFormatting(observed),
strWithNatLangFormatting(pval), similaryStatClassName))
core.divEnd()
core.divEnd()
core.end()
print str(core)
def getOptionsBoxMcfdrDepth(prevChoices):
analysisSpec = AnalysisDefHandler(REPLACE_TEMPLATES['$MCFDR$'])
return analysisSpec.getOptionsAsText().values()[0]
def prepareQ2(cls, choices, similarityStatClassName, trackTitles, randStrat, intensityTrackNameAsList=[]):
mcfdrDepth = choices.mcfdrDepth if choices.mcfdrDepth else \
AnalysisDefHandler(REPLACE_TEMPLATES['$MCFDR$']).getOptionsAsText().values()[0][0]
analysisDefString = REPLACE_TEMPLATES[
'$MCFDR$'] + ' -> GSuiteSimilarityToQueryTrackRankingsAndPValuesWrapperStat'
analysisSpec = AnalysisDefHandler(analysisDefString)
analysisSpec.setChoice('MCFDR sampling depth', mcfdrDepth)
analysisSpec.addParameter('assumptions', randStrat)
analysisSpec.addParameter('rawStatistic',
GSuiteStatUtils.PAIRWISE_STAT_LABEL_TO_CLASS_MAPPING[similarityStatClassName])
analysisSpec.addParameter('pairwiseStatistic', GSuiteStatUtils.PAIRWISE_STAT_LABEL_TO_CLASS_MAPPING[
similarityStatClassName]) # needed for call of non randomized stat for assertion
analysisSpec.addParameter('tail', 'more')
analysisSpec.addParameter('trackTitles', trackTitles)
analysisSpec.addParameter('queryTracksNum', str(1))
if intensityTrackNameAsList:
analysisSpec.addParameter('trackNameIntensity', '|'.join(intensityTrackNameAsList))
return analysisSpec
def prepareQ3(cls, choices, similarityStatClassName, summaryFunc, randStrat):
mcfdrDepth = choices.mcfdrDepth if choices.mcfdrDepth else \
AnalysisDefHandler(REPLACE_TEMPLATES['$MCFDR$']).getOptionsAsText().values()[0][0]
analysisDefString = REPLACE_TEMPLATES[
'$MCFDRv3$'] + ' -> TrackSimilarityToCollectionHypothesisWrapperStat'
analysisSpec = AnalysisDefHandler(analysisDefString)
analysisSpec.setChoice('MCFDR sampling depth', mcfdrDepth)
analysisSpec.addParameter('assumptions', randStrat)
analysisSpec.addParameter('rawStatistic', 'SummarizedInteractionWithOtherTracksV2Stat')
analysisSpec.addParameter('pairwiseStatistic',
GSuiteStatUtils.PAIRWISE_STAT_LABEL_TO_CLASS_MAPPING[
similarityStatClassName]) # needed for call of non randomized stat for assertion
analysisSpec.addParameter('summaryFunc',
GSuiteStatUtils.SUMMARY_FUNCTIONS_MAPPER[summaryFunc])
analysisSpec.addParameter('tail', 'right-tail')
return analysisSpec