Python机器学习 第三章第4题(答案代码及中文解释)
据美国疾病预防中心数据显示,美国约1/7的成年人患有糖尿病…
#!/usr/bin/python
# -*- coding: utf-8 -*-
import matplotlib.pyplot as plt
import numpy as np
from sklearn import datasets, linear_model
# Load the diabetes dataset 加载糖尿病的数据集
diabetes = datasets.load_diabetes()
# print(diabetes)
# Use only one feature 仅使用其中的一个特性 怀孕次数
diabetes_X = diabetes.data[:, np.newaxis, 2]
print(diabetes_X)
# Split the data into training / testing sets
# 20个为测试数据集,其余为训练集
diabetes_X_train = diabetes_X[:-20]
diabetes_X_test = diabetes_X[-20:]
# Split the targets into training / testing sets
diabetes_Y_train = diabetes_X[:-20]
diabetes_Y_test = diabetes_X[-20:]
# Create linear regression object 创建一个线性回归模型
regr = linear_model.LinearRegression()
# Train the model using the train sets 用线性回归模型和训练集训练
regr.fit(diabetes_X_train, diabetes_Y_train)
# The coefficients 线性模型的斜率
print("Coefficients: \n", regr.coef_)
# The mean square error 平方差
print("Residual sum of squares: %.2f" % np.mean(regr.predict(diabetes_X_test) - diabetes_Y_test) ** 2)
# Explained variance score : 1 is perfect prediction 平方差结果越接近1,模拟度就越好
print("Variance score :%.2f" % regr.score(diabetes_X_test,diabetes_Y_test))
# Plot outputs
plt.scatter(diabetes_X_test, diabetes_Y_test, color="black")
plt.plot(diabetes_X_test, regr.predict(diabetes_X_test), color='blue', linewidth=3)
plt.xticks(())
plt.yticks(())
plt.show()