def plot():
high_mcc_total, low_mcc_total = sampleMcc()
import plot
x1 = high_mcc_total[:, 0]
y1 = high_mcc_total[:, 1]
x2 = low_mcc_total[:, 0]
y2 = low_mcc_total[:, 1]
plot.two_scatter(x1, y1, x2, y2, ofn="lnr_high_low_scatter.pdf")
def plotScatter():
import plot
import plot_mcc_ener
high_mcc_diff = h5py.File('/ddnB/work/jaydy/dat/output/linr_out/08ff_all_decoy.h5')['/gaussian_nb/high/high_mcc_diff'][()]
low_mcc_diff = h5py.File('/ddnB/work/jaydy/dat/output/linr_out/08ff_all_decoy.h5')['/gaussian_nb/low/low_mcc_diff'][()]
high_mcc_diff = plot_mcc_ener.sampleMccTotalByMcc(high_mcc_diff, 2000)
low_mcc_diff = plot_mcc_ener.sampleMccTotalByMcc(low_mcc_diff, 2000)
x1 = high_mcc_diff[:, 0]
y1 = high_mcc_diff[:, 1]
x2 = low_mcc_diff[:, 0]
y2 = low_mcc_diff[:, 1]
plot.two_scatter(x1, y1, x2, y2, ofn='08_gaussian_mcc_diff_scatter.pdf')
