def test_get_best_exon_from_log2_fold_change(self):
# First example
exon_list = ExonList("A1BG")
exon_list.add_exon(Exon("PSR19025918.hg.1","A1BG",log2_fold_change=-0.056323333))
exon_list.add_exon(Exon("PSR19025921.hg.1","A1BG",log2_fold_change=0.075113333))
exon_list.add_exon(Exon("PSR19025922.hg.1","A1BG",log2_fold_change=0.037316667))
exon_list.add_exon(Exon("PSR19025925.hg.1","A1BG",log2_fold_change=-0.10211))
exon_list.add_exon(Exon("PSR19025926.hg.1","A1BG",log2_fold_change=0.059556667))
exon_list.add_exon(Exon("PSR19025927.hg.1","A1BG",log2_fold_change=0.119746667))
exon_list.add_exon(Exon("PSR19025928.hg.1","A1BG",log2_fold_change=-0.090296667))
exon_list.add_exon(Exon("PSR19025929.hg.1","A1BG",log2_fold_change=-0.02433))
exon_list.add_exon(Exon("PSR19025930.hg.1","A1BG",log2_fold_change=0.02569))
best_exon = exon_list.best_exon('log2_fold_change')
self.assertEqual(best_exon.name,"PSR19025927.hg.1")
# Second example
exon_list = ExonList("A1BG")
exon_list.add_exon(Exon("PSR19025918.hg.1","A1BG",log2_fold_change=-0.056323333))
exon_list.add_exon(Exon("PSR19025921.hg.1","A1BG",log2_fold_change=0.075113333))
exon_list.add_exon(Exon("PSR19025922.hg.1","A1BG",log2_fold_change=0.037316667))
exon_list.add_exon(Exon("PSR19025925.hg.1","A1BG",log2_fold_change=-0.10211))
exon_list.add_exon(Exon("PSR19025926.hg.1","A1BG",log2_fold_change=0.059556667))
exon_list.add_exon(Exon("PSR19025928.hg.1","A1BG",log2_fold_change=-0.090296667))
exon_list.add_exon(Exon("PSR19025929.hg.1","A1BG",log2_fold_change=-0.02433))
exon_list.add_exon(Exon("PSR19025930.hg.1","A1BG",log2_fold_change=0.02569))
best_exon = exon_list.best_exon('log2_fold_change')
self.assertEqual(best_exon.name,"PSR19025925.hg.1")
def test_get_best_exon_from_log2_fold_change(self):
# First example
exon_list = ExonList("A1BG")
exon_list.add_exon(
Exon("PSR19025918.hg.1", "A1BG", log2_fold_change=-0.056323333))
exon_list.add_exon(
Exon("PSR19025921.hg.1", "A1BG", log2_fold_change=0.075113333))
exon_list.add_exon(
Exon("PSR19025922.hg.1", "A1BG", log2_fold_change=0.037316667))
exon_list.add_exon(
Exon("PSR19025925.hg.1", "A1BG", log2_fold_change=-0.10211))
exon_list.add_exon(
Exon("PSR19025926.hg.1", "A1BG", log2_fold_change=0.059556667))
exon_list.add_exon(
Exon("PSR19025927.hg.1", "A1BG", log2_fold_change=0.119746667))
exon_list.add_exon(
Exon("PSR19025928.hg.1", "A1BG", log2_fold_change=-0.090296667))
exon_list.add_exon(
Exon("PSR19025929.hg.1", "A1BG", log2_fold_change=-0.02433))
exon_list.add_exon(
Exon("PSR19025930.hg.1", "A1BG", log2_fold_change=0.02569))
best_exon = exon_list.best_exon('log2_fold_change')
self.assertEqual(best_exon.name, "PSR19025927.hg.1")
# Second example
exon_list = ExonList("A1BG")
exon_list.add_exon(
Exon("PSR19025918.hg.1", "A1BG", log2_fold_change=-0.056323333))
exon_list.add_exon(
Exon("PSR19025921.hg.1", "A1BG", log2_fold_change=0.075113333))
exon_list.add_exon(
Exon("PSR19025922.hg.1", "A1BG", log2_fold_change=0.037316667))
exon_list.add_exon(
Exon("PSR19025925.hg.1", "A1BG", log2_fold_change=-0.10211))
exon_list.add_exon(
Exon("PSR19025926.hg.1", "A1BG", log2_fold_change=0.059556667))
exon_list.add_exon(
Exon("PSR19025928.hg.1", "A1BG", log2_fold_change=-0.090296667))
exon_list.add_exon(
Exon("PSR19025929.hg.1", "A1BG", log2_fold_change=-0.02433))
exon_list.add_exon(
Exon("PSR19025930.hg.1", "A1BG", log2_fold_change=0.02569))
best_exon = exon_list.best_exon('log2_fold_change')
self.assertEqual(best_exon.name, "PSR19025925.hg.1")
def test_get_best_exon_from_p_value(self):
exon_list = ExonList("A1BG")
exon_list.add_exon(Exon("PSR19025918.hg.1", "A1BG", p_value=0.5347865))
exon_list.add_exon(Exon("PSR19025921.hg.1", "A1BG", p_value=0.5820691))
exon_list.add_exon(Exon("PSR19025922.hg.1", "A1BG", p_value=0.7582407))
exon_list.add_exon(Exon("PSR19025925.hg.1", "A1BG", p_value=0.4111732))
exon_list.add_exon(Exon("PSR19025926.hg.1", "A1BG", p_value=0.6550312))
exon_list.add_exon(Exon("PSR19025927.hg.1", "A1BG", p_value=0.5002532))
exon_list.add_exon(Exon("PSR19025928.hg.1", "A1BG", p_value=0.3521274))
exon_list.add_exon(Exon("PSR19025929.hg.1", "A1BG", p_value=0.7716908))
exon_list.add_exon(Exon("PSR19025930.hg.1", "A1BG", p_value=0.7720515))
best_exon = exon_list.best_exon('p_value')
self.assertEqual(best_exon.name, "PSR19025928.hg.1")
def test_get_best_exon_from_p_value(self):
exon_list = ExonList("A1BG")
exon_list.add_exon(Exon("PSR19025918.hg.1","A1BG",p_value=0.5347865))
exon_list.add_exon(Exon("PSR19025921.hg.1","A1BG",p_value=0.5820691))
exon_list.add_exon(Exon("PSR19025922.hg.1","A1BG",p_value=0.7582407))
exon_list.add_exon(Exon("PSR19025925.hg.1","A1BG",p_value=0.4111732))
exon_list.add_exon(Exon("PSR19025926.hg.1","A1BG",p_value=0.6550312))
exon_list.add_exon(Exon("PSR19025927.hg.1","A1BG",p_value=0.5002532))
exon_list.add_exon(Exon("PSR19025928.hg.1","A1BG",p_value=0.3521274))
exon_list.add_exon(Exon("PSR19025929.hg.1","A1BG",p_value=0.7716908))
exon_list.add_exon(Exon("PSR19025930.hg.1","A1BG",p_value=0.7720515))
best_exon = exon_list.best_exon('p_value')
self.assertEqual(best_exon.name,"PSR19025928.hg.1")
