QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3978|回复: 0
打印 上一主题 下一主题

[代码资源] 线性回归实例-鸢尾花数据集

[复制链接]
字体大小: 正常 放大

1178

主题

15

听众

1万

积分

  • TA的每日心情
    开心
    2023-7-31 10:17
  • 签到天数: 198 天

    [LV.7]常住居民III

    自我介绍
    数学中国浅夏
    跳转到指定楼层
    1#
    发表于 2021-10-23 20:42 |只看该作者 |正序浏览
    |招呼Ta 关注Ta
                                           文章目录# v0 O* d' I7 g4 ^9 m! y
    一、具体实现步骤+ m/ I5 u% u2 E4 m! y- g3 G5 g
    1. 导入Iris鸢尾花数据集, D" j8 P! @2 Y. T! x# T4 t1 K
    2. 提取花瓣数据
    / F, B$ s% }- z1 R/ V# D2 b% D" Q3. 拆分数据
    / [8 t+ {7 ~+ H1 B4. 训练模型
    9 O& P' ~/ E2 b& G二、可视化结果展示9 P* s# d& S  v+ B- h4 E/ }7 x
    1. 训练集! B, b# D% M2 N7 _
    2. 测试集& m) u! X% u4 i- A
    三、相关知识点讲解5 B! |' u! I* k; A4 J% ^
    1. train_test_split()函数; @2 ~8 p; T; ~% i& h3 v
    2. LinearRegression()函数  g$ T+ ?3 S& K
    3. 散点图与折线统计图的绘制& P* }/ F6 G2 W0 ?) Z+ Y
    这篇文章中,我们要通过鸢尾花的花瓣长度预测花瓣宽度  Y# d) V3 d3 i& K/ @9 z
    ) ~9 m6 L  `2 P. q  a2 A

    / U8 C; d. I0 r* x% V环境:Python3.6.5
    8 M) p! W. C' W$ H7 |  K编译器:jupyter notebook
    , f0 I* I& G- E0 Q" W1 R" T' }8 j8 r6 y
    一、具体实现步骤
    8 a# a2 D" p$ {$ L" e" C1. 导入Iris鸢尾花数据集
    / g- t' Y- c% _4 y& zIris鸢尾花数据集共有150条记录,分别是:) b+ l4 }3 E" }. F4 L
    50条山鸢尾 (Iris-setosa)' Y& h* k& K$ e
    50条变色鸢尾(Iris-versicolor)$ t$ m" H, Y- ^* ]  Z& d" j
    50条维吉尼亚鸢尾(Iris-virginica)7 ?( f5 [8 Z9 b! M+ ]
    url = "https://archive.ics.uci.edu/ml/machine-learning-databases/iris/iris.data"  0 g8 X% m3 y8 g  D: R5 j9 T: {
    names = ['花萼-length', '花萼-width', '花瓣-length', '花瓣-width', 'class']  7 X7 q) h% K: N: e4 \1 v9 }
    dataset = pd.read_csv(url, names=names)# C; E- ]6 W2 o0 M3 W
    , D: d9 B/ L" C- a
    2. 提取花瓣数据* u1 \* i3 O3 r4 e$ O0 \
    下面我们提取数据集中花瓣宽度与花瓣长度数据,将花瓣数据分为训练数据与测试数据,训练数据用于训练线性回归模型,测试数据用于检测我们的模型的准确率。9 v9 B/ L1 v- h7 [

    + s/ M! l1 D0 a; N# M
    ; H& {* v, g" {6 _1 A
    最终我们要达到的效果是:输入花瓣宽度,通过模型预测花瓣宽度。
    8 u1 L  V+ h+ [! _) I' S2 }3 L
    3 V- W6 K0 J, e; {4 m: W
      |6 ]) n, q. ]
    X = dataset["花瓣-length"]) _4 `/ [6 {! F8 L6 Q% x9 s* p
    Y = dataset["花瓣-width"]
    8 p- q: k' }  `X = X.reshape(len(X),1)
    6 S  _7 e8 D! B5 `+ L; FY = Y.reshape(len(Y),1)! F# t2 N: |3 U+ a+ U* E! ]
    ' }: u2 |8 u! p0 N& b: \$ S6 A
    3. 拆分数据4 Z, n& M. n# Z1 {! ^
    将数据集拆分数据集成训练集、测试集5 A0 e+ Q; _2 i  n

    . H  ]$ c9 M* B: L2 L$ \
    * |% l/ V2 S9 C7 Q4 t
    from sklearn.model_selection import train_test_split: }( {& T' Q( C* r
    X_train, X_test, Y_train, Y_test = train_test_split(X, Y, test_size=0.2, random_state=0)0 V! b7 P5 w. E2 M6 ?9 N: |: P

    & O$ e2 n" m0 C3 i! d4. 训练模型7 ]9 E( T0 ?: p3 Z( N, T
    这里我们需要将我们的训练数据喂给模型进行训练。
    8 n' q( N1 \! [5 M; C! L7 I" Q3 K) C# C; l
    2 R$ A% N3 A$ j
    from sklearn.linear_model import LinearRegression1 B# k' {4 ~' v5 D  z( S
    regressor = LinearRegression()
    $ x# o& n* x6 X( _regressor = regressor.fit(X_train, Y_train)( D- B. J8 \5 s/ V8 _
    $ o# f; t: I; {; M- G1 U6 a  z
    二、可视化结果展示
    / s/ |& y7 E- M. V" f1. 训练集' T' u/ h4 _' y" R: Q3 h( `
    将训练集中每一朵花的花瓣数据与线性回归模型预测的结果放到同一张统计图中。
      e; X! [) H- {" R1 |: p
    1 _# {* ~: Y( E3 J/ h% K2 [
      X9 Z' W# d# i9 D
    import matplotlib.pyplot as plt& p, [5 M  a$ I5 n3 P7 h& Z3 T

    $ S9 k# u. ]. u0 x- e  o
      d% ]. u$ ]' C
    plt.scatter(X_train, Y_train, color='red')5 z  l7 t2 d+ F/ I8 M- h
    plt.plot(X_train, regressor.predict(X_train), color='green')4 |3 M9 T* E- m0 Y$ i
    plt.xlabel("Iris-length")
    & o8 E; a5 |  r5 l  J0 p4 Iplt.ylabel("Iris-width")
    ; n) ?' g2 ^0 R3 \! v- wplt.title("This is train dataset-kzb")
    8 y- k, T' ?$ ?* O, ]* yplt.show()
    ( H/ t) i: u% e) v5 A+ f; `/ P* @6 J! W9 J
    红色的点是训练数据集中的花瓣数据,我们不难看出花瓣长度与宽度是一个线性关系,绿色的线是我们模型拟合的结果。
    7 F4 L9 v6 Y0 ^
    / P. N% n' T* o! T9 x/ i6 A: m! V1 e7 F$ L+ g& p

    - {/ T6 u% k! b& N5 u: y2. 测试集- z" I7 w1 j: j  a% X/ x. a
    将测试集中每一朵花的花瓣数据与线性回归模型预测的结果放到同一张统计图中。0 H# u7 C' V, m7 n% B3 }

    / A+ W, S) n* o2 X+ U2 v- d! Q

    / l3 Y+ P3 }8 U$ w4 ?, @7 z7 vplt.scatter(X_test, Y_test, color='blue')
    : w: [# S; V+ m' iplt.plot(X_train, regressor.predict(X_train), color='green')
    $ M8 b9 d" y# t; z" [plt.xlabel("Iris-length") * Z) o  Z: {# v! K  D8 ~4 E+ v1 B  J
    plt.ylabel("Iris-width")
    , m, o! F0 G' g+ d* ~( W+ @plt.title("This is test dataset-kzb")
    + ^* k* K  g. z; ?7 h0 G# H2 zplt.show()' a: i+ f# E, n7 k  |- ~
    0 B, c1 ^1 X  J. g" m
    绿色的点是测试数据集中的花瓣数据,我们可以看出这部分数据也是符合线性关系的,随着集的增大,线性关系会更加明显。8 x6 W1 m8 w. q% |* o

    # n; r5 _) Q1 k( x1 b
    # I% m- P9 Z, j# e, M
    * B& R- P" F1 C  i/ t) N
    三、相关知识点讲解
    & ^) b- O' z+ L6 Z1. train_test_split()函数0 K# }" |9 z+ R* Q5 X8 n
    train_test_split():将数据集划分为测试集与训练集。
    9 u9 d! e. @" w! N
    , {4 S8 L* d. a' Z

    ' F+ K$ T& g5 _8 \8 _0 lX:所要划分的整体数据的特征集;# Q& {& z3 u+ w4 C7 a
    Y:所要划分的整体数据的结果;
    $ U: z4 q: d) K1 w/ {+ D; f- n; Gtest_size:测试集数据量在整体数据量中的占比(可以理解为X_test与X的比值);2 G: Y1 d7 [: ]
    random_state:①若不填或者填0,每次生成的数据都是随机,可能不一样。②若为整数,每次生成的数据都相同;5 E$ ~0 F6 V) h# O6 l: T2 c
    from sklearn.model_selection import train_test_split( g4 m5 A: e- Q( [
    X_train, X_test, Y_train, Y_test = train_test_split(X, Y, test_size=0.2, random_state=0)
    # H' \' _, w! {3 E, |% K& d. T4 T; O; ~
    2. LinearRegression()函数: N0 T# X0 n6 K3 W
    sklearn.linear_model包实现了广义线性模型,包括线性回归、Ridge回归、Bayesian回归等。LinearRegression是其中较为简单的线性回归模型。
    % N% i' ^) Z( F
    ; s% J+ I, S* K! ?, y/ w

    & x! r; ~  E/ U解释一下什么是回归:回归最简单的定义是,给出一个点集D,用一个函数去拟合这个点集,并且使得点集与拟合函数间的误差最小,如果这个函数曲线是一条直线,那就被称为线性回归,如果曲线是一条二次曲线,就被称为二次回归。
    ! k. b2 I$ x9 q3 b6 ]! g/ X# n$ n8 P1 K( x& I; E; Q! ^
    # [, Q' ]5 V7 w  o
    3. 散点图与折线统计图的绘制
    8 ?5 {5 Y# Y; O& `2 k1 e' kplt.scatter():绘画出数据的散点图
    , y1 z# q5 I1 Q: L6 k! }* Bplt.plot():绘画出依据模型(LinearRegression的线性回归模型)生成的直线
    . [- i3 G# S( n; r有问题请加我QQ32437105601 l2 |7 [$ O! `) O0 n9 [5 f7 H6 _

    " D8 G3 z3 {1 `! Y* t: D, k$ z% E7 F$ |2 [; a/ e$ v: K
    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-7-27 00:23 , Processed in 0.399304 second(s), 52 queries .

    回顶部