建模培训讲座第四讲(聚类分析及SAS).ppt

《建模培训讲座第四讲(聚类分析及SAS).ppt》由会员分享，可在线阅读，更多相关《建模培训讲座第四讲(聚类分析及SAS).ppt（118页珍藏版）》请在三一办公上搜索。

1、聚类分析,第一节 引言,第二节 相似性的量度,第三节 系统聚类分析法,第四节 K均值聚类分析,第五节 有序样品的聚类分析法,第六节 实例分析与计算机实现,第一节 引言,“物以类聚，人以群分”。对事物进行分类，是人们认识事物的出发点，也是人们认识世界的一种重要方法。因此，分类学已成为人们认识世界的一门基础科学。在生物、经济、社会、人口等领域的研究中，存在着大量量化分类研究。例如：在生物学中，为了研究生物的演变，生物学家需要根据各种生物不同的特征对生物进行分类。在经济研究中，为了研究不同地区城镇居民生活中的收入和消费情况，往往需要划分不同的类型去研究。在地质学中，为了研究矿物勘探，需要根据各种矿石

2、的化学和物理性质和所含化学成分把它们归于不同的矿石类。在人口学研究中，需要构造人口生育分类模式、人口死亡分类状况，以此来研究人口的生育和死亡规律。,但历史上这些分类方法多半是人们主要依靠经验作定性分类，致使许多分类带有主观性和任意性，不能很好地揭示客观事物内在的本质差别与联系；特别是对于多因素、多指标的分类问题，定性分类的准确性不好把握。为了克服定性分类存在的不足，人们把数学方法引入分类中，形成了数值分类学。后来随着多元统计分析的发展，从数值分类学中逐渐分离出了聚类分析方法。随着计算机技术的不断发展，利用数学方法研究分类不仅非常必要而且完全可能，因此近年来，聚类分析的理论和应用得到了迅速的发展

3、。聚类分析就是分析如何对样品（或变量）进行量化分类的问题。通常聚类分析分为Q型聚类和R型聚类。Q型聚类是对样品进行分类处理，R型聚类是对变量进行分类处理。,第二节 相似性的量度,一 样品相似性的度量,二 变量相似性的度量,一、样品相似性的度量,在聚类之前，要首先分析样品间的相似性。Q型聚类分析，常用距离来测度样品之间的相似程度。每个样品有p个指标（变量）从不同方面描述其性质，形成一个p维的向量。如果把n个样品看成p维空间中的n个点，则两个样品间相似程度就可用p维空间中的两点距离公式来度量。两点距离公式可以从不同角度进行定义，令dij 表示样品Xi与Xj的距离，存在以下的距离公式：1明考夫斯基距

4、离(5.1)明考夫斯基距离简称明氏距离，按的取值不同又可分成：,欧氏距离是常用的距离，大家都比较熟悉，但是前面已经提到，在解决多元数据的分析问题时，欧氏距离就显示出了它的不足之处。一是它没有考虑到总体的变异对“距离”远近的影响，显然一个变异程度大的总体可能与更多样品近些，既使它们的欧氏距离不一定最近；另外，欧氏距离受变量的量纲影响，这对多元数据的处理是不利的。为了克服这方面的不足，可用“马氏距离”的概念。,2马氏距离 设Xi与Xj是来自均值向量为，协方差为=（0）的总体G中的p维样品，则两个样品间的马氏距离为(5.5)马氏距离又称为广义欧氏距离。显然，马氏距离与上述各种距离的主要不同就是它考虑

5、了观测变量之间的相关性。如果各变量之间相互独立，即观测变量的协方差矩阵是对角矩阵，则马氏距离就退化为用各个观测指标的标准差的倒数作为权数的加权欧氏距离。马氏距离还考虑了观测变量之间的变异性，不再受各指标量纲的影响。将原始数据作线性变换后，马氏距离不变。,3兰氏距离(5.6)它仅适用于一切Xij0的情况，这个距离也可以克服各个指标之间量纲的影响。这是一个自身标准化的量，由于它对大的奇异值不敏感，它特别适合于高度偏倚的数据。虽然这个距离有助于克服明氏距离的第一个缺点，但它也没有考虑指标之间的相关性。,4距离选择的原则 一般说来，同一批数据采用不同的距离公式，会得到不同的分类结果。产生不同结果的原因

