在sklearn中提供一些数据包可以用来进行简单的机器学习,今天我们以糖尿病数据为例,来进行线性回归分析
1.导所需要的包
import numpy as np
from sklearn.datasets import load_diabetes
from sklearn.utils import shuffle
from sklearn.model_selection import train_test_split
import matplotlib.pyplot as plt
from sklearn import linear_model
from sklearn.metrics import mean_squared_error, r2_score
2.导入数据,对数据集进行划分
在sklearn中自带一些数据集供用户使用,例如其中diabetes就是部分糖尿病人的数据,数据中存在的键值有dict_keys(['data', 'target', 'frame', 'DESCR', 'feature_names', 'data_filename', 'target_filename'])
这些键解释分别是
-
data:特征值 (数组)
-
target:标签值 (数组)
-
target_names:标签 (列表)
-
DESCR:数据集描述
-
feature_names:特征 (列表)
-
filename:iris.csv 文件路径
下面我们导入数据集并进行划分
diabetes = load_diabetes()
data = diabetes.data
target = diabetes.target
X, y = shuffle(data, target, random_state=13)
X= X.astype(np.float32)
y = y.reshape((-1, 1))
X_train,X_test, y_train, y_test = train_test_split(X, y, test_size=0.25, random_state=42)
print(X_train.shape, y_train.shape, X_test.shape, y_test.shape)
3.建立数据模型并训练数据
reg = linear_model.LinearRegression()
reg.fit(X_train, y_train)
4.对模型进行评估
y_pred = reg.predict(X_test)
print('Coefficients: n',reg.coef_)
print("sMean squared error: %.2f"
%mean_squared_error(y_test, y_pred))
print('Variance score: %.2f'%r2_score(y_test, y_pred))
print(r2_score(y_test, y_pred))
plt.scatter(range(X_test.shape[0]), y_test, color='red')
plt.plot(range(X_test.shape[0]), y_pred, color='green', linewidth=3)
plt.xticks(())
plt.yticks(())
plt.show();