多元统计分析PPT优秀PPT.ppt

《多元统计分析PPT优秀PPT.ppt》由会员分享，可在线阅读，更多相关《多元统计分析PPT优秀PPT.ppt（45页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、1问题引入问题引入2思路点思路点拨拨3判判别别分析方法分析方法4DNA序列分序列分类问题类问题的求解的求解5.参参考文考文献献目目 录录 首先，我首先，我们来们来考考虑虑一下一下2000年年“网网易杯易杯”全全国国高校生高校生数学数学建模建模竞赛竞赛的的A题题是是关关于于“DNA序列分序列分类类”的的问题问题 1问题问题引入引入 人类基因组中的人类基因组中的DNA全序列是由全序列是由4个碱基个碱基A，T，C，G按确定依次排成的长约按确定依次排成的长约30亿的序列，毫无亿的序列，毫无疑问，这是一本记录着人类自身生老病死及遗传疑问，这是一本记录着人类自身生老病死及遗传进化的全部信息的进化的全部信息

2、的“天书天书”。但是，除了这四种。但是，除了这四种碱基外，人们对它所包含的内容知之甚少，如何碱基外，人们对它所包含的内容知之甚少，如何破译这部破译这部“天书天书”是二十一世纪最重要的任务之是二十一世纪最重要的任务之一。在这个目标中，探讨一。在这个目标中，探讨DNA全序列具有什么结全序列具有什么结构，由这构，由这4个字符排成的看似随机的序列中隐藏着个字符排成的看似随机的序列中隐藏着什么规律，又是解读这部天书的基础，是生物信什么规律，又是解读这部天书的基础，是生物信息学（息学（Bioinformatics）最重要的课题之一。）最重要的课题之一。虽然人类对这部虽然人类对这部“天书天书”知之甚少，但也

3、知之甚少，但也发觉了发觉了DNA序列中的一些规律性和结构。例序列中的一些规律性和结构。例如，在全序列中有一些是用于编码蛋白质如，在全序列中有一些是用于编码蛋白质的序列片段，即由这的序列片段，即由这4个字符组成的个字符组成的64种不种不同的同的3字符串，其中大多数用于编码构成蛋字符串，其中大多数用于编码构成蛋白质的白质的20种氨基酸。又例如，在不用于编种氨基酸。又例如，在不用于编码蛋白质的序列片段中，码蛋白质的序列片段中，A和和T的含量特殊的含量特殊多些，于是以某些碱基特殊丰富作为特征多些，于是以某些碱基特殊丰富作为特征去探讨去探讨DNA序列的结构也取得了一些结果。序列的结构也取得了一些结果。此

4、外，利用统计的方法还发觉序列的某些此外，利用统计的方法还发觉序列的某些片段之间具有相关性，等等。这些发觉让片段之间具有相关性，等等。这些发觉让人们信任，人们信任，DNA序列中存在着局部的和全局序列中存在着局部的和全局性的结构，充分发掘序列的结构对理解性的结构，充分发掘序列的结构对理解DNA全序列是特别有意义的。全序列是特别有意义的。作为探讨作为探讨DNA序列的结构的尝试，试对以下序列序列的结构的尝试，试对以下序列进行分类：进行分类：问题一：下面有问题一：下面有20个已知类别的人工制造的序列个已知类别的人工制造的序列（见附件（见附件1），其中序列标号），其中序列标号110 为为A类，类，11-2

5、0为为B类。请从中提取特征，构造分类方法，并用这类。请从中提取特征，构造分类方法，并用这些已知类别的序列，衡量你的方法是否足够好。些已知类别的序列，衡量你的方法是否足够好。然后用你认为满足的方法，对另外然后用你认为满足的方法，对另外20个未标明类个未标明类别的人工序列（标号别的人工序列（标号2140）进行分类，把结果）进行分类，把结果用序号（按从小到大的依次）标明它们的类别用序号（按从小到大的依次）标明它们的类别（无法分类的不写入）：（无法分类的不写入）：A类类 ；B类类 。问题二：请对问题二：请对 182个自然个自然DNA序序列（列（:/mcm.edu/mcm00/problems.htm）