6、，主要是由于不同的距离公式的侧重点和实际意义都有不同。因此我们在进行聚类分析时，应注意距离公式的选择。通常选择距离公式应注意遵循以下的基本原则：（1）要考虑所选择的距离公式在实际应用中有明确的意义。如欧氏距离就有非常明确的空间距离概念。马氏距离有消除量纲影响的作用。（2）要综合考虑对样本观测数据的预处理和将要采用的聚类分析方法。如在进行聚类分析之前已经对变量作了标准化处理，则通常就可采用欧氏距离。（3）要考虑研究对象的特点和计算量的大小。样品间距离公式的选择是一个比较复杂且带有一定主观性的问题，我们应根据研究对象的特点不同做出具体分折。实际中，聚类分析前不妨试探性地多选择几个距离公式分别进行聚

7、类，然后对聚类分析的结果进行对比分析，以确定最合适的距离测度方法。,二、变量相似性的度量,多元数据中的变量表现为向量形式，在几何上可用多维空间中的一个有向线段表示。在对多元数据进行分析时，相对于数据的大小，我们更多地对变量的变化趋势或方向感兴趣。因此，变量间的相似性，我们可以从它们的方向趋同性或“相关性”进行考察，从而得到“夹角余弦法”和“相关系数”两种度量方法。1、夹角余弦两变量Xi与Xj看作p维空间的两个向量，这两个向量间的夹角余弦可用下式进行计算(5.7)显然，cos ij 1。,2相关系数相关系数经常用来度量变量间的相似性。变量Xi与Xj的相关系数定义为(5.8)显然也有，rij 1。

8、,无论是夹角余弦还是相关系数，它们的绝对值都小于1，作为变量近似性的度量工具，我们把它们统记为cij。当cij=1时，说明变量Xi与Xj完全相似；当cij近似于1时，说明变量Xi与Xj非常密切；当cij=0时，说明变量Xi与Xj完全不一样；当cij近似于0时，说明变量Xi与Xj差别很大。据此，我们把比较相似的变量聚为一类，把不太相似的变量归到不同的类内。在实际聚类过程中，为了计算方便，我们把变量间相似性的度量公式作一个变换为 dij=1 cij(5.9)或者 dij2=1 cij2(5.10)用表示变量间的距离远近，小则与先聚成一类，这比较符合人们的一般思维习惯。,第三节 系统聚类分析法,一

9、系统聚类的基本思想,二 类间距离与系统聚类法,三 类间距离的统一性,一、系统聚类的基本思想,系统聚类的基本思想是：距离相近的样品（或变量）先聚成类，距离相远的后聚成类，过程一直进行下去，每个样品（或变量）总能聚到合适的类中。系统聚类过程是：假设总共有n个样品（或变量），第一步将每个样品（或变量）独自聚成一类，共有n类；第二步根据所确定的样品（或变量）“距离”公式，把距离较近的两个样品（或变量）聚合为一类，其它的样品（或变量）仍各自聚为一类，共聚成n 1类；第三步将“距离”最近的两个类进一步聚成一类，共聚成n 2类；，以上步骤一直进行下去，最后将所有的样品（或变量）全聚成一类。为了直观地反映以上

