some minor fixes after the publication

Author: vlewir

backbone/backbone_colored_circle_just_selected.py

@@ -6,15 +6,15 @@
 
 #usage: python backbone_colored_circle_just_selected.py
 
-sorting_index = np.loadtxt('../MaxAbs_edds_sorting_index_ltof.txt', dtype=int)
+sorting_index = np.loadtxt('../MaxAbs_ewds_sorting_index_ltof.txt', dtype=int)
 
 chList = np.loadtxt('../ChList.txt', dtype=str)
 
 chList_color = np.loadtxt('../ChList_colored.txt', dtype=str)
 
 source, target = np.loadtxt('../edgelist.txt', unpack=True, dtype=int)
 
-data = np.loadtxt('../MaxAbs_edds_just_3.txt', skiprows=3)
+data = np.loadtxt('../MaxAbs_ewds_just_3.txt', skiprows=3)
 
 data = data[:3570,:]
 

backbone/plot_edge_colors_selected.py

@@ -7,11 +7,11 @@
 
 edge_colors = np.loadtxt('edge_colors_selected.txt')
 
-data = np.loadtxt('../MaxAbs_edds_just_3.txt', skiprows=3)
+data = np.loadtxt('../MaxAbs_ewds_just_3.txt', skiprows=3)
 
 data = data[:3570,:]
 
-sorting_index = np.loadtxt('../MaxAbs_edds_sorting_index_ltof.txt', dtype=int)
+sorting_index = np.loadtxt('../MaxAbs_ewds_sorting_index_ltof.txt', dtype=int)
 
 fig, axes = plt.subplots(nrows=1, ncols=2, figsize=(30,40))
 

backbone/plot_edge_colors_selected_both.py

@@ -7,11 +7,11 @@
 
 edge_colors = np.loadtxt('edge_colors_selected.txt')
 
-data = np.loadtxt('../MaxAbs_edds_just_3.txt', skiprows=3)
+data = np.loadtxt('../MaxAbs_ewds_just_3.txt', skiprows=3)
 
 data = data[:3570,:]
 
-sorting_index = np.loadtxt('../MaxAbs_edds_sorting_index_ltof.txt', dtype=int)
+sorting_index = np.loadtxt('../MaxAbs_ewds_sorting_index_ltof.txt', dtype=int)
 
 new_data = []
 new_edge_colors = []

cliffs/ndd_by_nodes_Cliffs_individual_together.py

@@ -1,11 +1,7 @@
 import numpy as np
-import matplotlib.pyplot as plt
-from matplotlib import cm
-import matplotlib
 import sys
-import os
 
-#usage: python ndd_by_nodes_Cliffs_individual_together.py MaxAbs ChList.txt MaxAbs_ndds_ln_by_each_raw_data.txt
+#usage: python ndd_by_nodes_Cliffs_individual_together.py MaxAbs
 
 def ED(A,B):
 
@@ -28,32 +24,9 @@ def ED(A,B):
 
 name = sys.argv[1]
 
-cl = sys.argv[2]
-
-bins = sys.argv[3]
-
-ticks = np.loadtxt(cl, dtype=str)
-
-lines = [line.rstrip('\n') for line in open(bins)]
-
-bin_min, bin_max = lines[0].split(" ")
-
-bin_min = np.round(float(bin_min),2)
-bin_max = np.round(float(bin_max),2)
-
-fontsize = 36
-
-matplotlib.rcParams.update({'font.size': fontsize})
-
 names = ["Control", "EtOH"]
-tnames = {}
-tnames["Control"] = "G1 - control"
-tnames["EtOH"] = "G1 - EtOH"
 
 #names = ["TP1", "TP2"]
-#tnames = {}
-#tnames["TP1"] = "G2 - tp1"
-#tnames["TP2"] = "G2 - tp2"
 
 # -----
 
@@ -66,7 +39,6 @@ def ED(A,B):
 
 		if (subplot_index == 0):
 			A = data
-			plt.subplots_adjust(right=0.9)
 		else:
 			B = data
 

cliffs/newd_by_nodes_Cliffs_individual_together.py

@@ -1,11 +1,7 @@
 import numpy as np
-import matplotlib.pyplot as plt
-from matplotlib import cm
-import matplotlib
 import sys
-import os
 
-#usage: python newd_by_nodes_Cliffs_individual_together.py MaxAbs ChList.txt MaxAbs_newds_by_each_raw_data.txt
+#usage: python newd_by_nodes_Cliffs_individual_together.py MaxAbs
 
 def ED(A,B):
 
@@ -52,33 +48,12 @@ def ED(A,B):
 
 name = sys.argv[1]
 