6、进行）进行分类。它们都较长。用你的分类分类。它们都较长。用你的分类方法对它们进行分类，并给出分方法对它们进行分类，并给出分类结果。类结果。看了这道题，我们应当从何看了这道题，我们应当从何处入手呢，我们应当怎样进行分处入手呢，我们应当怎样进行分析呢析呢2思路点拨思路点拨细读全题细读全题对未知事物进行分类对未知事物进行分类 问题的问题的本本 质质对另外对另外20个未标明类别个未标明类别的的DNA序列进行分类序列进行分类 依据依据所给所给的的20个已个已知类知类别的别的DNA序列序列所供所供应的应的信息信息 对对182个自然个自然DNA序列进行分类序列进行分类 假如将每一个假如将每一个DNADNA序

7、列都看作样本，那序列都看作样本，那么该问题就进一步提炼成一个纯粹的数学么该问题就进一步提炼成一个纯粹的数学问题：设有两个总体（类）问题：设有两个总体（类）和和 ，其分布特征（来自各个总体的样本）已知，其分布特征（来自各个总体的样本）已知，对给定的新品对给定的新品 ，我们须要推断其属，我们须要推断其属于哪个总体（类）。于哪个总体（类）。对于上面的数学问题，可以用很多成熟对于上面的数学问题，可以用很多成熟的方法来解决，例如：的方法来解决，例如：（1 1）BPBP神经网络；神经网络；（2 2）聚类分析；）聚类分析；（3 3）判别分析；等等。）判别分析；等等。如何选取方法是建模过程中须要解决的另外一个

8、如何选取方法是建模过程中须要解决的另外一个问题：问题：BPBP神经网络是人工神经网络的一种，它通过对神经网络是人工神经网络的一种，它通过对训练样本的学习，提取样本的隐含信息，进而对新样训练样本的学习，提取样本的隐含信息，进而对新样本的类别进行预料。本的类别进行预料。BPBP神经网络可以用以解决上面的神经网络可以用以解决上面的DNADNA序列分类问题，但是，如何提取特征、如何提高网序列分类问题，但是，如何提取特征、如何提高网络的训练效率、如何提高网络的容错实力、如何建立络的训练效率、如何提高网络的容错实力、如何建立网络结构是能否成功解决网络结构是能否成功解决DNADNA序列分类问题的关键所在；序

9、列分类问题的关键所在；聚类分析和判别分析都是多元统计分析中的经典方法，聚类分析和判别分析都是多元统计分析中的经典方法，都可以用来将对象（或观测值）分成不同的集合或类都可以用来将对象（或观测值）分成不同的集合或类别，但是，聚类分析更侧重于别，但是，聚类分析更侧重于“探究探究”对象（或观测对象（或观测值）的自然分组方式，而判别分析则侧重于将未知类值）的自然分组方式，而判别分析则侧重于将未知类别的对象（或观测值）别的对象（或观测值）“归结归结”（或者说，安排）到（或者说，安排）到已知类别中。明显，判别分析更适合用来解决上面的已知类别中。明显，判别分析更适合用来解决上面的DNADNA序列分类问题。序列

10、分类问题。3判判别别分析方法分析方法 判别分析是用于判别样品所属类别的一种判别分析是用于判别样品所属类别的一种多元统计分析方法。判别分析问题都可以这样多元统计分析方法。判别分析问题都可以这样描述：设有描述：设有 个个 维的总体维的总体 ，其分布，其分布特征已知（如已知分布函数分别为特征已知（如已知分布函数分别为 或者已知来自各个总体的样本），对给定的一或者已知来自各个总体的样本），对给定的一个新样品个新样品 ，我们须要推断其属于哪个总体。，我们须要推断其属于哪个总体。一般来说，依据判别规则的不同，可以得到不一般来说，依据判别规则的不同，可以得到不同的判别方法同的判别方法，例如，距离判别、贝叶斯

11、，例如，距离判别、贝叶斯（Bayes）判别、费希尔（）判别、费希尔（Fisher）判别、逐步）判别、逐步判别、序贯判别等。这里，我们简洁介绍三个判别、序贯判别等。这里，我们简洁介绍三个常用的判别方法：距离判别、贝叶斯（常用的判别方法：距离判别、贝叶斯（Bayes）判别和费希尔（判别和费希尔（Fisher）判别。）判别。判判 别别 分分 析析 方方 法法1.1.距离判距离判别别2.2.贝贝叶斯（叶斯（BayesBayes）判）判别别3.3.费费希尔（希尔（FisherFisher）判）判别别4.4.判判别别分析模型的分析模型的 显显著性著性检验检验 3.1 距离判距离判别别 距离判别的基本思想：