10、的系统聚类过程，可以把整个分类系统画成一张谱系图。所以有时系统聚类也称为谱系分析。除系统聚类法外，还有有序聚类法、动态聚类法、图论聚类法、模糊聚类法等，限于篇幅，我们只介绍系统聚类方法。,二、类间距离与系统聚类法,在进行系统聚类之前，我们首先要定义类与类之间的距离，由类间距离定义的不同产生了不同的系统聚类法。常用的类间距离定义有8种之多，与之相应的系统聚类法也有8种，分别为最短距离法、最长距离法、中间距离法、重心法、类平均法、可变类平均法、可变法和离差平方和法。它们的归类步骤基本上是一致的，主要差异是类间距离的计算方法不同。以下用dij表示样品Xi与Xj之间距离，用Dij表示类Gi与Gj之间的

11、距离。,1.最短距离法定义类与之间的距离为两类最近样品的距离，即为(5.11)设类与合并成一个新类记为，则任一类与的距离为(5.12),最短距离法进行聚类分析的步骤如下：（1）定义样品之间距离，计算样品的两两距离，得一距离 阵记为D（0），开始每个样品自成一类，显然这时Dij=dij。（2）找出距离最小元素，设为Dpq，则将Gp和Gq合并成一个 新类，记为Gr，即Gr=Gp，Gq。（3）按（5.12）计算新类与其它类的距离。（4）重复（2）、（3）两步，直到所有元素。并成一类为 止。如果某一步距离最小的元素不止一个，则对应这些 最小元素的类可以同时合并。,【例5.1】设有六个样品，每个只测量一

12、个指标，分别是1，2，5，7，9，10，试用最短距离法将它们分类。（1）样品采用绝对值距离，计算样品间的距离阵D（0），见表5.1,表5.1,（2）D（0）中最小的元素是D12D561，于是将G1和G2合并成G7，G5和G6合并成G8，并利用（5.12）式计算新类与其它类的距离D（1），见表5.2,表5.2,（3）在D（1）中最小值是D34D482，由于G4与G3合并，又与G8合并，因此G3、G4、G8合并成一个新类G9，其与其它类的距离D（2），见表5.3,表5.3,（4）最后将G7和G9合并成G10，这时所有的六个样品聚为一类，其过程终止。上述聚类的可视化过程见图5.1所示，横坐标的刻度表

13、示并类的距离。这里我们应该注意，聚类的个数要以实际情况所定，其详细内容将在后面讨论。,图5.1 最短距离聚类法的过程,再找距离最小两类并类，直至所有的样品全归为一类为止。可以看出最长距离法与最短距离法只有两点不同：一是类与类之间的距离定义不同；另一是计算新类与其它类的距离所用的公式不同。,3.中间距离法最短、最长距离定义表示都是极端情况，我们定义类间距离可以既不采用两类之间最近的距离也不采用两类之间最远的距离，而是采用介于两者之间的距离，称为中间距离法。中间距离将类Gp与Gq类合并为类Gr，则任意的类Gk和Gr的距离公式为(14 0)(5.15)设DkqDkp，如果采用最短距离法，则Dkr=D

14、kp，如果采用最长距离法，则Dkr=Dkq。如图5.2所示，(5.15)式就是取它们（最长距离与最短距离）的中间一点作为计算Dkr的根据。,特别当=14，它表示取中间点算距离，公式为(5.16),图5.2 中间距离法,【例5.2】针对例5.1的数据，试用重心法将它们聚类。（1）样品采用欧氏距离，计算样品间的平方距离阵D2（0），见表5.4所示。,表5.4,（2）D2（0）中最小的元素是D212D2561，于是将G1和G2合并成G7，G5和G6合并成G8，并利用（5.18）式计算新类与其它类的距离得到距离阵D2（1），见表5.5：其中，其它结果类似可以求得,（3）在D2（1）中最小值是D2344

15、，那么G3与G4合并一个新类G9，其与与其它类的距离D2（2），见表5.6：,表5.6,（4）在中最小值是12.5，那么与合并一个新类，其与与其它类的距离，见表5.7：,表5.7,（5）最后将G7和G10合并成G11，这时所有的六个样品聚为一类，其过程终止。上述重心法聚类的可视化过程见图5.3所示，横坐标的刻度表示并类的距离。,图5.3 重心聚类法的过程,6.可变类平均法由于类平均法中没有反映出Gp和Gq之间的距离Dpq的影响，因此将类平均法进一步推广，如果将Gp和Gq合并为新类Gr，类Gk与新并类Gr的距离公式为：（5.22）其中是可变的且 1，称这种系统聚类法为可变类平均法。,8.离差平方

