数学建模社区-数学中国

标题: 线性回归实例-鸢尾花数据集 [打印本页]
作者: 1047521767    时间: 2021-10-23 20:42
标题: 线性回归实例-鸢尾花数据集
                                       文章目录
- B0 |& P7 N8 F; P5 x5 u一、具体实现步骤; Z* E" w  [' j" {$ ?+ d3 @! X
1. 导入Iris鸢尾花数据集- v* O3 {2 M/ `
2. 提取花瓣数据! l& N" x' h( k, n2 R
3. 拆分数据0 x3 X; X# i1 }8 R- R: G6 R6 p: k
4. 训练模型6 |: O3 V$ e1 X/ @& p4 k
二、可视化结果展示
3 p. B8 X" {1 o2 k* o1. 训练集" A. N4 |- {8 W1 O2 w/ f/ a$ I
2. 测试集/ s9 u- m* ?& q: `6 [0 y7 T; |
三、相关知识点讲解; U/ j" ^0 h0 }% {$ z
1. train_test_split()函数
* W+ z+ \, p8 C; C2. LinearRegression()函数2 \$ l/ `- Y2 @7 ?2 }- }5 p& @
3. 散点图与折线统计图的绘制: b4 Q* v7 M/ z* z. A
这篇文章中，我们要通过鸢尾花的花瓣长度预测花瓣宽度( i! |7 Q' e) G6 i) Z

$ A) U9 ]8 u) l  w* L
* e3 s" W$ Z( @环境：Python3.6.5/ t& i- {0 j, r+ q  t5 F3 U1 n
编译器：jupyter notebook
' B1 M0 ^7 _; A, h/ M8 H0 [: z) M9 i' b& h& N  }/ U4 E
一、具体实现步骤
5 f, g' U0 m" E; ?0 M- a( r1. 导入Iris鸢尾花数据集
& L3 U1 r" m# L, TIris鸢尾花数据集共有150条记录，分别是：
$ X0 x8 ^- Z' E; Q/ H# B50条山鸢尾 (Iris-setosa)
2 ]! [& Z( q/ \. a; x. o50条变色鸢尾(Iris-versicolor)
4 W+ J4 D# L) B  R7 Y! V  _( g50条维吉尼亚鸢尾(Iris-virginica)% L; a4 G- }6 m+ z5 Q( ]/ Q2 b
url = "https://archive.ics.uci.edu/ml/machine-learning-databases/iris/iris.data"  * s+ _8 r+ _( u% N
names = ['花萼-length', '花萼-width', '花瓣-length', '花瓣-width', 'class']  
. R9 U: n+ Q% ~' l$ f1 X4 c4 H, x/ `dataset = pd.read_csv(url, names=names)9 \$ b' W2 p. N/ C

6 c0 b* `! x0 O9 K: S& T2. 提取花瓣数据
9 \0 N. }1 E. |9 N* X5 J& |下面我们提取数据集中花瓣宽度与花瓣长度数据，将花瓣数据分为训练数据与测试数据，训练数据用于训练线性回归模型，测试数据用于检测我们的模型的准确率。0 c' e' [! m  U: c/ ^/ {0 v6 ^

5 U9 k% m7 p6 ?- E5 X3 E% }# n! i$ m# A. X
最终我们要达到的效果是：输入花瓣宽度，通过模型预测花瓣宽度。: E! z% @, K. N! t- C$ @6 x" ^
% Z1 F3 i2 u( G# v. B
# e* b( y. d/ a& q
X = dataset["花瓣-length"]
4 X5 B8 p, g* a9 D+ ?& X( U% q( N8 [Y = dataset["花瓣-width"]
: r3 h. @. B: Z+ _' {  C2 I" r/ KX = X.reshape(len(X),1)8 d$ t9 i4 x* S! E8 `4 F' H# F3 ^
Y = Y.reshape(len(Y),1)
/ V& p) ]: H/ d7 x
- o- e  t+ U& N  g2 `3. 拆分数据8 q# w9 W' x) x4 e) n( J
将数据集拆分数据集成训练集、测试集
8 t( @3 b7 D" @4 a0 J. d6 L1 ]! w/ |0 u- G5 \3 L1 K9 n
; [/ u$ v: h2 w$ g/ c: Y* W
from sklearn.model_selection import train_test_split
' |; s# p; x1 h* h& p& SX_train, X_test, Y_train, Y_test = train_test_split(X, Y, test_size=0.2, random_state=0)
5 s" u3 \  g0 v$ c- v
' G/ Z: o1 T# }9 z: D4. 训练模型
- u% d3 ]9 m( }8 c; y/ n$ _' @这里我们需要将我们的训练数据喂给模型进行训练。& M) |! d3 n6 g

; m# Q! R1 l0 n; e' m. Z% z& e; [4 v* K4 S  u; s
from sklearn.linear_model import LinearRegression
9 ^; L# a7 U/ L  h8 ~$ {2 jregressor = LinearRegression(). H7 ]4 w7 F1 A/ W6 G
regressor = regressor.fit(X_train, Y_train)
( a7 ^5 W$ i" [+ h+ D' `! F+ I& e" u* x0 B) b
二、可视化结果展示
. i* ]/ m3 _8 x/ d- }( p1. 训练集1 x9 c" K% h6 ~0 @) R
将训练集中每一朵花的花瓣数据与线性回归模型预测的结果放到同一张统计图中。1 r; u, c6 ~$ |# {

7 A& @9 y& d9 ?' s9 y" Z. o
0 z; d3 m+ a% n- w; oimport matplotlib.pyplot as plt
7 K1 z; b% W+ S1 K4 F) ~6 j' X2 }5 G  L
) b' i7 c; r+ W# ]& h' w& R: i
plt.scatter(X_train, Y_train, color='red')
, W+ g( s8 ?( X) z5 c. n/ i2 jplt.plot(X_train, regressor.predict(X_train), color='green')
1 J* n0 q0 r/ j6 @+ P3 ]- [" T( Splt.xlabel("Iris-length")
4 s- Y# b4 j5 h, cplt.ylabel("Iris-width") - K) ?9 @: j  r; V2 X( U
plt.title("This is train dataset-kzb"), M2 O5 D# C: v2 f
plt.show()
, L! }: u3 ]& Q; c& {- H; J" i9 U$ h2 y. R) D3 j* E+ a
红色的点是训练数据集中的花瓣数据，我们不难看出花瓣长度与宽度是一个线性关系，绿色的线是我们模型拟合的结果。
2 h  e  V/ R! H( ?
/ f# Q; E" J8 U5 k
' \0 r' i, H/ F# {# e. v( C& x
7 W8 H& L, v8 G: d- w7 S# P2. 测试集5 U; l- K8 R+ ~- }( @
将测试集中每一朵花的花瓣数据与线性回归模型预测的结果放到同一张统计图中。% I3 s' d; C& v1 o" q4 f
/ f, O  l* b' g; J/ A
( Y2 l1 Q  ^9 a/ R( a. I, K
plt.scatter(X_test, Y_test, color='blue')+ E( z9 t8 q. l/ O* d; m
plt.plot(X_train, regressor.predict(X_train), color='green')
, J, c+ W& b, c' ~; Q  Bplt.xlabel("Iris-length")
/ C: S. }+ `% _. I9 u" [: vplt.ylabel("Iris-width") & C. n4 g3 ^# u- L
plt.title("This is test dataset-kzb")0 {8 C$ m* x4 h+ k
plt.show(): f  m$ V; h6 f' t. G, E
# g( E- {) p( z/ J/ K0 Z
绿色的点是测试数据集中的花瓣数据，我们可以看出这部分数据也是符合线性关系的，随着集的增大，线性关系会更加明显。: c+ A6 G3 ?' h" ?! W

! S. w) U' F( k; e. {' l
. Q( f  P8 N5 i4 G( {. }! E" C. O" G0 \. X
三、相关知识点讲解
4 l/ B+ o! s, p* S% a8 n1. train_test_split()函数
( U  I% d3 O! {" L3 Vtrain_test_split()：将数据集划分为测试集与训练集。
6 {" ]% ]2 N& _, v! p  e" L
# \3 B2 ]+ _: q  Y( j+ C  Y* X* V. J/ p6 z, q# p' A# N* i5 G
X：所要划分的整体数据的特征集；2 k+ c. Q5 o0 }# U7 c
Y：所要划分的整体数据的结果；
1 O' Q7 h* |# ^# m3 M3 O8 B; U# v6 ptest_size：测试集数据量在整体数据量中的占比（可以理解为X_test与X的比值）；8 U2 f# d, ]$ ^6 y% E
random_state：①若不填或者填0，每次生成的数据都是随机，可能不一样。②若为整数，每次生成的数据都相同；& {: P) R1 ~6 T" o+ [1 g0 }
from sklearn.model_selection import train_test_split. G0 c# q& S2 X  y
X_train, X_test, Y_train, Y_test = train_test_split(X, Y, test_size=0.2, random_state=0)( j1 ?& k  C: I: a5 ~6 c

2 u4 X9 P2 a+ h; U# f2. LinearRegression()函数
/ P% f3 Z  w8 osklearn.linear_model包实现了广义线性模型，包括线性回归、Ridge回归、Bayesian回归等。LinearRegression是其中较为简单的线性回归模型。- I; d, ]- T. u3 J/ {
( c4 e: D. N  c3 W' Z% m

1 t) P% ?0 U1 R# ^0 t1 w解释一下什么是回归：回归最简单的定义是，给出一个点集D，用一个函数去拟合这个点集，并且使得点集与拟合函数间的误差最小，如果这个函数曲线是一条直线，那就被称为线性回归，如果曲线是一条二次曲线，就被称为二次回归。
9 ?, `$ _2 Y- K' L
$ g" I* M( v9 w  l. q' N6 y$ ~5 @0 @9 \
3. 散点图与折线统计图的绘制
& z' W3 |; v4 Xplt.scatter()：绘画出数据的散点图
4 a  l# \4 p1 `/ {! Q; x3 x* Tplt.plot()：绘画出依据模型（LinearRegression的线性回归模型）生成的直线
9 f- i: Y( |) k; a& {有问题请加我QQ3243710560
! O- ^8 s% q2 c. n' |8 O  p1 T% ]4 d0 |  S; w) l! c5 e
6 ?% |, h: M/ t% i0 o8 C




欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5