예제 #1
0
#Make custom header (if there is a new line, remove)
newHeader = data.header
newHeader[-1] = re.sub(r"\n", "",newHeader[-1])
newHeader.append("A")
newHeader.append("B")
newHeader.append("dGf")
newHeader.append("sd dGf")
newHeader.append("m")
newHeader.append("sd m")
newHeader.append("C 1/2")
newHeader.append("sd C 1/2")
newHeader.append("b")
newHeader.append("sd b")
newHeader[-1] = newHeader[-1]+"\n"
#Write header to file
fw.writeList(data.header)

def makeModel(x, dGf, m):
    kox = 0.013 #1/s
    t = 180 #s
    RT = 0.592154 #kcal/mol
    Kfold = 1+np.exp(-(dGf + np.multiply(m, x))/RT)
    return B+(A-B)*np.exp(-(kox*t)/Kfold)

def makeChalfModel(x, chalf, b):
    return A + (B-A)/(1+np.exp(-(x-chalf)/(b)))

allattempts = []
failedOnRuntime = []
failedOnCovMatrix = []
numSuccess = 0