paramsearchGA_DopDep_nonlinear.py 文件源码

def diffD1D2(Flags):
    # To check if the difference between D1 and D2 amplifies downstream
    # First decide which model to use
    delay = 1
    if Flags == "Allsym":
        (d1,d2,fsi,ti,ta,stn,gpi,ipctx,A,B,params) = calcRates(Flags,delay)
        # D = Direct pathway, ID = Indirect pathway, HD = Hyperdirect pathway
        # Reducing a full recurrent matrix leads to postive and negative contributions in ID and HD instead of pure just positive contributions
        D = params['gpid1']
        print "Direct",D

        de1 = 1. + (params['d1d1']*params['fsifsi']) - (params['stnti']*params['stnstn'])
        Ex1 = (params['stnti']*params['d1d1']*params['fsifsi']*params['gpistn'])/de1        
        Ex2 = (params['gpid1']*params['stnstn']*params['d1ta']*params['fsifsi'])/de1        
        IDpos = params['stnti']*params['tid2']*(params['d1ta']*params['gpid1']*params['tistn']+ Ex1 + Ex2)

        print "IDpos",IDpos

        Ex3 = (params['d1ta']+params['d1ta']*params['tata']+((params['stnta']*params['fsiti']*params['d1fsi'])/de1))
        IDneg = params['gpid1']+params['gpiti']*params['tid2']+params['gpid1']*params['stnstn']*params['tid2']*Ex3  
        print "IDneg",IDneg

        HDpos = (params['jstnctx']*params['gpid1']*params['stnstn']*params['fsiti']*params['d1fsi'])/de1
        print "HDpos",HDpos
        Ex4 = params['jstnctx']*params['d1d1']*params['fsifsi']*params['stnti']*params['gpistn']
        Ex5 = params['jstnctx']*params['gpid1']*params['stnstn']*params['d1ta']*params['fsifsi']
        HDneg = (Ex4 + Ex5)/de1 
        print "HDneg",HDneg

        d1fix = np.mean(d1[:-10])
        d2fix = np.mean(d2[:-10])
        DelMSN = d1fix - d2fix
        DelGpi = (D*d1fix) + ((IDpos - IDneg)*d2fix)

    return (DelMSN,DelGpi)