贝叶斯实验报告

1、精选优质文档-倾情为你奉上HUNAN UNIVERSITY人工智能实验报告题 目 实验三：分类算法实验 学生姓名 匿名 学生学号 xx 专业班级 智能科学与技术1302班 指导老师 袁进 一实验目的1. 了解朴素贝叶斯算法的基本原理；2. 能够使用朴素贝叶斯算法对数据进行分类3. 了解最小错误概率贝叶斯分类器和最小风险概率贝叶斯分类器4. 学会对于分类器的性能评估方法二、实验的硬件、软件平台硬件：计算机软件：操作系统：WINDOWS 10 应用软件：C,Java或者Matlab相关知识点:贝叶斯定理：表示事件B已经发生的前提下，事件A发生的概率，叫做事件B发生下事件A的条件

2、概率，其基本求解公式为：贝叶斯定理打通了从P(A|B)获得P(B|A)的道路。直接给出贝叶斯定理：  朴素贝叶斯分类是一种十分简单的分类算法，叫它朴素贝叶斯分类是因为这种方法的思想真的很朴素，朴素贝叶斯的思想基础是这样的：对于给出的待分类项，求解在此项出现的条件下各个类别出现的概率，哪个最大，就认为此待分类项属于哪个类别。  朴素贝叶斯分类的正式定义如下：   1、设为一个待分类项，而每个a为x的一个特征属性。      2、有类别集合。      3、计

3、算。      4、如果，则。      那么现在的关键就是如何计算第3步中的各个条件概率。我们可以这么做：      1、找到一个已知分类的待分类项集合，这个集合叫做训练样本集。      2、统计得到在各类别下各个特征属性的条件概率估计。即      3、如果各个特征属性是条件独立的，则根据贝叶斯定理有如下推导：  

4、0;         因为分母对于所有类别为常数，因为我们只要将分子最大化皆可。又因为各特征属性是条件独立的，所以有：整个朴素贝叶斯分类分为三个阶段：      第一阶段: 准备工作阶段，这个阶段的任务是为朴素贝叶斯分类做必要的准备，主要工作是根据具体情况确定特征属性，并对每个特征属性进行适当划分，然后由人工对一部分待分类项进行分类，形成训练样本集合。这一阶段的输入是所有待分类数据，输出是特征属性和训练样本。这一阶段是整个朴素贝叶斯分类中唯一需要人工完成的阶段，其质量对整个过

5、程将有重要影响，分类器的质量很大程度上由特征属性、特征属性划分及训练样本质量决定。      第二阶段: 分类器训练阶段，这个阶段的任务就是生成分类器，主要工作是计算每个类别在训练样本中的出现频率及每个特征属性划分对每个类别的条件概率估计，并将结果记录。其输入是特征属性和训练样本，输出是分类器。这一阶段是机械性阶段，根据前面讨论的公式可以由程序自动计算完成。      第三阶段: 应用阶段。这个阶段的任务是使用分类器对待分类项进行分类，其输入是分类器和待分类项，输出是待分类项与类别的映射关系。这

6、一阶段也是机械性阶段，由程序完成。三、实验内容及步骤实验内容：A.利用贝叶斯算法进行数据分类操作，并统计其预测正确率,数据集：汽车评估数据集（learn作为学习集，test作为测试集合）B.随机产生10000组正样本和20000负样本高斯分布的数据集合（维数设为二维），要求正样本：均值为1;3，方差为2 0;0 2；负样本：均值为10;20，方差为10 0;0 10.先验概率按样本量设定为1/3和2/3.分别利用最小错误概率贝叶斯分类器和最小风险概率贝叶斯分类器对其分类。（假设风险程度正样本分错风险系数为0.6，负样本分错风险为0.4，该设定仅用于最小风险分析）相关概念:1.贝叶斯法则机器学习

7、的任务：在给定训练数据D时，确定假设空间H中的最佳假设。最佳假设：一种方法是把它定义为在给定数据D以及H中不同假设的先验概率的有关知识下的最可能假设。贝叶斯理论提供了一种计算假设概率的方法，基于假设的先验概率、给定假设下观察到不同数据的概率以及观察到的数据本身。2.先验概率和后验概率用P(h)表示在没有训练数据前假设h拥有的初始概率。P(h)被称为h的先验概率。先验概率反映了关于h是一正确假设的机会的背景知识如果没有这一先验知识，可以简单地将每一候选假设赋予相同的先验概率。类似地，P(D)表示训练数据D的先验概率，P(D|h)表示假设h成立时D的概率。机器学习中，我们关心的是P(h|D)，即给

