+averaging example, removed empty averaging sec from docs

2 files changed, 42 insertions, 5 deletions

diff --git a/Readme.rst b/Readme.rst
index c6ded74..8228d7f 100644
--- a/Readme.rst
+++ b/Readme.rst
@@ -249,11 +249,6 @@ The following example computes and plots the ISI- and SPIKE-distance matrix, whe
     plt.show()
 
 
-Time Averages
--------------
-
-
-
 
 .. _ISI: http://www.scholarpedia.org/article/Measures_of_spike_train_synchrony#ISI-distance
 .. _SPIKE: http://www.scholarpedia.org/article/SPIKE-distance
diff --git a/examples/averages.py b/examples/averages.py
new file mode 100644
index 0000000..a06871a
--- /dev/null
+++ b/examples/averages.py
@@ -0,0 +1,42 @@
+""" averages.py
+
+Simple example showing how to compute averages of distance profiles
+
+Copyright 2014, Mario Mulansky <mario.mulansky@gmx.net>
+
+Distributed under the BSD License
+"""
+
+from __future__ import print_function
+
+import pyspike as spk
+
+spike_trains = spk.load_spike_trains_from_txt("PySpike_testdata.txt",
+                                              time_interval=(0, 4000))
+
+f = spk.isi_profile(spike_trains[0], spike_trains[1])
+
+print("ISI-distance: %.8f" % f.avrg())
+
+isi1 = f.avrg(interval=(0, 1000))
+isi2 = f.avrg(interval=(1000, 2000))
+isi3 = f.avrg(interval=(2000, 3000))
+isi4 = f.avrg(interval=(3000, 4000))
+
+print("ISI-distance (0-1000):    %.8f" % isi1)
+print("ISI-distance (1000-2000): %.8f" % isi2)
+print("ISI-distance (2000-3000): %.8f" % isi3)
+print("ISI-distance (3000-4000): %.8f" % isi4)
+print()
+
+print("SPIKE-distance: %.8f" % f.avrg())
+
+spike1 = f.avrg(interval=(0, 1000))
+spike2 = f.avrg(interval=(1000, 2000))
+spike3 = f.avrg(interval=(2000, 3000))
+spike4 = f.avrg(interval=(3000, 4000))
+
+print("SPIKE-distance (0-1000):    %.8f" % spike1)
+print("SPIKE-distance (1000-2000): %.8f" % spike2)
+print("SPIKE-distance (2000-3000): %.8f" % spike3)
+print("SPIKE-distance (3000-4000): %.8f" % spike4)