multreelib.py 文件源码

def get_conflicts(self):
        """Find irreconcilable connected components of leg."""
        conflicts = set()  # connected components with conflict
        for cc in nx.connected_components(self.leg):
            # key = species, val = set of loci in species for this cc
            loci_dct = collections.defaultdict(set)

            for label in cc:
                species, locus = label
                loci_dct[species].add(locus)

            for sp, loci in loci_dct.iteritems():
                # conflict if a species has more than one loci in this cc
                if len(loci) >= 2:
                    conflicts.add(tuple(cc))
                    break
        return conflicts