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 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 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 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)
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)
