Java自学者论坛

 找回密码
 立即注册

手机号码,快捷登录

恭喜Java自学者论坛(https://www.javazxz.com)已经为数万Java学习者服务超过8年了!积累会员资料超过10000G+
成为本站VIP会员,下载本站10000G+会员资源,会员资料板块,购买链接:点击进入购买VIP会员

JAVA高级面试进阶训练营视频教程

Java架构师系统进阶VIP课程

分布式高可用全栈开发微服务教程Go语言视频零基础入门到精通Java架构师3期(课件+源码)
Java开发全终端实战租房项目视频教程SpringBoot2.X入门到高级使用教程大数据培训第六期全套视频教程深度学习(CNN RNN GAN)算法原理Java亿级流量电商系统视频教程
互联网架构师视频教程年薪50万Spark2.0从入门到精通年薪50万!人工智能学习路线教程年薪50万大数据入门到精通学习路线年薪50万机器学习入门到精通教程
仿小米商城类app和小程序视频教程深度学习数据分析基础到实战最新黑马javaEE2.1就业课程从 0到JVM实战高手教程MySQL入门到精通教程
查看: 521|回复: 0

[Python]基于K-Nearest Neighbors[K-NN]算法的鸢尾花分类问题解决方案

[复制链接]
  • TA的每日心情
    奋斗
    2024-4-6 11:05
  • 签到天数: 748 天

    [LV.9]以坛为家II

    2034

    主题

    2092

    帖子

    70万

    积分

    管理员

    Rank: 9Rank: 9Rank: 9

    积分
    705612
    发表于 2021-7-13 16:18:14 | 显示全部楼层 |阅读模式

      看了原理,总觉得需要用具体问题实现一下机器学习算法的模型,才算学习深刻。而写此博文的目的是,网上关于K-NN解决此问题的博文很多,但大都是调用Python高级库实现,尤其不利于初级学习者本人对模型的理解和工程实践能力的提升也不利于Python初学者实现该模型。

      本博文的特点:

        一 全面性地总结K-NN模型的特征、用途

        二  基于Python的内置模块,不调用任何第三方库实现

      博文主要分为四部分:

        基本模型(便于理清概念、回顾模型)

        对待解决问题的重述

        模型(算法)和评价(一来,以便了解模型特点,为以后举一反三地应用作铺垫;二来,有利于以后快速复习)、

        编程实现(Code)。

      特别声明:

        1.劳动成果开源,未经同意博主(千千寰宇:http://cnblogs.com/johnnyzen),不得以任何形式转载、复制。

        2.如有纰漏或者其他看法,欢迎共同探讨~

    零 基本模型

      (本部分内容,均来源于引用[1],其原理讲解十分通俗易懂)

      ①K-近邻算法,即K-Nearest Neighbor algorithm,简称K-NN算法。单从名字来猜想,可以简单粗暴的认为是:K个最近的邻居,当K=1时,算法便成了最近邻算法,即寻找最近的那个邻居。

      ②所谓K-NN算法,即是给定一个训练数据集,对新的输入实例,在训练数据集中找到与该实例最邻近的K个实例(也就是K个邻居), 这K个实例的多数属于某个类,就把该输入实例分类到这个类中。

      ③实例

        猜猜看:有一个未知形状(绿色圆点),如何判断其是什么形状?

        问题:给这个绿色的圆分类?

        对噪声数据过于敏感。为了解决这个问题,我们可以把位置样本周边的多个最近样本计算在内,扩大参与决策的样本量,以避免个别数据直接决定决策结果。

        有两类不同的样本数据,分别用蓝色的小正方形和红色的小三角形表示,而图正中间的那个绿色的圆所标示的数据则是待分类的数据。

        如果K=3,判定绿色的这个待分类点属于红色的三角形一类。

        如果K=5,判定绿色的这个待分类点属于蓝色的正方形一类。

    一 问题

      题目:ML之k-NN:k-NN实现对150朵共三种花的实例的萼片长度、宽,花瓣长、宽数据统计,根据一朵新花的四个特征来预测其种类
      数据源:https://archive.ics.uci.edu/ml/machine-learning-databases/iris/iris.data
      数据源说明:https://archive.ics.uci.edu/ml/machine-learning-databases/iris/iris.names

    二 解决过程及模型评价

       一 审题/准备数据集
         1.1 明确问题基本模型,及涉及要素(特征值、有无标记、可考虑的基本算法模型):K-NN、分类、监督学习、有标记、鸢尾花分类
         1.2 准备数据集及其处理方法
          数据源:https://archive.ics.uci.edu/ml/machine-learning-databases/iris/iris.data
          数据源说明:https://archive.ics.uci.edu/ml/machine-learning-databases/iris/iris.names
       二 分析并设计解决思路(算法步骤)
         2.0 零 样本及样本空间、测试空间的程序表示
            设计样本的class类:Iris(鸢尾花类):
              便于将样本空间以数组化表示 + 便于后期训练对各样本的取值、标记、比对等操作
         2.1 一 读取数据集,装载样本集
            提取文件数据至计算机进程(程序)中,并恰当地表示样本、样本集、测试集
          文件I/O操作:
            open(filepath+filename,mode)、
            file.readline()、
            file.readlines()、
            file.next()、
          字符串处理:
            str.split():切割
            float(str)、int(str):数值转换
          注意:对数据集文件中无关数据的处理(跳行、文件转进制、跳列等过滤)
         2.2 二 训练测试样本
          算法描述:[以鸢尾花分类为例]
            遍历测试样本
              计算测试样本与已标记的样本的欧式距离
            对各欧式距离升序排序
            选择前K项的已样本作为一子集( 即 选择最近的K项邻居作为参照标准)
              遍历统计,已标记子集的花朵种类何种花朵数目种类最多
            设置当前测试样本的预测花朵种类为该种
            结束。
          注释:花的种类分别为:Iris-setosa、Iris-versicolor、Iris-virginica;共计3种。
         2.3 三 计算预测准确率
          rate = Count(测试验证集中,已标记值 == 预测值 的样本) / 测试验证集总数
         2.4 四 模型检验(未做)
       三 编程实现
        (见第三部分示范代码)
       四 总结模型[K-NN]
         针对本问题
          0 样本数据集的样本顺序排列,选拔初始的参照样本空间时,需要先将数据集乱序排序,再做选拔;否则,分类准确率往往趋于0
            原因:对未纳入训练样本对应种类的待测样本,几乎无真实参照集合,所以将导致预测结果差得离谱(准确率趋于0)。
         涉及因素
          0 特征向量Features
            影响算法的处理效率:特征向量数越多,计算量越大,处理时间越长。
            影响算法的结果正确性:K值越多,可参照指标越多,结果越趋于正确。
          1 K:K个最近邻居元素的控制
            控制待预测样本的参照边界
            影响算法的结果正确性:
            K=1时,选拔最近的元素的标记值作为预测值
          2 Weight权值:对参照邻居的选拔标准
             标准:欧式距离、切比雪夫距离、马氏距离(这个,目前不清楚)、最近者赋予更大权值(增加其重要性)
            控制参照样本对待测样本的重要性
            有利于选拔最贴近真实预测结果、最符合的样本
          3 已预测并标记的样本是否加入参照样本集,影响后续待测样本预测?
          4 样本数量与种类
             标准:种类全面、数量按比例丰富;在上前提下,样本越多越好
              否则,尤其是种类覆盖不全面时,可能严重影响对偏未标记种类的样本预测效果
            样本种类↑,数量↑,预测准确率↑,预测处理效率↓
         总结性评价
          1 监督学习、分类问题
          2 基于实例数据的非参数学习算法
            Input:已标记数据集中的K个最近的训练样本组成
            Output:判别待测样本之类型
          3 可用于非线性分类
          4 样本种类↑,数量↑,特征向量↑:预测准确率↑,预测处理效率↓
          5 性能与完美的准确性不能绝对两全
            总体处理的复杂度:O(n)
          6 对异常值不敏感
          7 总体来说,准确度高,无人为因素。
          8 K值大小对结果的准确性要视具体情况而定,人为调参优化
          9 样本数量与种类:
             标准:种类全面、数量按比例丰富;在上前提下,样本越多越好
            否则,尤其是种类覆盖不全面时,可能严重影响对偏未标记种类的样本预测效果
          10    已预测的样本是否加入参照样本空间,将影响后续待测样本对预测结果起多大作用是一个可考虑的问题
          11 其他:
            需要提前知道样本空间之所有种类
            预测样本时,计算量大,[不利于实时预测],偏胖服务器模式
          12 模型常见实践:
            手写数字识别系统
            鸢尾花分类
            爱情片与动作片分类
            约会网站匹配
            对新贷款用户的还款情况预测
            ========================
            文本分类
    三 编程实现(For Python)
      工程文件分为三部分:
         __init__.py【main()启动函数、核心算法】
         Iris.py【设计数据结构(类)、模块(职责分离)】
         file_handle.py【数据提取、文件处理】

     

        由于前面第二部分已经详细叙述,且代码中注释已经十分详细,便不在对代码进行解释,阅读注释便容易懂。

        __init__.py

    import random;
    import math;
    import Iris; # 自定义
    import file_handle; # 自定义
    
    #if __name__ == '__main__': #__name__ == '__main__'是Python的main函数入口
    def main(print_test=False,print_samples=False):
        follows = []; # 样本集空间(前sampleAmount项)+测试集空间(后 sampleAll-sampleAmount项)
        data = "";    # 样本数据(前sampleAmount项)+测试集空间(后 sampleAll-sampleAmount项)
        sampleAll = 150;
        sampleAmount = 100; # 标记样本集数目(剩余的便作为测试集)
        k = 5;
        test_print = print_test; 
        
        ########## 一 读取数据集,装载样本集
        ##### 1.1 加载数据集数据
        data = file_handle.read("./dataset/data.txt",1,'r');#1:忽略第一行
        # print(data);
        list = data.split('\n');
        i = 0;
        for line in list: # 如:line = "5.1,3.5,1.4,0.2,Iris-setosa"
            item = line.split(',');# 如:item = [5.1,3.5,1.4,0.2,Iris-setosa]
            label_species = item.pop();#移除最后一项:标记种类
            #print("[test] item:", item,"\tlabel_species:", label_species); # test
            follows.append(Iris.Iris(item,label_species));
            #print("[ ",i," ] ",follows.toString());
            i += 1;
            pass;
        random.shuffle(follows); # 【千万注意!!!】由于原数据集是有序的,如果不做乱序处理,预测结果会及其不理想(准确率,趋近于0),当然,这也是这一模型的缺陷之一
        ##### 1.2 选择前100项 作为已标记样本集
        #i = 0;
        #for i in range(sampleAmount):
        #    follows.setPredictSpecies(follows.label_species);
        #    pass;
        
        ########## 二 训练测试样本
        ##### 2.1 对101 - 150 项的测试集进行训练/预测
        ## 算法描述:
        ##     遍历测试样本
        ##        计算测试样本与已标记的样本的欧式距离
        ##     对各欧式距离升序排序
        ##     选择前K项的已样本作为一子集( 即 选择最近的K项邻居作为参照标准)
        ##        遍历统计,已标记子集的花朵种类何种花朵数目种类最多
        ##        设置当前测试样本的预测花朵种类为该种
        ##     结束。  
        ## 注释:花的种类分别为:Iris-setosa、Iris-versicolor、Iris-virginica;共计3种。
        offset = 0; # 测试空间偏移量:目的是为了将通过偏移量,增大原已标记样本空间的样本数量 即 使已预测的测试样本加入参照样本空间。
        for x in range(sampleAmount,sampleAll):# x:测试样本下标
            weights = [];# 对各欧式距离(权值)的升序排序列表
            for y in range(0,sampleAmount+offset):
                result = (math.sqrt( + \
                                   math.pow(follows[y].features[0] - follows[x].features[0],2) + \
                                   math.pow(follows[y].features[1] - follows[x].features[1],2) + \
                                   math.pow(follows[y].features[2] - follows[x].features[2],2) + \
                                   math.pow(follows[y].features[3] - follows[x].features[3],2)), y);# 存储x,方便排序后定位花朵
                #print("[test] weights[x]:", result);
                weights.append(result);
                pass;
            weights.sort(key = lambda item:item[0]); # 以各元组内第一首项[欧氏距离]为键,默认升序排序
            if test_print:
                for m in range(len(weights)): # 输出预测权重
                    print("[test] weights[",m,"]:",weights[m],"\tweights[",m,"][1] > ",weights[m][1],":",follows[weights[m][1]].toString());
            kinds_count = {"Iris-setosa":0,"Iris-versicolor":0,"Iris-virginica":0}; # 对已标记样本空间中各种花的数目统计作初始化
            for z in range(0,k): # 选择前K项的已样本作为一子集( 即 选择最近的K项邻居作为参照标准)
                if test_print:         
                    print("[test] 排名前",z+1,"项 follows[",z,"]:",follows[weights[z][1]].toString());
                label_species = follows[weights[z][1]].label_species;
                if(label_species == 'Iris-setosa'):
                    kinds_count["Iris-setosa"] += 1;
                elif label_species == 'Iris-versicolor':
                    kinds_count["Iris-versicolor"] += 1;
                elif label_species == 'Iris-virginica':
                    kinds_count['Iris-virginica'] += 1;
                else:
                    print("[ERROR:Unknown Species] follows[",weight[z][1],"]:",follows[weight[z][1]]);
                pass;
            result =  max(kinds_count.items(), key = lambda item:item[1]); # 取统计花类数字典中最大值对应的序列
            follows[x].predict_species = result[0]; # 标记预测种类
            if test_print:
                print("[test] 预测结果",result, " [follows[",x,"].predict_species]:", follows[x].predict_species); # test
            offset += 1;
            #for test in range(len(weights)): # 测试-输出距离权值结果
            #    print("[",test,"] weights:",weights[test][0],"\t",follows[weights[test][1]].toString());
            #    pass;
            pass;
            
        ########## 三 计算预测准确率
        rate = 0.0;
        i = 0;
        for i in range(sampleAmount,len(follows)):
            if(follows.label_species == follows.predict_species):
                rate += 1;
            else:
                print("[预测错误样本] follow[",i,"]:",follows.toString());
            pass;
        pass;
        rate = rate / (sampleAll - sampleAmount);
        print("预测准确率:",rate);
         
        if print_samples:
            for i in range(0,len(follows)):
                print(follows.toString());
                pass;
        pass;
    
    main(False,True);

        Iris.py

    'Iris module [class] '
    
    __author__ = 'Johnny Zen'
    
    class Iris:
        """
        Iris花(类)
        
        [Demo]
        iris = Iris([5.1,3.5,1.4,0.2],"Iris-setosa");
        print(iris.toString());
        iris.setPredictSpecies('Iris-setosa');
        print(iris.toString());
        print(iris.label_species);
        =======================
        [features][5.1, 3.5, 1.4, 0.2]    [label-species]Iris-setosa    [predict-species]None
        [features][5.1, 3.5, 1.4, 0.2]    [label-species]Iris-setosa    [predict-species]Iris-setosa
        Iris-setosa
        """
        features = [];
        label_species = None; # 标记种类
        predict_species = None; # 预测种类
        def __init__(self,features,label_species=None):
            if type(features).__name__ == 'list':
                self.features = features;
            else:
                self.features = list(features); # 此list方法对list对象执行将产生错误
                pass;
            for x in range(len(self.features)): # 列表内元素字符串转实数
                self.features[x] = float(self.features[x]);
            self.label_species = label_species;
            pass;
        def setPredictSpecies(self,predict_species=None):#设置预测种类
            self.predict_species = predict_species;
            pass;
        def toString(self):#与一般函数定义不同,类方法必须包含参数 self[第一个参数]
            return "[features]" + str(self.features) + "\t[label]" + str(self.label_species + "\t[predict]" + str(self.predict_species));
            pass;
        pass;

        file_handle.py

    "file_handle module [function]:read(filepath,ignore=0,mode='r')"
    
    def read(filepath,ignore=0,mode='r'):
        try:
            file = open(filepath,mode);
            ## file_content = file.read();
            file_content = '';
            i = 0;
            for i in range(0,ignore):
                file.readline();
                ##print(i);
                ##print(file.readline()); 
            for line in file.readlines():
                file_content += line;
        finally:
            if file:
                file.close();
            #print(file_content);
            return file_content;
        pass;

    四 参考文献

      [1] K-NN和K-Means算法

    哎...今天够累的,签到来了1...
    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 立即注册

    本版积分规则

    QQ|手机版|小黑屋|Java自学者论坛 ( 声明:本站文章及资料整理自互联网,用于Java自学者交流学习使用,对资料版权不负任何法律责任,若有侵权请及时联系客服屏蔽删除 )

    GMT+8, 2024-5-14 07:00 , Processed in 0.070047 second(s), 29 queries .

    Powered by Discuz! X3.4

    Copyright © 2001-2021, Tencent Cloud.

    快速回复 返回顶部 返回列表