doc/oldcode/indexing/standard.py
def proto_01_13_steps025dual(abf=exampleABF):
    """IC steps. See how hyperpol. step affects things."""
    swhlab.ap.detect(abf)
    standard_groupingForInj(abf,200)

    for feature in ['freq','downslope']:
        swhlab.ap.plot_values(abf,feature,continuous=False) #plot AP info
        swhlab.plot.save(abf,tag='A_'+feature)

    f1=swhlab.ap.getAvgBySweep(abf,'freq',None,1)
    f2=swhlab.ap.getAvgBySweep(abf,'freq',1,None)
    f1=np.nan_to_num(f1)
    f2=np.nan_to_num(f2)
    Xs=abf.clampValues(abf.dataX[int(abf.protoSeqX[1]+.01)])
    swhlab.plot.new(abf,title="gain function",xlabel="command current (pA)",
                    ylabel="average inst. freq. (Hz)")
    pylab.plot(Xs,f1,'.-',ms=20,alpha=.5,label="step 1",color='b')
    pylab.plot(Xs,f2,'.-',ms=20,alpha=.5,label="step 2",color='r')
    pylab.legend(loc='upper left')
    pylab.axis([Xs[0],Xs[-1],None,None])
    swhlab.plot.save(abf,tag='gain')
def proto_01_13_steps025dual(abf=exampleABF):
    """IC steps. See how hyperpol. step affects things."""
    swhlab.ap.detect(abf)
    standard_groupingForInj(abf,200)

    for feature in ['freq','downslope']:
        swhlab.ap.plot_values(abf,feature,continuous=False) #plot AP info
        swhlab.plot.save(abf,tag='A_'+feature)

    f1=swhlab.ap.getAvgBySweep(abf,'freq',None,1)
    f2=swhlab.ap.getAvgBySweep(abf,'freq',1,None)
    f1=np.nan_to_num(f1)
    f2=np.nan_to_num(f2)
    Xs=abf.clampValues(abf.dataX[int(abf.protoSeqX[1]+.01)])
    swhlab.plot.new(abf,title="gain function",xlabel="command current (pA)",
                    ylabel="average inst. freq. (Hz)")
    pylab.plot(Xs,f1,'.-',ms=20,alpha=.5,label="step 1",color='b')
    pylab.plot(Xs,f2,'.-',ms=20,alpha=.5,label="step 2",color='r')
    pylab.legend(loc='upper left')
    pylab.axis([Xs[0],Xs[-1],None,None])
    swhlab.plot.save(abf,tag='gain')
def proto_01_13_steps025dual ( abf = exampleABF ) : swhlab . ap . detect ( abf ) standard_groupingForInj ( abf , 200 ) for feature in [ 'freq' , 'downslope' ] : swhlab . ap . plot_values ( abf , feature , continuous = False ) #plot AP info swhlab . plot . save ( abf , tag = 'A_' + feature ) f1 = swhlab . ap . getAvgBySweep ( abf , 'freq' , None , 1 ) f2 = swhlab . ap . getAvgBySweep ( abf , 'freq' , 1 , None ) f1 = np . nan_to_num ( f1 ) f2 = np . nan_to_num ( f2 ) Xs = abf . clampValues ( abf . dataX [ int ( abf . protoSeqX [ 1 ] + .01 ) ] ) swhlab . plot . new ( abf , title = "gain function" , xlabel = "command current (pA)" , ylabel = "average inst. freq. (Hz)" ) pylab . plot ( Xs , f1 , '.-' , ms = 20 , alpha = .5 , label = "step 1" , color = 'b' ) pylab . plot ( Xs , f2 , '.-' , ms = 20 , alpha = .5 , label = "step 2" , color = 'r' ) pylab . legend ( loc = 'upper left' ) pylab . axis ( [ Xs [ 0 ] , Xs [ - 1 ] , None , None ] ) swhlab . plot . save ( abf , tag = 'gain' )
IC steps. See how hyperpol. step affects things.
IC steps . See how hyperpol . step affects things .