f=1-count[1]/(2.*p*q)
if f<0:
continue
alpha_=a+d*(q-p)
alpha[0]=-alpha_*2*p
alpha[1]=alpha_*(q-p)
alpha[2]=alpha_*2*q
dhwe=2*d*f*p*q
delta[0]=-2*p**2 * d + dhwe
delta[1]=2*p*q * d + dhwe
delta[2]=-2*q**2 * d + dhwe
for x in range(0, N):
vcf.set_sample(alpha, delta, iota, g, x, genotype[x] )
namefile=mappy.open(sys.argv[1])
name=mappy.read(namefile).split('\t')[1:]
N=len(name)
SIZE=int(sys.argv[2])
LOCI=int(sys.argv[3])
#sys.stderr.write("size: "+str(SIZE)+'\n')
SIZE=SIZE*32
#sys.stderr.write("size: "+str(SIZE)+'\n')
inb=open(sys.argv[4])
a=float(sys.argv[5])
d=float(sys.argv[6])
e=float(sys.argv[7])
F=[]
#for x in range(0, N):
import mappy
File=mappy.open("names-file.txt")
name=mappy.read(File)
print name