diff --git a/explicability/compute_shap_inter_vals.py b/explicability/compute_shap_inter_vals.py
index 828a1bf0a0d646cf2ccf08dbae25f002c859e460..c2f51001aabc3484ab8e3268d5ccfb836c7abf36 100644
--- a/explicability/compute_shap_inter_vals.py
+++ b/explicability/compute_shap_inter_vals.py
@@ -60,9 +60,10 @@ if __name__ == "__main__":
         X_test = data_dic['X_test_' + group]
         y_test = data_dic['y_test_' + group]
         for j, method in enumerate(['', '', 'over_', 'under_']):
-            if j != 1:
-                print('Skip')
-                continue
+            # Remove (used to isolate RF)
+            # if j != 1:
+            #     print('Skip')
+            #     continue
             print(f"{group}-{method_names[j]}")
             method_name = method_names[j]
             model_name = model_choices[method_name]
diff --git a/explicability/shap_plots.py b/explicability/shap_plots.py
new file mode 100644
index 0000000000000000000000000000000000000000..b8f5134115718aea03aa45aea7a891ec25a88838
--- /dev/null
+++ b/explicability/shap_plots.py
@@ -0,0 +1,54 @@
+# Libraries
+# --------------------------------------------------------------------------------------------------------
+import pandas as pd
+import numpy as np
+import shap
+# --------------------------------------------------------------------------------------------------------
+
+# Reading test data
+# --------------------------------------------------------------------------------------------------------
+def read_test_data(attribute_names):
+    # Load test data
+    X_test_pre = np.load('../gen_train_data/data/output/pre/X_test_pre.npy', allow_pickle=True)
+    y_test_pre = np.load('../gen_train_data/data/output/pre/y_test_pre.npy', allow_pickle=True)
+    X_test_post = np.load('../gen_train_data/data/output/post/X_test_post.npy', allow_pickle=True)
+    y_test_post = np.load('../gen_train_data/data/output/post/y_test_post.npy', allow_pickle=True)
+
+    # Type conversion needed    
+    data_dic = {
+        "X_test_pre": pd.DataFrame(X_test_pre, columns=attribute_names).convert_dtypes(),
+        "y_test_pre": y_test_pre,
+        "X_test_post": pd.DataFrame(X_test_post, columns=attribute_names).convert_dtypes(),
+        "y_test_post": y_test_post,
+    }
+    return data_dic
+# --------------------------------------------------------------------------------------------------------
+
+if __name__ == "__main__":
+    # Setup
+    # --------------------------------------------------------------------------------------------------------
+    # Retrieve attribute names in order
+    attribute_names = list(np.load('../gen_train_data/data/output/attributes.npy', allow_pickle=True))
+    # Reading data
+    data_dic = read_test_data(attribute_names)
+    method_names = {
+        0: "ORIG",
+        1: "ORIG_CW",
+        2: "OVER",
+        3: "UNDER"
+    }
+    # --------------------------------------------------------------------------------------------------------
+
+    # Plot generation
+    # --------------------------------------------------------------------------------------------------------
+    for i, group in enumerate(['pre', 'post']):
+        # Get test dataset based on group, add column names
+        X_test = data_dic['X_test_' + group]
+        y_test = data_dic['y_test_' + group]
+        for j, method in enumerate(['', '', 'over_', 'under_']):
+            print(f"{group}-{method_names[j]}")
+            method_name = method_names[j]
+            shap_vals = np.load(f'./output/shap_values/{group}_{method_name}.npy')
+            print(f'Loaded SHAP values. Shape: {shap_vals.shape}')
+            shap_inter_vals = np.load(f'./output/shap_inter_values/{group}_{method_name}.npy')
+            print(f'Loaded SHAP INTER values. Shape: {shap_inter_vals.shape}')