#!/usr/bin/python from socal_caller import * from cvxopt import matrix from optparse import OptionParser from random import shuffle from socal_trainer import * import sys import math import time c_aa = matrix([0,0],(2,1)) E_aa = matrix([1,0,0,1],(2,2)) c_ab = matrix([1,1],(2,1)) E_ab = matrix([1,0,0,1],(2,2)) c_bb = matrix([2,2],(2,1)) E_bb = matrix([1,0,0,1],(2,2)) sc = socal_caller(c_aa,E_aa,c_ab,E_ab,c_bb,E_ab,32) xs = [matrix([0,0],(2,1)),matrix([1,1],(2,1)),matrix([2,2],(2,1))] for x in xs: print sc.mindist(x)
c_aa = None E_aa = None c_ab = None E_ab = None c_bb = None E_bb = None if(e_aa != None): c_aa = e_aa['c'] E_aa = e_aa['E'] if(e_ab != None): c_ab = e_ab['c'] E_ab = e_ab['E'] if(e_bb != None): c_bb = e_bb['c'] E_bb = e_bb['E'] sc = socal_caller(c_aa,E_aa,c_ab,E_ab,c_bb,E_bb,c4,e_aa,e_ab,e_bb) # run the caller for i in xrange(len(indv_list)): info = cols[i+3].split(':') geno = info[1] sensea_int = float(info[0].split('/')[0]) senseb_int = float(info[0].split('/')[1]) # do calling x = matrix([sensea_int, senseb_int],(2,1)) result = sc.mindist(x) call = result[0] score = result[1] # write out to file