16、和法该方法是Ward提出来的，所以又称为Ward法。该方法的基本思想来自于方差分析，如果分类正确，同类样品的离差平方和应当较小，类与类的离差平方和较大。具体做法是先将n个样品各自成一类，然后每次缩小一类，每缩小一类，离差平方和就要增大，选择使方差增加最小的两类合并，直到所有的样品归为一类为止。设将n个样品分成k类G1，G2，Gk，用Xit表示Gt中的第I个样品，nt表示Gt中样品的个数，是Gt的重心，则Gt的样品离差平方和为,这种系统聚类法称为离差平方和法或Ward方法。下面论证离差平方和法的距离递推（5.26）式。,由于,三、类间距离的统一性,上述八种系统聚类法的步骤完全一样，只是距离的递推

17、公式不同。兰斯（Lance）和威廉姆斯（Williams）于1967年给出了一个统一的公式。(5.28)其中ap、aq、是参数，不同的系统聚类法，它们取不同的数，详见表5.8。这里应该注意，不同的聚类方法结果不一定完全相同，一般只是大致相似。如果有很大的差异，则应该仔细考查，找到问题所在；另外，可将聚类结果与实际问题对照，看哪一个结果更符合经验。,表5.8 系统聚类法参数表,第四节 K均值聚类分析,系统聚类法需要计算出不同样品或变量的距离，还要在聚类的每一步都要计算“类间距离”，相应的计算量自然比较大；特别是当样本的容量很大时，需要占据非常大的计算机内存空间，这给应用带来一定的困难。而K均值法

18、是一种快速聚类法，采用该方法得到的结果比较简单易懂，对计算机的性能要求不高，因此应用也比较广泛。K均值法是麦奎因（MacQueen，1967）提出的，这种算法的基本思想是将每一个样品分配给最近中心（均值）的类中，具体的算法至少包括以下三个步骤：1将所有的样品分成K个初始类；2通过欧氏距离将某个样品划入离中心最近的类中，并对获得样品与失去样品的类，重新计算中心坐标；3重复步骤2，直到所有的样品都不能再分配时为止。,K均值法和系统聚类法一样，都是以距离的远近亲疏为标准进行聚类的，但是两者的不同之处也是明显的：系统聚类对不同的类数产生一系列的聚类结果，而K均值法只能产生指定类数的聚类结果。具体类数的

19、确定，离不开实践经验的积累；有时也可以借助系统聚类法以一部分样品为对象进行聚类，其结果作为K均值法确定类数的参考。下面通过一个具体问题说明K均值法的计算过程。,【例5.3】假定我们对A、B、C、D四个样品分别测量两个变量和得到结果见表5.9。试将以上的样品聚成两类。,表5.9 样品测量结果,第一步：按要求取K=2，为了实施均值法聚类，我们将这些样品随意分成两类，比如（A、B）和（C、D），然后计算这两个聚类的中心坐标，见表5.10所示。表5.10中的中心坐标是通过原始数据计算得来的，比如（A、B）类的，等等。,表5.10 中心坐标,第二步：计算某个样品到各类中心的欧氏平方距离，然后将该样品分配

20、给最近的一类。对于样品有变动的类，重新计算它们的中心坐标，为下一步聚类做准备。先计算A到两个类的平方距离：由于A到（A、B）的距离小于到（C、D）的距离，因此A不用重新分配。计算B到两类的平方距离：,由于B到（A、B）的距离大于到（C、D）的距离，因此B要分配给（C、D）类，得到新的聚类是（A）和（B、C、D）。更新中心坐标如表5.11所示。,表5.11 更新后的中心坐标,第三步：再次检查每个样品，以决定是否需要重新分类。计算各样品到各中心的距离平方，得结果见表5.12。到现在为止，每个样品都已经分配给距离中心最近的类，因此聚类过程到此结束。最终得到K=2的聚类结果是A独自成一类，B、C、D聚