-cl = sys.argv[2]
-
-bins = sys.argv[3]
-
-ticks = np.loadtxt(cl, dtype=str)
-
-lines = [line.rstrip('\n') for line in open(bins)]
-
-bin_min, bin_max = lines[0].split(" ")
-
-bin_min = np.round(float(bin_min),2)
-bin_max = np.round(float(bin_max),2)
-
-fontsize = 36
-
-matplotlib.rcParams.update({'font.size': fontsize})
 
 names = ["Control", "EtOH"]
-tnames = {}
-tnames["Control"] = "G1 - control"
-tnames["EtOH"] = "G1 - EtOH"
 
 # names = ["TP1", "TP2"]
-# tnames = {}
-# tnames["TP1"] = "G2 - tp1"
-# tnames["TP2"] = "G2 - tp2"
 
+# -----
 
 for i in range(1,19):
 
@@ -89,7 +64,6 @@ def ED(A,B):
 
 		if (subplot_index == 0):
 			A = data
-			plt.subplots_adjust(right=0.9)
 		else:
 			B = data
 

plot/ewd_by_nodes_sorting_just_3_fig.py

@@ -5,7 +5,7 @@
 import sys
 import os
 
-#usage: python ewds_by_nodes_sorting_just_3_fig.py MaxAbs 30 MaxAbs_edds_just_3.txt MaxAbs_edds_sorting_index_ltof.txt
+#usage: python ewds_by_nodes_sorting_just_3_fig.py MaxAbs 30 MaxAbs_ewds_just_3.txt MaxAbs_ewds_sorting_index_ltof.txt
 
 name = sys.argv[1]
 

plot/ndd_by_nodes_sorting_just_2_Cliffs_average_FIG_sml.py

@@ -17,11 +17,11 @@
 
 sorting = sys.argv[5]
 
-ED_average = sys.argv[6]
+Cliffs_average = sys.argv[6]
 
-ED_values = np.loadtxt(ED_average)
+Cliffs_values = np.loadtxt(Cliffs_average)
 
-ED_values *= -1.0 # switch order
+Cliffs_values *= -1.0 # switch order
 
 data = np.loadtxt(filename, skiprows=4)
 
@@ -63,33 +63,13 @@
 
 # names = ["G2 - tp1\n", "G2 - tp2\n"]
 
-# names = ["G1 - control\n", "G2 - tp1\n"]
-
-# names = ["G2 - tp2\n", "G1 - EtOH\n"]
-
-# names = ["G1 - control\n", "G2 - tp2\n"]
-
-# names = ["G1 - EtOH\n", "G1 - control\n"]
-
-# names = ["G1 - EtOH\n", "G2 - tp1\n"]
-
-# names = ["G1 - EtOH\n", "G2 - tp2\n"]
-
-# names = ["G2 - tp1\n", "G1 - control\n"]
-
-# names = ["G2 - tp1\n", "G1 - EtOH\n"]
-
-# names = ["G2 - tp2\n", "G1 - control\n"]
-
-# names = ["G2 - tp2\n", "G2 - tp1\n"]
-
 ax = plt.subplot(144)
 
 ax.set_position([0.89, 0.0925, 0.04, 0.839])
 
 plt.title("$\delta$\n")
 
-tmp = ED_values.copy()
+tmp = Cliffs_values.copy()
 
 tmp = np.abs(tmp)
 
@@ -123,11 +103,11 @@
 
 plt.title("$\delta$\n")
 
-ED_values = np.expand_dims(ED_values, axis=1)
+Cliffs_values = np.expand_dims(Cliffs_values, axis=1)
 
-print(np.amin(ED_values), np.amax(ED_values))
+print(np.amin(Cliffs_values), np.amax(Cliffs_values))
 
-im = plt.imshow(ED_values[sorting_index], cmap=cm.coolwarm, interpolation='nearest', aspect='auto', vmin=-0.5, vmax=0.5)
+im = plt.imshow(Cliffs_values[sorting_index], cmap=cm.coolwarm, interpolation='nearest', aspect='auto', vmin=-0.5, vmax=0.5)
 
 plt.yticks([])
 

plot/ndd_by_nodes_sorting_just_2_Marozzi_Cliffs_average_FIG.py

@@ -17,13 +17,13 @@
 
 sorting = sys.argv[5]
 
-ED_average = sys.argv[6]
+Cliffs_average = sys.argv[6]
 
 Marozzi = sys.argv[7]
 
-ED_values = np.loadtxt(ED_average)
+Ciffs_values = np.loadtxt(Cliffs_average)
 
