【发布时间】:2016-09-30 00:50:02
【问题描述】:
我试图用 scipy curve_fit 对我的数据进行多元对数回归,结果期望得到一条线,但得到一条曲线。 这是我使用的代码:
Quercetin=[23,195,6,262,272,158,79,65,136,198]
Naringenin=[11,4,8,6,6,7,6,9,7,9]
Rutin=[178,165,93,239,202,3325,4427,7607,3499,1762]
TEAC=[23,189,37,265,290,267,362,388,364,321]
import matplotlib.pyplot as plt
import scipy
from scipy.optimize import curve_fit
import numpy as np
def func(x, a, b, c, d,e):
m=np.log(a*x[0]+b*x[1]+c*x[2])
return d*(m)+e
x=scipy.array([Quercetin, Naringenin,Rutin])
y=scipy.array(TEAC)
popt, pcov = curve_fit(func, x ,y)
print (popt)
plt.plot(func(x,*popt),y,'ro-')
plt.show()
虽然我想得到这样的东西:
谁能告诉我我做错了什么? 如果这很重要,我会在 Windows 10 上使用来自 Anaconda 的 Python 3.5。
【问题讨论】:
标签: python scipy curve-fitting