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