Changeset 97 for trunk/examples/retina/benchmark_noise.py

Timestamp:

01/14/08 16:15:03 (4 years ago)

Author:

LaurentPerrinet

Message:

adapting the retina in the examples folder to work with the new framework

Files:

• trunk/examples/retina/benchmark_noise.py (modified) (4 diffs)

trunk/examples/retina/benchmark_noise.py

@@ -5 +5 @@
 ==================
 
-Just studying how different background noise current is integrated by IAF_SFA neurons on one fiber.
+Just studying how different background noise current is integrated by the neurons on one fiber.
 
 Laurent Perrinet, INCM, CNRS
@@ -13 +13 @@
 """
 
-import os, datetime, numpy, pylab
+import os, datetime, numpy
 import retina as model
 import NeuroTools.benchmark as benchmark
@@ -29 +29 @@
     """
 
-    pylab.rcParams.update(pylab_params(fig_width_pt =497.9/2., ratio = 1))
-    pylab.figure(num = 1, dpi=200, facecolor='w', edgecolor='k')
-
-    Lmargin, Rmargin, dmargin, umargin = 0.15, 0.05, 0.15,  0.05
-    pylab.axes([Lmargin, dmargin , 1.0 - Rmargin- Lmargin,1.0-umargin-dmargin]) # [left, bottom, width, height]
 
     params = benchmark.get('params')
@@ -40 +35 @@
     noise_std, CRF = numpy.zeros(len(experiments)), numpy.zeros(len(experiments))
     for count, experiment in enumerate(experiments):
-        out_ON_DATA = benchmark.get('out',experiment)['out_ON_DATA'].as_list_id_list_time()
+        out_ON_DATA = benchmark.get('out',experiment)['out_ON_DATA']#.as_list_id_list_time()
         noise_std[count] = benchmark.get('experiments')[experiment]['noise_std']
-        CRF[count] = len(out_ON_DATA[0]) / N / N / simtime * 1000.0
+        CRF[count] = out_ON_DATA.mean_rate()#[0]) / N / N / simtime * 1000.0
 
 
+    import pylab
+
+    #pylab.rcParams.update(pylab_params(fig_width_pt =497.9/5, ratio = 1/1.2))
+    pylab.figure(num = 1)#, dpi=200, facecolor='w', edgecolor='k')
+
+    #Lmargin, Rmargin, dmargin, umargin = 0.15, 0.05, 0.15,  0.05
+    #pylab.axes([Lmargin, dmargin , 1.0 - Rmargin- Lmargin,1.0-umargin-dmargin]) # [left, bottom, width, height]
+
     #pylab.subplot(111)
-    pylab.plot(noise_std,CRF,'go')
+    idx_ = numpy.argsort(noise_std) # extract lines
+    pylab.plot(noise_std[idx_],CRF[idx_],'go-', label='line 1', linewidth=2)
     pylab.ylabel('Firing Frequency (Hz)')
     pylab.xlabel('Noise amplitude')
-    pylab.axis('tight')
-    pylab.savefig(benchmark.name + '/benchmark_noise.png')#, dpi=300) #
+    #pylab.axis('tight')
+    pylab.savefig(benchmark.filename + '/benchmark_noise.pdf')#, dpi=300) #
     #command = 'epstopdf  "%s" "%s"'% (pdffile, psfile)
     #os.system(command)