def generateMccTotalEner4NormedData(h5_path, dt_path, weight_fn):
from myh5 import myh5
f = myh5(h5_path)
dt = f.loadDt(dt_path)
weight = loadWeights(weight_fn)
mcc = dt[:, 0]
ener = dt[:, 1:]
total_ener = np.dot(ener, weight)
mcc_total = np.column_stack((mcc, total_ener))
return mcc_total
def sampleMcc():
from myh5 import myh5
h5_path = "/ddnB/work/jaydy/dat/output/linr_out/08_lnr.h5"
dt_path = "weighted_high"
f = myh5(h5_path)
high_mcc_total = f.loadDt(dt_path)
dt_path = "weighted_low"
low_mcc_total = f.loadDt(dt_path)
sample_sz = 2000
high_mcc_total = sampleMccTotalByMcc(high_mcc_total, sample_sz)
low_mcc_total = sampleMccTotalByMcc(low_mcc_total, sample_sz)
return high_mcc_total, low_mcc_total
