def get():
l = []
solution = get_solution()
for r in numpy.arange(0, 1200, 200):
r = numpy.float(r)
d = {}
misc.dict_update(d, solution['mul'])
l += [pl(d, '*', **{'name': ''})]
d = {}
misc.dict_update(d, solution['equal'])
dd = {
'conn': {
'GAp_GA_gaba': {
'fan_in': 5
},
'GIp_GA_gaba': {
'fan_in': 25
}
},
'node': {
'EAp': {
'rate': r
}
}
}
misc.dict_update(d, dd)
l[-1] += pl(d, '=', **{'name': 'GAr_{0}'.format(r)})
return l
def get():
l = []
solution = get_solution()
d = {}
misc.dict_update(d, solution['mul'])
l += [pl(d, '*', **{'name': ''})]
d = {}
misc.dict_update(d, solution['equal'])
dd = {
'conn': {
'GAp_GA_gaba': {
'fan_in': 5
},
'GIp_GA_gaba': {
'fan_in': 25
}
}
}
misc.dict_update(d, dd)
l[-1] += pl(d, '=', **{'name': 'control_sim'})
return l
def get():
l = []
solution = get_solution()
d = {}
misc.dict_update(d, solution['mul'])
l += [pl(d, '*', **{'name': ''})]
d = {}
misc.dict_update(d, solution['equal'])
dd = {
'conn': {
'GAp_GA_gaba': {
'fan_in': 5
},
'GIp_GA_gaba': {
'fan_in': 25
}
},
'nest': {
'ST_GA_ampa': {
'weight': 0.259
}
}
} #estimate from beta opt
misc.dict_update(d, dd)
l[-1] += pl(d, '=', **{'name': 'control_sim'})
return l
def get():
l = []
solution = get_solution()
'''
From
script_0000_rEI_1600.0_rEA_0.0_rCS_250.0_wGA_1.0-amp_0.08_1.0_stn_3.0_0
in new_beginning beta slow (1.77)
'''
# misc.dict_update(solution, {'equal':{'node':{'M2p':{'rate':2.9}}}})
l = []
solution = get_solution()
d = {}
misc.dict_update(d, solution['mul'])
l += [pl(d, '*', **{'name': ''})]
d = {}
misc.dict_update(d, solution['equal'])
l[-1] += pl(d, '=', **{'name': 'control_sim'})
return l
def get():
l = []
solution = get_solution()
d = {}
misc.dict_update(d, solution['mul'])
l += [pl(d, '*', **{'name': ''})]
d = {}
misc.dict_update(d, solution['equal'])
l[-1] += pl(d, '=', **{'name': 'control_sim'})
r_EI = [1000, 1600] #range(1000, 1200, 50)+range(1550, 1900, 50)
r_EA = range(100, 1500, 100)
r_CS = [180, 250]
w_STGA = numpy.linspace(1, 0.0, 15)
agg = [[r0, r1, r2, w] for r0, r2 in zip(r_EI, r_CS) for r1 in r_EA
for w in w_STGA]
for rEI, rEA, rCS, w in agg:
rEI = float(rEI)
rEA = float(rEA)
rCS = float(rCS)
w = float(w)
d = {}
misc.dict_update(d, solution['mul'])
d['nest']['ST_GA_ampa']['weight'] *= w
l += [pl(d, '*', **{'name': ''})]
d = {}
misc.dict_update(d, solution['equal'])
d['node']['EI']['rate'] = rEI
d['node']['EA']['rate'] = rEA
misc.dict_update(d, {'node': {'CS': {'rate': rCS}}})
s = 'rEI_{0}_rEA_{1}_rCS_{2}_STGAw_{3}'.format(rEI, rEA, rCS,
round(w, 2))
l[-1] += pl(d, '=', **{'name': s})
return l
def get():
l = []
solution = get_solution()
'''
From
script_0000_rEI_1600.0_rEA_0.0_rCS_250.0_wGA_1.0-amp_0.08_1.0_stn_3.0_0
in new_beginning beta slow (1.77)
'''
# misc.dict_update(solution, {'equal':{'node':{'M2p':{'rate':2.9}}}})
misc.dict_update(solution, {'equal': {'node': {'M2p': {'rate': 2.7}}}})
d = {}
misc.dict_update(d, solution['mul'])
l += [pl(d, '*', **{'name': ''})]
d = {}
misc.dict_update(d, solution['equal'])
l[-1] += pl(d, '=', **{'name': 'control_sim'})
d = {}
misc.dict_update(d, solution['mul'])
dd = {
'nest': {
'ST_GA_ampa': {
'weight': 1.
},
'GA_GA_gaba': {
'weight': 1.
},
'GI_GA_gaba': {
'weight': 1.
}
}
}
misc.dict_update(d, dd)
l += [pl(d, '*', **{'name': ''})]
d = {}
misc.dict_update(d, solution['equal'])
l[-1] += pl(d, '=', **{'name': 'GAw_*1.'})
return l
def get():
l = []
solution = get_solution()
rEI = 1700.0
rEA = 200.0
rCS = 250.0
d = {}
misc.dict_update(d, solution['mul'])
l += [pl(d, '*', **{'name': ''})]
d = {}
misc.dict_update(d, solution['equal'])
d['node']['EI']['rate'] = rEI
d['node']['EA']['rate'] = rEA
misc.dict_update(d, {'node': {'CS': {'rate': rCS}}})
s = 'rEI_{0}_rEA_{1}_rCS_{2}'.format(rEI, rEA, rCS)
l[-1] += pl(d, '=', **{'name': s})
rEI = 800.0
rEA = 100.0
rCS = 170.0
d = {}
misc.dict_update(d, solution['mul'])
l += [pl(d, '*', **{'name': ''})]
d = {}
misc.dict_update(d, solution['equal'])
d['node']['EI']['rate'] = rEI
d['node']['EA']['rate'] = rEA
misc.dict_update(d, {'node': {'CS': {'rate': rCS}}})
s = 'rEI_{0}_rEA_{1}_rCS_{2}'.format(rEI, rEA, rCS)
l[-1] += pl(d, '=', **{'name': s})
return l
def get():
l = []
solution = get_solution()
'''
From
script_0000_rEI_1600.0_rEA_0.0_rCS_250.0_wGA_1.0-amp_0.08_1.0_stn_3.0_0
in new_beginning beta slow (1.77)
'''
# misc.dict_update(solution, {'equal':{'node':{'M2p':{'rate':2.9}}}})
'''
From estimate single GA, ST, subtract since single not estimate overestimated
exitation needed. IN network inhibition onto ST from GP is weaker at the same
frequency. And ST_GA_ampa need to be stronger in full network than in single
mode to optain similar firing rates.
'''
misc.dict_update(solution, {'equal': {'node': {'CS': {'rate': 190 - 10}}}})
misc.dict_update(
solution, {'equal': {
'nest': {
'ST_GA_ampa': {
'weight': 0.236 * 1.1
}
}
}})
misc.dict_update(solution, {'equal': {'node': {'EA': {'rate': 210.}}}})
misc.dict_update(solution, {'equal': {'node': {'EI': {'rate': 1200.0}}}})
d = {}
misc.dict_update(d, solution['mul'])
l += [pl(d, '*', **{'name': ''})]
d = {}
misc.dict_update(d, solution['equal'])
l[-1] += pl(d, '=', **{'name': 'control_sim'})
return l
def get():
l = []
solution = get_solution()
rEA = 200.0
rES = 3200.0
for rCS, rM2, rEI in [
[250.0, 700.0, 1700.0],
[200.0, 740.0, 1700.0],
[250.0, 700.0, 1400.0],
]:
d = {}
misc.dict_update(d, solution['mul'])
l += [pl(d, '*', **{'name': ''})]
d = {}
misc.dict_update(d, solution['equal'])
d['node']['EI']['rate'] = rEI
d['node']['EA']['rate'] = rEA
d['node']['C2']['rate'] = rM2
misc.dict_update(d, {'node': {'CS': {'rate': rCS}}})
misc.dict_update(d, {'node': {'ES': {'rate': rES}}})
s = 'rEI_{0}_rEA_{1}_rCS_{2}_rES_{3}_rM2_{4}'.format(
rEI, rEA, rCS, rES, rM2)
l[-1] += pl(d, '=', **{'name': s})
rEI = 1700.0
rEA = 200.0
rCS = 250.0
rES = 1800.0
for wGISN in [0.8, 0.6, 0.4, 0.2]:
d = {}
misc.dict_update(d, solution['mul'])
misc.dict_update(d, {'nest': {'GI_SN_gaba': {'weight': wGISN}}})
l += [pl(d, '*', **{'name': ''})]
d = {}
misc.dict_update(d, solution['equal'])
d['node']['EI']['rate'] = rEI
d['node']['EA']['rate'] = rEA
d['node']['C2']['rate'] = rM2
misc.dict_update(d, {'node': {'CS': {'rate': rCS}}})
misc.dict_update(d, {'node': {'ES': {'rate': rES}}})
s = 'rEI_{0}_rEA_{1}_rCS_{2}_rES_{3}_rM2_{4}_wGISN_{5}'.format(
rEI, rEA, rCS, rES, rM2, wGISN)
l[-1] += pl(d, '=', **{'name': s})
rEA = 200.0
rES = 3200.0
for rCS, rM2, rEI in [
[250.0, 740.0, 1700.0],
]:
d = {}
misc.dict_update(d, solution['mul'])
l += [pl(d, '*', **{'name': ''})]
d = {}
misc.dict_update(d, solution['equal'])
d['node']['EI']['rate'] = rEI
d['node']['EA']['rate'] = rEA
d['node']['C2']['rate'] = rM2
misc.dict_update(d, {'node': {'CS': {'rate': rCS}}})
misc.dict_update(d, {'node': {'ES': {'rate': rES}}})
s = 'rEI_{0}_rEA_{1}_rCS_{2}_rES_{3}_rM2_{4}'.format(
rEI, rEA, rCS, rES, rM2)
l[-1] += pl(d, '=', **{'name': s})
return l
def get():
l = []
solution = get_solution()
# d={}
# misc.dict_update(d, solution['mul'])
# l+=[pl(d, '*', **{'name':''})]
#
# d={}
# misc.dict_update(d, solution['equal'])
# l[-1]+=pl(d, '=', **{'name':'control_sim'})
r_EI = range(800, 1000, 50)
r_EA = range(0, 500, 100)
r_CS = range(140, 190, 10)
agg = [[r0, r1, r2] for r0 in r_EI for r1 in r_EA for r2 in r_CS]
for rEI, rEA, rCS in agg:
rEI = float(rEI)
rEA = float(rEA)
rCS = float(rCS)
d = {}
misc.dict_update(d, solution['mul'])
l += [pl(d, '*', **{'name': ''})]
d = {}
misc.dict_update(d, solution['equal'])
d['node']['EI']['rate'] = rEI
d['node']['EA']['rate'] = rEA
misc.dict_update(d, {'node': {'CS': {'rate': rCS}}})
s = 'rEI_{0}_rEA_{1}_rCS_{2}'.format(rEI, rEA, rCS)
l[-1] += pl(d, '=', **{'name': s})
r_EI = range(1550, 1750, 50)
r_EA = range(0, 500, 100)
r_CS = range(250, 300, 10)
agg = [[r0, r1, r2] for r0 in r_EI for r1 in r_EA for r2 in r_CS]
for rEI, rEA, rCS in agg:
rEI = float(rEI)
rEA = float(rEA)
rCS = float(rCS)
d = {}
misc.dict_update(d, solution['mul'])
l += [pl(d, '*', **{'name': ''})]
d = {}
misc.dict_update(d, solution['equal'])
d['node']['EI']['rate'] = rEI
d['node']['EA']['rate'] = rEA
misc.dict_update(d, {'node': {'CS': {'rate': rCS}}})
s = 'rEI_{0}_rEA_{1}_rCS_{2}'.format(rEI, rEA, rCS)
l[-1] += pl(d, '=', **{'name': s})
return l
def get():
l = []
solution = get_solution()
# misc.dict_update(solution, {'mul':d})
d = {}
misc.dict_update(d, solution['mul'])
l += [pl(d, '*', **{'name': ''})]
d = {}
misc.dict_update(d, solution['equal'])
dd = {
'conn': {
'GAp_GA_gaba': {
'fan_in': 5
},
'GIp_GA_gaba': {
'fan_in': 25
}
}
}
misc.dict_update(d, dd)
l[-1] += pl(d, '=', **{'name': 'control_sim'})
d = {}
misc.dict_update(d, solution['mul'])
l += [pl(d, '*', **{'name': ''})]
d = {}
misc.dict_update(d, solution['equal'])
dd = {
'conn': {
'GAp_GA_gaba': {
'fan_in': 5
},
'GIp_GA_gaba': {
'fan_in': 25
}
},
'node': {
'EAp': {
'rate': 400.0
}
}
}
misc.dict_update(d, dd)
l[-1] += pl(d, '=', **{'name': 'control_sim'})
d = {}
misc.dict_update(d, solution['mul'])
l += [pl(d, '*', **{'name': ''})]
d = {}
misc.dict_update(d, solution['equal'])
dd = {
'conn': {
'GAp_GA_gaba': {
'fan_in': 5
},
'GIp_GA_gaba': {
'fan_in': 25
}
},
'node': {
'EAp': {
'rate': 800.0
}
}
}
misc.dict_update(d, dd)
l[-1] += pl(d, '=', **{'name': 'EA800.'})
d = {}
misc.dict_update(d, solution['mul'])
l += [pl(d, '*', **{'name': ''})]
d = {}
misc.dict_update(d, solution['equal'])
dd = {
'conn': {
'GAp_GA_gaba': {
'fan_in': 5
},
'GIp_GA_gaba': {
'fan_in': 25
}
},
'node': {
'EAp': {
'rate': 1200.0
}
}
}
misc.dict_update(d, dd)
l[-1] += pl(d, '=', **{'name': 'EA800.'})
# d={}
# misc.dict_update(d, solution['mul'])
# l+=[pl(d, '*', **{'name':''})]
#
# d={}
# # GI predominently connect to to GA
#
# misc.dict_update(d, solution['equal'])
# dd={'conn': {'GAp_GA_gaba':{'fan_in': 15},
# 'GIp_GA_gaba':{'fan_in': 15 }}}
# misc.dict_update(d, dd)
# l[-1]+=pl(d, '=', **{'name':'fanin1515.'})
#
#
#
#
#
# d={}
# misc.dict_update(d, solution['mul'])
# dd={'nest':{'ST_GA_ampa':{'weight': 1.},
# 'GA_GA_gaba':{'weight': 1.},
# 'GI_GA_gaba':{'weight': 1.}}}
#
# misc.dict_update(d, dd)
# l+=[pl(d, '*', **{'name':''})]
#
# d={}
# misc.dict_update(d, solution['equal'])
# dd={'conn': {'GAp_GA_gaba':{'fan_in': 5},
# 'GIp_GA_gaba':{'fan_in': 25}}}
# misc.dict_update(d, dd)
#
# l[-1]+=pl(d, '=', **{'name':'GAw_*1.'})
#
# d={}
# misc.dict_update(d, solution['mul'])
# x=0.1
# dd={'nest':{'ST_GA_ampa':{'weight':x},
# 'GA_GA_gaba':{'weight':x},
# 'GI_GA_gaba':{'weight':x}}}
#
# misc.dict_update(d, dd)
# l+=[pl(d, '*', **{'name':''})]
#
# d={}
# misc.dict_update(d, solution['equal'])
# l[-1]+=pl(d, '=', **{'name':'GAw_*.1'})
#
#
# d={}
# misc.dict_update(d, solution['mul'])
# dd={'node':{'EAp':{'rate':0.0}}}
# misc.dict_update(d, dd)
# l+=[pl(d, '*', **{'name':''})]
#
# d={}
# misc.dict_update(d, solution['equal'])
# l[-1]+=pl(d, '=', **{'name':'EA0.'})
#
#
# d={}
# misc.dict_update(d, solution['mul'])
#
# misc.dict_update(d, dd)
# l+=[pl(d, '*', **{'name':''})]
#
# d={}
# misc.dict_update(d, solution['equal'])
# dd={'node':{'EAp':{'rate':1200*0.25}}}
# l[-1]+=pl(d, '=', **{'name':'EA+.'})
return l