gin commit from UX410UAK

New files: 2

python_code/read_data_pkl.py

@@ -0,0 +1,19 @@
+import pandas as pd
+import matplotlib.pylab as plt
+
+filename = '../Data/interictal/01ref/pkl/activity_SU_4ms.pkl'
+
+data = pd.read_pickle(filename)
+print(data.columns)
+
+existing_binnings = data['binning'].tolist()
+
+for binning in existing_binnings:
+    activity = data[data['binning'] == binning]['activity'].tolist()[0]
+    plt.figure()
+    plt.plot(activity)
+    plt.xlabel('Time [4 ms]')
+    plt.ylabel('Activity [# recorded spikes]')
+    plt.title('Activity in the {} hemisphere'.format(binning))
+
+plt.show()

python_code/read_data_txt.py

@@ -0,0 +1,17 @@
+import matplotlib.pylab as plt
+
+filename = '../Data/interictal/01ref/txt/activity_SU_4ms_right.txt'
+
+with open(filename, 'r') as f:
+    lines = f.readlines()
+    activity = []
+    for line in lines:
+        if line.startswith('#'):
+            continue
+        activity.append(line.split(' ')[1])
+
+plt.figure()
+plt.plot(activity)
+plt.xlabel('Time [4 ms]')
+plt.ylabel('Activity [# recorded spikes]')
+plt.show()