def get():
l = []
solution, s_mul, s_equal = get_solution_slow_GP_striatum_2()
for y in numpy.arange(1.0, -0.1, -0.25):
for z in numpy.arange(1.0, -0.1, -0.25):
x = 2.5
d = {}
for keys in s_mul:
update(solution, d, keys)
misc.dict_update(d, {'nest': {'GI_GA_gaba': {'weight': y}}})
misc.dict_update(d, {'nest': {'ST_GA_ampa': {'weight': z}}})
l += [pl(d, '*', **{'name': ''})]
d = {}
for keys in s_equal:
update(solution, d, keys)
l[-1] += pl(d, '=',
**{'name': 'GIGA_' + str(y) + '_STGA_' + str(z)})
return l
def get():
l=[]
solution, s_mul, s_equal=get_solution()
d0=0.8
f_beta_rm=lambda f: (1-f)/(d0+f*(1-d0))
for nodes in [['M1', 'M2'],
['M1', 'M2', 'FS'],
['M1', 'M2', 'FS','ST', 'GA', 'GI','SN']]:
x=2.5
d={}
for keys in s_mul: update(solution, d, keys)
l+=[pl(d, '*', **{'name':''})]
d={}
for keys in s_equal: update(solution, d, keys)
misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}})
misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}})
d['node']['EA']['rate']*=0.7
for node in nodes:
misc.dict_update(d,{'node':{node:{'edge_wrap':False}}})
l[-1]+=pl(d, '=', **{'name':'mod_ST_beta_'+str(x)})
return l
def get():
l=[]
solution, s_mul, s_equal=get_solution()
# Decrease/increase C1,C2,CF,CS, EI and EA. Test idea that at
# slow wave cortical and thalamic input is decreased.
mod=numpy.arange(0.5, 1.5, 0.05)
for x in mod:
d={}
for keys in s_mul: update(solution, d, keys)
misc.dict_update(d, {'node':{'CF':{'rate':x}}} )
misc.dict_update(d, {'node':{'CS':{'rate':x}}} )
l+=[pl(d, '*', **{'name':''})]
d={}
for keys in s_equal: update(solution, d, keys)
d['node']['C1']['rate']*=x
d['node']['C2']['rate']*=x
d['node']['EI']['rate']*=x
d['node']['EA']['rate']*=x
l[-1]+=pl(d, '=', **{'name':'mod_C1_C2_CF_CS_EI_EA_'+str(x)})
return l
def get():
l = []
solution, s_mul, s_equal = get_solution()
d0 = 0.8
f_beta_rm = lambda f: (1 - f) / (d0 + f * (1 - d0))
v0 = numpy.arange(1.25, 3.75, 0.5)
v1 = numpy.arange(2.5, 5., 0.5)
v = [[x, y] for x in v0 for y in v1]
for x, y in v:
d = {}
for keys in s_mul:
update(solution, d, keys)
l += [pl(d, '*', **{'name': ''})]
d = {}
for keys in s_equal:
update(solution, d, keys)
misc.dict_update(d, {'nest': {'ST': {'beta_I_GABAA_1': f_beta_rm(x)}}})
misc.dict_update(d, {'nest': {'ST': {'beta_I_AMPA_1': f_beta_rm(y)}}})
misc.dict_update(d, {'nest': {'ST': {'beta_I_NMDA_1': f_beta_rm(y)}}})
d['node']['EA']['rate'] *= 0.7
l[-1] += pl(d, '=',
**{'name': ('mod_ST_beta_' + str(x) + '_' + str(y))})
return l
def get():
l = []
solution, s_mul, s_equal = get_solution_slow_GP_striatum()
d0 = 0.8
f_beta_rm = lambda f: (1 - f) / (d0 + f * (1 - d0))
v = STN_ampa_gaba_input_magnitude()
for _, x, y in v:
# x=2.5-1
d = {}
for keys in s_mul:
update(solution, d, keys)
l += [pl(d, '*', **{'name': ''})]
d = {}
for keys in s_equal:
update(solution, d, keys)
# misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}})
# misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}})
# d['node']['EA']['rate']*=0.7
misc.dict_update(d, {'node': {'CS': {'rate': x}}})
misc.dict_update(d, {'nest': {'GI_ST_gaba': {'weight': y}}})
l[-1] += pl(d, '=', **{'name': 'mod_ST_inp_' + str(x) + '_' + str(y)})
return l
def get():
l=[]
solution, s_mul, s_equal=get_solution()
d0=0.8
f_beta_rm=lambda f: (1-f)/(d0+f*(1-d0))
x=2.5-1.25
d={}
for keys in s_mul: update(solution, d, keys)
l+=[pl(d, '*', **{'name':''})]
d={}
for keys in s_equal: update(solution, d, keys)
misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}})
misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}})
d['node']['EA']['rate']*=0.7
misc.dict_update(d,{'node':{'CS':{'rate': 830.}}})
misc.dict_update(d,{'nest':{'GI_ST_gaba':{'weight':0.35 }}})
l[-1]+=pl(d, '=', **{'name':'mod_ST_beta_'+str(x)})
return l
def get():
l=[]
solution, s_mul, s_equal=get_solution_slow_GP_striatum()
d0=0.8
f_beta_rm=lambda f: (1-f)/(d0+f*(1-d0))
x=2.5
d={}
for keys in s_mul: update(solution, d, keys)
l+=[pl(d, '*', **{'name':''})]
d={}
for keys in s_equal: update(solution, d, keys)
misc.dict_update(d,{'nest':{'GI':{'beta_I_GABAA_1': f_beta_rm(2.)}}})
misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}})
misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}})
d['node']['EA']['rate']*=0.7
l[-1]+=pl(d, '=', **{'name':'mod_ST_beta_'+str(x)})
return l
def get():
l = []
solution, s_mul, s_equal = get_solution()
d0 = 0.8
f_beta_rm = lambda f: (1 - f) / (d0 + f * (1 - d0))
x = 2.5 - 1.25
d = {}
for keys in s_mul:
update(solution, d, keys)
l += [pl(d, '*', **{'name': ''})]
d = {}
for keys in s_equal:
update(solution, d, keys)
misc.dict_update(d, {'nest': {'ST': {'beta_I_AMPA_1': f_beta_rm(x)}}})
misc.dict_update(d, {'nest': {'ST': {'beta_I_NMDA_1': f_beta_rm(x)}}})
d['node']['EA']['rate'] *= 0.7
misc.dict_update(d, {'node': {'CS': {'rate': 830.}}})
misc.dict_update(d, {'nest': {'GI_ST_gaba': {'weight': 0.35}}})
l[-1] += pl(d, '=', **{'name': 'mod_ST_beta_' + str(x)})
return l
def get():
l=[]
solution, s_mul, s_equal=get_solution()
# d0=0.8
# f_beta_rm=lambda f: (1-f)/(d0+f*(1-d0))
# x=2.5
for x in numpy.arange(0.7, 1.5 ,0.2):
d={}
for keys in s_mul: update(solution, d, keys)
l+=[pl(d, '*', **{'name':''})]
d={}
for keys in s_equal: update(solution, d, keys)
# misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}})
# misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}})
d['node']['EA']['rate']*=x
l[-1]+=pl(d, '=', **{'name':'mod_EA_'+str(x)})
return l
def get():
l = []
solution, s_mul, s_equal = get_solution_slow_GP_striatum()
d0 = 0.8
f_beta_rm = lambda f: (1 - f) / (d0 + f * (1 - d0))
for y in numpy.arange(0.9, 0.0, -0.2):
x = 2.5
d = {}
for keys in s_mul:
update(solution, d, keys)
misc.dict_update(d, {'conn': {'nest': {'GI_GA_gaba': {'weight': y}}}})
l += [pl(d, '*', **{'name': ''})]
d = {}
for keys in s_equal:
update(solution, d, keys)
misc.dict_update(d, {'nest': {'ST': {'beta_I_AMPA_1': f_beta_rm(x)}}})
misc.dict_update(d, {'nest': {'ST': {'beta_I_NMDA_1': f_beta_rm(x)}}})
d['node']['EA']['rate'] *= 0.7
l[-1] += pl(d, '=', **{'name': 'mod_ST_beta_' + str(x)})
return l
def get():
l=[]
solution, s_mul, s_equal=get_solution()
d0=0.8
f_beta_rm=lambda f: (1-f)/(d0+f*(1-d0))
v0=numpy.arange(1.25,3.75,0.5)
v1=numpy.arange(2.5, 5.,0.5)
v=[[x,y] for x in v0 for y in v1]
for x, y in v:
d={}
for keys in s_mul: update(solution, d, keys)
l+=[pl(d, '*', **{'name':''})]
d={}
for keys in s_equal: update(solution, d, keys)
misc.dict_update(d,{'nest':{'ST':{'beta_I_GABAA_1': f_beta_rm(x)}}})
misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(y)}}})
misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(y)}}})
d['node']['EA']['rate']*=0.7
l[-1]+=pl(d, '=', **{'name':('mod_ST_beta_'
+str(x)+'_'+str(y))
})
return l
def get():
l = []
solution, s_mul, s_equal = get_solution()
d0 = 0.8
f_beta_rm = lambda f: (1 - f) / (d0 + f * (1 - d0))
d = {}
for x in numpy.arange(3.4, 6, 0.4):
# x=3.4
y = 2.25
for keys in s_mul:
update(solution, d, keys)
l += [pl(d, '*', **{'name': ''})]
d = {}
for keys in s_equal:
update(solution, d, keys)
print f_beta_rm(x)
misc.dict_update(d, {'nest': {'ST': {'beta_I_AMPA_1': f_beta_rm(x)}}})
misc.dict_update(d, {'nest': {'ST': {'beta_I_NMDA_1': f_beta_rm(x)}}})
misc.dict_update(d, {'nest': {'ST': {'beta_I_GABAA_1': f_beta_rm(y)}}})
# d['node']['EA']['rate']*=0.7
l[-1] += pl(
d, '=',
**{'name': 'mod_CS_STbeta_' + str(x) + '_GA_STbeta_' + str(y)})
return l
def get():
l=[]
solution, s_mul, s_equal=get_solution_slow_GP_striatum()
d0=0.8
f_beta_rm=lambda f: (1-f)/(d0+f*(1-d0))
x=2.5
d={}
for keys in s_mul: update(solution, d, keys)
l+=[pl(d, '*', **{'name':''})]
d={}
for keys in s_equal: update(solution, d, keys)
misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}})
misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}})
d['node']['EA']['rate']*=0.7
l[-1]+=pl(d, '=', **{'name':'mod_ST_beta_'+str(x)})
return l
def get():
l = []
solution, s_mul, s_equal = get_solution()
# Decrease/increase C1,C2,CF,CS, EI and EA. Test idea that at
# slow wave cortical and thalamic input is decreased.
mod = numpy.arange(0.5, 1.5, 0.05)
for x in mod:
d = {}
for keys in s_mul:
update(solution, d, keys)
misc.dict_update(d, {'node': {'CF': {'rate': x}}})
misc.dict_update(d, {'node': {'CS': {'rate': x}}})
l += [pl(d, '*', **{'name': ''})]
d = {}
for keys in s_equal:
update(solution, d, keys)
d['node']['C1']['rate'] *= x
d['node']['C2']['rate'] *= x
d['node']['EI']['rate'] *= x
d['node']['EA']['rate'] *= x
l[-1] += pl(d, '=', **{'name': 'mod_C1_C2_CF_CS_EI_EA_' + str(x)})
return l
def get():
l=[]
solution, s_mul, s_equal=get_solution_slow_GP_striatum()
d0=0.8
f_beta_rm=lambda f: (1-f)/(d0+f*(1-d0))
v=STN_ampa_gaba_input_magnitude()
for _, x, y in v:
# x=2.5-1
d={}
for keys in s_mul: update(solution, d, keys)
l+=[pl(d, '*', **{'name':''})]
d={}
for keys in s_equal: update(solution, d, keys)
# misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}})
# misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}})
# d['node']['EA']['rate']*=0.7
misc.dict_update(d,{'node':{'CS':{'rate': x}}})
misc.dict_update(d,{'nest':{'GI_ST_gaba':{'weight':y }}})
l[-1]+=pl(d, '=', **{'name':'mod_ST_inp_'+str(x)+'_'+str(y)})
return l
def get():
l=[]
solution, s_mul, s_equal=get_solution_slow_GP_striatum()
d0=0.8
f_beta_rm=lambda f: (1-f)/(d0+f*(1-d0))
d={}
for x in numpy.arange(3.4, 6, 0.4):
# x=3.4
y=2.25
for keys in s_mul: update(solution, d, keys)
l+=[pl(d, '*', **{'name':''})]
d={}
for keys in s_equal: update(solution, d, keys)
misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}})
misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}})
misc.dict_update(d,{'nest':{'ST':{'beta_I_GABAA_1': f_beta_rm(y)}}})
d['node']['EA']['rate']*=0.7
l[-1]+=pl(d, '=', **{'name':'mod_CS_STbeta_'+str(x)+'_GA_STbeta_'+str(y)})
return l
def get():
l=[]
solution, s_mul, s_equal=get_solution()
d0=0.8
f_beta_rm=lambda f: (1-f)/(d0+f*(1-d0))
y=2.5
v=numpy.arange(0.1,1.2,0.1)
for x in v:
d={}
for keys in s_mul: update(solution, d, keys)
l+=[pl(d, '*', **{'name':''})]
d={}
for keys in s_equal: update(solution, d, keys)
misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(y)}}})
misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(y)}}})
misc.dict_update(d,{'nest':{'GA_M1_gaba':{'weight': x}}})
misc.dict_update(d,{'nest':{'GA_M2_gaba':{'weight': x}}})
misc.dict_update(d,{'nest':{'GA_FS_gaba':{'weight': x}}})
d['node']['EA']['rate']*=0.7
l[-1]+=pl(d, '=', **{'name':('mod_GA_striatum_1' + str(x))
})
return l
def get():
l = []
solution, s_mul, s_equal = get_solution()
# d0=0.8
# f_beta_rm=lambda f: (1-f)/(d0+f*(1-d0))
# x=2.5
for x in numpy.arange(0.7, 1.5, 0.2):
d = {}
for keys in s_mul:
update(solution, d, keys)
l += [pl(d, '*', **{'name': ''})]
d = {}
for keys in s_equal:
update(solution, d, keys)
# misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}})
# misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}})
d['node']['EA']['rate'] *= x
l[-1] += pl(d, '=', **{'name': 'mod_EA_' + str(x)})
return l
def get():
l = []
solution, s_mul, s_equal = get_solution()
d0 = 0.8
f_beta_rm = lambda f: (1 - f) / (d0 + f * (1 - d0))
y = 2.5
v = numpy.arange(0.1, 1.2, 0.1)
for x in v:
d = {}
for keys in s_mul:
update(solution, d, keys)
l += [pl(d, '*', **{'name': ''})]
d = {}
for keys in s_equal:
update(solution, d, keys)
misc.dict_update(d, {'nest': {'ST': {'beta_I_AMPA_1': f_beta_rm(y)}}})
misc.dict_update(d, {'nest': {'ST': {'beta_I_NMDA_1': f_beta_rm(y)}}})
misc.dict_update(d, {'nest': {'GA_M1_gaba': {'weight': x}}})
misc.dict_update(d, {'nest': {'GA_M2_gaba': {'weight': x}}})
misc.dict_update(d, {'nest': {'GA_FS_gaba': {'weight': x}}})
d['node']['EA']['rate'] *= 0.7
l[-1] += pl(d, '=', **{'name': ('mod_GA_striatum_1' + str(x))})
return l
def get():
l = []
solution, s_mul, s_equal = get_solution()
d0 = 0.8
f_beta_rm = lambda f: (1 - f) / (d0 + f * (1 - d0))
# v=STN_ampa_gaba_input_magnitude()[0:5]
for _, u, w in v:
v0 = numpy.arange(1.25, 3.75, 1.)
v1 = numpy.arange(1.5, 4, 1.)
v = [[x, y] for x in v0 for y in v1]
for x, y in v:
d = {}
for keys in s_mul:
update(solution, d, keys)
l += [pl(d, '*', **{'name': ''})]
d = {}
for keys in s_equal:
update(solution, d, keys)
misc.dict_update(
d, {'nest': {
'ST': {
'beta_I_GABAA_1': f_beta_rm(x)
}
}})
misc.dict_update(d,
{'nest': {
'ST': {
'beta_I_AMPA_1': f_beta_rm(y)
}
}})
misc.dict_update(d,
{'nest': {
'ST': {
'beta_I_NMDA_1': f_beta_rm(y)
}
}})
misc.dict_update(d, {'node': {'CS': {'rate': u}}})
misc.dict_update(d, {'nest': {'GI_ST_gaba': {'weight': w}}})
d['node']['EA']['rate'] *= 0.7
l[-1] += pl(
d, '=', **{
'name': ('mod_ST_beta_' + str(u) + '_' + str(w) + '_' +
str(x) + '_' + str(y))
})
return l
def get():
l = []
solution, s_mul, s_equal = get_solution()
d0 = 0.8
f_beta_rm = lambda f: (1 - f) / (d0 + f * (1 - d0))
# betas MS-MS
v = [0.1, 0.25, 0.5, 0.75]
for x, y in [[v1, v2] for v1 in v for v2 in v]:
d = {}
for keys in s_mul:
update(solution, d, keys)
l += [pl(d, '*', **{'name': ''})]
d = {}
for keys in s_equal:
update(solution, d, keys)
misc.dict_update(d,
{'nest': {
'ST': {
'beta_I_AMPA_1': f_beta_rm(2.5)
}
}})
misc.dict_update(d,
{'nest': {
'ST': {
'beta_I_NMDA_1': f_beta_rm(2.5)
}
}})
d['node']['EA']['rate'] *= 0.7
misc.dict_update(d, {'nest': {'MS': {'beta_I_GABAA_2': f_beta_rm(x)}}})
for conn in ['M1_M1_gaba', 'M1_M2_gaba', 'M2_M1_gaba', 'M2_M2_gaba']:
misc.dict_update(d,
{'conn': {
conn: {
'beta_fan_in': f_beta_rm(y)
}
}})
l[-1] += pl(d, '=', **{'name': 'mod_MSg_' + str(x) + '_MSc_' + str(y)})
#beta CTX-STN
return l
def get():
l=[]
solution, s_mul, s_equal=get_solution_slow_GP_striatum()
d0=0.8
f_beta_rm=lambda f: (1-f)/(d0+f*(1-d0))
x=2.5
d={}
for y in [12]+ list(numpy.arange(5,55,5)):
# x=3.4
# y=2.25
for keys in s_mul: update(solution, d, keys)
l+=[pl(d, '*', **{'name':''})]
d={}
for keys in s_equal: update(solution, d, keys)
ratio=12./y
dd={'nest':{'GA_M1_gaba':{'weight':0.8*ratio},
'GA_M2_gaba':{'weight':0.8*ratio}}}
misc.dict_update(d,dd)
# Decreasing from 2 leads to ...
# Increasing from 2 leads to ...
dd={'nest':{'GA_FS_gaba':{'weight':2.*ratio}}}
misc.dict_update(d,dd)
# Just assumed to be 12 ms
dd={'nest':{'M1_low':{'GABAA_3_Tau_decay':12./ratio},
'M2_low':{'GABAA_3_Tau_decay':12./ratio},
'FS_low':{'GABAA_2_Tau_decay':12./ratio},
}}
misc.dict_update(d,dd)
misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}})
misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}})
d['node']['EA']['rate']*=0.7
l[-1]+=pl(d, '=', **{'name':'mod_GAtau_'+str(y)})
return l
def get():
l = []
solution, s_mul, s_equal = get_solution()
d0 = 0.8
f_beta_rm = lambda f: (1 - f) / (d0 + f * (1 - d0))
y = 2.5
w = numpy.arange(0.5, 1.2, 0.1)
v = [0.1, 0.25, 0.5]
ll = [[v1, v2] for v1 in v for v2 in w]
pp(l)
for x, z in ll:
print x, z
d = {}
for keys in s_mul:
update(solution, d, keys)
l += [pl(d, '*', **{'name': ''})]
d = {}
for keys in s_equal:
update(solution, d, keys)
misc.dict_update(d, {'nest': {'ST': {'beta_I_AMPA_1': f_beta_rm(y)}}})
misc.dict_update(d, {'nest': {'ST': {'beta_I_NMDA_1': f_beta_rm(y)}}})
misc.dict_update(d, {'nest': {'MS': {'beta_I_GABAA_2': f_beta_rm(x)}}})
for conn in ['M1_M1_gaba', 'M1_M2_gaba', 'M2_M1_gaba', 'M2_M2_gaba']:
misc.dict_update(d,
{'conn': {
conn: {
'beta_fan_in': f_beta_rm(x)
}
}})
misc.dict_update(d, {'nest': {'GA_M1_gaba': {'weight': z}}})
misc.dict_update(d, {'nest': {'GA_M2_gaba': {'weight': z}}})
misc.dict_update(d, {'nest': {'GA_FS_gaba': {'weight': z}}})
d['node']['EA']['rate'] *= 0.7
l[-1] += pl(d, '=',
**{'name': ('mod_GA_str_' + str(z) + '_MS_' + str(x))})
return l
def get():
l = []
solution, s_mul, s_equal = get_solution_slow_GP_striatum_2()
d = {}
for keys in s_mul:
update(solution, d, keys)
l += [pl(d, '*', **{'name': ''})]
d = {}
for keys in s_equal:
update(solution, d, keys)
l[-1] += pl(d, '=', **{'name': ''})
return l
def get():
l=[]
solution, s_mul, s_equal=get_solution_slow_GP_striatum_2()
d={}
for keys in s_mul: update(solution, d, keys)
l+=[pl(d, '*', **{'name':''})]
d={}
for keys in s_equal: update(solution, d, keys)
l[-1]+=pl(d, '=', **{'name':''})
return l
def get():
l=[]
solution, s_mul, s_equal=get_solution()
d0=0.8
f_beta_rm=lambda f: (1-f)/(d0+f*(1-d0))
y=2.5
w=numpy.arange(0.5, 1.2, 0.1)
v=[0.1,0.25, 0.5]
ll=[[v1,v2] for v1 in v for v2 in w]
pp(l)
for x, z in ll:
print x,z
d={}
for keys in s_mul: update(solution, d, keys)
l+=[pl(d, '*', **{'name':''})]
d={}
for keys in s_equal: update(solution, d, keys)
misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(y)}}})
misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(y)}}})
misc.dict_update(d,{'nest':{'MS':{'beta_I_GABAA_2': f_beta_rm(x)}}})
for conn in ['M1_M1_gaba', 'M1_M2_gaba','M2_M1_gaba','M2_M2_gaba']:
misc.dict_update(d,{'conn':{conn:{'beta_fan_in': f_beta_rm(x)}}})
misc.dict_update(d,{'nest':{'GA_M1_gaba':{'weight': z}}})
misc.dict_update(d,{'nest':{'GA_M2_gaba':{'weight': z}}})
misc.dict_update(d,{'nest':{'GA_FS_gaba':{'weight': z}}})
d['node']['EA']['rate']*=0.7
l[-1]+=pl(d, '=', **{'name':('mod_GA_str_' + str(z)+'_MS_'+str(x))
})
return l
def get():
l=[]
solution, s_mul, s_equal=get_solution()
d0=0.8
f_beta_rm=lambda f: (1-f)/(d0+f*(1-d0))
# v=STN_ampa_gaba_input_magnitude()[0:5]
for _, u, w in v:
v0=numpy.arange(1.25, 3.75, 1.)
v1=numpy.arange(1.5, 4, 1.)
v=[[x,y] for x in v0 for y in v1]
for x, y in v:
d={}
for keys in s_mul: update(solution, d, keys)
l+=[pl(d, '*', **{'name':''})]
d={}
for keys in s_equal: update(solution, d, keys)
misc.dict_update(d,{'nest':{'ST':{'beta_I_GABAA_1': f_beta_rm(x)}}})
misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(y)}}})
misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(y)}}})
misc.dict_update(d,{'node':{'CS':{'rate': u}}})
misc.dict_update(d,{'nest':{'GI_ST_gaba':{'weight':w }}})
d['node']['EA']['rate']*=0.7
l[-1]+=pl(d, '=', **{'name':('mod_ST_beta_'
+str(u)+'_'+str(w)+'_'
+str(x)+'_'+str(y))
})
return l
def get():
l = []
solution, s_mul, s_equal = get_solution()
d0 = 0.8
f_beta_rm = lambda f: (1 - f) / (d0 + f * (1 - d0))
y = 2.5
v = [1. / 10., 3. / 10, 5. / 10, 7. / 10, 9. / 10]
for x in v:
d = {}
for keys in s_mul:
update(solution, d, keys)
misc.dict_update(d, {'conn': {'GA_M1_gaba': {'fan_in0': x}}})
misc.dict_update(d, {'conn': {'GA_M2_gaba': {'fan_in0': x}}})
misc.dict_update(d, {'nest': {'GA_M1_gaba': {'weight': 1. / x}}})
misc.dict_update(d, {'nest': {'GA_M2_gaba': {'weight': 1. / x}}})
l += [pl(d, '*', **{'name': ('mod_GA_MS_' + str(x))})]
d = {}
for keys in s_equal:
update(solution, d, keys)
misc.dict_update(d, {'nest': {'ST': {'beta_I_AMPA_1': f_beta_rm(y)}}})
misc.dict_update(d, {'nest': {'ST': {'beta_I_NMDA_1': f_beta_rm(y)}}})
d['node']['EA']['rate'] *= 0.7
l[-1] += pl(d, '=', **{'name': ''})
for x in v:
d = {}
for keys in s_mul:
update(solution, d, keys)
misc.dict_update(d, {'conn': {'GA_FS_gaba': {'fan_in0': x}}})
misc.dict_update(d, {'nest': {'GA_FS_gaba': {'weight': 1. / x}}})
l += [pl(d, '*', **{'name': ('mod_GA_FS_' + str(x))})]
d = {}
for keys in s_equal:
update(solution, d, keys)
misc.dict_update(d, {'nest': {'ST': {'beta_I_AMPA_1': f_beta_rm(y)}}})
misc.dict_update(d, {'nest': {'ST': {'beta_I_NMDA_1': f_beta_rm(y)}}})
d['node']['EA']['rate'] *= 0.7
l[-1] += pl(d, '=', **{'name': ''})
for x in v:
d = {}
for keys in s_mul:
update(solution, d, keys)
misc.dict_update(d, {'conn': {'GA_FS_gaba': {'fan_in0': x}}})
misc.dict_update(d, {'conn': {'GA_M1_gaba': {'fan_in0': x}}})
misc.dict_update(d, {'conn': {'GA_M2_gaba': {'fan_in0': x}}})
misc.dict_update(d, {'nest': {'GA_M1_gaba': {'weight': 1. / x}}})
misc.dict_update(d, {'nest': {'GA_M2_gaba': {'weight': 1. / x}}})
misc.dict_update(d, {'nest': {'GA_FS_gaba': {'weight': 1. / x}}})
l += [pl(d, '*', **{'name': ('mod_GA_str_' + str(x))})]
d = {}
for keys in s_equal:
update(solution, d, keys)
misc.dict_update(d, {'nest': {'ST': {'beta_I_AMPA_1': f_beta_rm(y)}}})
misc.dict_update(d, {'nest': {'ST': {'beta_I_NMDA_1': f_beta_rm(y)}}})
d['node']['EA']['rate'] *= 0.7
l[-1] += pl(d, '=', **{'name': ''})
return l
def get():
l = []
solution, s_mul, s_equal = get_solution()
d0 = 0.8
f_beta_rm = lambda f: (1 - f) / (d0 + f * (1 - d0))
# betas MS-MS
v = numpy.arange(1.5, 3, 0.25)
for x in v:
d = {}
for keys in s_mul:
update(solution, d, keys)
l += [pl(d, '*', **{'name': ''})]
d = {}
for keys in s_equal:
update(solution, d, keys)
misc.dict_update(d, {'nest': {'ST': {'beta_I_AMPA_1': f_beta_rm(x)}}})
misc.dict_update(d, {'nest': {'ST': {'beta_I_NMDA_1': f_beta_rm(x)}}})
l[-1] += pl(d, '=', **{'name': 'mod_ST_beta_' + str(x)})
v1 = numpy.arange(0.1, 1, 0.2)
for x in v1:
d = {}
for keys in s_mul:
update(solution, d, keys)
l += [pl(d, '*', **{'name': ''})]
d = {}
for keys in s_equal:
update(solution, d, keys)
d['node']['EA']['rate'] *= x
l[-1] += pl(d, '=', **{'name': 'mod_EA_' + str(x)})
v2 = numpy.arange(1.1, 2, 0.2)
for x in v2:
d = {}
for keys in s_mul:
update(solution, d, keys)
l += [pl(d, '*', **{'name': ''})]
d = {}
for keys in s_equal:
update(solution, d, keys)
d['node']['EI']['rate'] *= x
l[-1] += pl(d, '=', **{'name': 'mod_EI_' + str(x)})
return l
def get():
l=[]
solution, s_mul, s_equal=get_solution()
d0=0.8
f_beta_rm=lambda f: (1-f)/(d0+f*(1-d0))
d={}
x=2.5
for keys in s_mul: update(solution, d, keys)
l+=[pl(d, '*', **{'name':''})]
d={}
for keys in s_equal: update(solution, d, keys)
misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}})
misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}})
d['node']['EA']['rate']*=0.7
l[-1]+=pl(d, '=', **{'name':'mod_ST_beta_'+str(x)})
for y in numpy.arange(1,11,1):
d={}
for keys in s_mul: update(solution, d, keys)
l+=[pl(d, '*', **{'name':''})]
d={}
for keys in s_equal: update(solution, d, keys)
misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}})
misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}})
d['node']['EA']['rate']*=0.7
misc.dict_update(d,{'conn':{'GI_FS_gaba':{'lesion': False}}})
misc.dict_update(d,{'conn':{'GI_FS_gaba':{'fan_in0': y}}})
misc.dict_update(d,{'conn':{'GI_M1_gaba':{'lesion': False}}})
misc.dict_update(d,{'conn':{'GI_M1_gaba':{'fan_in0': y}}})
misc.dict_update(d,{'conn':{'GI_M2_gaba':{'lesion': False}}})
misc.dict_update(d,{'conn':{'GI_M2_gaba':{'fan_in0': y}}})
l[-1]+=pl(d, '=', **{'name':'mod_GI_M2_'+str(y)})
for y, EA_rate in zip(numpy.arange(25,125,25),
numpy.arange(0.9,1.3,0.1)):
d={}
for keys in s_mul: update(solution, d, keys)
l+=[pl(d, '*', **{'name':''})]
d={}
for keys in s_equal: update(solution, d, keys)
misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}})
misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}})
d['node']['EA']['rate']*=EA_rate
misc.dict_update(d,{'conn':{'M2_GA_gaba':{'lesion': False}}})
misc.dict_update(d,{'conn':{'M2_GA_gaba':{'fan_in0': y}}})
l[-1]+=pl(d, '=', **{'name':'M2GA_'+str(y)+'_EArate_'+str(EA_rate)})
return l
def get():
l=[]
solution, s_mul, s_equal=get_solution()
d0=0.8
f_beta_rm=lambda f: (1-f)/(d0+f*(1-d0))
# betas MS-MS
v=numpy.arange(1.5,3,0.25)
for x in v:
d={}
for keys in s_mul: update(solution, d, keys)
l+=[pl(d, '*', **{'name':''})]
d={}
for keys in s_equal: update(solution, d, keys)
misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}})
misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}})
l[-1]+=pl(d, '=', **{'name':'mod_ST_beta_'+str(x)})
v1=numpy.arange(0.1, 1, 0.2)
for x in v1:
d={}
for keys in s_mul: update(solution, d, keys)
l+=[pl(d, '*', **{'name':''})]
d={}
for keys in s_equal: update(solution, d, keys)
d['node']['EA']['rate']*=x
l[-1]+=pl(d, '=', **{'name':'mod_EA_'+str(x)})
v2=numpy.arange(1.1, 2, 0.2)
for x in v2:
d={}
for keys in s_mul: update(solution, d, keys)
l+=[pl(d, '*', **{'name':''})]
d={}
for keys in s_equal: update(solution, d, keys)
d['node']['EI']['rate']*=x
l[-1]+=pl(d, '=', **{'name':'mod_EI_'+str(x)})
return l
def get():
l = []
solution, s_mul, s_equal = get_solution_slow_GP_striatum()
d0 = 0.8
f_beta_rm = lambda f: (1 - f) / (d0 + f * (1 - d0))
for y in numpy.arange(1.0, -0.1, -0.5):
for z in numpy.arange(1.0, -0.1, -0.5):
x = 2.5
d = {}
for keys in s_mul:
update(solution, d, keys)
misc.dict_update(d, {'nest': {'GI_GA_gaba': {'weight': y}}})
misc.dict_update(d, {'nest': {'ST_GA_ampa': {'weight': z}}})
l += [pl(d, '*', **{'name': ''})]
d = {}
for keys in s_equal:
update(solution, d, keys)
misc.dict_update(d,
{'nest': {
'ST': {
'beta_I_AMPA_1': f_beta_rm(x)
}
}})
misc.dict_update(d,
{'nest': {
'ST': {
'beta_I_NMDA_1': f_beta_rm(x)
}
}})
d['node']['EA']['rate'] *= 0.7
l[-1] += pl(d, '=',
**{'name': 'GIGA_' + str(y) + '_STGA_' + str(z)})
for y in numpy.arange(1., 6.):
for z in numpy.arange(1., 6.):
x = 2.5
d = {}
for keys in s_mul:
update(solution, d, keys)
l += [pl(d, '*', **{'name': ''})]
d = {}
for keys in s_equal:
update(solution, d, keys)
misc.dict_update(d, {'nest': {'ST_GA_ampa': {'delay': y}}})
misc.dict_update(d, {'nest': {'ST_GI_ampa': {'delay': y}}})
misc.dict_update(d, {'nest': {'GI_ST_gaba': {'delay': z}}})
misc.dict_update(d,
{'nest': {
'ST': {
'beta_I_AMPA_1': f_beta_rm(x)
}
}})
misc.dict_update(d,
{'nest': {
'ST': {
'beta_I_NMDA_1': f_beta_rm(x)
}
}})
d['node']['EA']['rate'] *= 0.7
l[-1] += pl(
d, '=',
**{'name': 'STGPdelay_' + str(y) + '_GPSTdelay_' + str(z)})
for y in numpy.arange(2.5, 21., 2.5):
x = 2.5
d = {}
for keys in s_mul:
update(solution, d, keys)
l += [pl(d, '*', **{'name': ''})]
d = {}
for keys in s_equal:
update(solution, d, keys)
misc.dict_update(d, {'nest': {'CS_ST_ampa': {'delay': y}}})
misc.dict_update(d, {'nest': {'CS_ST_nmda': {'delay': y}}})
misc.dict_update(d, {'nest': {'ST': {'beta_I_AMPA_1': f_beta_rm(x)}}})
misc.dict_update(d, {'nest': {'ST': {'beta_I_NMDA_1': f_beta_rm(x)}}})
d['node']['EA']['rate'] *= 0.7
l[-1] += pl(d, '=', **{'name': 'CSSTdelay_' + str(y)})
for y in numpy.arange(2.5, 21., 2.5):
x = 2.5
d = {}
for keys in s_mul:
update(solution, d, keys)
l += [pl(d, '*', **{'name': ''})]
d = {}
for keys in s_equal:
update(solution, d, keys)
misc.dict_update(d, {'nest': {'C1_M1_ampa': {'delay': y}}})
misc.dict_update(d, {'nest': {'C1_M1_nmda': {'delay': y}}})
misc.dict_update(d, {'nest': {'C2_M2_ampa': {'delay': y}}})
misc.dict_update(d, {'nest': {'C2_M2_nmda': {'delay': y}}})
misc.dict_update(d, {'nest': {'CF_FS_ampa': {'delay': y}}})
misc.dict_update(d, {'nest': {'ST': {'beta_I_AMPA_1': f_beta_rm(x)}}})
misc.dict_update(d, {'nest': {'ST': {'beta_I_NMDA_1': f_beta_rm(x)}}})
d['node']['EA']['rate'] *= 0.7
l[-1] += pl(d, '=', **{'name': 'CSTRdelay_' + str(y)})
for y in numpy.arange(1, 8., 1.):
x = 2.5
d = {}
for keys in s_mul:
update(solution, d, keys)
l += [pl(d, '*', **{'name': ''})]
d = {}
for keys in s_equal:
update(solution, d, keys)
misc.dict_update(d, {'nest': {'GA_M1_gaba': {'delay': y}}})
misc.dict_update(d, {'nest': {'GA_M2_gaba': {'delay': y}}})
misc.dict_update(d, {'nest': {'GA_FS_gaba': {'delay': y}}})
d['node']['EA']['rate'] *= 0.7
l[-1] += pl(d, '=', **{'name': 'GASTRdelay_' + str(y)})
return l
def get():
l=[]
solution, s_mul, s_equal=get_solution_slow_GP_striatum()
d0=0.8
f_beta_rm=lambda f: (1-f)/(d0+f*(1-d0))
for y in numpy.arange(1.0,-0.1,-0.5):
for z in numpy.arange(1.0,-0.1,-0.5):
x=2.5
d={}
for keys in s_mul: update(solution, d, keys)
misc.dict_update(d,{'nest':{'GI_GA_gaba':{'weight':y}}})
misc.dict_update(d,{'nest':{'ST_GA_ampa':{'weight':z}}})
l+=[pl(d, '*', **{'name':''})]
d={}
for keys in s_equal: update(solution, d, keys)
misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}})
misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}})
d['node']['EA']['rate']*=0.7
l[-1]+=pl(d, '=', **{'name':'GIGA_'+str(y)+'_STGA_'+str(z)})
for y in numpy.arange(1.,6.):
for z in numpy.arange(1.,6.):
x=2.5
d={}
for keys in s_mul: update(solution, d, keys)
l+=[pl(d, '*', **{'name':''})]
d={}
for keys in s_equal: update(solution, d, keys)
misc.dict_update(d,{'nest':{'ST_GA_ampa':{'delay':y}}})
misc.dict_update(d,{'nest':{'ST_GI_ampa':{'delay':y}}})
misc.dict_update(d,{'nest':{'GI_ST_gaba':{'delay':z}}})
misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}})
misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}})
d['node']['EA']['rate']*=0.7
l[-1]+=pl(d, '=', **{'name':'STGPdelay_'+str(y)+'_GPSTdelay_'+str(z)})
for y in numpy.arange(2.5, 21., 2.5):
x=2.5
d={}
for keys in s_mul: update(solution, d, keys)
l+=[pl(d, '*', **{'name':''})]
d={}
for keys in s_equal: update(solution, d, keys)
misc.dict_update(d,{'nest':{'CS_ST_ampa':{'delay':y}}})
misc.dict_update(d,{'nest':{'CS_ST_nmda':{'delay':y}}})
misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}})
misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}})
d['node']['EA']['rate']*=0.7
l[-1]+=pl(d, '=', **{'name':'CSSTdelay_'+str(y)})
for y in numpy.arange(2.5, 21., 2.5):
x=2.5
d={}
for keys in s_mul: update(solution, d, keys)
l+=[pl(d, '*', **{'name':''})]
d={}
for keys in s_equal: update(solution, d, keys)
misc.dict_update(d,{'nest':{'C1_M1_ampa':{'delay':y}}})
misc.dict_update(d,{'nest':{'C1_M1_nmda':{'delay':y}}})
misc.dict_update(d,{'nest':{'C2_M2_ampa':{'delay':y}}})
misc.dict_update(d,{'nest':{'C2_M2_nmda':{'delay':y}}})
misc.dict_update(d,{'nest':{'CF_FS_ampa':{'delay':y}}})
misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(x)}}})
misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(x)}}})
d['node']['EA']['rate']*=0.7
l[-1]+=pl(d, '=', **{'name':'CSTRdelay_'+str(y)})
for y in numpy.arange(1, 8., 1.):
x=2.5
d={}
for keys in s_mul: update(solution, d, keys)
l+=[pl(d, '*', **{'name':''})]
d={}
for keys in s_equal: update(solution, d, keys)
misc.dict_update(d,{'nest':{'GA_M1_gaba':{'delay':y}}})
misc.dict_update(d,{'nest':{'GA_M2_gaba':{'delay':y}}})
misc.dict_update(d,{'nest':{'GA_FS_gaba':{'delay':y}}})
d['node']['EA']['rate']*=0.7
l[-1]+=pl(d, '=', **{'name':'GASTRdelay_'+str(y)})
return l
def get():
l = []
solution, s_mul, s_equal = get_solution()
d0 = 0.8
f_beta_rm = lambda f: (1 - f) / (d0 + f * (1 - d0))
d = {}
x = 2.5
for keys in s_mul:
update(solution, d, keys)
l += [pl(d, '*', **{'name': ''})]
d = {}
for keys in s_equal:
update(solution, d, keys)
misc.dict_update(d, {'nest': {'ST': {'beta_I_AMPA_1': f_beta_rm(x)}}})
misc.dict_update(d, {'nest': {'ST': {'beta_I_NMDA_1': f_beta_rm(x)}}})
d['node']['EA']['rate'] *= 0.7
l[-1] += pl(d, '=', **{'name': 'mod_ST_beta_' + str(x)})
for y in numpy.arange(1, 11, 1):
d = {}
for keys in s_mul:
update(solution, d, keys)
l += [pl(d, '*', **{'name': ''})]
d = {}
for keys in s_equal:
update(solution, d, keys)
misc.dict_update(d, {'nest': {'ST': {'beta_I_AMPA_1': f_beta_rm(x)}}})
misc.dict_update(d, {'nest': {'ST': {'beta_I_NMDA_1': f_beta_rm(x)}}})
d['node']['EA']['rate'] *= 0.7
misc.dict_update(d, {'conn': {'GI_FS_gaba': {'lesion': False}}})
misc.dict_update(d, {'conn': {'GI_FS_gaba': {'fan_in0': y}}})
misc.dict_update(d, {'conn': {'GI_M1_gaba': {'lesion': False}}})
misc.dict_update(d, {'conn': {'GI_M1_gaba': {'fan_in0': y}}})
misc.dict_update(d, {'conn': {'GI_M2_gaba': {'lesion': False}}})
misc.dict_update(d, {'conn': {'GI_M2_gaba': {'fan_in0': y}}})
l[-1] += pl(d, '=', **{'name': 'mod_GI_M2_' + str(y)})
for y, EA_rate in zip(numpy.arange(25, 125, 25),
numpy.arange(0.9, 1.3, 0.1)):
d = {}
for keys in s_mul:
update(solution, d, keys)
l += [pl(d, '*', **{'name': ''})]
d = {}
for keys in s_equal:
update(solution, d, keys)
misc.dict_update(d, {'nest': {'ST': {'beta_I_AMPA_1': f_beta_rm(x)}}})
misc.dict_update(d, {'nest': {'ST': {'beta_I_NMDA_1': f_beta_rm(x)}}})
d['node']['EA']['rate'] *= EA_rate
misc.dict_update(d, {'conn': {'M2_GA_gaba': {'lesion': False}}})
misc.dict_update(d, {'conn': {'M2_GA_gaba': {'fan_in0': y}}})
l[-1] += pl(d, '=',
**{'name': 'M2GA_' + str(y) + '_EArate_' + str(EA_rate)})
return l
def get():
l = []
solution, s_mul, s_equal = get_solution()
d0 = 0.8
f_beta_rm = lambda f: (1 - f) / (d0 + f * (1 - d0))
for z0 in numpy.arange(1, 6., 1):
for z1 in numpy.arange(1, 6., 1):
y = 2.5
x = 2.5
d = {}
for keys in s_mul:
update(solution, d, keys)
l += [pl(d, '*', **{'name': ''})]
d = {}
for keys in s_equal:
update(solution, d, keys)
misc.dict_update(d, {'nest': {'ST_GA_ampa': {'delay': z0}}})
misc.dict_update(d, {'nest': {'ST_GI_ampa': {'delay': z0}}})
misc.dict_update(d, {'nest': {'GI_ST_gaba': {'delay': z1}}})
misc.dict_update(d, {'nest': {'C1_M1_ampa': {'delay': y}}})
misc.dict_update(d, {'nest': {'C1_M1_nmda': {'delay': y}}})
misc.dict_update(d, {'nest': {'C2_M2_ampa': {'delay': y}}})
misc.dict_update(d, {'nest': {'C2_M2_nmda': {'delay': y}}})
misc.dict_update(d, {'nest': {'CF_FS_ampa': {'delay': y}}})
misc.dict_update(d,
{'nest': {
'ST': {
'beta_I_AMPA_1': f_beta_rm(x)
}
}})
misc.dict_update(d,
{'nest': {
'ST': {
'beta_I_NMDA_1': f_beta_rm(x)
}
}})
d['node']['EA']['rate'] *= 0.7
l[-1] += pl(d, '=', **{'name': 'CSTRdelay_' + str(y)})
for y in numpy.arange(1, 8., 1.):
z = 2.5 # ms
x = 2.5
d = {}
for keys in s_mul:
update(solution, d, keys)
l += [pl(d, '*', **{'name': ''})]
d = {}
for keys in s_equal:
update(solution, d, keys)
misc.dict_update(d, {'nest': {'C1_M1_ampa': {'delay': z}}})
misc.dict_update(d, {'nest': {'C1_M1_nmda': {'delay': z}}})
misc.dict_update(d, {'nest': {'C2_M2_ampa': {'delay': z}}})
misc.dict_update(d, {'nest': {'C2_M2_nmda': {'delay': z}}})
misc.dict_update(d, {'nest': {'CF_FS_ampa': {'delay': z}}})
misc.dict_update(d, {'nest': {'GA_M1_gaba': {'delay': y}}})
misc.dict_update(d, {'nest': {'GA_M2_gaba': {'delay': y}}})
misc.dict_update(d, {'nest': {'GA_FS_gaba': {'delay': y}}})
d['node']['EA']['rate'] *= 0.7
l[-1] += pl(d, '=', **{'name': 'GASTRdelay_' + str(y)})
for y in numpy.arange(1, 11., 2):
for z in numpy.arange(1, 11., 2):
x = 2.5
d = {}
for keys in s_mul:
update(solution, d, keys)
l += [pl(d, '*', **{'name': ''})]
d = {}
for keys in s_equal:
update(solution, d, keys)
misc.dict_update(d, {'nest': {'C1_M1_ampa': {'delay': y}}})
misc.dict_update(d, {'nest': {'C1_M1_nmda': {'delay': y}}})
misc.dict_update(d, {'nest': {'C2_M2_ampa': {'delay': y}}})
misc.dict_update(d, {'nest': {'C2_M2_nmda': {'delay': y}}})
misc.dict_update(d, {'nest': {'CF_FS_ampa': {'delay': y}}})
misc.dict_update(d, {'nest': {'CS_ST_ampa': {'delay': z}}})
misc.dict_update(d, {'nest': {'CS_ST_nmda': {'delay': z}}})
misc.dict_update(d,
{'nest': {
'ST': {
'beta_I_AMPA_1': f_beta_rm(x)
}
}})
misc.dict_update(d,
{'nest': {
'ST': {
'beta_I_NMDA_1': f_beta_rm(x)
}
}})
d['node']['EA']['rate'] *= 0.7
l[-1] += pl(
d, '=',
**{'name': 'CSTRdelay_' + str(y) + '_CSTdelay_' + str(z)})
return l
def get():
l=[]
solution, s_mul, s_equal=get_solution()
d0=0.8
f_beta_rm=lambda f: (1-f)/(d0+f*(1-d0))
y=2.5
v=[1./10., 3./10, 5./10,7./10, 9./10]
for x in v:
d={}
for keys in s_mul: update(solution, d, keys)
misc.dict_update(d,{'conn':{'GA_M1_gaba':{'fan_in0': x}}})
misc.dict_update(d,{'conn':{'GA_M2_gaba':{'fan_in0': x}}})
misc.dict_update(d,{'nest':{'GA_M1_gaba':{'weight': 1./x}}})
misc.dict_update(d,{'nest':{'GA_M2_gaba':{'weight': 1./x}}})
l+=[pl(d, '*', **{'name':('mod_GA_MS_' + str(x))})]
d={}
for keys in s_equal: update(solution, d, keys)
misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(y)}}})
misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(y)}}})
d['node']['EA']['rate']*=0.7
l[-1]+=pl(d, '=', **{'name':'' })
for x in v:
d={}
for keys in s_mul: update(solution, d, keys)
misc.dict_update(d,{'conn':{'GA_FS_gaba':{'fan_in0': x}}})
misc.dict_update(d,{'nest':{'GA_FS_gaba':{'weight': 1./x}}})
l+=[pl(d, '*', **{'name':('mod_GA_FS_' + str(x))})]
d={}
for keys in s_equal: update(solution, d, keys)
misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(y)}}})
misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(y)}}})
d['node']['EA']['rate']*=0.7
l[-1]+=pl(d, '=', **{'name':'' })
for x in v:
d={}
for keys in s_mul: update(solution, d, keys)
misc.dict_update(d,{'conn':{'GA_FS_gaba':{'fan_in0': x}}})
misc.dict_update(d,{'conn':{'GA_M1_gaba':{'fan_in0': x}}})
misc.dict_update(d,{'conn':{'GA_M2_gaba':{'fan_in0': x}}})
misc.dict_update(d,{'nest':{'GA_M1_gaba':{'weight': 1./x}}})
misc.dict_update(d,{'nest':{'GA_M2_gaba':{'weight': 1./x}}})
misc.dict_update(d,{'nest':{'GA_FS_gaba':{'weight': 1./x}}})
l+=[pl(d, '*', **{'name':('mod_GA_str_' + str(x))})]
d={}
for keys in s_equal: update(solution, d, keys)
misc.dict_update(d,{'nest':{'ST':{'beta_I_AMPA_1': f_beta_rm(y)}}})
misc.dict_update(d,{'nest':{'ST':{'beta_I_NMDA_1': f_beta_rm(y)}}})
d['node']['EA']['rate']*=0.7
l[-1]+=pl(d, '=', **{'name':'' })
return l
