-1

グローバル フィット問題 (schild 分析) に lmfit を使用しようとしています。いくつかの共有パラメーターと、これらの共有に基づいて計算されたものがあります。ある時点で、関数は負の数のログに遭遇し、[nan] リストをスローして失敗します。どうすればそれを防ぐことができますか? ありがとうございました。

def g1(params,xdata,ydata):

    hillSlope = params['hillSlope'].value
    schildSlope = params['SchildSlope'].value
    top = params['top'].value
    bottom = params['bottom'].value
    pA2 = params['pA2'].value


    EC50_1 = params['ec50_2'].value
    B_1 = params['B_2'].value
    Antag_1 = 1+(B_1/(10**(-1*pA2)))**schildSlope
    LogEC_1=np.log10(EC50_1*Antag_1)
    y_model_1 = y_model_1 = bottom + (top-bottom)/(1+10**((LogEC_1-xdata)*hillSlope))
4

2 に答える 2