21、成一类。,表5.12 样品聚类结果,第五节 有序样品的聚类分析法,一 有序样品可能的分类数目,二 费希尔最优求解法,三 一个典型例子,以上的系统聚类和K均值聚类中，样品的地位是彼此独立的，没有考虑样品的次序。但在实际应用中，有时样品的次序是不能变动的，这就产生了有序样品的聚类分析问题。例如对动植物按生长的年龄段进行分类，年龄的顺序是不能改变的，否则就没有实际意义了；又例如在地质勘探中，需要通过岩心了解地层结构，此时按深度顺序取样，样品的次序也不能打乱。如果用X（1），X（2），X（n）表示n个有序的样品，则每一类必须是这样的形式，即X（i），X（i+1)，X（j），其中1 r n，且j n，简

22、记为Gi=i，i+1，j。在同一类中的样品是次序相邻的。这类问题称为有序样品的聚类分析。,一、有序样品可能的分类数目,n个有序样品分成k类，则一切可能的分法有 种。实际上，n个有序样品共有（n 1）个间隔，分成k类相当于在这（n 1）个间隔中插入k 1根“棍子”。由于不考虑棍子的插入顺序，是一个组合问题，共有 种插法。图5.4 有序样品的分类法这就是n个有序样品分成k类的一切可能分法。因此，对于有限的n和k，有序样品的所有可能分类结果是有限的，可以在某种损失函数意义下，求得最优解。所以有序样品聚类分析又称为最优分割，该算法是费希尔（Fisher）最先提出来的，故也称之为费希尔最优求解法。,二、

23、费希尔最优求解法,这里需要注意，若要寻找将n个样品分为k类的最优分割，则对于任意的j（k j n）,先将前面j 1个样品最优分割为k 1类，得到p（j 1，k 1），否则从j到n这最后一类就不可能构成k类的最优分割，参见图5.6。再考虑使Lb（n，k）最小的j，得到p（n，k）。因此我们得到费希尔最优求解法的递推公式为（5.23）,图5.6 最优分割,三、一个典型例子,【例5.4】为了了解儿童的生长发育规律，今随机抽样统计了男孩从出生到11岁每年平均增长的重量数据表5.13，试问男孩发育可分为几个阶段？在分析这是一个有序样品的聚类问题时，我们通过图形可以看到男孩增重随年龄顺序变化的规律，从图5

24、.6中发现男孩发育确实可以分为几个阶段。,表5.13 111岁儿童每年平均增长的重量,图5.7 儿童成长阶段分析,下面通过有序样品的聚类分析确定男孩发育分成几个阶段较合适。步骤如下：,表5.14 直径 D（i，j）,（3）分类个数的确定。如果能从生理角度事先确定k当然最好；有时不能事先确定k时，可以从Lp（l，k）随k的变化趋势图中找到拐点处，作为确定k的根据。当曲线拐点很平缓时，可选择的k很多，这时需要用其它的办法来确定，比如均方比和特征根法，限于篇幅此略，有兴趣的读者可以查看其它资料。本例从表5.15中的最后一行可以看出k=3，4处有拐点，即分成3类或4类都是较合适的，从图5.8中可以更明

25、显看出这一点。,第六节 实例分析与计算机实现,一 在SAS中利用系统聚类法进行聚类分析,二 在SAS中利用K均值法进行聚类分析,例,为了研究辽宁等5省1991年城镇居民生活消费情况的分布规律，根据调查资料做类型分类，用最短距离做类间分类。数据如下：,将每一个省区视为一个样本，先计算5个省区之间的欧式距离，用D0表示距离矩阵（对称阵，故给出下三角阵）因此将3.4合并为一类，为类6，替代了3、4两类类6与剩余的1、2、5之间的距离分别为：d(3,4)1=min(d31,d41)=min(13.80,13.12)=13.12 d(3,4)2=min(d32,d42)=min(24.63,24.06)

