File size: 1,616 Bytes
db2fca9 | 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 | """final2.py seeds(comma) name : L8 u12 650 steps, variants sum3 / sum3b. saves out/graph2{v}_{name}_c{23,25,26}.npy (float32, v21seq val order), eval arrays, prints table"""
import sys,glob,numpy as np; sys.dont_write_bytecode=True
from evalrho import *; from ana import *; import ana
sys.path.insert(0,'/workspace/huawei/v21seq')
from evalseq import tologit,yv
seeds=[int(s) for s in sys.argv[1].split(',')]; name=sys.argv[2]; WSX=(.05,.1,.15,.2,.25,.3)
V={'s3':(1,1,1,0),'s3b':(2,0,1,0)}
y26=yv(26); fp=sorted(glob.glob('/workspace/huawei/v21seq/ens/ctl6_A_c26_s*.npy')); M26=nr(np.mean([tologit(np.load(f)[:,-1]) for f in fp],0)); a26=AUC(y26,M26)
Z={c:(lambda A:np.stack([zs(A[:,j]) for j in range(4)],1))(avg('u12',c,seeds,650)) for c in (23,25,26,27)}
print(name,'seeds',seeds)
for v,w in V.items():
w=np.array(w,float)
for c in (23,25,26):
G=Z[c]@w; np.save(O+f'graph2{v}_{name}_c{c}.npy',G.astype(np.float32))
if c==26:
Gn=nr(G); print(f'C26 {v} auc {AUC(y26,G):.5f} rhoM {np.corrcoef(Gn,M26)[0,1]:.3f} vsM '+' '.join('w%.2f %+.1f'%(x,(AUC(y26,(1-x)*M26+x*Gn)-a26)*1e5) for x in WSX),flush=True)
else:
a,r,o=best(c,G,WSX); print(f'C{c} {v} auc {a:.5f} rhoM {r:.3f} '+' | '.join(f'{k}: best {q[0]:+.1f}@{q[1]} {q[2]}' for k,q in o.items()),flush=True)
G=Z[27]@w; np.save(O+f'graph2{v}_{name}_eval_v21seqorder.npy',G.astype(np.float32)); gm=to_merged(G); np.save(O+f'graph2{v}_{name}_eval_merged_rank.npy',nr(gm))
print(f'eval {v}: rho EN {rho(gm):.4f} rho G7 {rho(gm,nr(G7)):.4f}',flush=True)
np.save(O+'graph2_eval_rowid.npy',np.asarray(rid))
|