| @@ -0,0 +1,64 @@ | |||||
| #!/usr/bin/env python3 | |||||
| # -*- coding: utf-8 -*- | |||||
| """ | |||||
| Compare gram matrices | |||||
| Created on Mon Sep 24 10:52:25 2018 | |||||
| @author: ljia | |||||
| """ | |||||
| import numpy as np | |||||
| import matplotlib.pyplot as plt | |||||
| N = 7 | |||||
| egmin = [-3.425001366427846e-15, | |||||
| -5.513191435356332e-15, | |||||
| -1.1563146193980238e-15, | |||||
| -1.3999833987183273e-15, | |||||
| -5.811474553224136e-15, | |||||
| 2.196833029054622e-10, | |||||
| 0.0001002381061317695] | |||||
| egmin = np.abs(egmin) | |||||
| egmin2 = [-9.433792343294819e-15, | |||||
| np.NaN, | |||||
| -7.502900269338164e-16, | |||||
| -1.3999833987183273e-15, | |||||
| -8.73626400337456e-15, | |||||
| np.NaN, | |||||
| -4.04460628433013e-14] | |||||
| egmin2 = np.abs(egmin2) | |||||
| egmax = [142.86649135778595, | |||||
| 140.08307372708344, | |||||
| 64.31844814063015, | |||||
| 92.38382991977493, | |||||
| 160.72585558445357, | |||||
| 943.9175660197347, | |||||
| 299.17895175532897] | |||||
| egmax2 = [172.4203026547106, | |||||
| np.NaN, | |||||
| 65.53092059526354, | |||||
| 92.38382991977493, | |||||
| 180.374192331094, | |||||
| np.NaN, | |||||
| 529.3691973508182] | |||||
| fig, ax = plt.subplots() | |||||
| ind = np.arange(N) # the x locations for the groups | |||||
| width = 0.20 # the width of the bars: can also be len(x) sequence | |||||
| p1 = ax.bar(ind, egmin, width) | |||||
| p2 = ax.bar(ind, egmax, width, bottom=egmin) | |||||
| p3 = ax.bar(ind + width, egmin2, width) | |||||
| p4 = ax.bar(ind + width, egmax2, width, bottom=egmin2) | |||||
| ax.set_yscale('log', nonposy='clip') | |||||
| ax.set_xlabel('datasets') | |||||
| ax.set_ylabel('absolute eigen values') | |||||
| ax.set_title('Absolute eigen values of gram matrices on all datasets') | |||||
| plt.xticks(ind + width / 2, ('Acyclic', 'Alkane', 'MAO', 'PAH', 'MUTAG', 'Letter-med', 'ENZYMES')) | |||||
| #ax.set_yticks(np.logspace(-16, -3, num=20, base=10)) | |||||
| ax.set_ylim(bottom=1e-15) | |||||
| ax.legend((p1[0], p2[0], p3[0], p4[0]), ('min1', 'max1', 'min2', 'max2'), loc='upper right') | |||||
| plt.savefig('../check_gm/compare_eigen_values.eps', format='eps', dpi=300) | |||||
| plt.show() | |||||