26、=24.06 d(3,4)5=min(d35,d45)=min(3.51,2.21)=2.21,得到新矩阵合并类6和类5，得到新类7类7与剩余的1、2之间的距离分别为：d(5,6)1=min(d51,d61)=min(12.80,13.12)=12.80 d(5,6)2=min(d52,d62)=min(23.54,24.06)=23.54,得到新矩阵合并类1和类2，得到新类8此时，我们有两个不同的类：类7和类8。它们的最近距离d(7,8)=min(d71,d72)=min(12.80,23.54)=12.80,得到矩阵最后合并为一个大类。这就是按最短距离定义类间距离的系统聚类方法。最长距离法

27、类似！,3.重心法(CENtroid method),4.类平均法(AVErage method),中间距离,5.离差平方和法(WARD),基本思想来源于方差分析。它认为：如果分类正确，同类间的类差平方和应较小，类与类之间的离差平方和应较大.具体做法是，先将n个样本分成一类,然后每次缩小一类,每缩小一类离差平方和就要增大.,离差平方和法(WARD),系统聚类方法的统一,系统聚类法参数表,类的个数的确定,由适当的阈值确定；根据数据点的散布直观地确定类的个数；根据统计量确定分类个数；,类的个数的确定,根据谱系图确定分类个数的准则：各类重心间的距离必须很大；类中保包含的元素不要太多；类的个数必须符合

28、实际应用；如果采用几种不同的聚类方法处理，则在各种聚类图中应该发现相同的类。,四、系统聚类的参数选择聚类类别：统计图：树型谱系图 冰柱谱系图聚类方法1.Between-groups linkage 类间平均法 两类距离为两类元素两两之间平均平方距离2.Within-groups linkage 类内平均法两类距离为合并后类中可能元素两两之间平均平方距离3.Nearest neighbor 最短距离法 4.Furthest neighbor 最长距离法5.Centroid clustering 重心法(欧式距离)6.Median clustering 中间距离法(欧式距离)7.Ward Meth

29、od 离差平方法(欧式距离),五、系统聚类法之例：地区按经济效益分类某年全国28个省区经济效益指标表，演示,类间平均法,类内平均法,最短距离法,最长距离法,重心法,中间距离法,离差平方法,六、快速聚类法,快速聚类法的聚类数由用户指定，分类是唯一的。1.分类数2.聚类方法:method:iterate and clussify clussify only3.聚类中心:centers4.迭代次数:iterate5.保存分类结果:save,8.2 判别分析8.2.1 判别分析的基本概念8.2.2 SAS中作判别分析的过程8.2.3 分析实例,8.2.1 判别分析的基本概念1.距离判别法 距离判别的基

30、本思想是：样品和哪个总体的距离最近，就判断它属于哪个总体。(1)两总体情况 设有两个总体G1和G2，若定义样品x到G1和G2的距离分别为d2(x,G1)和d2(x,G2)，则距离判别的判别规则是：即未知所属总体的样品x离哪个总体较近，就判x属于哪个总体。,当总体G1和G2为正态总体，距离选用马氏距离，公式为：式中，1，2，1，2分别为总体G1和G2的均值和协差阵。当1=2=时，可以证明：d2(x,G1)d2(x,G2)=2x(1+2)/2 1(12),令于是判别规则可表示为 称W(x)为判别函数。由于它是x的线性函数，因此又称为线性判别函数。线性判别的应用最为广泛。如果1与2不等，则判别函数W

