def test_span_gaps_nomerge_neg(self):
p = PlotLog('x', 20, 2, 0.05, 0.05, 1, 4, 1, True, None)
l = p.span_block_gaps([
0, 0, 0, -9, -9, -9, 0, 0, 0, 0, 0, -8, -8, -8, -8, 0, 0, 0, 0, 0
])
self.assertEqual([
0, 0, 0, -9, -9, -9, 0, 0, 0, 0, 0, -8, -8, -8, -8, 0, 0, 0, 0, 0
], l)
def test_span_gaps_block_at_end_neg(self):
p = PlotLog('x', 20, 2, 0.05, 0.05, 1, 4, 1, True, None)
l = p.span_block_gaps([
0, 0, 0, -9, -9, -9, 0, 0, 0, 0, 0, -8, -8, -8, -8, 0, 0, -7, -7,
-7
])
self.assertEqual([
0, 0, 0, -9, -9, -9, 0, 0, 0, 0, 0, -8, -8, -8, -8, -2, -2, -7, -7,
-7
], l)
def generate_logfc_plot(self, analysis_type, essentiality_files):
files = [
getattr(essentiality_files[plotfile],
analysis_type).tradis_essentiality_filename
for plotfile in self.plotfiles
]
only_ess_files = [
getattr(essentiality_files[plotfile],
analysis_type).only_essential_filename
for plotfile in self.plotfiles
]
annotation_files = [
essentiality_files[plotfile].embl_filename
for plotfile in self.plotfiles
]
mid = int(len(files) / 2)
t = TradisComparison(files[:mid], files[mid:], self.verbose,
self.minimum_block, only_ess_files[:mid],
only_ess_files[mid:])
t.run()
p = PlotLog(t.output_filename, self.genome_length, self.minimum_logfc,
self.pvalue, self.qvalue, self.minimum_logcpm,
self.window_size, self.span_gaps,
self.report_decreased_insertions, annotation_files[0])
p.construct_plot_file()
renamed_csv_file = os.path.join(self.prefix, analysis_type + ".csv")
renamed_plot_file = os.path.join(self.prefix, analysis_type + ".plot")
shutil.copy(t.output_filename, renamed_csv_file)
shutil.copy(p.output_filename, renamed_plot_file)
os.remove(t.output_filename)
os.remove(p.output_filename)
if self.verbose:
print("Comprison:\t" + renamed_csv_file)
print("Plot log:\t" + renamed_plot_file)
return renamed_plot_file
def test_span_gaps_block_at_end(self):
p = PlotLog('x', 20, 2, 0.05, 0.05, 1, 4, 1, True, None)
l = p.span_block_gaps(
[0, 0, 0, 9, 9, 9, 0, 0, 0, 0, 0, 8, 8, 8, 8, 0, 0, 7, 7, 7])
self.assertEqual(
[0, 0, 0, 9, 9, 9, 0, 0, 0, 0, 0, 8, 8, 8, 8, 2, 2, 7, 7, 7], l)
def test_span_gaps_merge(self):
p = PlotLog('x', 20, 2, 0.05, 0.05, 1, 4, 1, True, None)
l = p.span_block_gaps(
[0, 0, 0, 9, 9, 9, 9, 0, 0, 0, 0, 8, 8, 8, 8, 0, 0, 0, 0, 0])
self.assertEqual(
[0, 0, 0, 9, 9, 9, 9, 2, 2, 2, 2, 8, 8, 8, 8, 0, 0, 0, 0, 0], l)