Octave Scripts and Results


In [1]:
setenv GNUTERM svg

Description

This repository contains scripts that demonstrates presence of high-order cell assemblies from wide-field calcium imaging data. This is a collabo.pngrative project with Prof. Spencer Smith and Dr. Yiyi Yu from University of North Carolina, Chapel Hill. Part of the results has been presented at the 3rd International Conference on Mathematical Neuroscience. The conference submission and presentation can be found in icmns_submission

Data preprocessing

Run the script combine_datasets.m to pre-process the data. The data should be stored in the data/ directory. Each directory contains a different experiment. For each experiment, there are simultaneous recording from the primary visual cortex (V1, Combo3_V1.mat) and a higher visual cortex (area AL, Combo3_AL.mat).

These files can be opened via Matlab, or using the mat4py library.


In [2]:
source('combine_datasets.m')

Decrease sparsity over trials

Run the script sparsity_over_trial.m to see the decreasing sparsity across repeated trials. You can change the area variable to see the result for V1 or AL. The summary plots will be saved in figures/sparsity/. Below is an example result for area AL.


In [2]:
source('sparsity_over_trial.m')


Number of cells: 352 
Number of cells: 375 
Number of cells: 147 
Number of cells: 226 
Number of cells: 375 
Number of cells: 590 
Number of cells: 189 
Gnuplot Produced by GNUPLOT 5.2 patchlevel 2 0.005 0.0055 0.006 0.0065 0.007 0.0075 0.008 0.0085 0 5 10 15 20 trial # gnuplot_plot_1a
Gnuplot Produced by GNUPLOT 5.2 patchlevel 2 0.004 0.005 0.006 0.007 0.008 0.009 0 5 10 15 20 Sparseness Trial # Area V1 Stimulus 2 n=2254 gnuplot_plot_1a

Presence of high-order cell assemblies

