医学统计学课件-主成分分析与因子分析(第20章)

主成分分析与因子分析PrincipalComponentsAnalysis&FactorAnalysis2020/12/191医学统计学讲课内容：第一节主成分分析第二节因子分析2020/12/192医学统计学第一节主成分分析PrincipalComponentsAnalysis2020/12/193医学统计学一、基本思想数据的降维、数据的解释将原来众多具有一定相关性的指标，组合成一组新的相互无关的综合指标。

从中选取几个较少的综合指标尽可能多的反映原来众多指标的信息。这种既减少了指标的数目又抓住了主要矛盾的做法有利于问题的分析和处理。2020/12/194医学统计学2020/12/195医学统计学

如何利用这些指标对每一儿童的生长发育作出正确评价？仅用单一指标：结论片面；没有充分利用原有数据信息。利用所有指标：各指标评价的结论可能不一致，使综合评价困难；工作量大。2020/12/196医学统计学找出几个综合指标(长度、围度、特体)，这些综合指标是原始指标的线性组合，既保留了原始指标的信息，且互不相关。各综合指标提供的“信息”量大小用其方差来衡量。衡量一个指标的好坏除了正确性与精确性外，还必须能充分反映个体间的变异，一项指标在个体间的变异越大，提供的信息量越多。2020/12/197医学统计学二、数学模型及几何意义2020/12/198医学统计学Z=AX2020/12/199医学统计学第一主成分在所有Zi中最大2020/12/1910医学统计学第二主成分……理论上主成分个数最多为m个(指标个数)实际工作中确定的主成分个数总是小于m个在所有Zi中为第2大。无关，互相垂直：2020/12/1911医学统计学X1X2112-2-2-1-120相关变异2020/12/1912医学统计学X1X2Z1Z2112-2-2-2-211-1-1-1-122202020/12/1913医学统计学Z1Z2-2-211-1-1220相关变异2020/12/1914医学统计学三、主成分的求法及性质2020/12/1915医学统计学（一）主成分的求法

1.对各原始指标值进行标准化为了方便，仍用Xij表示Xij’。2020/12/1916医学统计学标准化后的数据矩阵X=2020/12/1917医学统计学2.

求出X1,X2,…,Xm

的相关矩阵RR=Cov(X)

=2020/12/1918医学统计学Pearson相关系数标准化后的协方差协方差2020/12/1919医学统计学3.

求出矩阵R的全部特征值(eigenvalue)i,

第i个主成分的组合系数ai1,ai2,,aim满足方程组:(r11－i)ai1+r12ai2++r1maim=0r21ai1+(r22－i)ai2++r2maim=0

rm1ai1+rm2ai2++(rmm－i)aim=0

2020/12/1920医学统计学(r11－i)ai1+r12ai2++r1maim=0r21ai1+(r22－i)ai2++r2maim=0

rm1ai1+rm2ai2++(rmm－i)aim=0i为矩阵R的第i个特征值，共有m个非负特征值，由大到小的顺序排列为：1≥2≥≥m≥0i=Var(Zi)2020/12/1921医学统计学4.由以上方程组，求出相应于特征值i的特征向量(eigenvector)(ai1,ai2,,aim)’2020/12/1922医学统计学（二）主成分的性质

1.各主成分互不相关

2020/12/1923医学统计学2.主成分的贡献率与累积贡献率

(原始指标值标准化)(指标个数)贡献率累积贡献率2020/12/1924医学统计学3.主成分个数的选取（1）前k个主成分的累积贡献率>70%。（2）主成分Zi的特征值i≥

1。4.因子载荷（第i主成分Zi与第j原始指标Xi间相关系数）2020/12/1925医学统计学5.样品的主成分得分2020/12/1926医学统计学四、实例2020/12/1927医学统计学2020/12/1928医学统计学2020/12/1929医学统计学2020/12/1930医学统计学1.主成分个数的选取

3很接近于1；

3

与2的贡献率相差不大，为25%左右，若舍去3不合理。取前三个主成分。2020/12/1931医学统计学2.列出主成分表达式Z1为急性炎症成分（X1转氨酶、X2肝大指数）Z2为慢性炎症成分（X3硫酸锌浊度）Z3为癌变成分（X4甲胎球蛋白）2020/12/1932医学统计学3.求出因子载荷阵2020/12/1933医学统计学4.主成分得分2020/12/1934医学统计学标准化指标主成分还原为原始指标主成分2020/12/1935医学统计学

