在python中使回归线适合指数函数

时间:2016-09-19 03:49:42

标签: python python-3.x numpy machine-learning linear-regression

我正在尝试学习如何解释用Python创建的指数函数的线性回归模型。我通过首先通过获取自然对数将指数Y数据转换为直线来创建模型。然后我创建一个线性模型并记下斜率和截距。最后,我尝试使用斜率和截距计算样本值。具体来说,我尝试在X = 1.1时计算Y. Y应该是~2.14但我的模型解释产生的Y值为3.78。

问题1:在解释模型时我做错了什么。

问题2:我必须重塑X阵列或者在regr.fit中出错。为什么我必须重塑X阵列。

代码如下:

import matplotlib.pyplot as plt
import numpy as np
from sklearn import datasets, linear_model

# create some exponential data
X = np.arange(1, 10, 0.1)
print(X)
Y = np.power(2, X)
print(Y)

# transform the exponential Y data to make it a straight line
ln_Y = np.log(Y)

# show the exponential plot
plt.scatter(X, Y)
plt.show()

# Create linear regression object
regr = linear_model.LinearRegression()

# reshape the X to avoid regr.fit errors
X = np.reshape(X, (X.size, 1))

# Train the model using the training sets
regr.fit(X,ln_Y)

# The coefficients
print('Slope: \n', regr.coef_)
print('Intercept: \n', regr.intercept_)

# predict Y when X = 1.1 (should be approximately 2.14354693)
# equation = e^(0.00632309*1.1) + 2.7772517886
print("predicted val = ", np.exp(0.00632309*1.1) + 2.7772517886)

1 个答案:

答案 0 :(得分:1)

确保您已获得最新版本的scikit;我有不同的系数给你:

Slope:
 [ 0.69314718]
Intercept:
 4.4408920985e-16

你需要获取整个表达式的exp,而不仅仅是x术语:

In [17]: np.exp(0.69314718*1.1 + 4.4408920985e-16)
Out[17]: 2.1435469237522917