def test_exon_list_sort_on_log2_fold_change(self): 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)) exon_list.sort('log2_fold_change') last_log2_fold_change = None for exon in exon_list.exons: if last_log2_fold_change is not None: self.assertTrue(last_log2_fold_change >= exon.log2_fold_change) last_log2_fold_change = exon.log2_fold_change
def test_exon_list_sort_on_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)) exon_list.sort('p_value') last_p_value = None for exon in exon_list.exons: if last_p_value is not None: self.assertTrue(last_p_value <= exon.p_value) last_p_value = exon.p_value
def test_exon_list_sort_on_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)) exon_list.sort('p_value') last_p_value = None for exon in exon_list.exons: if last_p_value is not None: self.assertTrue(last_p_value <= exon.p_value) last_p_value = exon.p_value
def test_exon_list_sort_on_log2_fold_change(self): 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)) exon_list.sort('log2_fold_change') last_log2_fold_change = None for exon in exon_list.exons: if last_log2_fold_change is not None: self.assertTrue(last_log2_fold_change >= exon.log2_fold_change) last_log2_fold_change = exon.log2_fold_change