将该肝病患者的四项肝功能指标代入原始指标主成分表达式：

Z1=2.50865Z2=－1.06626Z3=－1.22943

该肝病患者可能为急性炎症。2020/12/1936医学统计学五、主成分分析的应用

1.对原始指标进行综合以互不相关的较少个综合指标反应众多原始指标提供的信息。

主成分回归(解决多元共线问题)。

2.进行综合评价

2020/12/1937医学统计学3.进行探索性分析利用因子载荷阵，找出影响各综合指标的主要原始指标。

4.对样品进行分类利用主成分得分对样品进行分类：

Z1为急性炎症成分

Z2为慢性炎症成分

Z3为癌变成分2020/12/1938医学统计学第二节因子分析FactorAnalysis2020/12/1939医学统计学一、因子分析基本思想从分析多个可观测的原始指标的相关关系入手，找到支配这种相关关系的有限个不可观测的潜在变量。是多元分析中处理降维的一种统计方法。如：脑部疾病患者的意识清醒状态可由语言能力、辩识能力、记忆能力、理解能力与思维逻辑能力等可观测的指标反映。2020/12/1940医学统计学二、因子分析数学模型

X1:收缩压X2:舒张压X3:心跳间隔X4:呼吸间隔X5:舌下温度F1:交感神经F2:副交感神经commonfactor2020/12/1941医学统计学specificfactorcommonfactor2020/12/1942医学统计学Xi:观测指标(标准化数据)Fi:公因子ei:特殊因子aij:因子载荷(计算关键项)2020/12/1943医学统计学X=AF+e2020/12/1944医学统计学2020/12/1945医学统计学2020/12/1946医学统计学三、因子模型的性质

矩阵A的统计意义1.公共度(共性方差)2020/12/1947医学统计学因子的共性方差2020/12/1948医学统计学2.因子贡献与因子贡献率矩阵A第j列元素反映了第j个公因子Fj对所有原始指标的影响;数据标准化后全部原始指标的总方差为指标个数m。Fj对原始指标的方差贡献率2020/12/1949医学统计学各因子的贡献2020/12/1950医学统计学3.因子载荷及因子载荷阵A2020/12/1951医学统计学四、因子载荷阵的求解及计算步骤

1.

收集原始数据并整理为下表2020/12/1952医学统计学2.对各指标进行标准化3.求指标间的相关系数矩阵RX4.求指标间的约相关系数矩阵R*

（1）R*的非对角线元素与相关矩阵RX的非对角线元素相等

（2）R*的对角线元素为共性方差2020/12/1953医学统计学5.求出约关系数矩阵R*所有大于零的特征值及相应的特征向量6.写出因子载荷阵A，得出原始指标X的公因子表达式2020/12/1954医学统计学要求：1.

保留公因子个数q小于指标个数m，原则：

j≥1

前k个公因子累积贡献率≥70%2.

各共性方差接近于1。3.

各原始指标在同一公因子Fj上的因子载荷之间的差别应尽可能大。2020/12/1955医学统计学五、实例2020/12/1956医学统计学1.主成分解2020/12/1957医学统计学2020/12/1958医学统计学2020/12/1959医学统计学主成分解：除因子1可初步认定为综合因子外，其余3个因子的专业意义不明显。2.主因子解：除因子1可初步认定为综合因子外，其余3个因子的专业意义不明显。2020/12/1960医学统计学六、因子旋转当各公因子的专业意义难以解释时，可以通过因子旋转来解决。如求得的因子载荷阵A不甚理想，可右乘一个正交阵T，使AT有更好的实际意义，

使各原始指标在同一公因子上之间差别尽可能增大。称因子正交旋转。正交旋转可保持各指标的共性方差不变；各公因子互不相关。常用方差最大旋转法等。2020/12/1961医学统计学2020/12/1962医学统计学2020/12/1963医学统计学七、几点注意

1.因子分析的解不唯一（1）同一问题可以有不同的因子分析解：主成分解、主因子解、极大似然解（2）进行因子旋转以获得更为满意的解。2.因子得分不能直接进行计算，但可以估计。2020/12/1964医学统计学3.主成分分析与因子分析间的关系（1）两者的分析重点不一致

Z=AX主成分为原始变量线性组合，重点在综合原始变量信息。

X=AF+e原始变量为公因子与特殊因子线性组合，公因子重点反映支配原始变量的不可观测的潜在因素。重