view inbreeding_and_kinship.py @ 39:e56023008e36 default tip

Changed revision of package_fisher_0_1_4 to be2fc454d121 Changed revision of package_matplotlib_1_2 to a03ee94316b5
author miller-lab
date Mon, 06 Jul 2015 10:32:24 -0400
parents a631c2f6d913
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)