Page 56 - 《中国药房》2021年21期

2.6 麸炒北苍术饮片炮制过程样品的多元统计分析
本研究以麸炒北苍术饮片炮制过程样品的色度值 200
以及质量指标等33个变量为指标，采用Simca 14.1软件
150
进 行 多 元 统 计 分 析 [ 主 成 分 分 析（PCA）、聚 类 分 析
（HCA）、偏最小二乘法判别分析（PLS-DA）] 。 平方欧氏距离 100
[13]
2.6.1 PCA 和 HCA 对麸炒北苍术饮片炮制过程样品
50
进行PCA时，共提取了4个特征值＞1的主成分因子，其
特征值和方差贡献率见表7，PCA得分图见图5。由表7 0 S0 S1 S2 S3 S4 S5 S6 S7 S8 S9 S10 S11 S12 S13 S14 S15
和图5可知，这4个主成分的作用显著，可用于反映麸炒 图6 麸炒北苍术饮片炮制过程样品的HCA图
北苍术饮片炮制过程样品90％以上的特征信息。 Fig 6 HCA diagram in sample of A. chinensis decoc-
表7 主成分因子的特征值和方差贡献率 tion piece stir-fired with bran during processing
Tab 7 Characteristic value and variance contribution
rate of principal component factors
4
主成分因子 特征值 方差贡献率，％ 累积方差贡献率，％
2
1 23.10 70.00 70.00
2 3.99 12.08 82.09 0
3 1.79 5.43 87.51 t[2]
－2
4 1.14 3.44 90.96
－4
－6
4 －15 －10 －5 0 5 10
t[1]
2
图7 麸炒北苍术饮片炮制过程样品的PLS-DA得分图
0
t[2] Fig 7 PLS-DA score diagram in sample of A. chinen-
－2
sis decoction piece stir-fired with bran during
－4
processing
－6
－15 －10 －5 0 5 10
t[1]
2.5
图5 麸炒北苍术饮片炮制过程样品的PCA得分图 1.5 2
Fig 5 PCA score diagram in sample of A. chinensis de- VIP值 0.5 1
coction piece stir-fired with bran during pro- 0
－0.5
cessing a * L * b *
以平方欧氏距离为测度，进行 HCA，结果见图 6。 峰2相对峰面积 峰19相对峰面积 峰9相对峰面积 峰7相对峰面积 峰8相对峰面积 峰10相对峰面积 峰13相对峰面积 峰1相对峰面积 峰4相对峰面积 峰21相对峰面积 醇溶性浸出物含量 峰20相对峰面积 峰14相对峰面积 峰25相对峰面积 峰3相对峰面积 峰22相对峰面积 峰12相对峰面积 峰16相对峰面积 峰17相对峰面积 挥发油含量 水溶性浸出物含量 峰11相对峰面积 峰23相对峰面积 峰24相对峰面积 峰15相对峰面积 峰6相对峰面积 峰18相对峰面积 峰26相对峰面积 峰27相对峰面积 峰5相对峰面积
由图6可知，麸炒北苍术饮片炮制过程样品可聚为3类：
S0～S5 样品聚为一类，S6～S11 样品聚为一类，S12～ 图8 麸炒北苍术饮片炮制过程样品的VIP值图
S15 样品聚为一类。基于此，笔者将与生品（S0）接近的 Fig 8 VIP value chart in sample of A. chinensis decoc-
样品拟定为炮制轻度样品（S0～S5），与 S15 样品（该样 tion piece stir-fired with bran during processing
品中峰15丢失，其原因可能为炮制过度）接近的样品拟
2.7 麸炒北苍术饮片的炮制时间研究
定为炮制过度样品（S12～S15），其余为炮制适度样品
结合“2.6”项下结果可知，麸炒北苍术饮片炮制适度
（S6～S11，炮制时间为18～33 min）。
的时间为 18～33 min，但该时间范围较宽，无法较好地
2.6.2 PLS-DA 对麸炒北苍术饮片炮制过程样品进行
确定麸炒北苍术饮片的炮制时间。基于此，本研究参考
PLS-DA 后，其模型得分图见图 7。由图 7 可知，麸炒北
[15]
*
灰色关联度法 ，将“2.6”项下的重要质量指标a 和峰2、
苍术饮片炮制过程样品可聚为3类，聚类结果与HCA结
果一致。PLS-DA 模型中 R Y＝0.86，Q ＝0.69，表明该 19、9、7、8、10、13、1、4、21 的 相 对 峰 面 积 以 及 峰 3
2
2
[16]
模型拟合程度好、预测能力强 。进一步对该模型中33 （5-HMF，其具有一定毒性 ）的相对峰面积作为评价指
[14]
个指标的变量投影重要性（VIP）进行分析，结果见图 8。 标，进行麸炒北苍术饮片炮制过程样品的综合质量评
选择VIP值大于1的指标作为麸炒北苍术饮片炮制过程 价，以确定其炮制时间，具体步骤如下。
中的重要质量指标 ，得到 a 和峰 2、19、9、7、8、10、13、 第一，确定分析序列。将麸炒北苍术饮片炮制过程
[15]
*
1、4、21 的相对峰面积是该饮片炮制过程中的重要质量 的 16 批样品组成 16 个评价单元序列 xi （k）（i＝1，2，
指标。 3……m；k＝1，2，3……n。本研究中 m＝16，n＝12），并
·2610 · China Pharmacy 2021 Vol. 32 No. 21 中国药房 2021年第32卷第21期