def pos_neg(mat):
pos = []
neg = []
for i in range(len(mat)):
for j in range(i+1, len(mat)):
if (i==0 and j==1) or (i==3 and j==4):
pos.append(mat[i][j])
else: neg.append(mat[i][j])
print "pos", pos
print "neg", neg
datafile = "datasimulationout-proposedmcmcp_exclusive_len50sig15vari1.5noise8.0.out"
m = []
f = open(datafile, 'r')
for line in f:
list = line.replace('\n', '').split(' ')
m.append([int(x) for x in list])
f.close()
poi_pmat_p, poi_pmat_s = poi.poisson_p(m)
print "pppppppppppppppppppppppppppppp"
pos_neg(poi_pmat_p)
print "ssssssssssssssssssssssssssssss"
pos_neg(poi_pmat_s)
#!/usr/bin/env python
import transfer_concretedata as tc
import poiss_binomial as poi
datafile = "datasimulationout-proposedmcmcp_exclusive_len50sig15vari1.5noise8.0.out"
m = []
f = open(datafile, 'r')
for line in f:
list = line.replace('\n', '').split(' ')
m.append([int(x) for x in list])
f.close()
poi_pmat_p, poi_pmat_s = poi.poisson_p(m)
def prplot(fisher, proposed, existing, poisson_p, poisson_s, rn, leng, signal, variance, noise, num, roop_n):
#pid = os.getpid()
print "prplot starts!"
sys.stdout.flush()
fpos_fisher = "pos_"+fisher+"_"+str(rn)+".out"
fneg_fisher = "neg_"+fisher+"_"+str(rn)+".out"
pos_fisher = []
neg_fisher = []
fpos_poi_p = "pos_"+poisson_p+"_"+str(rn)+".out"
fneg_poi_p = "neg_"+poisson_p+"_"+str(rn)+".out"
pos_poi_p = []
neg_poi_p = []
fpos_poi_s = "pos_"+poisson_s+"_"+str(rn)+".out"
fneg_poi_s = "neg_"+poisson_s+"_"+str(rn)+".out"
pos_poi_s = []
neg_poi_s = []
for i in range(5):
print "iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii=%d", i
sys.stdout.flush()
m = ins.normal_dist(leng,leng,signal, 0.0, variance, noise, 1, "b")
datafile = "data"+proposed+".out"
bm = tc.read_and_writedata_sim_withoutfile(m, 0.5, datafile)
ft = fisher_test2_exclusive.fisher_test(m)
poi_pmat_p, poi_pmat_s = poi.poisson_p(m)
#poi_pmat_p = [[0 for i in range(50)] for j in range(50)]
#poi_pmat_s = [[0 for i in range(50)] for j in range(50)]
inproposed1 = proposed+"_num"+str(num*5)
print "file1", inproposed1
sys.stdout.flush()
cmd="sh child_poissonandproposed_exclusive.sh "+str(leng)+" "+str(num)+" "+datafile+" "+inproposed1+""
os.system(cmd)
print "inproposed1"
inproposed2= proposed+"_num"+str(num*5*10)
print "file2", inproposed2
sys.stdout.flush()
cmd="sh child_poissonandproposed_exclusive.sh "+str(leng)+" "+str(num*10)+" "+datafile+" "+inproposed2+""
os.system(cmd)
print "inproposed2"
inexisting1 = existing+"_num"+str(num*5)
print "inexisting1", inexisting1
cmd="sh child_poissonandexisting_exclusive.sh "+str(leng)+" "+str(num)+" "+datafile+" "+inexisting1+""
os.system(cmd)
inexisting2 = existing+"_num"+str(num*5*10)
print "inexisting2", inexisting2
cmd="sh child_poissonandexisting_exclusive.sh "+str(leng)+" "+str(num*10)+" "+datafile+" "+inexisting2+""
os.system(cmd)
"""
tmp = 1
while tmp >0:
time.sleep(10)
tmp = int(commands.getoutput("qstat| grep "+"ch_"+str(pid)+" |wc").split()[0])
"""
for k in range(leng):
for l in range(k+1, leng):
if(k == 0 and l == 1) or (k == 3 and l == 4):
pos_fisher.append(str(ft[k][l])+" ")
pos_poi_p.append(str(poi_pmat_p[k][l])+" ")
pos_poi_s.append(str(poi_pmat_s[k][l])+" ")
else:
neg_fisher.append(str(ft[k][l])+" ")
neg_poi_p.append(str(poi_pmat_p[k][l])+" ")
neg_poi_s.append(str(poi_pmat_s[k][l])+" ")
ffp = open(fpos_fisher, "w")
ffn = open(fneg_fisher, "w")
ffp.writelines(pos_fisher)
ffn.writelines(neg_fisher)
ffp.close()
ffn.close()
ffp = open(fpos_poi_p, "w")
ffn = open(fneg_poi_p, "w")
ffp.writelines(pos_poi_p)
ffn.writelines(neg_poi_p)
ffp.close()
ffn.close()
ffp = open(fpos_poi_s, "w")
ffn = open(fneg_poi_s, "w")
ffp.writelines(pos_poi_s)
ffn.writelines(neg_poi_s)
ffp.close()
ffn.close()
return 0
def prplot(fisher, proposed, existing, poisson_p, poisson_s, rn, leng, signal, variance, noise, num, roop_n):
#pid = os.getpid()
print "prplot starts!"
sys.stdout.flush()
fpos_fisher = "pos_"+fisher+"_"+str(rn)+".out"
fneg_fisher = "neg_"+fisher+"_"+str(rn)+".out"
pos_fisher = []
neg_fisher = []
fpos_poi_p = "pos_"+poisson_p+"_"+str(rn)+".out"
fneg_poi_p = "neg_"+poisson_p+"_"+str(rn)+".out"
pos_poi_p = []
neg_poi_p = []
fpos_poi_s = "pos_"+poisson_s+"_"+str(rn)+".out"
fneg_poi_s = "neg_"+poisson_s+"_"+str(rn)+".out"
pos_poi_s = []
neg_poi_s = []
for i in range(5):
print "iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii=%d", i
sys.stdout.flush()
m = ins.normal_dist(leng,leng,signal, 0.0, variance, noise, 1, "b")
datafile = "data"+proposed+".out"
bm = tc.read_and_writedata_sim_withoutfile(m, 0.5, datafile)
ft = fisher_test2_together.fisher_test(m)
poi_pmat_p, poi_pmat_s = poi.poisson_p(m)
#poi_pmat_p = [[0 for i in range(50)] for j in range(50)]
#poi_pmat_s = [[0 for i in range(50)] for j in range(50)]
inproposed1 = proposed+"_num"+str(num*5)
print "file1", inproposed1
sys.stdout.flush()
cmd="sh child_poissonandproposed_together.sh "+str(leng)+" "+str(num)+" "+datafile+" "+inproposed1+""
os.system(cmd)
print "inproposed1"
inproposed2= proposed+"_num"+str(num*5*10)
print "file2", inproposed2
sys.stdout.flush()
cmd="sh child_poissonandproposed_together.sh "+str(leng)+" "+str(num*10)+" "+datafile+" "+inproposed2+""
os.system(cmd)
print "inproposed2"
inexisting1 = existing+"_num"+str(num*5)
print "inexisting1", inexisting1
cmd="sh child_poissonandexisting_together.sh "+str(leng)+" "+str(num)+" "+datafile+" "+inexisting1+""
os.system(cmd)
inexisting2 = existing+"_num"+str(num*5*10)
print "inexisting2", inexisting2
cmd="sh child_poissonandexisting_together.sh "+str(leng)+" "+str(num*10)+" "+datafile+" "+inexisting2+""
os.system(cmd)
"""
tmp = 1
while tmp >0:
time.sleep(10)
tmp = int(commands.getoutput("qstat| grep "+"ch_"+str(pid)+" |wc").split()[0])
"""
for k in range(leng):
for l in range(k+1, leng):
if(k == 0 and l == 1) or (k == 3 and l == 4):
pos_fisher.append(str(ft[k][l])+" ")
pos_poi_p.append(str(poi_pmat_p[k][l])+" ")
pos_poi_s.append(str(poi_pmat_s[k][l])+" ")
else:
neg_fisher.append(str(ft[k][l])+" ")
neg_poi_p.append(str(poi_pmat_p[k][l])+" ")
neg_poi_s.append(str(poi_pmat_s[k][l])+" ")
ffp = open(fpos_fisher, "w")
ffn = open(fneg_fisher, "w")
ffp.writelines(pos_fisher)
ffn.writelines(neg_fisher)
ffp.close()
ffn.close()
ffp = open(fpos_poi_p, "w")
ffn = open(fneg_poi_p, "w")
ffp.writelines(pos_poi_p)
ffn.writelines(neg_poi_p)
ffp.close()
ffn.close()
ffp = open(fpos_poi_s, "w")
ffn = open(fneg_poi_s, "w")
ffp.writelines(pos_poi_s)
ffn.writelines(neg_poi_s)
ffp.close()
ffn.close()
return 0
