文档介绍:主成分分析与主成分回归
第1页,共48页,编辑于2022年,星期二
1. Introduction
Chemometrics
Necessary Knowledge
第2页,共48页,编辑于2022年,星期12页,共48页,编辑于2022年,星期二
Grade dik received by student i from professor k is
矩阵: 一组相同大小的矢量组合经典例子: 教授给学生打分
True score of student i
Relative loading(importance) given by professor k
j: factors (i,e., subjects)
chem., physics, math., etc.
第13页,共48页,编辑于2022年,星期二
Four students three professortwo subject: Chemistry and English
Professors
1 2 3
Students
1
2
3
4
Students
1
2
3
4
Professors
1 2 3
factors
1 2
1
2
Factors
3教授给4学生写留学推荐信
[S] is the matrix of true scores, called the score matrix
[L] is the matrix of importance, called the loading matrix
得分矩阵
载荷矩阵
第14页,共48页,编辑于2022年,星期二
矩阵的秩:对于A(m×n), 其秩是A中
最大线性无关的行数(或列数)。
秩=组分数?
秩为几?三种组分,吸收光谱各不相同(s1, s2 ,s3)
6组溶液,各组分浓度不同
吸光度矩阵A(20×6)
Rank =Number of Eigenvalue
秩=不为0的特征值的数目
矩阵: 一组不同浓度组合的混合溶液测得的光谱集合
矢量: 一条光谱
第15页,共48页,编辑于2022年,星期二
Eigenvalue 特征值
奇异值分解法:Y=USVt
S: 对角矩阵,收集了Y的特征值
U: 标准列正交矩阵(Scores Matrix)
Vt:标准行正交矩阵(Loadings Matrix)
用Matlab 很方便!一句话!
BACK
第16页,共48页,编辑于2022年,星期二
2. PCA 主成分分析 Principal Component Analysis
目的1
基本步骤2
应用实例3
第17页,共48页,编辑于2022年,星期二
主成分分析(PCA)的目的
BACK
现代仪器获得两维数据(矩阵)
矩阵处理
确定秩为多少
确定复杂分析体系中的物种数
PCA的目的-定性
有几种物种species
定性
第18页,共48页,编辑于2022年,星期二
PCA的步骤
BACK
矩阵分解
真实误差法
收集特征值
特征值比值法
Y=USVt
在S中
比较RSD与RE
Max
第19页,共48页,编辑于2022年,星期二
BACK
NIPALS分解
矩阵分解
奖金10000元
=
10000×1
5000×2
100×100
1×10000
Y=TP
奇异值(SVD)分解
Single Value Decomposition
Y=USVt
S: 对角矩阵,收集了Y的特征值
U: 标准列正交矩阵(Scores Matrix)
Vt:标准行正交矩阵(Loadings Matrix)
用Matlab 很方便!一句话!
怎么分解?
看了头大!
分解成正交矩阵的乘积
第20页,共48页,编辑于2022年,星期二
Y(m×n)有d个主成分
真实误差法-确定主成分数d
+表示来自主因子
0表示来误差
=
真实误差RE (Real Error,可以知道)
RE=RSD (剩余标准偏差)
Residual Standard Deviation
第21页,共48页,编辑于2022年,星期二
确定或设定RE
d=1…n-1计算RSD(d)
d=1
RSD(d)≤RE
YES
此时d即为主成分数
No
d