Merge pull request #415 from Neuroblox/MSN_G_M_parameter

change default parameters in Adam's model and allow changing G_M
Neuroblox · Sep 26, 2024 · b7637b8 · b7637b8
2 parents e24ecab + 7d2f66a
commit b7637b8
Show file tree

Hide file tree

Showing 2 changed files with 22 additions and 19 deletions.
diff --git a/src/blox/DBS_Model_Blox_Adam_Brown.jl b/src/blox/DBS_Model_Blox_Adam_Brown.jl
@@ -64,13 +64,14 @@ struct Striatum_MSN_Adam <: CompositeBlox
         namespace = nothing,
         N_inhib = 100,
         E_syn_inhib=-80,
-        I_bg=1.19*ones(N_inhib),
+        I_bg=1.172*ones(N_inhib),
         freq=zeros(N_inhib),
         phase=zeros(N_inhib),
         τ_inhib=13,
-        σ=4.0,
+        σ=0.11,
         density=0.3,
         weight=0.1,
+        G_M=1.3,
         connection_matrix=nothing
     )
         n_inh = [
@@ -82,7 +83,8 @@ struct Striatum_MSN_Adam <: CompositeBlox
                     I_bg = I_bg[i],
                     freq = freq[i],
                     phase = phase[i],
-                    σ=σ
+                    σ=σ,
+                    G_M=G_M
             ) 
             for i in Base.OneTo(N_inhib)
         ]
@@ -132,12 +134,12 @@ struct Striatum_FSI_Adam  <: CompositeBlox
         namespace = nothing,
         N_inhib = 50,
         E_syn_inhib=-80,
-        I_bg=5.5*ones(N_inhib),
+        I_bg=6.2*ones(N_inhib),
         freq=zeros(N_inhib),
         phase=zeros(N_inhib),
         τ_inhib=11,
         τ_inhib_s=6.5,
-        σ=6.0,
+        σ=1.2,
         density=0.58,
         g_density=0.33,
         weight=0.6,
@@ -215,11 +217,11 @@ struct GPe_Adam <: CompositeBlox
         namespace = nothing,
         N_inhib = 80,
         E_syn_inhib=-80,
-        I_bg=3*ones(N_inhib),
+        I_bg=3.4*ones(N_inhib),
         freq=zeros(N_inhib),
         phase=zeros(N_inhib),
         τ_inhib=10,
-        σ=8.0,
+        σ=1.7,
         density=0.0,
         weight=0.0,
         connection_matrix=nothing
@@ -284,11 +286,11 @@ struct STN_Adam <: CompositeBlox
         namespace = nothing,
         N_exci = 40,
         E_syn_exci=0.0,
-        I_bg=1.9*ones(N_exci),
+        I_bg=1.8*ones(N_exci),
         freq=zeros(N_exci),
         phase=zeros(N_exci),
         τ_exci=2,
-        σ=8.0,
+        σ=1.7,
         density=0.0,
         weight=0.0,
         connection_matrix=nothing

diff --git a/src/blox/neuron_models.jl b/src/blox/neuron_models.jl
@@ -175,15 +175,16 @@ struct HHNeuronInhib_MSN_Adam_Blox <: AbstractInhNeuronBlox
         name, 
         namespace=nothing,
         E_syn=-80.0, 
-        I_bg=1.19,
+        I_bg=1.172,
         freq=0,
         phase=0,
         τ=13,
         Cₘ=1.0,
-		σ=4.0,
+		σ=0.11,
 		a=2,
 		b=4,
-		T=37
+		T=37,
+		G_M=1.3
     )
 		sts = @variables begin 
 			V(t)=-63.83 
@@ -206,7 +207,7 @@ struct HHNeuronInhib_MSN_Adam_Blox <: AbstractInhNeuronBlox
 			G_Na = 100 
 			G_K  = 80 
 			G_L = 0.1 
-			G_M = 1.3
+			G_M = G_M
 			E_Na = 50 
 			E_K = -100 
 			E_L = -67 
@@ -262,13 +263,13 @@ struct HHNeuronInhib_FSI_Adam_Blox <: AbstractInhNeuronBlox
         name, 
         namespace=nothing,
         E_syn=-80.0, 
-        I_bg=5.5,
+        I_bg=6.2,
         freq=0,
         phase=0,
         τ=11,
 		τₛ=6.5,
         Cₘ=1.0,
-		σ=60.0,
+		σ=1.2,
 		a=4,
 		b=10,
 		T=37
@@ -354,12 +355,12 @@ struct HHNeuronExci_STN_Adam_Blox <: AbstractExciNeuronBlox
         name, 
         namespace=nothing,
         E_syn=0.0, 
-        I_bg=1.9,
+        I_bg=1.8,
         freq=0,
         phase=0,
         τ=2,
         Cₘ=1.0,
-		σ=80.0,
+		σ=1.7,
 		a=5,
 		b=4
     )
@@ -432,12 +433,12 @@ struct HHNeuronInhib_GPe_Adam_Blox <: AbstractInhNeuronBlox
         name, 
         namespace=nothing,
         E_syn=-80.0, 
-        I_bg=3.0,
+        I_bg=3.4,
         freq=0,
         phase=0,
         τ=10,
         Cₘ=1.0,
-		σ=80.0,
+		σ=1.7,
 		a=2,
 		b=4,
 		T=37