Run the script analyze_data_populaiton.m to see evidence of high-order cell assemblies. The script will generate the following result for each experiment session

  • sparsity over trial (figures/sparseness_over_timestim*, Fig. 1C)
  • subset index calculation (figures/subset_indexstim* Fig. 1D)
  • Number of unique high-order cell assemblies (cellAssembly_* Fig. 1E)
  • Accuracy of cell assembly occurance time (cellAssemblyjitter*, Fig. 1F)
  • Example raster plots (figures/exampleRasterStim* slide #27)

In [4]:
pkg load statistics
pkg load dataframe
source('analyze_data_populaiton.m')


Number of cells: 352 
 num binary words: 351
 mean frequency : 51.0171
 max frequency : 276
 num binary words: 39119
 mean frequency : 1.0464
 max frequency : 6
 num binary words: 19497
 mean frequency : 1.024
 max frequency : 4
 num binary words: 39119
 mean frequency : 1.0464
 max frequency : 6
 num binary words: 18668
 mean frequency : 1.0213
 max frequency : 4
 num binary words: 12790
 mean frequency : 1.0879
 max frequency : 4
 num binary words: 19450
 mean frequency : 1.021
 max frequency : 4
 num binary words: 12678
 mean frequency : 1.0312
 max frequency : 4
 num binary words: 18952
 mean frequency : 1.026
 max frequency : 5
 num binary words: 13210
 mean frequency : 1.0442
 max frequency : 3
 1.0464
 1.0260
 1.0442
 num binary words: 52207
 mean frequency : 1.0046
 max frequency : 3
 num binary words: 16169
 mean frequency : 1.0037
 max frequency : 4
 num binary words: 6882
 mean frequency : 1.0003
 max frequency : 2
 num binary words: 16131
 mean frequency : 1.0027
 max frequency : 2
 num binary words: 6394
 mean frequency : 1.0013
 max frequency : 2
 num binary words: 15377
 mean frequency : 1.0028
 max frequency : 2
 num binary words: 6579
 mean frequency : 1.0003
 max frequency : 2
 1.0046
 1.0028
 1.0003
 num binary words: 60537
 mean frequency : 1.0003
 max frequency : 2
 num binary words: 10576
 mean frequency : 1.0003
 max frequency : 2
 num binary words: 3119
 mean frequency : 1
 max frequency : 1
 num binary words: 12670
 mean frequency : 1.0002
 max frequency : 2
 num binary words: 2753
 mean frequency : 1
 max frequency : 1
 num binary words: 11357
 mean frequency : 1.0002
 max frequency : 2
 num binary words: 2223
 mean frequency : 1
 max frequency : 1
 1.0003
 1.0002
 1
 num binary words: 61826
 mean frequency : 1
 max frequency : 1
 num binary words: 4814
 mean frequency : 1
 max frequency : 1
 num binary words: 1074
 mean frequency : 1
 max frequency : 1
 num binary words: 9939
 mean frequency : 1.0002
 max frequency : 2
 num binary words: 1115
 mean frequency : 1
 max frequency : 1
 num binary words: 6278
 mean frequency : 1
 max frequency : 1
 num binary words: 895
 mean frequency : 1
 max frequency : 1
 1
 1
 1
Gnuplot Produced by GNUPLOT 5.2 patchlevel 2 0 5 10 15 20 0 5 10 15 20 25 30 Trial Time (s) Single Cell Spikes gnuplot_plot_1a gnuplot_plot_2a gnuplot_plot_3a gnuplot_plot_4a gnuplot_plot_5a gnuplot_plot_6a gnuplot_plot_7a gnuplot_plot_8a gnuplot_plot_9a gnuplot_plot_10a gnuplot_plot_11a gnuplot_plot_12a gnuplot_plot_13a gnuplot_plot_14a gnuplot_plot_15a gnuplot_plot_16a gnuplot_plot_17a gnuplot_plot_18a gnuplot_plot_19a gnuplot_plot_20a gnuplot_plot_21a gnuplot_plot_22a gnuplot_plot_23a gnuplot_plot_24a gnuplot_plot_25a gnuplot_plot_26a gnuplot_plot_27a gnuplot_plot_28a gnuplot_plot_29a gnuplot_plot_30a gnuplot_plot_31a gnuplot_plot_32a gnuplot_plot_33a gnuplot_plot_34a gnuplot_plot_35a gnuplot_plot_36a gnuplot_plot_37a gnuplot_plot_38a gnuplot_plot_39a gnuplot_plot_40a gnuplot_plot_41a gnuplot_plot_42a gnuplot_plot_43a gnuplot_plot_44a gnuplot_plot_45a gnuplot_plot_46a gnuplot_plot_47a gnuplot_plot_48a gnuplot_plot_49a gnuplot_plot_50a gnuplot_plot_51a gnuplot_plot_52a gnuplot_plot_53a gnuplot_plot_54a gnuplot_plot_55a gnuplot_plot_56a gnuplot_plot_57a gnuplot_plot_58a gnuplot_plot_59a gnuplot_plot_60a gnuplot_plot_61a gnuplot_plot_62a gnuplot_plot_63a gnuplot_plot_64a gnuplot_plot_65a gnuplot_plot_66a gnuplot_plot_67a gnuplot_plot_68a gnuplot_plot_69a gnuplot_plot_70a gnuplot_plot_71a gnuplot_plot_72a gnuplot_plot_73a gnuplot_plot_74a gnuplot_plot_75a gnuplot_plot_76a gnuplot_plot_77a gnuplot_plot_78a gnuplot_plot_79a gnuplot_plot_80a gnuplot_plot_81a gnuplot_plot_82a gnuplot_plot_83a gnuplot_plot_84a gnuplot_plot_85a gnuplot_plot_86a gnuplot_plot_87a gnuplot_plot_88a gnuplot_plot_89a gnuplot_plot_90a gnuplot_plot_91a gnuplot_plot_92a gnuplot_plot_93a gnuplot_plot_94a gnuplot_plot_95a gnuplot_plot_96a gnuplot_plot_97a gnuplot_plot_98a gnuplot_plot_99a gnuplot_plot_100a gnuplot_plot_101a gnuplot_plot_102a gnuplot_plot_103a gnuplot_plot_104a gnuplot_plot_105a gnuplot_plot_106a gnuplot_plot_107a gnuplot_plot_108a gnuplot_plot_109a gnuplot_plot_110a gnuplot_plot_111a gnuplot_plot_112a gnuplot_plot_113a gnuplot_plot_114a gnuplot_plot_115a gnuplot_plot_116a gnuplot_plot_117a gnuplot_plot_118a gnuplot_plot_119a gnuplot_plot_120a gnuplot_plot_121a gnuplot_plot_122a gnuplot_plot_123a gnuplot_plot_124a gnuplot_plot_125a gnuplot_plot_126a gnuplot_plot_127a gnuplot_plot_128a gnuplot_plot_129a gnuplot_plot_130a gnuplot_plot_131a gnuplot_plot_132a gnuplot_plot_133a gnuplot_plot_134a gnuplot_plot_135a gnuplot_plot_136a gnuplot_plot_137a gnuplot_plot_138a gnuplot_plot_139a gnuplot_plot_140a gnuplot_plot_141a gnuplot_plot_142a gnuplot_plot_143a gnuplot_plot_144a gnuplot_plot_145a gnuplot_plot_146a gnuplot_plot_147a gnuplot_plot_148a gnuplot_plot_149a gnuplot_plot_150a gnuplot_plot_151a gnuplot_plot_152a gnuplot_plot_153a gnuplot_plot_154a gnuplot_plot_155a gnuplot_plot_156a gnuplot_plot_157a gnuplot_plot_158a gnuplot_plot_159a gnuplot_plot_160a gnuplot_plot_161a gnuplot_plot_162a gnuplot_plot_163a gnuplot_plot_164a gnuplot_plot_165a gnuplot_plot_166a gnuplot_plot_167a gnuplot_plot_168a gnuplot_plot_169a gnuplot_plot_170a gnuplot_plot_171a gnuplot_plot_172a gnuplot_plot_173a gnuplot_plot_174a gnuplot_plot_175a gnuplot_plot_176a gnuplot_plot_177a gnuplot_plot_178a gnuplot_plot_179a gnuplot_plot_180a gnuplot_plot_181a gnuplot_plot_182a gnuplot_plot_183a gnuplot_plot_184a gnuplot_plot_185a gnuplot_plot_186a gnuplot_plot_187a gnuplot_plot_188a gnuplot_plot_189a gnuplot_plot_190a gnuplot_plot_191a gnuplot_plot_192a gnuplot_plot_193a gnuplot_plot_194a gnuplot_plot_195a gnuplot_plot_196a gnuplot_plot_197a gnuplot_plot_198a gnuplot_plot_199a gnuplot_plot_200a gnuplot_plot_201a gnuplot_plot_202a gnuplot_plot_203a gnuplot_plot_204a gnuplot_plot_205a gnuplot_plot_206a gnuplot_plot_207a gnuplot_plot_208a gnuplot_plot_209a gnuplot_plot_210a gnuplot_plot_211a gnuplot_plot_212a gnuplot_plot_213a gnuplot_plot_214a gnuplot_plot_215a gnuplot_plot_216a gnuplot_plot_217a gnuplot_plot_218a gnuplot_plot_219a gnuplot_plot_220a gnuplot_plot_221a gnuplot_plot_222a gnuplot_plot_223a gnuplot_plot_224a gnuplot_plot_225a gnuplot_plot_226a gnuplot_plot_227a gnuplot_plot_228a gnuplot_plot_229a gnuplot_plot_230a gnuplot_plot_231a gnuplot_plot_232a gnuplot_plot_233a gnuplot_plot_234a gnuplot_plot_235a gnuplot_plot_236a gnuplot_plot_237a gnuplot_plot_238a gnuplot_plot_239a gnuplot_plot_240a gnuplot_plot_241a gnuplot_plot_242a gnuplot_plot_243a gnuplot_plot_244a gnuplot_plot_245a gnuplot_plot_246a gnuplot_plot_247a gnuplot_plot_248a gnuplot_plot_249a gnuplot_plot_250a gnuplot_plot_251a gnuplot_plot_252a gnuplot_plot_253a gnuplot_plot_254a gnuplot_plot_255a gnuplot_plot_256a gnuplot_plot_257a gnuplot_plot_258a gnuplot_plot_259a gnuplot_plot_260a gnuplot_plot_261a gnuplot_plot_262a gnuplot_plot_263a gnuplot_plot_264a gnuplot_plot_265a gnuplot_plot_266a gnuplot_plot_267a gnuplot_plot_268a gnuplot_plot_269a gnuplot_plot_270a gnuplot_plot_271a gnuplot_plot_272a gnuplot_plot_273a gnuplot_plot_274a gnuplot_plot_275a gnuplot_plot_276a 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0 5 10 15 20 25 30 Probability Time (s) gnuplot_plot_1b gnuplot_plot_2b gnuplot_plot_3b 0 5 10 15 20 0 5 10 15 20 25 30 Trial Time (s) Shuffled Spikes gnuplot_plot_1c gnuplot_plot_2c gnuplot_plot_3c gnuplot_plot_4c gnuplot_plot_5c gnuplot_plot_6c gnuplot_plot_7c gnuplot_plot_8c gnuplot_plot_9c gnuplot_plot_10c gnuplot_plot_11c gnuplot_plot_12c gnuplot_plot_13c gnuplot_plot_14c gnuplot_plot_15c gnuplot_plot_16c gnuplot_plot_17c gnuplot_plot_18c gnuplot_plot_19c gnuplot_plot_20c gnuplot_plot_21c gnuplot_plot_22c gnuplot_plot_23c gnuplot_plot_24c gnuplot_plot_25c gnuplot_plot_26c gnuplot_plot_27c gnuplot_plot_28c gnuplot_plot_29c gnuplot_plot_30c gnuplot_plot_31c gnuplot_plot_32c gnuplot_plot_33c gnuplot_plot_34c gnuplot_plot_35c gnuplot_plot_36c gnuplot_plot_37c gnuplot_plot_38c gnuplot_plot_39c gnuplot_plot_40c gnuplot_plot_41c gnuplot_plot_42c gnuplot_plot_43c gnuplot_plot_44c gnuplot_plot_45c gnuplot_plot_46c gnuplot_plot_47c gnuplot_plot_48c gnuplot_plot_49c gnuplot_plot_50c gnuplot_plot_51c gnuplot_plot_52c gnuplot_plot_53c gnuplot_plot_54c gnuplot_plot_55c gnuplot_plot_56c gnuplot_plot_57c gnuplot_plot_58c gnuplot_plot_59c gnuplot_plot_60c gnuplot_plot_61c gnuplot_plot_62c gnuplot_plot_63c gnuplot_plot_64c gnuplot_plot_65c gnuplot_plot_66c gnuplot_plot_67c gnuplot_plot_68c gnuplot_plot_69c gnuplot_plot_70c gnuplot_plot_71c gnuplot_plot_72c gnuplot_plot_73c gnuplot_plot_74c gnuplot_plot_75c gnuplot_plot_76c gnuplot_plot_77c gnuplot_plot_78c gnuplot_plot_79c gnuplot_plot_80c gnuplot_plot_81c gnuplot_plot_82c gnuplot_plot_83c gnuplot_plot_84c gnuplot_plot_85c gnuplot_plot_86c gnuplot_plot_87c gnuplot_plot_88c gnuplot_plot_89c gnuplot_plot_90c gnuplot_plot_91c gnuplot_plot_92c gnuplot_plot_93c gnuplot_plot_94c gnuplot_plot_95c gnuplot_plot_96c gnuplot_plot_97c gnuplot_plot_98c gnuplot_plot_99c gnuplot_plot_100c gnuplot_plot_101c gnuplot_plot_102c gnuplot_plot_103c gnuplot_plot_104c gnuplot_plot_105c gnuplot_plot_106c gnuplot_plot_107c gnuplot_plot_108c gnuplot_plot_109c gnuplot_plot_110c gnuplot_plot_111c gnuplot_plot_112c gnuplot_plot_113c gnuplot_plot_114c gnuplot_plot_115c gnuplot_plot_116c gnuplot_plot_117c gnuplot_plot_118c gnuplot_plot_119c gnuplot_plot_120c gnuplot_plot_121c gnuplot_plot_122c gnuplot_plot_123c gnuplot_plot_124c gnuplot_plot_125c gnuplot_plot_126c gnuplot_plot_127c gnuplot_plot_128c gnuplot_plot_129c gnuplot_plot_130c gnuplot_plot_131c gnuplot_plot_132c gnuplot_plot_133c gnuplot_plot_134c gnuplot_plot_135c gnuplot_plot_136c gnuplot_plot_137c gnuplot_plot_138c gnuplot_plot_139c gnuplot_plot_140c gnuplot_plot_141c gnuplot_plot_142c gnuplot_plot_143c gnuplot_plot_144c gnuplot_plot_145c gnuplot_plot_146c gnuplot_plot_147c gnuplot_plot_148c gnuplot_plot_149c gnuplot_plot_150c gnuplot_plot_151c gnuplot_plot_152c gnuplot_plot_153c gnuplot_plot_154c gnuplot_plot_155c gnuplot_plot_156c gnuplot_plot_157c gnuplot_plot_158c gnuplot_plot_159c gnuplot_plot_160c gnuplot_plot_161c gnuplot_plot_162c gnuplot_plot_163c gnuplot_plot_164c gnuplot_plot_165c gnuplot_plot_166c gnuplot_plot_167c gnuplot_plot_168c gnuplot_plot_169c gnuplot_plot_170c gnuplot_plot_171c gnuplot_plot_172c gnuplot_plot_173c gnuplot_plot_174c gnuplot_plot_175c gnuplot_plot_176c gnuplot_plot_177c gnuplot_plot_178c gnuplot_plot_179c gnuplot_plot_180c gnuplot_plot_181c gnuplot_plot_182c gnuplot_plot_183c gnuplot_plot_184c gnuplot_plot_185c gnuplot_plot_186c gnuplot_plot_187c gnuplot_plot_188c gnuplot_plot_189c gnuplot_plot_190c gnuplot_plot_191c gnuplot_plot_192c gnuplot_plot_193c gnuplot_plot_194c gnuplot_plot_195c gnuplot_plot_196c gnuplot_plot_197c gnuplot_plot_198c gnuplot_plot_199c gnuplot_plot_200c gnuplot_plot_201c gnuplot_plot_202c gnuplot_plot_203c gnuplot_plot_204c gnuplot_plot_205c gnuplot_plot_206c gnuplot_plot_207c gnuplot_plot_208c gnuplot_plot_209c gnuplot_plot_210c gnuplot_plot_211c gnuplot_plot_212c gnuplot_plot_213c gnuplot_plot_214c gnuplot_plot_215c gnuplot_plot_216c gnuplot_plot_217c gnuplot_plot_218c gnuplot_plot_219c gnuplot_plot_220c gnuplot_plot_221c gnuplot_plot_222c gnuplot_plot_223c gnuplot_plot_224c gnuplot_plot_225c gnuplot_plot_226c gnuplot_plot_227c gnuplot_plot_228c gnuplot_plot_229c gnuplot_plot_230c gnuplot_plot_231c gnuplot_plot_232c gnuplot_plot_233c gnuplot_plot_234c gnuplot_plot_235c gnuplot_plot_236c gnuplot_plot_237c gnuplot_plot_238c gnuplot_plot_239c gnuplot_plot_240c gnuplot_plot_241c gnuplot_plot_242c gnuplot_plot_243c gnuplot_plot_244c gnuplot_plot_245c gnuplot_plot_246c gnuplot_plot_247c gnuplot_plot_248c gnuplot_plot_249c gnuplot_plot_250c gnuplot_plot_251c gnuplot_plot_252c gnuplot_plot_253c gnuplot_plot_254c gnuplot_plot_255c gnuplot_plot_256c gnuplot_plot_257c gnuplot_plot_258c gnuplot_plot_259c gnuplot_plot_260c gnuplot_plot_261c gnuplot_plot_262c gnuplot_plot_263c gnuplot_plot_264c gnuplot_plot_265c gnuplot_plot_266c gnuplot_plot_267c gnuplot_plot_268c gnuplot_plot_269c gnuplot_plot_270c gnuplot_plot_271c gnuplot_plot_272c gnuplot_plot_273c gnuplot_plot_274c gnuplot_plot_275c gnuplot_plot_276c 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0 5 10 15 20 25 30 Probability Time (s) gnuplot_plot_1d gnuplot_plot_2d gnuplot_plot_3d 0 5 10 15 20 0 5 10 15 20 25 30 Trial Time (s) Cell Assembly gnuplot_plot_1e gnuplot_plot_2e gnuplot_plot_3e gnuplot_plot_4e gnuplot_plot_5e gnuplot_plot_6e 0 5 10 15 20 0 5 10 15 20 25 30 Trial Time (s) Cell Assembly with shuffled spikes gnuplot_plot_1f 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0 5 10 15 20 25 30 Probability Time (s) gnuplot_plot_1g 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0 5 10 15 20 25 30 Probability Time (s) gnuplot_plot_1h
Gnuplot Produced by GNUPLOT 5.2 patchlevel 2 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 Data Shuffled Poisson Subset index V1 gnuplot_plot_1a gnuplot_plot_2a gnuplot_plot_3a gnuplot_plot_4a gnuplot_plot_5a gnuplot_plot_6a gnuplot_plot_7a
Gnuplot Produced by GNUPLOT 5.2 patchlevel 2 0.002 0.004 0.006 0.008 0.01 0.012 0 5 10 15 20 sparseness trial # area V1 stim 2 gnuplot_plot_1a
Gnuplot Produced by GNUPLOT 5.2 patchlevel 2 0 5000 10000 15000 20000 25000 30000 35000 40000 Data Shuffled Poisson # 3-cell assemblies gnuplot_plot_1a gnuplot_plot_2a gnuplot_plot_3a gnuplot_plot_4a gnuplot_plot_5a gnuplot_plot_6a gnuplot_plot_7a gnuplot_plot_8a 0 500 1000 1500 2000 Data Shuffled Poisson # repeatable 3-cell assemblies gnuplot_plot_1b gnuplot_plot_2b gnuplot_plot_3b gnuplot_plot_4b gnuplot_plot_5b gnuplot_plot_6b gnuplot_plot_7b gnuplot_plot_8b 0 10000 20000 30000 40000 50000 60000 Data Shuffled Poisson # 4-cell assemblies gnuplot_plot_1c gnuplot_plot_2c gnuplot_plot_3c gnuplot_plot_4c gnuplot_plot_5c gnuplot_plot_6c gnuplot_plot_7c gnuplot_plot_8c 0 50 100 150 200 250 Data Shuffled Poisson # repeatable 4-cell assemblies gnuplot_plot_1d gnuplot_plot_2d gnuplot_plot_3d gnuplot_plot_4d gnuplot_plot_5d gnuplot_plot_6d gnuplot_plot_7d gnuplot_plot_8d 0 10000 20000 30000 40000 50000 60000 70000 Data Shuffled Poisson # 5-cell assemblies gnuplot_plot_1e gnuplot_plot_2e gnuplot_plot_3e gnuplot_plot_4e gnuplot_plot_5e gnuplot_plot_6e gnuplot_plot_7e gnuplot_plot_8e 0 5 10 15 20 Data Shuffled Poisson # repeatable 5-cell assemblies gnuplot_plot_1f gnuplot_plot_2f gnuplot_plot_3f gnuplot_plot_4f gnuplot_plot_5f gnuplot_plot_6f gnuplot_plot_7f gnuplot_plot_8f 0 10000 20000 30000 40000 50000 60000 70000 Data Shuffled Poisson # 6-cell assemblies gnuplot_plot_1g gnuplot_plot_2g gnuplot_plot_3g gnuplot_plot_4g gnuplot_plot_5g gnuplot_plot_6g gnuplot_plot_7g gnuplot_plot_8g 0 0.5 1 1.5 2 Data Shuffled Poisson # repeatable 6-cell assemblies gnuplot_plot_1h gnuplot_plot_2h gnuplot_plot_3h gnuplot_plot_4h gnuplot_plot_5h gnuplot_plot_6h gnuplot_plot_7h gnuplot_plot_8h
Gnuplot Produced by GNUPLOT 5.2 patchlevel 2 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 -1 -0.5 0 0.5 1 Prob. of observing repeated cell assembly Time jitter (sec) 3-order cell assembly 3-order cell assembly single cell single cell gnuplot_plot_3a
Gnuplot Produced by GNUPLOT 5.2 patchlevel 2 0 10000 20000 30000 40000 50000 60000 70000 3 4 5 6 # unique cell assembly Cell Assembly Order Data Data gnuplot_plot_2a gnuplot_plot_3a gnuplot_plot_4a gnuplot_plot_5a gnuplot_plot_6a gnuplot_plot_7a gnuplot_plot_8a gnuplot_plot_9a Shuffled Shuffled gnuplot_plot_11a gnuplot_plot_12a gnuplot_plot_13a gnuplot_plot_14a gnuplot_plot_15a gnuplot_plot_16a gnuplot_plot_17a Poisson Poisson gnuplot_plot_19a gnuplot_plot_20a gnuplot_plot_21a gnuplot_plot_22a gnuplot_plot_23a gnuplot_plot_24a gnuplot_plot_25a gnuplot_plot_26a gnuplot_plot_27a gnuplot_plot_28a