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): # F.append(0 )
import mappy File=mappy.open("names-file.txt") name=mappy.read(File) print name