Mercurial > repos > miller-lab > genome_diversity

comparison inbreeding_and_kinship.py @ 31:a631c2f6d913

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
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
comparison
equal deleted inserted replaced
30:4188853b940b 31:a631c2f6d913
1 #!/usr/bin/env python
2
3 import sys
4 import gd_util
5
6 ################################################################################
7
8 if len(sys.argv) != 6:
9 gd_util.die('Usage')
10
11 ped_input, ind_input, computed_value, output, kinship_input = sys.argv[1:]
12
13 ################################################################################
14
15 prog = 'inbreed'
16
17 args = [ prog ]
18 args.append(ped_input) # pedigree
19 args.append(ind_input) # specified individuals (e.g.,,potential breeding population)
20 args.append(kinship_input) # kinships of founders
21 args.append(computed_value) # 0 = inbreedng coefficients, 1 = kinships, 2 = mean kinships
22
23 with open(output, 'w') as fh:
24 gd_util.run_program(prog, args, stdout=fh)
25
26 ################################################################################
27
28 sys.exit(0)
29