def main(): loaddir = "../../../data/" probdir = "prob_trans.inputs.sigma2_u1.100.0.sigma2_u2.1.0.sigma2_v.1.0.m.20.nofu1s.5.N.25000" utils.setsizes() params, data = dataload(dir=loaddir+"lms/"+probdir, paramfile="stepsizes.txt") ymax = data[params==0.001] plot(params, data/ymax, label="LMS") #params, data = dataload(dir=loaddir+"bmp/"+probdir, paramfile="metastepsizes.txt") #plot(params, data/ymax, label="K-Y") params, data = dataload(dir=loaddir+"idbd/"+probdir, paramfile="metastepsizes.txt") plot(params, data/ymax, label="IDBD") params, data = dataload(dir=loaddir+"autostep/"+probdir, paramfile="metastepsizes.txt") plot(params, data/ymax, label="Autostep") params, data = dataload(dir=loaddir+"rls/"+probdir, paramfile="one_minus_gammas.txt") plot(params, data/ymax, label="RLS") legend(loc='lower left') xlim([10**-4, 1]) ylim([0, 1.2]) xscale("log") xlabel("tunable parameter", fontsize=utils.getxlabelsize(), labelpad=utils.getxlabelpadsize()) ylabel("$S$ (MSE relative to best LMS)", fontsize=utils.getxlabelsize(), labelpad=utils.getxlabelpadsize())
def main(): loaddir = "../../../data/" probdir = "prob.s2u.1.0.a.0.9998.s2v.1.0.s2o1.100.0.s2o1.1.0.m.20.o1s.5.N.50000" utils.setsizes() params, data = dataload(dir=loaddir+"lms/"+probdir, paramfile="stepsizes.txt") ymax = 1 plot(params, data/ymax, label="LMS") # params, data = dataload(dir=loaddir+"bmp/"+probdir, paramfile="metastepsizes.txt") # plot(params, data/ymax, label="BMP") params, data = dataload(dir=loaddir+"idbd/"+probdir, paramfile="metastepsizes.txt") plot(params, data/ymax, label="IDBD") params, data = dataload(dir=loaddir+"autostep/"+probdir, paramfile="metastepsizes.txt") plot(params, data/ymax, label="Autostep") params, data = dataload(dir=loaddir+"rls/"+probdir, paramfile="one_minus_gammas.txt") plot(params, data/ymax, label="RLS") legend(loc='lower left') xlim([10**-11, 1]) ylim([0, 140]) xscale("log") xlabel("tunable parameter", fontsize=utils.getxlabelsize(), labelpad=utils.getxlabelpadsize()) ylabel("RMSE", fontsize=utils.getxlabelsize(), labelpad=utils.getxlabelpadsize())
def main(): utils.setsizes() mpl.rcParams['lines.markeredgewidth'] = utils.geterrorbarlinewidth() loaddir = "../../../data/" probdir1 = "prob.s2u.1.0.a.0.9998.s2v.1.0.s2o1.0.1.s2o1.0.0.m.20.o1s.5.N.50000" probdir2 = "prob.s2u.1.0.a.0.9998.s2v.1.0.s2o1.100.0.s2o1.0.0.m.20.o1s.5.N.50000" params, data1, data2 = jointdataload(dir1=loaddir+"lms/"+probdir1, dir2=loaddir+"lms/"+probdir2, paramfile="stepsizes.txt") ymax1 = data1[params==0.05] ymax2 = data2[params==0.05] #plot(params, (data1/ymax1 + data2/ymax2)/2.0) params, data1, data2 = jointdataload(dir1=loaddir+"lms/"+probdir1, dir2=loaddir+"lms/"+probdir2, paramfile="stepsizes.txt") plot(params, (data1/ymax1 + data2/ymax2)/2.0, label="LMS") params, data1, data2 = jointdataload(dir1=loaddir+"idbd/"+probdir1, dir2=loaddir+"idbd/"+probdir2, paramfile="metastepsizes.txt") plot(params, (data1/ymax1 + data2/ymax2)/2.0, label="IDBD") params, data1, data2 = jointdataload(dir1=loaddir+"autostep/"+probdir1, dir2=loaddir+"autostep/"+probdir2, paramfile="metastepsizes.txt") plot(params, (data1/ymax1 + data2/ymax2)/2.0, label="Autostep") #params, data1, data2 = jointdataload(dir1=loaddir+"bmp/"+probdir1, dir2=loaddir+"bmp/"+probdir2, paramfile="metastepsizes.txt") #plot(params, (data1/ymax1 + data2/ymax2)/2.0, label="K-Y") params, data1, data2 = jointdataload(dir1=loaddir+"rls/"+probdir1, dir2=loaddir+"rls/"+probdir2, paramfile="one_minus_gammas.txt") plot(params, (data1/ymax1 + data2/ymax2)/2.0, label="RLS") legend(loc='lower left') xlim([10**-9, 1]) ylim([0.6, 1.2]) xscale("log") xlabel("tunable parameter", fontsize=utils.getxlabelsize()) ylabel("$S$ (RMSE relative to best LMS)", fontsize=utils.getxlabelsize(), labelpad=utils.getxlabelpadsize())