31、(x)为：W(x)=d2(x，G1)d2(x，G2)=即W(x)是x的二次函数。,，,(2)多总体情况 设有m个总体：G1，G2，Gm（m 2），它们的均值、协差阵分别为i，i，(i=1，2，m)。对任意给定的样品x，要判断它来自哪个总体。多总体按距离最近的准则对x进行判别归类时，首先计算样品x到m个总体的马氏距离di2(x)（i=1,2,m），然后进行比较，把x判归距离最小的那个总体。即若dh2(x)=mindi2(x)|i=1,2,m，则x Gh。计算马氏距离di2(x)（i=1,2,m）时，类似地可考虑1=2=m或i不全相等的两种情况，并用样本统计量作为i和i的估计.,2.费歇尔判别（F

32、isher）下面以两个总体为例说明费歇尔判别的思想。设有两个总体G1、G2，其均值分别为1和2，协方差阵分别1和2，并假定1=2=，考虑线性组合：y=Lx。通过寻求合适的L向量，使得来自两个总体的数据间的距离较大，而来自同一个总体数据间的差异较小。为此，可以证明，当选L=c1(12)，其中c 0时，所得的投影即满足要求。从而称c=1时的线性函数：y=Lx=(1 2)1x为费歇尔线性判别函数。,其判别规则为：其中，m为两个总体均值在投影方向上的中点，即当1、2和未知时，可由总体G1和G2中分别抽出n1和n2个样品，计算相应的样本均值和协方差阵作为1、2和的估计。,3.贝叶斯判别（Bayes）(1

33、)贝叶斯判别法 设有m个总体G1,Gm，假定它们各自的分布密度分别为f1(x)，f2(x)，fm(x)，各自的先验概率（可以根据经验事先给出或估出）分别为q1,q2,qm，显然 贝叶斯判别的方法是：当抽取了一个未知总体的样品x，要判断它属于哪个总体，可用著名的贝叶斯公式计算x属于第j个总体的后验概率：,当 时，判断x属于第h个总体。或者计算按先验概率加权的误判平均损失：然后再比较这m个误判平均损失的h1(x)，h2(x)，hm(x)的大小，选取其中最小的，就可以判定样品x来自该总体。上式中C(j|i)为假定本来属于Gi的样品被判为属于Gj时造成的损失。当然C(i|i)=0，C(j|i)0（i，

34、j=1，m）。,在实际问题中，错判的损失可以给出定性的分析，但很难用数值来表示，但应用贝叶斯判别的方法时，要求定量给出C(j|i)，C(j|i)的赋值常用以下两种方法：根据经验人为赋值；假定各种错判的损失都相等。,(2)错判概率 当样品xGi，用判别法D判别时，把x判归Gj(ij)，出现错判。用P(j|I;D)（或简记为P(j|i)）表示实属Gi的样品错判为Gj的概率，错判概率的估计方法有下面几种：利用训练样本作为检验集，即用判别方法对已知样品进行回判，统计错判的个数以及错判的比率，作为错判率的估计。此法得出的估计一般较低。当训练样本的大小足够大时，可留出一些已知类别的样品不参加建立判别准则，

35、而是作为检验集，并把错判的比率作为错判率的估计。此法当检验集较小时估计的方差大。,设一法（或称为交叉确认法），每次留出一个已知类别的样品，而用其他n 1个样品建立判别准则，然后对留出的这一个已知类别的样品进行判别归类。对训练样本中n个样品逐个处理后把错判的比率作为错判率的估计。(3)广义平方距离判别法 在正态总体的假定下，按贝叶斯判别的思想，在错判造成的损失认为相等的情况下得到的判别函数其实就是马氏距离判别在考虑先验概率及协方差阵是否相等情况下的推广，故在SAS的DISCRIM过程中称为广义平方距离判别法。,4.逐步判别(1)各变量判别能力的检验 当检验k个类的均值向量是否全都相等(即检验H0

36、：1=2=k)时，否定了这个假设H0(即表明各总体的均值向量有显著性差异)，也并不能保证其各分量的均值有显著差异，若第i个分量间没有显著差异时，说明相应的变量Xi对判别分类不起作用，应该剔除。关于各变量判别能力的检验问题是筛选判别变量的理论基础，也是逐步判别的理论基础。,变量判别能力的度量通常采用删去该变量后考察判别能力的变化，即考察该变量对区分k个类是否能提供更多的附加信息，然后由附加信息构造F统计量进行检验。利用F统计量对假设H0(i)(第i个变量在k个总体中的均值相等)作统计检验。若否定H0(i)，表示变量Xi对区分k个总体的判别能力是显著的(在显著水平下)。否则，变量Xi对区分k个总体

