summaryrefslogtreecommitdiff
path: root/examples/distance_matrix.py
diff options
context:
space:
mode:
authorMario Mulansky <mario.mulansky@gmx.net>2014-10-21 19:05:39 +0200
committerMario Mulansky <mario.mulansky@gmx.net>2014-10-21 19:05:39 +0200
commit1b0421a207f5a1eb43b12bb18d5e783e753b739e (patch)
treef3ad4ed8075ff731c7790777e105f0ee41e64efa /examples/distance_matrix.py
parent89287c100c7a49d8e5ee5ca848a77550a4ea36b5 (diff)
doc: distance matrix, improved example
Diffstat (limited to 'examples/distance_matrix.py')
-rw-r--r--examples/distance_matrix.py33
1 files changed, 33 insertions, 0 deletions
diff --git a/examples/distance_matrix.py b/examples/distance_matrix.py
new file mode 100644
index 0000000..38bd9c8
--- /dev/null
+++ b/examples/distance_matrix.py
@@ -0,0 +1,33 @@
+""" distance_matrix.py
+
+Simple example showing how to compute the isi distance matrix of a set of spike
+trains.
+
+Copyright 2014, Mario Mulansky <mario.mulansky@gmx.net>
+
+Distributed under the BSD License
+"""
+
+
+from __future__ import print_function
+
+import matplotlib.pyplot as plt
+
+import pyspike as spk
+
+# first load the data, interval ending time = 4000, start=0 (default)
+spike_trains = spk.load_spike_trains_from_txt("PySpike_testdata.txt", 4000)
+
+print(len(spike_trains))
+
+plt.figure()
+isi_distance = spk.isi_distance_matrix(spike_trains)
+plt.imshow(isi_distance, interpolation='none')
+plt.title("ISI-distance")
+
+plt.figure()
+spike_distance = spk.spike_distance_matrix(spike_trains)
+plt.imshow(spike_distance, interpolation='none')
+plt.title("SPIKE-distance")
+
+plt.show()