fast / slow with top/bottom 20%

Author: EloyPB

batch_analysis.py

@@ -15,11 +15,11 @@ def analyze(identifier):
                     pickles_root_path=pickles_path)
 
     pf = PlaceFields.current_instance(config)
-    # pf.plot_activations(fig_size=(4, 4))
+    pf.plot_activations(fig_size=(4, 4))
     pf.sizes_vs_mean_speed(plot=True)
     pf.slow_and_fast_sizes()
-    # pf.density_vs_mean_speed()
-    # pf.field_peak_shifts(plot=True)  # only for NetworkIntDriven and NetworkExtDriven
+    pf.density_vs_mean_speed()
+    pf.field_peak_shifts(plot=True)  # only for NetworkIntDriven and NetworkExtDriven
     # pf.size_vs_induction_speed()  # only for NetworkIndep
     del pf
 
@@ -28,10 +28,10 @@ def analyze(identifier):
     pp.fast_and_slow_slopes()
     del pp
 
-    # sweeps = ThetaSweeps.current_instance(config)
-    # sweeps.length_vs_mean_speed(plot=True)
-    # sweeps.ahead_and_behind_vs_mean_speed(plot=True)
-    # sweeps.behind_length_vs_peak_shift(plot=True)    # only for NetworkIntDriven and NetworkExtDriven
+    sweeps = ThetaSweeps.current_instance(config)
+    sweeps.length_vs_mean_speed(plot=True)
+    sweeps.ahead_and_behind_vs_mean_speed(plot=True)
+    sweeps.behind_length_vs_peak_shift(plot=True)    # only for NetworkIntDriven and NetworkExtDriven
 
 
 # RUN ANALYSES

batch_config.py

@@ -8,7 +8,7 @@
 # choose network class in NetworkClass.py
 
 variants = {'NetworkIntDriven': 'IntDrivenLog80', 'NetworkExtDriven': 'ExtDrivenLog100', 'NetworkIndep': 'IndepLog80'}
-variants |= {'LinearTrack': 'ManyLaps'}  # for fast / slow analysis with top and bottom 20% percentiles of speed
+variants |= {'LinearTrack': 'ManyLaps'}  # better for fast/slow analysis with top and bottom 20% percentiles of speed
 # variants |= {'PlaceFields': 'HighThreshold'}  # comment this out only when using NetworkIndep
 
 

batch_plots.py

