数学建模|多元分析（一）

发表于 2022-01-13 更新于 2022-05-10 Valine：本文字数： 7.8k 估计阅读时长 ≈ 7 分钟

多元分析是多变量的统计分析方法。

聚类分析

聚类分析一般分为Q型聚类分析和R型聚类分析。

Q型聚类分析是指对样品进行聚类分析
R型聚类分析是指对变量进行聚类。

根据处理方法的不同聚类分析又分为系统聚类法、有序样品聚类法、动态聚类法、模糊聚类法、图论聚类法。

聚类的一般过程：

数据预处理（标准化）
构造关系矩阵（亲疏关系的描述）
聚类（根据不同方法进行分类）
确定最佳分类（类别数）

Q型聚类分析

样本相似性度量

要用数量化的方法对事物进行分类，就必须用数量化的方法描述事物之间的相似程度。一个事物常常需要用多个变量来刻画。如果对于一群有待分类的样本点需用个变量描述，则每个样本点可以看成是空间中的一个点。因此，很自然地想到可以用距离来度量样本点间的相似程度。

记是样本点集，距离是的一个函数，满足条件：

$当且仅当$

这一距离的定义是我们所熟知的，它满足正定性，对称性和三角不等式。 在聚类分析中，对于定量变量，常用的闵氏(Minkowski)距离、绝对值距离、欧氏距离、切比雪夫距离。

闵氏距离

在 Minkowski 距离中，常用的是欧氏距离，它的主要优点是当坐标轴进行正交旋转时，欧氏距离是保持不变的。因此，如果对原坐标系进行平移和旋转变换，则变换后样本点间的距离和变换前完全相同。
值得注意的是在采用 Minkowski距离时，一定要采用相同量纲的变量。如果变量的量纲不同，测量值变异范围相差悬殊时，建议首先进行数据的标准化处理，然后再计算距离。在采用 Minkowski 距离时，还应尽可能地避免变量的多重相关性(multicollinearity)。多重相关性所造成的信息重叠，会片面强调某些变量的重要性。由于 Minkowski 距离的这些缺点，一种改进的距离就是马氏距离，定义如下:

马氏距离

其中为来自维总体的样本观测值；为的协方差矩阵，实际中往往是不知道的，常常需要用样本协方差来估计。马氏距离对一切线性变换是不变的，故不受量纲的影响。

类与类间的相似性度量

如果有两个样本类和，可以用下面的一系列方法度量它们间的距离：

(1)最短距离法

(2)最长距离法

(3)重心法

(4)类平均法

(5)离差平方和法

事实上，若内部点与点距离很小，则它们能很好地各自聚为一类，并且这两类又能够充分分离（即很大），这时必然有很大。因此，按定义可以认为，两类之间的距离很大。离差平方和法初是由 Ward 在 1936 年提出，后经 Orloci 等人 1976 年发展起来的，故又称为 Ward 方法。

聚类图

Q型聚类结果可由一个聚类图展示出来：

生成聚类图

案例

Q型聚类例子参考博客： matlab Q型聚类分析例题 Q型聚类分析

R型聚类

变量相似性度量

对变量进行聚类分析，首先要确定变量的相似性度量，常用的变量相似性度量有两张：

夹角余弦

各种定义的相似度量均应具有以下两个性质：

变量聚类法

类似样本集合聚类分析中最常用的最短距离法、最长距离法，变量聚类法采用了与系统聚类法相同的思路和过程。在变量聚类问题中，常用的有最长距离法、最短距离法。

最长距离法

在最长距离法中，定义两类变量的距离为：

$，$

式子中：或者，这时， $，$ 与两类中相似性最小的两变量间的相似性度量值有关。

最短距离法

在最短距离法中，定义两类变量的距离为：

$，$

式子中：或者，这时， $，$ 与两类中相似性最大的两变量间的相似性度量值有关。

案例

R型聚类分析

聚类算法分类

主成分分析(PCA)

主成分分析的主要目的是希望用较少的变量去解释原来资料中的大部分变异，将我们手中许多相关性很高的变量转化成彼此相互独立或不相关的变量。从数学的角度来说，是一种降维处理技术，降维可以去除噪声和不重要的特征，可以提升数据处理的速度。

基本思想

如果用表示门课程，表示各门课程的权重，那么加权之和是：

我们希望选择合适的权重能更好的区分学生的成绩。每个学生都对应一个综合成绩，记为为学生人数，如果这些值很分散，则表明区分的很好，但是现实情况往往是不能很好的分散，所以需要找到一个加权方式来使得尽可能的分散。 有如下的统计定义： 设表示以为样本观测值的随机变量，如果能找到一组权重3，使得

的值达到最大，则由于方差反应了数据差异的程度，也就表明了我们抓住了个变量的最大变异。当然，该式子也必须加上某种限制，不然权值无限大就没有意义了，一般来说规定

在此约束条件下，求得的最优解是纬空间的一个单位向量，它所表示的方向，就是主成分的方向。实际运用中，一个主成分往往不能代表原来的变量，所以我们需要寻找第二个、第三个……，且新寻找的主成分不再包含旧主成分的信息。设表示第个成分，,则可以表示为：

且对每个，均有

基本方法

假设有n个样本，p个指标，则可以构成大小为的矩阵

1. 进行标准化处理

按列计算均值和标准差

