add ipsirepr_inhib arg to L6_model() and separate main_sim.ipynb into…

… three different components
rythorpe · Oct 18, 2024 · 02f9d28 · 02f9d28
1 parent d35de24
commit 02f9d28
Show file tree

Hide file tree

Showing 6 changed files with 13,838 additions and 9,072 deletions.
diff --git a/hnn_core/network_models.py b/hnn_core/network_models.py
@@ -308,7 +308,7 @@ def calcium_model(params=None, add_drives_from_params=False,
 
 def L6_model(params=None, add_drives_from_params=False,
              legacy_mode=False, layer_6_fb=True, grid_shape=(12, 12),
-             rng=None):
+             rng=None, ipsirepr_inhib=1.0):
     """Instantiate the updated calcium model with layer 6 cell types.
 
     Returns
@@ -536,7 +536,7 @@ def L6_model(params=None, add_drives_from_params=False,
                 prob_e_e_6 = prob_e_e
                 prob_i_e_6 = prob_i_e
                 prob_e_i_6 = prob_e_i
-                prob_i_e_cross = 1.0
+                prob_i_e_cross = ipsirepr_inhib
             else:
                 # between-group connection probabilities
                 prob_e_e = 0.00
@@ -548,7 +548,7 @@ def L6_model(params=None, add_drives_from_params=False,
                 prob_e_e_6 = prob_e_e
                 prob_i_e_6 = prob_i_e
                 prob_e_i_6 = prob_e_i
-                prob_i_e_cross = 0.0
+                prob_i_e_cross = (1 - ipsirepr_inhib)
 
             # layer2 Pyr -> layer2 Pyr
             for receptor in ['nmda', 'ampa']:

diff --git a/rs_dd_project/main_sim.ipynb b/rs_dd_project/main_sim.ipynb
diff --git a/rs_dd_project/main_sim_ipsirepr_inhib_0.5.ipynb b/rs_dd_project/main_sim_ipsirepr_inhib_0.5.ipynb
diff --git a/rs_dd_project/main_sim_ipsirepr_inhib_0.ipynb b/rs_dd_project/main_sim_ipsirepr_inhib_0.ipynb
diff --git a/rs_dd_project/main_sim_ipsirepr_inhib_1.ipynb b/rs_dd_project/main_sim_ipsirepr_inhib_1.ipynb
diff --git a/rs_dd_project/plot_simulate_rep_L6.py b/rs_dd_project/plot_simulate_rep_L6.py
@@ -31,8 +31,9 @@
 # emp_dpl = read_dipole('S1_SupraT.txt')
 
 
-def sim_dev_spiking(dev_magnitude=-1, reps=4, n_trials=1, burn_in_time=300.0,
-                    n_procs=10, record_vsec=False, rng=None):
+def sim_dev_spiking(dev_magnitude=-1, reps=4, ipsirepr_inhib=1.0, n_trials=1,
+                    burn_in_time=300.0, n_procs=10, record_vsec=False,
+                    rng=None):
 
     # Hyperparameters of repetitive drive sequence
     stim_interval = 100.  # in ms; 10 Hz
@@ -80,7 +81,8 @@ def sim_dev_spiking(dev_magnitude=-1, reps=4, n_trials=1, burn_in_time=300.0,
 
     ###########################################################################
     # Let us first create our network
-    net = L6_model(layer_6_fb=True, rng=rng, grid_shape=grid_shape)
+    net = L6_model(layer_6_fb=True, rng=rng, grid_shape=grid_shape,
+                   ipsirepr_inhib=ipsirepr_inhib)
     net.set_cell_positions(inplane_distance=300.0)
 
     # before we continue, sample a random integer to serve as seed for the