2022年主成分分析原理 .pdf

第七章主成分分析（一）教学目的通过本章的学习，对主成分分析从总体上有一个清晰地认识，理解主成分分析的基本思想和数学模型，掌握用主成分分析方法解决实际问题的能力。（二）基本要求了解主成分分析的基本思想，几何解释，理解主成分分析的数学模型，掌握主成分分析方法的主要步骤。（三）教学要点1、主成分分析基本思想，数学模型，几何解释2、主成分分析的计算步骤及应用（四）教学时数3 课时（五）教学内容1、主成分分析的原理及模型2、主成分的导出及主成分分析步骤在实际问题中，我们经常会遇到研究多个变量的问题，而且在多数情况下，多个变量之间常常存在一定的相关性。由于变量个数较多再加上变量之间的相关性，势必增加了分析问题的复杂性。如何从多个变量中综合为少数几个代表性变量，既能够代表原始变量的绝大多数信息，又互不相关，并且在新的综合变量基础上，可以进一步的统计分析，这时就需要进行主成分分析。第一节主成分分析的原理及模型一、主成分分析的基本思想与数学模型（一）主成分分析的基本思想主成分分析是采取一种数学降维的方法，找出几个综合变量来代替原来众多的变量，使这些综合变量能尽可能地代表原来变量的信息量，而且彼此之间互不相关。这种将把多个变量化为少数几个互相无关的综合变量的统计分析方法就叫做主成分分析或主分量分析。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 10 页 - - - - - - - - - 主成分分析所要做的就是设法将原来众多具有一定相关性的变量，重新组合为一组新的相互无关的综合变量来代替原来变量。通常，数学上的处理方法就是将原来的变量做线性组合，作为新的综合变量，但是这种组合如果不加以限制，则可以有很多，应该如何选择呢？如果将选取的第一个线性组合即第一个综合变量记为，自然希望它尽可能多地反1F映原来变量的信息，这里“信息 ” 用方差来测量，即希望越大，表示包含的信)(1FVar1F息越多。因此在所有的线性组合中所选取的应该是方差最大的，故称为第一主成分。1F1F如果第一主成分不足以代表原来个变量的信息，再考虑选取即第二个线性组合，为了p2F有效地反映原来信息，已有的信息就不需要再出现在中，用数学语言表达就是要求1F2F，称为第二主成分，依此类推可以构造出第三、四 第个主成分。0),(21FFCov2Fp（二）主成分分析的数学模型对于一个样本资料，观测个变量，个样品的数据资料阵为：ppxxxL,21nnpnnppxxxxxxxxxXLMMMMLL212222111211pxxxL,21其中：pjxxxxnjjjjLM, 2, 1,21主成分分析就是将个观测变量综合成为个新的变量（综合变量） ，即ppppppppppppxaxaxaFxaxaxaFxaxaxaFLLLL22112222121212121111简写为：pjpjjjxxxFL2211名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 10 页 - - - - - - - - - pj, 2, 1L要求模型满足以下条件：互不相关（，）jiFF ,jipji, 2, 1,L的方差大于的方差大于的方差，依次类推1F2F3F.,2 , 1122221pkaaakpkkLL于是，称为第一主成分，为第二主成分，依此类推，有第个主成分。主成分1F2Fp又叫主分量。这里我们称为主成分系数。ija上述模型可用矩阵表示为：，其中AXFpFFFFM21pxxxXM21pppppppaaaaaaaaaaaaAMLMMMMLL21212222111211称为主成分系数矩阵。A二、主成分分析的几何解释假设有个样品，每个样品有二个变量，即在二维空间中讨论主成分的几何意义。设n个样品在二维空间中的分布大致为一个椭园，如下图所示：n图 7.1 主成分几何解释图名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 10 页 - - - - - - - - - 将坐标系进行正交旋转一个角度，使其椭圆长轴方向取坐标，在椭圆短轴方向取1y坐标，旋转公式为2ycos)sin(sincos212211jjjjjjxxyxxynjL2, 1写成矩阵形式为：nnyyyyyyY2222111211LLXUxxxxxxnn2222111211cossinsincosLL其中为坐标旋转变换矩阵，它是正交矩阵，即有，即满足UIUUUU,1。1cossin22经过旋转变换后，得到下图的新坐标：图 7.2 主成分几何解释图新坐标有如下性质：21yy(1)个点的坐标和的相关几乎为零。n1y2y(2) 二维平面上的个点的方差大部分都归结为轴上，而轴上的方差较小。n1y2y和称为原始变量和的综合变量。由于个点在轴上的方差最大，因而将1y2y1x2xn1y二维空间的点用在轴上的一维综合变量来代替，所损失的信息量最小，由此称轴为第1y1y一主成分，轴与轴正交，有较小的方差，称它为第二主成分。2y1y名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 10 页 - - - - - - - - - 三、主成分分析的应用主成分概念首先是由Karl parson 在 1901 年引进，但当时只对非随机变量来讨论的。1933 年 Hotelling将这个概念推广到随机变量。特别是近年来，随着计算机软件的应用，使得主成分分析的应用也越来越广泛。其中，主成分分析可以用于系统评估。系统评估是指对系统营运状态做出评估，而评估一个系统的营运状态往往需要综合考察许多营运变量，例如对某一类企业的经济效益作评估，影响经济效益的变量很多，很难直接比较其优劣，所以解决评估问题的焦点是希望客观、科学地将一个多变量问题综合成一个单变量形式，也就是说只有在一维空间中才能使排序评估成为可能，这正符合主成分分析的基本思想。在经济统计研究中，除了经济效益的综合评价研究外，对不同地区经济发展水平的评价研究，不同地区经济发展竞争力的评价研究，人民生活水平、生活质量的评价研究，等等都可以用主成分分析方法进行研究。另外，主成分分析除了用于系统评估研究领域外，还可以与回归分析结合，进行主成分回归分析，以及利用主成分分析进行挑选变量，选择变量子集合的研究。第二节主成分的导出及主成分分析的步骤一、主成分的导出根据主成分分析的数学模型的定义，要进行主成分分析，就需要根据原始数据，以及模型的三个条件的要求，如何求出主成分系数，以便得到主成分模型。这就是导出主成分所要解决的问题。1、根据主成分数学模型的条件要求主成分之间互不相关，为此主成分之间的协差阵应该是一个对角阵。即，对于主成分，AXF其协差阵应为，AXAXAXAXAXVarFVar)()()()(=pO212、设原始数据的协方差阵为，如果原始数据进行了标准化处理后则协方差阵等于V相关矩阵，即有，名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 10 页 - - - - - - - - - XXRV3、再由主成分数学模型条件和正交矩阵的性质，若能够满足条件最好要求为A正交矩阵，即满足IAA于是，将原始数据的协方差代入主成分的协差阵公式得AARAXAXFVar)(AARAAR展开上式得pppppppppppppppppppaaaaaaaaaaaaaaaaaarrrrrrrrrOLMMMMLLLMMMMLLLMMMMLL21212221212111212221212111212222111211展开等式两边，根据矩阵相等的性质，这里只根据第一列得出的方程为：0)(0)(0111221111212122112111121211111ppppppppparararararararararLLLLLLL为了得到该齐次方程的解，要求其系数矩阵行列式为0，即0121212221112111pppppprrrrrrrrrLMMMMLL01IR显然，是相关系数矩阵的特征值，是相应的特征向量。1paaaa112111,L根据第二列、第三列等可以得到类似的方程，于是是方程i0IR名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页，共 10 页 - - - - - - - - - 的个根，为特征方程的特征根，是其特征向量的分量。pija4、下面再证明主成分的方差是依次递减设相关系数矩阵的个特征根为，相应的特征向量为RppL21japppppppaaaaaaaaaaaaAMLMMMMLL21212222111211相对于的方差为1F111111)(aRaaXXaFVar同样有：，即主成分的方差依次递减。并且协方差为：iiFVar)(jijiRaaXaXaCov),(jpiaaaa)(1jiaaaajip,0)(1综上所述，根据证明有，主成分分析中的主成分协方差应该是对角矩阵，其对角线上的元素恰好是原始数据相关矩阵的特征值，而主成分系数矩阵的元素则是原始数据相关A矩阵特征值相应的特征向量。矩阵是一个正交矩阵。A于是，变量经过变换后得到新的综合变量pxxxL,21ppppppppppxaxaxaFxaxaxaFxaxaxaFLLLL22112222121212121111新的随机变量彼此不相关，且方差依次递减。二、主成分分析的计算步骤样本观测数据矩阵为：名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 7 页，共 10 页 - - - - - - - - - npnnppxxxxxxxxxXLMMMMLL212222111211第一步：对原始数据进行标准化处理。)var(*jjijijxxxx),2 , 1;, 2, 1(pjniLL其中niijjxnx1121)(11)var(jniijjxxnx), 2, 1(pjL第二步：计算样本相关系数矩阵。pppppprrrrrrrrrRLMLMMLL212222111211为方便，假定原始数据标准化后仍用表示，则经标准化处理后的数据的相关系数为: Xtjnttiijxxnr111),2 , 1,(pjiL第三步：用雅克比方法求相关系数矩阵的特征值（）和相应的特征向量RpL21,。piaaaaipiiiLL2, 1,21第四步：选择重要的主成分，并写出主成分表达式。主成分分析可以得到个主成分，但是，由于各个主成分的方差是递减的，包含的信p息量也是递减的，所以实际分析时，一般不是选取个主成分，而是根据各个主成分累计p贡献率的大小选取前个主成分，这里贡献率就是指某个主成分的方差占全部方差的比重，k实际也就是某个特征值占全部特征值合计的比重。即名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 8 页，共 10 页 - - - - - - - - - 贡献率 =piii1贡献率越大，说明该主成分所包含的原始变量的信息越强。主成分个数的选取，主k要根据主成分的累积贡献率来决定，即一般要求累计贡献率达到85% 以上，这样才能保证综合变量能包括原始变量的绝大多数信息。另外，在实际应用中，选择了重要的主成分后，还要注意主成分实际含义解释。主成分分析中一个很关键的问题是如何给主成分赋予新的意义，给出合理的解释。一般而言，这个解释是根据主成分表达式的系数结合定性分析来进行的。主成分是原来变量的线性组合，在这个线性组合中个变量的系数有大有小，有正有负，有的大小相当，因而不能简单地认为这个主成分是某个原变量的属性的作用，线性组合中各变量系数的绝对值大者表明该主成分主要综合了绝对值大的变量，有几个变量系数大小相当时，应认为这一主成分是这几个变量的总和，这几个变量综合在一起应赋予怎样的实际意义，这要结合具体实际问题和专业，给出恰当的解释，进而才能达到深刻分析的目的。第五步：计算主成分得分。根据标准化的原始数据，按照各个样品，分别代入主成分表达式，就可以得到各主成分下的各个样品的新数据，即为主成分得分。具体形式可如下。nknnkkFFFFFFFFFLMMMMLL212222111211第六步：依据主成分得分的数据，则可以进行进一步的统计分析。其中，常见的应用有主成份回归，变量子集合的选择，综合评价等。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 9 页，共 10 页 - - - - - - - - - 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 10 页，共 10 页 - - - - - - - - -