From 74b523eba183ebbeee3b5086529a9e440f775504 Mon Sep 17 00:00:00 2001
From: efish <efish@hummingbird>
Date: Thu, 15 Nov 2018 10:14:32 +0100
Subject: [PATCH 1/2] vector_phase.py

---
 code/analysis_rs.py |  4 ++++
 code/eod_cl.py      |  0
 code/utility.py     | 12 +++++++++++-
 3 files changed, 15 insertions(+), 1 deletion(-)
 delete mode 100644 code/eod_cl.py

diff --git a/code/analysis_rs.py b/code/analysis_rs.py
index 3a3cdf8..f7973f1 100644
--- a/code/analysis_rs.py
+++ b/code/analysis_rs.py
@@ -52,5 +52,9 @@ for i, idx in enumerate(eod_idx):
 	#embed()
 	#exit()
 	plt.plot(time[int(idx):int(eod_idx[i+1])], eod[int(idx):int(eod_idx[i+1])])
+
+
+
+
 plt.show()
 
diff --git a/code/eod_cl.py b/code/eod_cl.py
deleted file mode 100644
index e69de29..0000000
diff --git a/code/utility.py b/code/utility.py
index 5a55632..86f2407 100644
--- a/code/utility.py
+++ b/code/utility.py
@@ -5,4 +5,14 @@ def zero_crossing(eod,time):
 	threshold = 0;
 	shift_eod = np.roll(eod, 1)
 	eod_times = time[(eod >= threshold) & (shift_eod < threshold)]
-	return eod_times
+	sampling_rate = 40000.0
+	eod_idx = eod_times*sampling_rate
+	return eod_idx
+
+
+def vector_strength(spike_times, eod_durations)
+	alphas = spike_times/ eod_durations
+	cs = (1/len(spike_times))*np.sum(np.cos(alphas))^2
+	sn = (1/len(spike_times))*np.sum(np.sin(alphas))^2
+	vs = np.sprt(cs+sn)
+	return vs

From 477fa15dc430b3d9c42ac3e40c59d67b3075c007 Mon Sep 17 00:00:00 2001
From: efish <efish@localhost.localdomain>
Date: Thu, 15 Nov 2018 10:30:29 +0100
Subject: [PATCH 2/2] liste erstellt

---
 code/liste.py | 43 +++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 43 insertions(+)
 create mode 100644 code/liste.py

diff --git a/code/liste.py b/code/liste.py
new file mode 100644
index 0000000..cba7d05
--- /dev/null
+++ b/code/liste.py
@@ -0,0 +1,43 @@
+# 9.11.18
+
+aa: quality: poor, depth: -1341
+ab: quality: poor, depth: -1341
+ac: quality: good, depth: -1341
+ad: quality: good, depth: -200
+ae: quality: good, depth: -200
+af: quality: good, depth: -200
+ag: no info.dat
+
+
+# 13.11.18
+
+aa: good, -30 µm, maybe no reaction
+ab: good, -309 µm
+ac: poor, -309 µm
+ad: fair, -360 µm
+ae: fair, -350 µm
+af: good, -440 µm, bursting
+ag: fair, -174 µm
+ah: good, -209 µm
+ai: good, -66.9 µm
+aj: good, -132 µm
+ak: good, -284 µm
+al: good, -286 µm
+
+
+# 14.11.18
+
+aa: good, -184 µm
+ab: fair, -279 µm, no reaction
+ac: fair, -60 µm
+ad: good, -357 µm
+ae: fair, -357 µm
+af: fair, -527 µm, only baseline
+ag: fair, -533 µm
+ah: poor, -505 µm 
+ai: good, -500 µm, still same cell 3x
+aj: poor, -314 µm,no modulation
+ak: good, -140 µm
+al: good, -280 µm
+am: good, -320 µm
+an: good, -434 µm