|
@@ -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() |