def testOneMissingPath(self): db1 = {} db2 = {} clusters = [] exons1 = [self.g1] exons2 = [self.a2, self.b2, self.c2, self.d2, self.e2, self.f2] comp.addExon(db1, clusters, exons1) comp.addExon(db2, clusters, exons2) missingPaths, alteredPaths = comp.findPathDiff(db1, db2) self.assertEqual(len(missingPaths), 1) self.assertEqual(len(alteredPaths), 0) missingPaths, alteredPaths = comp.findPathDiff(db2, db1) self.assertEqual(len(missingPaths), 1) self.assertEqual(len(alteredPaths), 0)
def testTwoPathDiff(self): db1 = {} db2 = {} clusters = [] exons1 = [self.a1, self.b1, self.c1, self.d1, self.f1] exons2 = [self.a2, self.b2, self.c2, self.d2, self.e2, self.f2] exons3 = [self.a3, self.b3, self.c3, self.d3, self.e3, self.g3] comp.addExon(db1, clusters, exons1) comp.addExon(db2, clusters, exons2) comp.addExon(db2, clusters, exons3) missingPaths, alteredPaths = comp.findPathDiff(db1, db2) self.assertEqual(len(missingPaths), 0) self.assertEqual(len(alteredPaths), 1) missingPaths, alteredPaths = comp.findPathDiff(db2, db1) self.assertEqual(len(missingPaths), 0) self.assertEqual(len(alteredPaths), 2)
def testOnePathDiff(self): db1 = {} db2 = {} clusters = [] exons1 = [self.a1, self.c1, self.d1, self.e1, self.f1] exons2 = [self.a2, self.c2, self.d2, self.e2, self.f2, self.g2] comp.addExon(db1, clusters, exons1) comp.addExon(db2, clusters, exons2) self.assertEqual(len(db1), 5) self.assertEqual(len(db2), 6) self.assertEqual(len(db1[self.a1.coord].cluster.nodes()), 5) self.assertEqual(len(db2[self.a1.coord].cluster.nodes()), 6) missingPaths, alteredPaths = comp.findPathDiff(db2, db1) self.assertEqual(len(missingPaths), 0) self.assertEqual(len(alteredPaths), 1)
def testLeftRightExtension(self): db1 = {} db2 = {} clusters = [] exons1 = [self.aa2, self.c1, self.d1, self.e1, self.f1, self.gg1] exons2 = [self.aa1, self.c2, self.d2, self.e2, self.f2, self.gg2] comp.addExon(db1, clusters, exons1) comp.addExon(db2, clusters, exons2) missingPaths, alteredPaths = comp.findPathDiff(db1, db2) self.assertEqual(len(missingPaths), 0) self.assertEqual(len(alteredPaths), 1) cluster1 = self.gg1.cluster cluster2 = self.gg2.cluster leftExtension, rightExtension = \ comp.findExtendedEnd(cluster1, cluster2, db1, db2) self.assertEqual(len(leftExtension), 1) self.assertEqual(len(rightExtension), 1)