12、样品距离判别的基本思想：样品 X X离哪个总体的距离最近，就推断离哪个总体的距离最近，就推断 X X 属于哪个总体。属于哪个总体。这里的这里的“距离距离”是通常意义下的是通常意义下的距离（欧几里得距离：在距离（欧几里得距离：在 m m 维欧几里维欧几里得空间得空间 R R 中，两点中，两点 与与 的欧几里得距离，也就的欧几里得距离，也就是通常我们所说的距离为是通常我们所说的距离为 ）吗？）吗？带着这个疑问，我们来考虑这样带着这个疑问，我们来考虑这样一个问题一个问题 ：设有两个正态总体设有两个正态总体 ，和和 ，现在有一个新的样品位于现在有一个新的样品位于 A A 处（参见图处（参见图1 1）从

13、图中不难看出：从图中不难看出：，是否，是否 A A 处处的样品属于总体的样品属于总体 呢？呢？图图 1 明显不是，因为从概率的角度来看，明显不是，因为从概率的角度来看，总体总体 的样本比较分散，而总体的样本比较分散，而总体 的样的样本则特别集中，因此本则特别集中，因此 处的样品属于总处的样品属于总体体 的概率明显大于属于总体的概率明显大于属于总体 的概率，的概率，也就是说，也就是说，处的样品属于总体处的样品属于总体 的的“可可能性能性”明显大于属于总体明显大于属于总体 的的“可能性可能性”！这也说明白用欧几里得距离来度量！这也说明白用欧几里得距离来度量样品到总体距离的局限性。因此，须要样品到总

14、体距离的局限性。因此，须要引入新的距离概念引入新的距离概念这就是下面给出这就是下面给出的马氏距离。的马氏距离。定义定义1 1（马氏距离）：设总体（马氏距离）：设总体 G 为为 m 维总体维总体 （m m 个因素或指标），其均个因素或指标），其均值向量为值向量为 （这里（这里 T T 表示表示转置），协方差阵为转置），协方差阵为 ，则样品，则样品 到总体到总体 G G 的的马氏距离定义为马氏距离定义为3.1.1 3.1.1 两总体的距离判别两总体的距离判别 先考虑两个总体（先考虑两个总体（）的状况。设有两个总体）的状况。设有两个总体 和和 ,和和 分别是分别是 和和 的协方差阵，的协方差阵，和和

15、 分别是分别是 和和 的均值。对于新的样品的均值。对于新的样品 ，须要推断它来自那个总体。，须要推断它来自那个总体。设来自设来自 （）的训）的训练样本为练样本为其中其中 表示来自哪个表示来自哪个总体，总体，表示来自表示来自总体总体 的样本量。的样本量。要推断新样品要推断新样品 来自哪个总体，一般的想法是分来自哪个总体，一般的想法是分别计算新样品到两个总体的马氏距离别计算新样品到两个总体的马氏距离 和和 ：假如假如 则判定则判定 ；反之，假如反之，假如 则判定则判定 ：即即 （1 1）A.A.时的判别方法时的判别方法其中其中 ，记记 为了得到更简洁的判别规则，我们下面为了得到更简洁的判别规则，我

16、们下面计算新样品到两个总体的马氏距离计算新样品到两个总体的马氏距离 和和 的差的差明显，判别规则（明显，判别规则（1 1）式等价于）式等价于 （2 2）通常，称通常，称 为判别系数向量称为判别系数向量称 为线性判别函数。为线性判别函数。留意判别准则（留意判别准则（1 1）式或者（）式或者（2 2）式将）式将 维空间维空间 划分成两部分：划分成两部分：和和 也即也即 。距离判别的实质就是：给出。距离判别的实质就是：给出空间空间 的一个划分的一个划分 和和 ，假如样品，假如样品 落入落入 之中，之中，则判定则判定 ；假如样品；假如样品 落入落入 之中，则判定之中，则判定 。当当 时，依据判别准则（

17、时，依据判别准则（1 1）式，我们同）式，我们同样的给出判别函数样的给出判别函数 为为相应的判别规则为相应的判别规则为 （3 3）B.B.时的判别方法时的判别方法 在实际应用中，总体的均值和协方差阵一般在实际应用中，总体的均值和协方差阵一般是未知的，我们所知道的仅仅是一组样本或者观是未知的，我们所知道的仅仅是一组样本或者观测值，在这种状况下，就须要利用数理统计的学测值，在这种状况下，就须要利用数理统计的学问，对问，对 进行估计。进行估计。利用已知样本，易得利用已知样本，易得 的无偏的无偏估计分别为估计分别为C.C.的估计的估计 对于多个总体的状况，可以类似于两个总体的对于多个总体的状况，可以类

