INNER CODE UNIT · Python

read_king

brentp/peddy · peddy/king.py:44

def read_king(king_file):
    pairs = {}
    for d in reader(king_file):
        pairs[(d['ID1'], d['ID2'])] = float(d['Kinship'])
        pairs[(d['ID2'], d['ID1'])] = float(d['Kinship'])
    return pairs


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)

View source record →