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