18、似于两个总体的处理过程，我们给出如下的步骤：处理过程，我们给出如下的步骤：第一步：计算样品第一步：计算样品 到每个总体的马氏距离到每个总体的马氏距离 ；其次步：比较其次步：比较 的大小，将样品的大小，将样品 判为距离最小的那个总体。判为距离最小的那个总体。假如均值为：假如均值为：和协方差：和协方差：未知，可以类似两个总体的情形运用训练样原来进未知，可以类似两个总体的情形运用训练样原来进行估计。这里不再赘述。行估计。这里不再赘述。3.1.2 多多总总体的距离判体的距离判别别3.1.3 3.1.3 距离判别的不足距离判别的不足 距离判别方法简洁好用，距离判别方法简洁好用，简洁实现，并且结论的意义明

19、简洁实现，并且结论的意义明确。但是，距离判别没有考虑：确。但是，距离判别没有考虑：（1 1）各总体本身出现的可）各总体本身出现的可能性在距离判别中没有考虑；能性在距离判别中没有考虑；（2 2）错判造成的损失在距）错判造成的损失在距离判别中也没有考虑。离判别中也没有考虑。在很多状况下，不考虑上在很多状况下，不考虑上面的两种因素是不合理的。贝面的两种因素是不合理的。贝叶斯（叶斯（BayesBayes）判别方法克服）判别方法克服了距离判别的不足。了距离判别的不足。与前面距离判别方法不同的是：所谓贝与前面距离判别方法不同的是：所谓贝叶斯（叶斯（BayesBayes）判别，就是在考虑各总体）判别，就是在

20、考虑各总体的先验概率和错判损失的状况下，给出的先验概率和错判损失的状况下，给出空间空间 的一个划分：的一个划分：，使得运用此划分来判别归类时，所带来使得运用此划分来判别归类时，所带来的平均错判损失最小。的平均错判损失最小。3.2 贝贝叶斯判叶斯判别别贝叶斯（贝叶斯（BayesBayes）判别问题的数学描述为：设有）判别问题的数学描述为：设有 个个 维的总体维的总体 ，其密度函数分别为，其密度函数分别为 ，若已知这若已知这 个总体各自出现的概率（先验概率）个总体各自出现的概率（先验概率）为为 （，且，且 ），假设已知将原来属于总体），假设已知将原来属于总体 的样品错判为总的样品错判为总体体 所造

21、成的损失为所造成的损失为 。在这样的情形下，对于新的样品在这样的情形下，对于新的样品 ，须要推，须要推断它来自那个总体。断它来自那个总体。为了给出贝叶斯（为了给出贝叶斯（BayesBayes）判别准则，我们）判别准则，我们从下面几个方面来探讨：从下面几个方面来探讨：先先验概验概率表示率表示对对各各总总体的先知相体的先知相识识，或者，或者说说，事先，事先对对所探所探讨讨的的问题问题所具有的相所具有的相识识。一般。一般来说来说，先，先验概验概率率并并不不简洁获简洁获得，其更多得，其更多的的来来自于自于长长期累期累积积的的阅历阅历。先。先验概验概率可以通率可以通过过下面几下面几种赋值种赋值方法得到：

22、方法得到：基于基于阅历阅历或者或者历历史史资资料料进进行估行估计计；利用利用训练样训练样本中各本中各种样种样品所占的比例品所占的比例 作作为为 的的值值。其中。其中 表示第表示第 类总类总体的体的样样品品数数，是是总训练样总训练样本本数数，该该方法要求方法要求训练样训练样本是通本是通过随过随机抽机抽样样得到的；得到的；3.假定假定 。A.A.确定先验概率确定先验概率 这里用这里用 表示将原来属于总体表示将原来属于总体 的样品错判为的样品错判为总体总体 的概率，即误判概率，明显，依据概率的定义，的概率，即误判概率，明显，依据概率的定义，易得：易得：（）（4 4）事实上，（事实上，（4 4）式的几

