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()
    
    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()
    
    # 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))
    
    
#     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()

    # 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))
     
    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
Пример #7
0
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()

    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
    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))
    
    
    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))
    

    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))

    # 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))
    
    
    # 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))

    #     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()
     
    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))
    
   
    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
Пример #19
0
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