再标准化数据，原始样本矩阵经过标准化为：

2. 计算标准化样本的协方差矩阵

其中

其实这两步可以合成一步，

3. 计算的特征值和特征向量

特征值：

(是半定矩阵，且)

特征向量：

Matlab中计算特征值和特征向量的函数：

4. 计算主成分贡献率以及累计贡献率

贡献率：

累计贡献率：

5.写出主成分

一般取累计贡献率超过的特征值所对应的第一、第二、……、第m个主成分,第个主成分：

6. 根据系数分析主成分代表的意义

简而言之，系数越大，影响越大

主成分分析回归

主成分可用于聚类分析、回归分析。

例如Hald水泥问题，参考这里：水泥问题及matlab代码在文章尾部

案例

数学建模常用模型08 ：主成分分析

数学建模--主成分分析

数学建模算法一简述（4）主成分分析（PCA）

因子分析(FA)

因子分析可以看做是主成分分析的推广，也是利用降维的思想，由研究原始变量相关矩阵或协方差矩阵的内部依赖关系出发，把一些具有错综复杂关系的多个变量归结为少数几个综合因子的一种多元统计分析方法。

基本思想

把每个研究变量分解为几个影响因素变量，将每个原始变量分解成两部分因素，一部分是由所有变量共同具有的少数几个公共因子组成的，另一部分是每个变量独自具有的因素，即特殊因子。

与主成分分析的区别

主成分分析模型是原始变量的线性组合，是将原始变量加以综合、归纳，仅仅是变量变换；而因子分析是将原始变量加以分解，描述原始变量协方差矩阵结构的模型；只有当提取的公因子个数等于原始变量个数时，因子分析才对应变量变换。
主成分分析中每个主成分对应的系数是唯一确定的；因子分析中每个因子的相应系数即因子载荷不是唯一的。
因子分析中因子载荷的不唯一性有利于对公因子进行有效解释；而主成分分析对提取的主成分的解释能力有限。

因子分析模型

因子载荷矩阵中的几个统计性质

因子载荷的统计意义 因子载荷是第个变量与个公共因子的相关系数，反映了第个变量与第个公共因子的相关重要性，绝对值越大，相关的密切程度越高。
变量共同度的统计意义 变量的共同度是因子载荷矩阵的第行的元素的平方和，记为,对两边求方差，有

即

可以看出所有的公共因子对特殊因子对变量的贡献为，如果非常接近，说明从原变量空间到公共因子空间的转化效果好。
公共因子方差贡献的统计意义 因子载荷矩阵中各列元素的平方和
则称对所有的的方差贡献和，用于衡量的相对重要性。

因子载荷矩阵的估计方法

因子分析的一个基本问题是如何估计因子载荷，即如何求解因子模型式，下面介绍常用的因子载荷矩阵的估计方法。

主成分分析法

设为样本相关系数矩阵的特征值，为相应的标准正交化特征向量，设，则因子的荷载矩阵为

实例参考（含matlab代码）：因子分析—建立载荷矩阵

主因子法

主因子法是对主成分法的修正， 首先对变量标准化 $。$ 记式中，为约相关系数矩阵，对角线上的元素是。在实际情况中对特殊因子的方差一般是未知的，可以通过一组样本来估计，估计的方法如下：

取，在这种情况下主因子解与主成分解等价
取,这意味着与其余的的简单相关系数的绝对值最大者。记作

直接求的前个特征值和对应的正交特征向量。得到如下的因子载荷矩阵:

极大似然估计

Matlab工具箱求因子载荷矩阵使用的是最大似然估计法，命令是

案例

数学建模常用模型14 ：因子分析

因子旋转（正交变换）

为什么要旋转因子？

建立了因子分析数学目的不仅仅要找出公共因子以及对变量进行分组，更重要的要知道每个公共因子的意义，以便进行进一步的分析，如果每个公共因子的含义不清，则不便于进行实际背景的解释。由于因子载荷阵是不唯一的，所以应该对因子载荷阵进行旋转。目的是使每个变量在尽可能少的因子上有比较高的载荷，让某个变量在某个因子上的载荷趋于1，而在其他因子上的载荷趋于0。即：使载荷矩阵每列或行的元素平方值向0和1两极分化。 旋转的方法有：正交旋转；斜交旋转

正交旋转

由初始载荷矩阵A左乘一正交矩阵得到；目的是新的载荷系数尽可能的接近于0或尽可能的远离0；只是在旋转后的新的公因子仍保持独立性。

方差最大法： 方差最大法从简化因子载荷矩阵的每一列出发，使和每个因子有关的载荷的平方的方差最大。当只有少数几个变量在某个因子上有较高的载荷时，对因子的解释最简单。方差最大的直观意义是希望通过因子旋转后，使每个因子上的载荷尽量拉开距离，一部分的载荷趋于 1，另一部分趋于0。
四次方最大旋转： 四次方最大旋转是从简化载荷矩阵的行出发，通过旋转初始因子，使每个变量只在一个因子上有较高的载荷，而在其它的因子上尽可能低的载荷。如果每个变量只在一个因子上有非零的载荷，这时的因子解释是最简单的。四次方最大法通过使因子载荷矩阵中每一行的因子载荷平方的方差达到最大。
等量最大法： 等量最大法把四次方最大法和方差最大法结合起来求行和列因子载荷平方的方差的加权平均最大。