23、何意义是很明显的，见图）式的几何意义是很明显的，见图2 2。B.B.确定错判概率确定错判概率 图图2 2表示的是两个正态总体的误判概率示意图。表示的是两个正态总体的误判概率示意图。误判概率的估计方法有以下几种：误判概率的估计方法有以下几种：（1 1）利用训练样本为检验集，用判别方法）利用训练样本为检验集，用判别方法对训练样本进行推断，统计误判的样本个数，对训练样本进行推断，统计误判的样本个数，计算误判样本占总样本的比例，并作为误判概计算误判样本占总样本的比例，并作为误判概率的估计值；率的估计值；（2 2）当训练样本足够大时，从训练样本中）当训练样本足够大时，从训练样本中预留一部分作为检验集，并

24、记录判错的比率，预留一部分作为检验集，并记录判错的比率，作为误判概率的估计值；作为误判概率的估计值；（3 3）运用舍一法：每次预留一个样原来检）运用舍一法：每次预留一个样原来检验，用剩下验，用剩下 的样本建立判别准则，循环检的样本建立判别准则，循环检验完全部训练样本，记录判错的比率，以此作验完全部训练样本，记录判错的比率，以此作为误判概率的估计值。为误判概率的估计值。C.C.确定错判损失确定错判损失 错判必定带来损失。错判必定带来损失。现实中，错判的损失一般来现实中，错判的损失一般来说很难定量给出。但是可以说很难定量给出。但是可以运用赋值法来确定：运用赋值法来确定：（1 1）依据阅历或者实依据

25、阅历或者实际问题的特征人为际问题的特征人为 确定；确定；（2 2）假设各种误判损假设各种误判损失都相等。失都相等。基于前面的探讨，运用概率学问：判别基于前面的探讨，运用概率学问：判别法法 将原来属于总体将原来属于总体 的样本错判给其它的样本错判给其它总体的平均损失为总体的平均损失为那么，关于先验概率的平均错判损失那么，关于先验概率的平均错判损失 为为 （5 5）假如能找到假如能找到 使得平均损失使得平均损失 达到最小，那么达到最小，那么 就称为贝叶斯（就称为贝叶斯（BayesBayes）判别的解。判别的解。D.D.确定平均损失确定平均损失定理定理1 1：设有：设有 个总体：个总体：，已知，已知

26、 的联合密度函数为的联合密度函数为 ，先验概率为，先验概率为 ，错判损失为，错判损失为 ，则贝，则贝叶斯（叶斯（BayesBayes）判别的解）判别的解为为 其中其中 （6 6）下面给出贝叶斯（下面给出贝叶斯（BayesBayes）判别的解的主要结论）判别的解的主要结论：3.3 3.3 费希尔判别费希尔判别 费希尔判别的基本思想：费希尔判别的基本思想：借助于方差分析的思想，利用借助于方差分析的思想，利用投影将投影将 元的数据投影到某一元的数据投影到某一个方向，使得投影后组与组之个方向，使得投影后组与组之间的差异尽可能的大，然后依间的差异尽可能的大，然后依据确定的判别规则对新样本的据确定的判别规

27、则对新样本的类别进行推断。类别进行推断。首先构造一个线性判别函首先构造一个线性判别函数数 （7 7）可见上面的函数将可见上面的函数将 元的数据元的数据投影到了一个方向，系数投影到了一个方向，系数 的确定原则是使得总体间的差的确定原则是使得总体间的差异最大，总体内部的离差最小。异最大，总体内部的离差最小。A.A.确定线性判别函数确定线性判别函数 设有设有 个个 元总体：元总体：，它们的均值为：，它们的均值为：；协方差为：协方差为：。在在 的条件下，有的条件下，有 令令这里这里确定确定 ，使得，使得 个总体间的差异个总体间的差异最大，总体内部的离差最小，则最大，总体内部的离差最小，则 应当达到最大

28、。应当达到最大。为了确保为了确保 的唯一性，不妨设的唯一性，不妨设 。因此，问题转。因此，问题转化为：在条件化为：在条件 约束下，求约束下，求 使得使得 式达式达到最大到最大这是大家特别熟悉的条件极值的问题。依据这是大家特别熟悉的条件极值的问题。依据拉格朗日乘子法：拉格朗日乘子法：求解得求解得 （8 8）由方程（由方程（8 8）第一式知，）第一式知，是是 的特征根，的特征根，是是相应的特征向量。可以证明相应的特征向量。可以证明 （的大小可的大小可以衡量判别函数以衡量判别函数 的判别效果，故称的判别效果，故称 为为判别效率）。设判别效率）。设 的非零特征根为的非零特征根为 ，相应的满足约束条件的

