def runAnalysis():
spec = {'results_dir': getHighDataDir() + '/microhomology/mh_freqs_by_len',
'dirname_to_result_fn': lambda x: x,
'result_to_dirname_fn': lambda x: x,
'py_func_load': loadAllMHLenData,
'py_funcs_per_result': [(plotK562PercScatterAnalysis,'RegrLines'), (passData, 'Data')],
'py_funcs_all_results': [compareMHK562lines, plotGCContent],
'reads_colname': 'Non-Null Reads',
'check_output_fn': lambda x: True,
'id_colname': 'Oligo ID',
'min_reads': MIN_READS,
'partitions': ['Non-Targeting'],
'samples': ['K562 New']
}
analyseResultsPerPartition( spec )
spec = {'results_dir': getHighDataDir() + '/microhomology/mh_freqs_by_len',
'dirname_to_result_fn': lambda x: x,
'result_to_dirname_fn': lambda x: x,
'py_func_load': loadAllMHLenData,
'py_funcs_per_result': [(plotPercScatterAnalysis,'RegrLines')],
'py_funcs_all_results': [compareMHlines],
'reads_colname': 'Non-Null Reads',
'check_output_fn': lambda x: True,
'id_colname': 'Oligo ID',
'min_reads': MIN_READS,
'partitions': ['Non-Targeting'],
'samples': ['DPI7']
}
analyseResultsPerPartition( spec )
def runAnalysis():
spec = {'results_dir':getHighDataDir() + '/microhomology_mismatch/mh_mismatch_indel_frequencies',
'dirname_to_result_fn': lambda x: '%s.txt' % x,
'result_to_dirname_fn': lambda x: x.split('/')[-1][:-4],
'py_func_load': loadData,
'py_funcs_per_result': [(passData,'Data')],
'py_funcs_all_results': [plotMicrohomologyMismatches],
'check_output_fn': lambda x: True,
'reads_colname': 'Orig Non-Null Reads',
'min_reads': MIN_READS,
'id_colname': 'Oligo ID',
'partitions': ['Non-Targeting'],
'samples': ['K562 New']
}
analyseResultsPerPartition( spec )
def runAnalysis():
spec = {
'results_dir': getHighDataDir() + '/i1/i1_summaries',
'dirname_to_result_fn': lambda x: '%s.txt' % x,
'result_to_dirname_fn': lambda x: x.split('/')[-1][:-4],
'py_func_load': defaultLoadData,
'py_funcs_per_result': [(mergeWithIndelData, 'i1IndelData')],
'py_funcs_all_results':
[plotMergedPieDataWithAmbig, plotMergedI1Repeats],
'check_output_fn': lambda x: True,
'reads_colname': 'Total reads',
'min_reads': MIN_READS,
'id_colname': 'Oligo Id',
'partitions': ['Non-Targeting'],
'samples': ['K562 New']
}
analyseResultsPerPartition(spec)
def runAnalysis():
spec = {
'results_dir':
getHighDataDir() + '/indel_details/indel_pie_summaries_per_oligo',
'dirname_to_result_fn':
lambda x: '%s.txt' % x,
'result_to_dirname_fn':
lambda x: x.split('/')[-1][:-4],
'py_func_load':
loadData,
'py_funcs_per_result': [(perOligoCounts, 'perOligoCounts'),
(perOligoMCI, 'perOligoMCI'),
(computePercentages, 'PercData')],
'py_funcs_all_results': [plotSumPie, plotMCIPie, plotPercCorrelations],
'check_output_fn':
lambda x: True,
'reads_colname':
'Total reads',
'min_reads':
MIN_READS,
'id_colname':
'Oligo Id',
'partitions': ['Real Guides'],
'samples': ['K562 New']
}
analyseResultsPerPartition(spec)
spec = {
'results_dir':
getHighDataDir() + '/indel_details/indel_pie_summaries_per_oligo',
'dirname_to_result_fn': lambda x: '%s.txt' % x,
'result_to_dirname_fn': lambda x: x.split('/')[-1][:-4],
'py_func_load': loadData,
'py_funcs_per_result': [(computePieData, 'PieData')],
'py_funcs_all_results': [plotBarSummaryPieIndels],
'check_output_fn': lambda x: True,
'reads_colname': 'Total reads',
'min_reads': MIN_READS,
'id_colname': 'Oligo Id',
'partitions': ['Real Guides'],
'samples': ['DPI7']
}
analyseResultsPerPartition(spec)