-ED_values *= -1.0 # switch order
+Ciffs_values *= -1.0 # switch order
 
 pvalues = np.loadtxt(Marozzi)
 
@@ -57,26 +57,6 @@
 
 # names = ["G2 - tp1\n", "G2 - tp2\n"]
 
-# names = ["G1 - control\n", "G2 - tp1\n"]
-
-# names = ["G2 - tp2\n", "G1 - EtOH\n"]
-
-# names = ["G1 - control\n", "G2 - tp2\n"]
-
-# names = ["G1 - EtOH\n", "G1 - control\n"]
-
-# names = ["G1 - EtOH\n", "G2 - tp1\n"]
-
-# names = ["G1 - EtOH\n", "G2 - tp2\n"]
-
-# names = ["G2 - tp1\n", "G1 - control\n"]
-
-# names = ["G2 - tp1\n", "G1 - EtOH\n"]
-
-# names = ["G2 - tp2\n", "G1 - control\n"]
-
-# names = ["G2 - tp2\n", "G2 - tp1\n"]
-
 ax = plt.subplot(144)
 
 ax.set_position([0.89, 0.0925, 0.04, 0.839])
@@ -105,11 +85,11 @@
 
 plt.title("$\delta$\n")
 
-ED_values = np.expand_dims(ED_values, axis=1)
+Ciffs_values = np.expand_dims(Ciffs_values, axis=1)
 
-print(np.amin(ED_values), np.amax(ED_values))
+print(np.amin(Ciffs_values), np.amax(Ciffs_values))
 
-im = plt.imshow(ED_values[sorting_index], cmap=cm.coolwarm, interpolation='nearest', aspect='auto', vmin=-0.5, vmax=0.5)
+im = plt.imshow(Ciffs_values[sorting_index], cmap=cm.coolwarm, interpolation='nearest', aspect='auto', vmin=-0.5, vmax=0.5)
 
 plt.yticks([])
 

plot/newd_by_nodes_sorting_just_2_Cliffs_average_orig_FIG_2delta_sml.py

@@ -17,11 +17,11 @@
 
 sorting = sys.argv[5]
 
-ED_average = sys.argv[6]
+Cliffs_average = sys.argv[6]
 
-ED_values = np.loadtxt(ED_average)
+Cliffs_values = np.loadtxt(Cliffs_average)
 
-ED_values *= -1.0 # switch order
+Cliffs_values *= -1.0 # switch order
 
 data = np.loadtxt(filename, skiprows=4)
 
@@ -58,39 +58,17 @@
 
 fig.tight_layout(rect=[0.13, 0.07, 0.88, 0.95], w_pad=3, h_pad=8.0)
 
-# names = ["G1 - control - part 1\n", "G1 - EtOH - part 2\n"]
-
 names = ["G1 - control\n", "G1 - EtOH\n"]
 
 # names = ["G2 - tp1\n", "G2 - tp2\n"]
 
-# names = ["G1 - control\n", "G2 - tp1\n"]
-
-# names = ["G2 - tp2\n", "G1 - EtOH\n"]
-
-# names = ["G1 - control\n", "G2 - tp2\n"]
-
-# names = ["G1 - EtOH\n", "G1 - control\n"]
-
-# names = ["G1 - EtOH\n", "G2 - tp1\n"]
-
-# names = ["G1 - EtOH\n", "G2 - tp2\n"]
-
-# names = ["G2 - tp1\n", "G1 - control\n"]
-
-# names = ["G2 - tp1\n", "G1 - EtOH\n"]
-
-# names = ["G2 - tp2\n", "G1 - control\n"]
-
-# names = ["G2 - tp2\n", "G2 - tp1\n"]
-
 ax = plt.subplot(144)
 
 ax.set_position([0.89, 0.0925, 0.06, 0.839])
 
 plt.title("$\delta_n \delta_p$\n")
 
-tmp = ED_values.copy()
+tmp = Cliffs_values.copy()
 
 tmp = np.abs(tmp)
 
@@ -131,9 +109,9 @@
 
 plt.title("$\delta_n \delta_p$\n")
 
-print(np.amin(ED_values), np.amax(ED_values))
+print(np.amin(Cliffs_values), np.amax(Cliffs_values))
 
-im = plt.imshow(ED_values[sorting_index], cmap=cm.coolwarm, interpolation='nearest', aspect='auto', vmin=-0.5, vmax=0.5)
+im = plt.imshow(Cliffs_values[sorting_index], cmap=cm.coolwarm, interpolation='nearest', aspect='auto', vmin=-0.5, vmax=0.5)
 
 plt.yticks([])