29、特征向量为相应的满足约束条件的特征向量为 ，明显，明显，取取 时时 达到最大。达到最大。B.B.确定判别规则确定判别规则 假设系数假设系数 已经求出，那么已经求出，那么线性判别函数线性判别函数 就完全确定下来，对于就完全确定下来，对于一个新的样品一个新的样品 ，可以构造下面的判别规则：，可以构造下面的判别规则：（9 9）3.4 3.4 判别分析模型的显著性检验判别分析模型的显著性检验 建立了判别分析模型以后还建立了判别分析模型以后还须要对模型进行评价，这就须要须要对模型进行评价，这就须要对判别分析模型的显著性进行检对判别分析模型的显著性进行检验，主要包括两个方面：判别效验，主要包括两个方面：判

30、别效果的检验和各变量判别实力的检果的检验和各变量判别实力的检验验11，22。所谓判别效果的检验。所谓判别效果的检验就是检验就是检验 k k 个总体的均值是否有个总体的均值是否有显著的差异，反映了接受判别分显著的差异，反映了接受判别分析模型的有效性问题；各变量判析模型的有效性问题；各变量判别实力的检验反映的是各指标别实力的检验反映的是各指标（因素）对判别分类所起的作用（因素）对判别分类所起的作用时候显著。具体的检验方法因为时候显著。具体的检验方法因为过于困难，这里不再赘述。过于困难，这里不再赘述。须要说明的是，作为多元统须要说明的是，作为多元统计分析中的经典方法，判别分析计分析中的经典方法，判别

31、分析在很多关于多元统计分析的教材在很多关于多元统计分析的教材中均有具体而深刻的论述。中均有具体而深刻的论述。4 4DNADNA序列分类问题的求解序列分类问题的求解 关于关于DNADNA序列分类问题的探讨序列分类问题的探讨和分析，我们在第和分析，我们在第1 1部分和第部分和第2 2部分部分已经作了具体的分析和探讨。这里，已经作了具体的分析和探讨。这里，我们将依据多元统计分析的学问建我们将依据多元统计分析的学问建立判别分析模型来求解立判别分析模型来求解DNADNA序列的序列的分类问题（这里只求解问题一）。分类问题（这里只求解问题一）。首先，须要提取每条序列的首先，须要提取每条序列的所蕴含的特征所蕴

32、含的特征因为将序列的全因为将序列的全部信息都作为指标（因素）来建立部信息都作为指标（因素）来建立判别模型是不行能的。判别模型是不行能的。A.A.特征的提取特征的提取 DNA DNA序列中所蕴含的信息是特别序列中所蕴含的信息是特别丰富的，因此，如何提取特征、提丰富的，因此，如何提取特征、提取什么特征是一个特别困难的问题取什么特征是一个特别困难的问题这个问题涉及到生物学的学问，这个问题涉及到生物学的学问，此处将不深化绽开探讨。我们考虑此处将不深化绽开探讨。我们考虑接受序列中接受序列中4 4个碱基个碱基A A，T T，C C，G G的的含量百分比作为含量百分比作为DNADNA序列的特征。序列的特征。

33、为了便于探讨，我们用为了便于探讨，我们用 表示碱表示碱基基A A在序列中所占的百分比；在序列中所占的百分比；表示表示碱基碱基T T在序列中所占的百分比；在序列中所占的百分比；表表示碱基示碱基C C在序列中所占的百分比；在序列中所占的百分比；表示碱基表示碱基G G在序列中所占的百分比。在序列中所占的百分比。因为，因为，+=1+=1，因此，因此 、和和 中只有三个变量是独立的，中只有三个变量是独立的，不失一般性，我们选取不失一般性，我们选取 、和和 为指标，以为指标，以 表示表示第第 个个DNADNA序列的特征向量，换句序列的特征向量，换句话说，话说，“完全完全”代表第代表第 个个DNADNA序列

34、。序列。正如前面第正如前面第2 2部分中所分析的：在此部分中所分析的：在此DNADNA序列分类序列分类问题中，共有两个已知的总体（类）问题中，共有两个已知的总体（类）（A A类）和类）和 （B B类），而且，训练样本共有类），而且，训练样本共有 个，其中，个，其中，前前1010个样本（记为个样本（记为 ）属于总体（类）属于总体（类）；后；后 10 10个样本（记为个样本（记为 ）属于总体）属于总体 （类）（类）。须要解决的问题是：对给定的新样品须要解决的问题是：对给定的新样品 ，我们，我们须要推断其属于哪个总体（类）。须要推断其属于哪个总体（类）。B.B.建立判别模型建立判别模型 这里分别建立