8、定D时h的成立的概率，称为h的后验概率。3.贝叶斯公式贝叶斯公式提供了从先验概率P(h)、P(D)和P(D|h)计算后验概率P(h|D)的方法p(h|D)=P(D|H)*P(H)/P(D)P(h|D)随着P(h)和P(D|h)的增长而增长，随着P(D)的增长而减少，即如果D独立于h时被观察到的可能性越大，那么D对h的支持度越小。4.极大后验假设学习器在候选假设集合H中寻找给定数据D时可能性最大的假设h，h被称为极大后验假设（MAP）确定MAP的方法是用贝叶斯公式计算每个候选假设的后验概率，计算式如下:h_map=argmax P(h|D)=argmax (P(D|h)*P(h)/P(D)=ar

9、gmax P(D|h)*p(h) (h属于集合H)C.编写一个贝叶斯分类器。输入为：均指向量、先验概率、协方差矩阵、输入学习数据X,测试数据类别XLABEL,测试数据Y.输出为Y对应的类别。（选做）。四、实验步骤：1仔细阅读并了解实验数据集；2使用任何一种熟悉的计算机语言(比如C,Java或者matlab)实现朴素贝叶斯算法；3利用朴素贝叶斯算法在训练数据上学习分类器,训练数据的大小分别设置为：前100个数据，前200个数据，前500个数据，前700个数据，前1000个数据，前1350个数据；4利用测试数据对学习的分类器进行性能评估；5统计分析实验结果并上交实验报告；A源代码:package

10、Bayes;import java.io.BufferedReader;import java.io.FileInputStream;import java.io.FileNotFoundException;import java.io.FileReader;import java.io.IOException;import java.math.BigDecimal;import java.util.Vector;import Bayes.NaiveBayesTool.Property;public class NaiveBayesTool /* * 申明全局变量 / 前面是自己的属性，后面是

11、value的属性 * */int testTotal = 0;/ 训练样本数量int predictTotal = 0;/ 测试样本的数据int predictSucess = 0;/ 预测成功的数量/存储数量public int buy=new int44;/vhigh,high,med,lowpublic int maint=new int44;/vhigh,high,med,lowpublic int door=new int44;/2,3,4,5morepublic int person=new int34;/2,4,morepublic int lug_boot=new int34;

12、/small ,med,bigpublic int safe=new int34;/low,med,highpublic int ClassValues=new int4;/unacc,acc, good,vgoodString ClassValueName = "unacc", "acc", "good", "vgood" ;/存储概率float ClassValue_gl = new float4;/ unacc-0 acc-1 good-2 vgood-3float buy_Vlaue_gl = new fl

13、oat44; /前面是自己的属性，后面是value的属性float maint_Value_gl = new float44;float door_Value_gl = new float44;float person_Value_gl = new float34;float lugboot_Value_gl = new float34;float safe_Value_gl = new float34;/* * 主函数 */public static void main(String args) throws IOException NaiveBayesTool NBayes=new Nai

14、veBayesTool(); NBayes.ReadFile("learn.txt");/获取训练样本NBayes.Calculated_probability();/计算概率NBayes.TestData();/导入测试样本数据NBayes.show();/输出结果/*汽车属性类 * */public class Property/汽车有6个属性，每个属性都有几种类别，根据这6个属性来判断汽车的性价比Classvalue如何，public String buying;/vhigh,high,med,lowpublic String maint;/vhigh,high,me

15、d,lowpublic String doors;/2,3,4,5morepublic String persons;/2,4,morepublic String lug_boot;/small ,med,bigpublic String safety;/ low,med,highpublic String ClassValues;/unacc,acc, good,vgoodpublic String PredictResult = new String5;/ 记录预测结果public Property(String b,String m,String d,String p,String l,

16、String s,String c)buying=b;maint=m;doors=d;persons=p;lug_boot=l;safety=s;ClassValues=c;Vector<Property> Data=new Vector();/存储数据Vector<Property> DataTest=new Vector();/存储测试数据/* * 文件读写 获取训练样本 * */public void ReadFile(String filename) throws IOException BufferedReader br=new BufferedReader(

17、new FileReader(filename); String temp=null; temp=br.readLine(); String str = null; Property TempClass = null; while(temp!=null) str=temp.split(","); TempClass=new Property(str0,str1,str2,str3,str4,str5,str6); Statistics(TempClass);/统计个数 testTotal+; temp=br.readLine(); br.close();/* * 统计每一项

18、的个数vhigh,vhigh,2,2,small,low,unacc * */public void Statistics(Property car)for(int i=0;i<4;i+)if(car.ClassValues.equals(ClassValueNamei) ClassValuesi+;/vhigh,high,med,lowif(car.buying.equals("vhigh") buy0i+;else if(car.buying.equals("high") buy1i+;else if(car.buying.equals(&qu

19、ot;med") buy2i+;else buy3i+;/vhigh,high,med,lowif(car.maint.equals("vhigh") maint0i+;else if(car.maint.equals("high") maint1i+;else if(car.maint.equals("med") maint2i+;else maint3i+;/2,3,4,5moreif(car.doors.equals("2") door0i+;else if(car.doors.equals(&qu

20、ot;3") door1i+;else if(car.doors.equals("4") door2i+;else door3i+;/2,4,moreif(car.persons.equals("2") person0i+;else if(car.persons.equals("4") person1i+;else person2i+;/small ,med,bigif(car.lug_boot.equals("small") lug_boot0i+;else if(car.lug_boot.equals

21、("med") lug_boot1i+;else lug_boot2i+;/ low,med,highif(car.safety.equals("low") safe0i+;else if(car.safety.equals("med") safe1i+;else safe2i+;/*计算概率*/public void Calculated_probability()for(int i=0;i<ClassValues.length;i+)ClassValue_gli=(float)ClassValuesi/testTotal;f

22、or(int i=0;i<buy_Vlaue_gl.length;i+)for(int j=0;j<buy_Vlaue_gl0.length;j+)buy_Vlaue_glij=(float)buyij/ClassValuesj;maint_Value_glij=(float)maintij/ClassValuesj;door_Value_glij=(float)doorij/ClassValuesj;for(int i=0;i<person_Value_gl.length;i+)for(int j=0;j<person_Value_gl0.length;j+)pers

23、on_Value_glij=(float)personij/ClassValuesj;lugboot_Value_glij=(float)lug_bootij/ClassValuesj;safe_Value_glij=(float)safeij/ClassValuesj;/* * 获取测试数据 * */public void TestData() throws IOExceptionBufferedReader br=new BufferedReader(new FileReader("test.txt");String temp;temp=br.readLine();St

24、ring str = null; Property Car = null;while(temp!=null)str=temp.split(",");Car=new Property(str0,str1,str2,str3,str4,str5,str6);predictTotal+;Data.addElement(Car);calculate(Car);temp=br.readLine();/* * 对分类器进行性能测试，判断其成功率为多少 * param car */public void calculate(Property car)/ unacc,acc,good,vg

25、ood, P(yi)-ClassValueTotal_gl 、 P(x|yi)=low,vhigh,4,2,small,low,unacc 第一条float itemGl;/ 每一条的概率int b, m, d, p, l, s;b = m = d = p = l = s = -1;float MaxGl = 0;if (car.buying.equals("vhigh") b = 0;else if (car.buying.equals("high") b = 1;else if (car.buying.equals("med")b

26、 = 2;else b = 3;if (car.maint.equals("vhigh") m = 0;else if (car.maint.equals("high") m = 1;else if (car.maint.equals("med")m = 2;else m = 3;if (car.doors.equals("2") d = 0;else if (car.doors.equals("3") d = 1;else if (car.doors.equals("4")

27、d = 2;else d = 3;if (car.persons.equals("2") p = 0;else if (car.persons.equals("4")p = 1;else p = 2;if (car.lug_boot.equals("small") l = 0;else if (car.lug_boot.equals("med") l = 1;else l = 2;if (car.safety.equals("low") s = 0;else if (car.safety.equ

28、als("med") s= 1;else s = 2;int t = 0;/ 记录最大概率的下标int i;for (i = 0; i < ClassValue_gl.length; i+) / 计算在unacc,acc,good,vgood下的概率itemGl = 0;BigDecimal bigDecimal = new BigDecimal(Float.toString(ClassValue_gli),new BigDecimal(Float.toString(buy_Vlaue_glbi),new BigDecimal(Float.toString(maint

29、_Value_glmi),new BigDecimal(Float.toString(door_Value_gldi),new BigDecimal(Float.toString(person_Value_glpi),new BigDecimal(Float.toString(lugboot_Value_glli),new BigDecimal(Float.toString(safe_Value_glsi),;for (int j = 1; j < bigDecimal.length; j+) /加:a.add(b); 除:a.divide(b,2);/2为精度取值 bigDecimal

30、0 = bigDecimal0.multiply(bigDecimalj);/multiply乘itemGl = bigDecimal0.floatValue();car.PredictResulti = itemGl + "t"if (MaxGl < itemGl) MaxGl = itemGl;t = i;/ 判断结果是否正确if (car.ClassValues.equals(ClassValueNamet) / 预测结果和开始给定的结果相等car.PredictResulti = "true"predictSucess+; else car

31、.PredictResulti = "false"public void show()for (int i =0; i < predictTotal; i+) Property c = Data.get(i);for (int j = 0; j < c.PredictResult.length; j+)System.out.print(c.PredictResultj + "t");System.out.println();/ 分类器的准确率float t = (float) predictSucess / predictTotal;t=t*10000/100;System.out.pr