37、的判别能力不能提供附加信息，这个变量应剔除。,(2)逐步判别法的基本思想 前面讨论了用全部m个变量：X1,X2,Xm来建立判别函数对样品进行判别归类的几种方法。在这m个变量中，有的变量对区分k个总体的判别能力可能很强，有的可能很弱。如果不加区别地把m个变量全部用来建立判别函数，则必增加大量的计算，还可能因为变量间的相关性引起计算上的困难(病态或退化等)及计算精度的降低。另一方面由于一些对区分k个总体的判别能力很小的变量的引入，产生干扰，致使建立的判别函数不稳定，反而影响判别效果，因此自然提出一个变量的选择问题。即如何从m个变量中挑选出对区分k个总体有显著判别能力的变量，来建立判别函数，用以判别

38、归类。,逐步判别的基本思想和逐步回归是类似的。逐个引入变量，每次把一个判别能力最强的变量引入判别式，每引入一个新变量，对判别式中的老变量逐个进行检验，如其判别能力因新变量的引入而变得不显著，应把它从判别式中剔除。这种通过逐步筛选变量使得建立的判别函数中仅保留判别能力显著的变量的方法，就是逐步判别法。,(3)逐步判别法的基本步骤 1)逐步筛选变量：根据各变量对区分k个总体的判别能力的大小，按基本思想所介绍的过程来筛选变量。SAS中的STEPDISC过程专用于筛选变量子集。该过程利用向前选入、向后剔除或逐步筛选的方法来选择区分k个总体的最佳变量子集。2)判别归类：对已选出变量子集，使用以上几节介绍

39、的判别方法对样品进行判别归类。,8.2.2 SAS中作判别分析的过程1.DISCRIM过程简介 DISCRIM过程根据一个分类变量和若干数值变量的数据计算出各种判别函数(判别准则)，根据这个判别函数，再将该批数据或其他数据中的观测分别归入已知类别中去。DISCRIM过程用以获得判别准则的数据称之为训练数据集(也称训练样本)。,实例分析：为反映员工心理因素方面的差异性，某研究者调研了一工厂36名员工的6 项心理因素，如下表所示。请根据这一分数对员工进行分类，并回答：员工分成几类比较理想？,练习题,1.为研究不同公司的运营特点，调查了 15 个公司的组织文化、组织氛围、领导角色和员工发展4个方面的

40、内容。请将这15个公司按照其各自的特点划分成4种类型。数据如下表所示。,2.请每个同学都参加“人性的哲学修订量表”问卷调查，然后分别对参加调查的同学和问卷项目进行类型划分。,人性的哲学修订量表 说明：本问卷由一系列关于态度的陈述组成。每一陈述只是代表人们的一种普遍观点，并没有对错之分。您可能同意某些观点却不同意另一些观点。我们感兴趣的是您同意或不同意的程度如何。仔细阅读每一陈述，然后针对每一陈述打一分数以表明你同意或不同意的程度。这些数值的意义如下：6完全同意 5部分同意 4略微同意 3略微不同意 2部分不同意 1完全不同意 回答这些问题时，往往第一印象最为恰当。阅读每一陈述，决定您是否同意以及程度如何。务请回答每一问题。如您发现这些数值均不能反映你的观点，请选择与你的观点最接近的数值。,问卷项目,参见SAS程序,