#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