@@ -151,6 +151,8 @@ def plot_speed_ratios(name, class_def: Type[SmartSim], x_label, rel_path_y, bin_
     ax.set_ylabel("Count")
     ax.set_xlabel(x_label)
 
+    ax.set_xlim((0.65, 1.35))
+
     ax.xaxis.set_major_locator(ticker.MultipleLocator(0.2))
 
     ax.spines.right.set_visible(False)
@@ -170,29 +172,29 @@ def plot_speed_ratios(name, class_def: Type[SmartSim], x_label, rel_path_y, bin_
      z_binned_average=True, z_bin_size=2, fig_size=(4.94*CM, 5*CM))
 plot_speed_ratios("size_increments", PlaceFields, "Fast / slow", "slow_and_fast_sizes",
                   bin_size=0.05, fig_size=(3.75*CM, 2.3*CM))
-# plot("densities", PlaceFields, x_label, "density/speeds", "Place field density (peaks/cm)", "density/densities",
-#      alpha=0.5, x_binned_average=True, x_bin_size=2, fig_size=(3.73*CM, 4.9*CM))
+plot("densities", PlaceFields, x_label, "density/speeds", "Place field density (peaks/cm)", "density/densities",
+     alpha=0.5, x_binned_average=True, x_bin_size=2, fig_size=(3.73*CM, 4.9*CM))
 
 plot("slopes", PhasePrecession, x_label, "speeds", "1 / phase precession slope (cm/deg)", "slopes",
      "Position (cm)", "positions", z_binned_average=True, z_bin_size=2, fig_size=(5.25*CM, 5*CM), multiple_locator=0.05)
 plot_speed_ratios("slope_increments", PhasePrecession, "Fast / slow", "slow_and_fast_slopes", fig_size=(3.75*CM, 2.3*CM))
 
-# plot("sweep lengths", ThetaSweeps, x_label, "speeds", "Theta sweep length (cm)", "lengths", "Position (cm)", "positions",
-#      z_binned_average=True, z_bin_size=2, s=5, format='png', fig_size=small_fig)
-# plot("ahead lengths", ThetaSweeps, x_label, "ahead_and_behind/ahead_speeds", "Look-ahead distance (cm)",
-#      "ahead_and_behind/ahead_lengths", "Position (cm)", "ahead_and_behind/ahead_real_pos",
-#      z_binned_average=True, z_bin_size=2, s=5, format='png', fig_size=small_fig)
-# plot("behind lengths", ThetaSweeps, x_label, "ahead_and_behind/behind_speeds", "Look-behind distance (cm)",
-#      "ahead_and_behind/behind_lengths", "Position (cm)", "ahead_and_behind/behind_real_pos",
-#      z_binned_average=True, z_bin_size=2, s=5, format='png', fig_size=small_fig)
-#
-# # only for NetworkIntDriven and NetworkExtDriven:
-# plot("shifts", PlaceFields, x_label, "shifts/speeds", "Place field shift (cm)", "shifts/shifts", "Position (cm)",
-#      "shifts/positions", z_binned_average=True, z_bin_size=2, fig_size=small_fig)
-# plot("matched shifts", ThetaSweeps, "Place field shift (cm)", "ahead_and_behind/shifts",
-#      "Look-behind distance (cm)", "ahead_and_behind/behind_lengths", "Position (cm)",
-#      "ahead_and_behind/behind_real_pos", z_binned_average=True, z_bin_size=2, s=5, format='png',
-#      extra_plotting=lambda ax: ax.plot((0, 20), (0, 20), linestyle='dashed', color='black'), fig_size=small_fig)
+plot("sweep lengths", ThetaSweeps, x_label, "speeds", "Theta sweep length (cm)", "lengths", "Position (cm)", "positions",
+     z_binned_average=True, z_bin_size=2, s=5, format='png', fig_size=small_fig)
+plot("ahead lengths", ThetaSweeps, x_label, "ahead_and_behind/ahead_speeds", "Look-ahead distance (cm)",
+     "ahead_and_behind/ahead_lengths", "Position (cm)", "ahead_and_behind/ahead_real_pos",
+     z_binned_average=True, z_bin_size=2, s=5, format='png', fig_size=small_fig)
+plot("behind lengths", ThetaSweeps, x_label, "ahead_and_behind/behind_speeds", "Look-behind distance (cm)",
+     "ahead_and_behind/behind_lengths", "Position (cm)", "ahead_and_behind/behind_real_pos",
+     z_binned_average=True, z_bin_size=2, s=5, format='png', fig_size=small_fig)
+
+# only for NetworkIntDriven and NetworkExtDriven:
+plot("shifts", PlaceFields, x_label, "shifts/speeds", "Place field shift (cm)", "shifts/shifts", "Position (cm)",
+     "shifts/positions", z_binned_average=True, z_bin_size=2, fig_size=small_fig)
+plot("matched shifts", ThetaSweeps, "Place field shift (cm)", "ahead_and_behind/shifts",
+     "Look-behind distance (cm)", "ahead_and_behind/behind_lengths", "Position (cm)",
+     "ahead_and_behind/behind_real_pos", z_binned_average=True, z_bin_size=2, s=5, format='png',
+     extra_plotting=lambda ax: ax.plot((0, 20), (0, 20), linestyle='dashed', color='black'), fig_size=small_fig)
 
 # # only for NetworkIndep:
 # plot("sizes_vs_induction_speed", PlaceFields, "Induction speed (cm/s)", "induction_speeds/induction_speeds",

parameters/PlaceFields.json

@@ -7,5 +7,5 @@
   "last_unit": 170,
   "dens_window_size": 10,
   "dens_window_stride": 2,
-  "speed_percentile": 25
+  "speed_percentile": 20
 }