INNER CODE UNIT · Python

kingped

brentp/peddy · peddy/king.py:52

def kingped(ped_obj, king_pairs, sib_pairs, parent_child, cutoff=0.13):

    seen = set()
    print("sample_a\tsample_b\terror\tkinship")
    fmt = "%s\t%s\t%s\t%.3f"

    for s in sib_pairs:
        try:
            king_pairs[s]
            seen.add(s)
        except KeyError:
            continue

        v = king_pairs[s]
        if v < cutoff:
            print(fmt % (s[0], s[1], "sibs low kingship", v))

    for kp in parent_child:

View source record →