| import sys,glob,numpy as np,polars as pl; sys.dont_write_bytecode=True |
| from final import * |
| sys.path.insert(0,'/workspace/huawei/v21seq') |
| from evalseq import tologit,yv |
| G26=np.load(O+'graph_c26.npy') if False else None |
| for c in (26,27): |
| A=score(c); G=zs(A['s_uv'])+zs(A['s_h'])+zs(A['cos_mean']); np.save(O+f'graph_c{c}.npy',G.astype(np.float32)); print(c,G.shape,np.isfinite(G).all()) |
| G=np.load(O+'graph_c26.npy').astype(np.float64); y=yv(26) |
| fp=sorted(glob.glob('/workspace/huawei/v21seq/ens/ctl6_A_c26_s*.npy')); M=nr(np.mean([tologit(np.load(f)[:,-1]) for f in fp],0)); Gr=nr(G) |
| print('C26 y rows',len(y),len(G),'graph auc %.5f M auc %.5f rho %.3f'%(AUC(y,G),AUC(y,M),np.corrcoef(Gr,M)[0,1])) |
| a0=AUC(y,M); print('C26 gain over M (rank): '+' '.join(f'w{w}:{(AUC(y,(1-w)*M+w*Gr)-a0)*1e5:+.0f}' for w in (.05,.1,.125,.15,.2))) |
| E=np.load(O+'graph_c27.npy'); np.save(O+'graph_eval.npy',E.astype(np.float32)) |
| T=np.load(O+'tab_eval.npz'); np.save(O+'graph_eval_rowid.npy',T['rowid']); print('eval',E.shape,T['rowid'][:3]) |
|
|
