def compareMeasureNG(edges,sims,ratio,names,wName=None,ciss=None,link_func=maxSim):
'''
NG means this experiment is not using the grouping
of cis-element
'''
precs=dict.fromkeys(names)
recalls=dict.fromkeys(names)
ss=dict.fromkeys(names)
ee=flatenDictList(edges)#ee is true edges
e=mdf.global_remove(edges,ratio)
for name in names:
ss[name]=flatenDict(netAlignInfer\
(link_func,e,sims[name]))
for name in names:
precs[name]=[]
recalls[name]=[]
for thresh in np.arange(0,1.0,0.1):
(p,r)=prEval(ee,ss[name],thresh)
precs[name].append(p)
recalls[name].append(r)
if not wName==None:
(precs[wName],recalls[wName])=doEntropyWeights(ee,e,ciss,\
link_func=link_func)
return (precs,recalls,ss)
def compareMeasureNG2(edges,sims,ratio,names,proms,wName=None,ciss=None,link_func=maxSim):
precs=dict.fromkeys(names)
recalls=dict.fromkeys(names)
ss=dict.fromkeys(names)
ee=flatenDictList(edges)#ee is true edges
e=mdf.global_remove(edges,ratio)
for name in names:
ss[name]=flatenDict(netAlignInfer\
(link_func,e,sims[name]))
tss=getTS2(edges,ciss,proms,ss[name],0.1)
factors=getAssoc2(*tss,min_c=0.9)
for name in names:
precs[name]=[]
recalls[name]=[]
for thresh in np.arange(0,1.0,0.1):
for key in ss[name].keys():
if key[1] in factors.keys():
if key[0]+str(ciss[key[1]].slen) == tss[1][key[1]]:
#print 'BEFORE',ss[name][key]
ss[name][key]=1
#print 'AFTER',ss[name][key]
(p,r)=prEval(ee,ss[name],thresh)
precs[name].append(p)
recalls[name].append(r)
if not wName==None:
(precs[wName],recalls[wName])=doEntropyWeights(ee,e,ciss,\
link_func=link_func)
return (precs,recalls,ss)
def doEntropyWeights(ee,e,ciss,link_func=maxSim,proms=None):
precs=[]
recalls=[]
weights=getEntropyWeights(e,ciss)
cis_keys=ciss.keys()
sims=dict()
for x in cis_keys:
sims[x]=dict()
for y in cis_keys:
sims[x][y]=dmSim(ciss[x].seq,ciss[y].seq\
,w1=weights[x],w2=weights[y]\
,useLength=True)
ss=flatenDict(netAlignInfer(link_func,e,sims))
for thresh in np.arange(0,1,0.1):
(p,r)=prEval(ee,ss,thresh)
precs.append(p)
recalls.append(r)
return (precs,recalls)
# Repeats of experiments
N=5
ratio=0.3
x=[]
y=[]
x2=[]
y2=[]
x3=[]
y3=[]
ee=flatenDictList(edges)
for exp_i in range(N):
e=mdf.global_remove(edges,0.3)
scores1=netAlignInfer(maxSim,e,nwSims)
scores2=netAlignInfer(maxSim,e,nwSimsL)
scores3=netAlignInfer(maxSim,e,dmSimsL)
ss1=flatenDict(scores1)
ss2=flatenDict(scores2)
ss3=flatenDict(scores3)
for thresh in np.arange(0,1,0.1):
result=prEval(ee,ss1,thresh)
result2=prEval(ee,ss2,thresh)
result3=prEval(ee,ss3,thresh)
x.append(result[0])
y.append(result[1])
x2.append(result2[0])
y2.append(result2[1])
x3.append(result3[0])
y3.append(result3[1])
plt.figure(1)