35、距这里分别建立距离判别和费希尔（离判别和费希尔（FisherFisher）判别模）判别模型。型。（1 1）距离判别模型）距离判别模型在距离判别模型中，首先须要估在距离判别模型中，首先须要估计总体（类）计总体（类）和和 的均值的均值 和协方差阵和协方差阵 ，运用前面，运用前面3.1.13.1.1小节的公式，可小节的公式，可以得到均值以得到均值 和协方差阵和协方差阵 的无偏估计。的无偏估计。然后，计算样本然后，计算样本 到总体（类）到总体（类）和和 的马氏距离的马氏距离 ；依据判别准则（依据判别准则（1 1），可得：），可得：若若 ，则判定，则判定 为为A A类；类；若若 ，则判定，则判定 为为B

36、 B类。类。（2 2）费希尔判别模型）费希尔判别模型 首先，确定判别函数首先，确定判别函数 利用费希尔利用费希尔(Fisher)(Fisher)判别准则，可以求出判别准则，可以求出 确定了判别函数，就确定了判别函数，就可以依据费希尔可以依据费希尔(Fisher)(Fisher)判别准则（判别准则（9 9）式，推断新）式，推断新的的DNADNA序列的类别问题了。序列的类别问题了。C.C.模型求解和软件实现模型求解和软件实现 判别分析的计算机实现是特别简洁的，运用判别分析的计算机实现是特别简洁的，运用MatlabMatlab和和SPSSSPSS等软件都可以进行判别分析。经过计等软件都可以进行判别分

37、析。经过计算，结果如下：（算，结果如下：（1 1）距离判别：运用上面算法对）距离判别：运用上面算法对已知类别的训练本进行判定，结果为：已知类别的训练本进行判定，结果为：A A类：类：1 1，2 2，3 3，5 5，6 6，7 7，8 8，9 9，1010B B类：类：4 4，1111，1212，1313，1414，1515，1616，1717，1818，1919，2020只有第只有第4 4个样本发生了误判，故正确率为个样本发生了误判，故正确率为95%95%。运用上面算法对运用上面算法对21-4021-40进行分类，结果为进行分类，结果为A A类：类：2222，2323，2525，2727，29

38、29，3030，3232，3333，3434，3535，36 36，3737B B类：类：2121，2424，2626，2828，3131，3838，3939，4040（2 2）费希尔判别）费希尔判别 运用上面算法对已知类别的训练本进行判定，结运用上面算法对已知类别的训练本进行判定，结果为：果为：A A类：类：1 1，2 2，3 3，5 5，6 6，7 7，8 8，9 9，1010 B B类：类：4 4，1111，1212，1313，1414，1515，1616，1717，18 18，1919，2020 只有第只有第4 4个样本发生了误判，故正确率为个样本发生了误判，故正确率为95%95%。运

39、用上面算法对运用上面算法对21-4021-40进行分类，结果为：进行分类，结果为：A A类：类：2222，2323，2525，2727，2929，34 34，3535，3636，3737B B类：类：2121，2424，2626，2828，3030，3131，3232，3333，38 38，39 39，4040参考文献参考文献1 1 高惠璇，应用多元统计高惠璇，应用多元统计分析，北京高校出版社，分析，北京高校出版社，2005 2005。2 2 张尧庭张尧庭,方开泰，多元方开泰，多元统计分析引论，科学出版社，统计分析引论，科学出版社，1982 1982。3 3 陈家鼎，孙山泽，李东陈家鼎，孙山泽，李东风，数理统计学讲义，高等风，数理统计学讲义，高等教教 育出版社，育出版社，19981998。4 4 梁之舜，邓集贤，杨维梁之舜，邓集贤，杨维权等，概率论与数理统计权等，概率论与数理统计（第（第 二版）（上），高等教二版）（上），高等教化出版社，化出版社，19981998。5 5 朱建平，应用多元统计朱建平，应用多元统计分析，科学出版社，分析，科学出版社，20062006。6 6 全国高校生数学建模竞全国高校生数学建模竞赛优秀论文汇编（赛优秀论文汇编（1999-1999-2000 2000），中国物价出版），中国物价出版社，社，20022002。