Page 103 - 《中国药房》2022年7期

60 的变化，从而影响其质量，且色度值对区分金樱子生品
A B A A
与盐炙品有一定影响。
50
表7 金樱子与盐金樱子的PCA载荷
40 指标 PC1 PC2 PC3 PC4 PC5
0.316
色度值 30 金樱子 1号峰 0.321 0.174 －0.175 0. 0.240 0.318
0.480
0.132
2号峰
0.
盐金樱子
3号峰 0. 0.331 0. 0.499 －0.260
20
4号峰 0.136 －0.113 －0.137 0.578 －0.366
5号峰 0.358 0. 0.422 0. 0.
10
6号峰 0.338 －0.107 0. 0. －0.380
7号峰 0.289 0. －0.233 －0.401 －0.396
0
L 值 a 值 b 值 E 值 8号峰 0. 0.361 －0.409 0. 0.333
*
*
*
*
9号峰 0.100 0.494 0. －0.170 0.
A：与生品比较，P＜0.05；B：与生品比较，P＞0.05
10号峰 0.246 0.147 0.497 0. 0.141
图5 金樱子与盐金樱子色度值差异分析柱状图
色度值L * 0.320 －0.234 －0.287 0. 0.340
2.2.8 多元统计分析 采用 CA、PCA 、OPLS-DA 对各 色度值a * －0.254 0.345 0.114 －0.256 －0.340
色度值b * 0.410 0.126 －0.161 0. 0.
批金樱子和盐金樱子的测定结果进行分析。 总色度值E * 0.355 0.115 －0.274 －0.286 0.
（1）CA：以金樱子和盐金樱子共有峰的相对峰面积 标准差 2.167 1.816 1.327 1.124 1.078
为变量，导入SPSS 26.0统计软件进行系统CA。将数据 方差贡献率/％ 33.54 23.57 12.57 9.02 8.30
累计方差贡献率/％ 33.54 57.11 69.68 78.70 87.00
标准化后采用瓦尔德法，以欧氏距离作为距离公式计
算，结果见图6。当欧氏距离为20时，金樱子样品分为2
类，即生品与盐炙品各聚为一类。 4
欧氏距离
23.57％ 0 金樱子
0 5 10 15 20 25
PC2： 盐金樱子
－4

－4 －2 0 2
PC1：33.54％
图7 金樱子与盐金樱子PCA得分图

（3）OPLS-DA：利用 R 语言中的 ropls 数据包对金樱
子、盐金樱子的共有峰峰面积和色度值进行有监督模式
的 OPLS-DA，结果见图 8A。模型中累积解释能力参数
图6 金樱子与盐金樱子CA树状图
（R Y）为 0.931，表明拟合度好，预测能力参数（Q Y）为
2
2
（2）PCA：利用R语言软件中princomp函数，以金樱
[12]
0.864，说明预测能力较好 ，可用于区别金樱子与盐金
子、盐金樱子共有峰峰面积与色度值为变量进行无监督
樱子。同时以各样本在投影平面内以及正交投影面的
模式的 PCA，分析结果见表 7。利用 R 语言中的 ggplot2
距离作图，结果见图8B。由图8B可知，没有样品超过观
数据包制作二维散点图，结果见图7。根据分析结果，选
测距离，说明样品间无较大差异的个例。对建立的模型
择标准差大于1的前5个主成分（PC1～PC5），可反映金
进行随机置换检验，结果见图 8C，可知实际模型的 Q Y
2
樱子与盐金樱子的指纹图谱和色度值 87.00％的信息。
值与验证模型的Q Y值（图上散点）均小于真实值（图上
2
*
其中 PC1 的方差贡献率为 33.54％，色度值 b（0.410）、5
横线），表明所构建的OPLS-DA模型未出现过拟合。提取
号峰（0.358）与总色度值 E（0.355）的载荷较大；PC2 的
*
OPLS-DA 模型中变量重要性投影（variable importance
方差贡献率为23.57％，2号峰（0.480）与9号峰（0.494）的
载荷较大；PC3 的方差贡献率为 12.57％，5 号峰（0.422） in the projection，VIP），结果见表8。筛选VIP值＞1的成
和 10 号峰（0.497）有较大载荷；PC4 的方差贡献率为 分作为区分金樱子生品、盐炙品的重要因素，分别为5号
9.02％，3 号峰（0.499）和 4 号峰（0.578）有较大载荷；PC5 峰（VIP 值为 1.537）、1 号峰（VIP 值为 1.445）、色度值 b *
*
的方差贡献率为8.30％，6号峰（－0.380）与7号峰（－0.396） （VIP 值为 1.376）、总色度值 E（VIP 值为 1.251）、10 号峰
有较大载荷。由图 7 可知，金樱子与盐金樱子沿 PC1 轴 （VIP 值为 1.195）、6 号峰（VIP 值为 1.161），这些因素可
各自聚为一类，说明金樱子盐炙后发生了多种化学成分 显著影响金樱子生品和盐炙品的分类。
中国药房 2022年第33卷第7期 China Pharmacy 2022 Vol. 33 No. 7 ·865 ·