Python get_solution_slow_GP_striatum Examples

Example #1

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

Example #2

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

Example #3

File: oscillation_perturbations13_slow.py Project: mickelindahl/bgmodel

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

Example #4

File: oscillation_perturbations22_slow.py Project: sunzhe839/bgmodel

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

Example #5

File: oscillation_perturbations13_slow_dopMSGP.py Project: sunzhe839/bgmodel

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

Example #6

File: oscillation_perturbations15_slow.py Project: mickelindahl/bgmodel

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

Example #7

File: oscillation_perturbations23.py Project: sunzhe839/bgmodel

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

Example #8

File: oscillation_perturbations27_slow.py Project: mickelindahl/bgmodel

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

Example #9

File: oscillation_perturbations27_slow.py Project: sunzhe839/bgmodel

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