view inbreeding_and_kinship.py @ 31:a631c2f6d913

Update to Miller Lab devshed revision 3c4110ffacc3
author Richard Burhans <burhans@bx.psu.edu>
date Fri, 20 Sep 2013 13:25:27 -0400
parents
children
line wrap: on
line source

#!/usr/bin/env python

import sys
import gd_util

################################################################################

if len(sys.argv) != 6:
    gd_util.die('Usage')

ped_input, ind_input, computed_value, output, kinship_input = sys.argv[1:]

################################################################################

prog = 'inbreed'

args = [ prog ]
args.append(ped_input)      # pedigree
args.append(ind_input)      # specified individuals (e.g.,,potential breeding population)
args.append(kinship_input)  # kinships of founders
args.append(computed_value) # 0 = inbreedng coefficients, 1 = kinships, 2 = mean kinships

with open(output, 'w') as fh:
    gd_util.run_program(prog, args, stdout=fh)

################################################################################

sys.exit(0)