Mercurial > repos > miller-lab > genome_diversity
diff 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 diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/inbreeding_and_kinship.py Fri Sep 20 13:25:27 2013 -0400 @@ -0,0 +1,29 @@ +#!/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) +