使用 sklearn 进行多项式回归的最简单方法?
Simplest way to make a polynomial regression with sklearn?
我有一些数据不适合线性回归:
其实应该拟合一个二次函数'exactly':
P = R*I**2
我正在做这个:
model = sklearn.linear_model.LinearRegression()
X = alambres[alambre]['mediciones'][x].reshape(-1, 1)
Y = alambres[alambre]['mediciones'][y].reshape(-1, 1)
model.fit(X,Y)
是否有机会通过以下方式解决它:
model.fit([X,X**2],Y)
你可以使用 numpy 的 polyfit.
import numpy as np
from matplotlib import pyplot as plt
X = np.linspace(0, 100, 50)
Y = 23.24 + 2.2*X + 0.24*(X**2) + 10*np.random.randn(50) #added some noise
coefs = np.polyfit(X, Y, 2)
print(coefs)
p = np.poly1d(coefs)
plt.plot(X, Y, "bo", markersize= 2)
plt.plot(X, p(X), "r-") #p(X) evaluates the polynomial at X
plt.show()
输出:
[ 0.24052058 2.1426103 25.59437789]
使用多项式特征。
import numpy as np
from sklearn.preprocessing import PolynomialFeatures
x = np.array([[1,],[2,],[3,]])
X = PolynomialFeatures(degree=2).fit_transform(x)
X
输出:
array([[1., 1., 1.],
[1., 2., 4.],
[1., 3., 9.]])
我有一些数据不适合线性回归:
其实应该拟合一个二次函数'exactly':
P = R*I**2
我正在做这个:
model = sklearn.linear_model.LinearRegression()
X = alambres[alambre]['mediciones'][x].reshape(-1, 1)
Y = alambres[alambre]['mediciones'][y].reshape(-1, 1)
model.fit(X,Y)
是否有机会通过以下方式解决它:
model.fit([X,X**2],Y)
你可以使用 numpy 的 polyfit.
import numpy as np
from matplotlib import pyplot as plt
X = np.linspace(0, 100, 50)
Y = 23.24 + 2.2*X + 0.24*(X**2) + 10*np.random.randn(50) #added some noise
coefs = np.polyfit(X, Y, 2)
print(coefs)
p = np.poly1d(coefs)
plt.plot(X, Y, "bo", markersize= 2)
plt.plot(X, p(X), "r-") #p(X) evaluates the polynomial at X
plt.show()
输出:
[ 0.24052058 2.1426103 25.59437789]
使用多项式特征。
import numpy as np
from sklearn.preprocessing import PolynomialFeatures
x = np.array([[1,],[2,],[3,]])
X = PolynomialFeatures(degree=2).fit_transform(x)
X
输出:
array([[1., 1., 1.],
[1., 2., 4.],
[1., 3., 9.]])