diff --git a/GRO_seq_TFSEE/box_plot_cluster_1.png b/GRO_seq_TFSEE/box_plot_cluster_1.png
index d7fcb41e968d761ab856b20d8aa7ad93335af9bb..7e7b6c801ca543f39565bddec369ca0bfaccd11a 100644
Binary files a/GRO_seq_TFSEE/box_plot_cluster_1.png and b/GRO_seq_TFSEE/box_plot_cluster_1.png differ
diff --git a/GRO_seq_TFSEE/box_plot_cluster_2.png b/GRO_seq_TFSEE/box_plot_cluster_2.png
index eb3d945cf08867214337728b2e2e3cf8df80816e..3fccabfbb497d8881379df81c9bc86b9d90d7bfc 100644
Binary files a/GRO_seq_TFSEE/box_plot_cluster_2.png and b/GRO_seq_TFSEE/box_plot_cluster_2.png differ
diff --git a/GRO_seq_TFSEE/box_plot_cluster_3.png b/GRO_seq_TFSEE/box_plot_cluster_3.png
index d3a847fb44059d052926279776eff2ed8a912819..e7eb9b31ef281a277cc4337394bc93a01707ecef 100644
Binary files a/GRO_seq_TFSEE/box_plot_cluster_3.png and b/GRO_seq_TFSEE/box_plot_cluster_3.png differ
diff --git a/GRO_seq_TFSEE/box_plot_cluster_3_tfs_fpkm.png b/GRO_seq_TFSEE/box_plot_cluster_3_tfs_fpkm.png
index f692823b4e0c9bb76c7f19e733f995499e40cdb0..0e0f72dbd9e64d5567a9e1cca81f8e74eb6512d5 100644
Binary files a/GRO_seq_TFSEE/box_plot_cluster_3_tfs_fpkm.png and b/GRO_seq_TFSEE/box_plot_cluster_3_tfs_fpkm.png differ
diff --git a/GRO_seq_TFSEE/box_plot_cluster_4.png b/GRO_seq_TFSEE/box_plot_cluster_4.png
index ea948ae27e5472d232b943fa134a66ed2e971866..e08838f8df26be9a6cddb6509a510479d080c48d 100644
Binary files a/GRO_seq_TFSEE/box_plot_cluster_4.png and b/GRO_seq_TFSEE/box_plot_cluster_4.png differ
diff --git a/GRO_seq_TFSEE/box_plot_cluster_4_tfs_fpkm.png b/GRO_seq_TFSEE/box_plot_cluster_4_tfs_fpkm.png
index 741510df362de2bdfb8a91ebf8df4da56325feb4..3e027b9ff4f69139f3e20f530339f6c2573247b5 100644
Binary files a/GRO_seq_TFSEE/box_plot_cluster_4_tfs_fpkm.png and b/GRO_seq_TFSEE/box_plot_cluster_4_tfs_fpkm.png differ
diff --git a/GRO_seq_TFSEE/matrix_analysis.py b/GRO_seq_TFSEE/matrix_analysis.py
index b1ebb3fda2ee97e3712a2faf1782ce0ae1a1fa59..7ada58ec0d62b62da5f72c4aae4309b3d7ab45af 100644
--- a/GRO_seq_TFSEE/matrix_analysis.py
+++ b/GRO_seq_TFSEE/matrix_analysis.py
@@ -339,7 +339,6 @@ needed_rows = [row for row in rpkm_scaled.index if row in list(motif_enhancers_s
 rpkm_robust_filtered= rpkm_scaled.loc[needed_rows]
 H3K27ac_robust_filtered= H3K27ac_scaled.loc[needed_rows]
 H3K4me1_robust_filtered= H3K4me1_scaled.loc[needed_rows]
-fd
 
 # 1. add H3K27ac and H3K4me1 signal
 H3K27ac_H3K4me1 = H3K27ac_scaled.add(H3K4me1_scaled)
@@ -453,8 +452,8 @@ colors = ["#FFD66F","#2E6A44","#862743", "#4FA6C7", "#3398CC"]
 medianprops = dict(linestyle='-', linewidth=2, color='black')
 
 box = cell_tf_values_std_cluster_1.boxplot(column=['ES_D0', 'ES_D2', 'ES_D5', 'ES_D7',  'ES_D10'],patch_artist=True,showfliers=False,manage_xticks=False,widths = 0.6, medianprops = medianprops)
-plt.setp(box['whiskers'], color='k', linestyle='-', linewidth = 3)
-plt.setp(box['boxes'], color='k', linestyle='-', linewidth = 3)
+plt.setp(box['whiskers'], color='k', linestyle='-', linewidth = 5)
+plt.setp(box['boxes'], color='k', linestyle='-', linewidth = 5)
 
 for patch, color in zip(box['boxes'], colors):
     patch.set_facecolor(color)
@@ -470,9 +469,8 @@ plt.clf()
 
 
 box = cell_tf_values_std_cluster_2.boxplot(column=['ES_D0', 'ES_D2', 'ES_D5', 'ES_D7',  'ES_D10'],patch_artist=True,showfliers=False,manage_xticks=False,widths = 0.6, medianprops = medianprops)
-plt.setp(box['whiskers'], color='k', linestyle='-', linewidth = 3)
-plt.setp(box['boxes'], color='k', linestyle='-', linewidth = 3)
-plt.setp(box['boxes'], color='k', linestyle='-', linewidth = 1.5)
+plt.setp(box['whiskers'], color='k', linestyle='-', linewidth = 5)
+plt.setp(box['boxes'], color='k', linestyle='-', linewidth = 5)
 
 for patch, color in zip(box['boxes'], colors):
     patch.set_facecolor(color)
@@ -487,8 +485,8 @@ plt.savefig('box_plot_cluster_2.png')
 plt.clf()
 
 box = cell_tf_values_std_cluster_3.boxplot(column=['ES_D0', 'ES_D2', 'ES_D5', 'ES_D7',  'ES_D10'],patch_artist=True,showfliers=False,manage_xticks=False,widths = 0.6, medianprops = medianprops)
-plt.setp(box['whiskers'], color='k', linestyle='-', linewidth = 3)
-plt.setp(box['boxes'], color='k', linestyle='-', linewidth = 3)
+plt.setp(box['whiskers'], color='k', linestyle='-', linewidth = 5)
+plt.setp(box['boxes'], color='k', linestyle='-', linewidth = 5)
 
 for patch, color in zip(box['boxes'], colors):
     patch.set_facecolor(color)
@@ -504,8 +502,8 @@ plt.clf()
 
 
 box = cell_tf_values_std_cluster_4.boxplot(column=['ES_D0', 'ES_D2', 'ES_D5', 'ES_D7',  'ES_D10'],patch_artist=True,showfliers=False,manage_xticks=False,widths = 0.6, medianprops = medianprops)
-plt.setp(box['whiskers'], color='k', linestyle='-', linewidth = 3)
-plt.setp(box['boxes'], color='k', linestyle='-', linewidth = 3)
+plt.setp(box['whiskers'], color='k', linestyle='-', linewidth = 5)
+plt.setp(box['boxes'], color='k', linestyle='-', linewidth = 5)
 
 for patch, color in zip(box['boxes'], colors):
     patch.set_facecolor(color)
@@ -592,8 +590,8 @@ cluster4_tfs = tf_cell_lines.loc[cell_tf_values_std_cluster_4.index.values]
 cluster4_tfs.to_csv("cluster_4_tfs.csv")
 
 box = cluster4_tfs.boxplot(column=['ES_D0', 'ES_D2', 'ES_D5', 'ES_D7',  'ES_D10'],patch_artist=True,showfliers=False,manage_xticks=False,widths = 0.6, medianprops = medianprops)
-plt.setp(box['whiskers'], color='k', linestyle='-', linewidth = 3)
-plt.setp(box['boxes'], color='k', linestyle='-', linewidth = 3)
+plt.setp(box['whiskers'], color='k', linestyle='-', linewidth = 5)
+plt.setp(box['boxes'], color='k', linestyle='-', linewidth = 5)
 
 for patch, color in zip(box['boxes'], colors):
     patch.set_facecolor(color)
@@ -671,8 +669,8 @@ cluster3_tfs = tf_cell_lines.loc[cell_tf_values_std_cluster_3.index.values]
 cluster3_tfs.to_csv("cluster_3_tfs.csv")
 
 box = cluster3_tfs.boxplot(column=['ES_D0', 'ES_D2', 'ES_D5', 'ES_D7',  'ES_D10'],patch_artist=True,showfliers=False,manage_xticks=False,widths = 0.6, medianprops = medianprops)
-plt.setp(box['whiskers'], color='k', linestyle='-', linewidth = 3)
-plt.setp(box['boxes'], color='k', linestyle='-', linewidth = 3)
+plt.setp(box['whiskers'], color='k', linestyle='-', linewidth = 5)
+plt.setp(box['boxes'], color='k', linestyle='-', linewidth = 5)
 
 for patch, color in zip(box['boxes'], colors):
     patch.set_facecolor(color)