附子的多元统计分析方法与流程

clustering)。
13.进一步的，所述高效液相色谱串联质谱联用检测中附子制品及煻灰炮附子需制备相应的附子供试品溶液；
14.所述附子供试品溶液是取附子制品/煻灰炮附子加入氨水溶液，再加入体积比为1∶1的异丙醇-乙酸乙酯混合溶液，经超声提取、定容、过滤，所得续滤液经浓缩后，再次加入体积比为1∶1的异丙醇-二氯甲烷混合溶液溶解，过滤制得。
15.进一步的，所述高效液相色谱串联质谱联用检测中附子制品及煻灰炮附子需制备相应的附子供试品溶液；
16.所述附子供试品溶液是取4重量份附子制品/煻灰炮附子加入3体积份浓度为10wt％的氨水溶液，再加入50体积份体积比为1∶1的异丙醇-乙酸乙酯混合溶液，经超声提取、定容、过滤，再取40体积份所得续滤液浓缩后，再次加入1体积份体积比为1∶1的异丙醇-二氯甲烷混合溶液溶解，过滤制得；
17.其中，重量份与体积份之间的对应关系为g：ml。
18.进一步的，所述高效液相色谱串联质谱联用检测中的质谱分析采用esi离子源，在正负离子模式下进行检测。
19.进一步的，所述高效液相色谱串联质谱联用检测中的质谱条件为：
20.第一个质量分析器(ms1)的检测范围为100～1000质荷比(m/z)，离子积累时间设定为30ms；
21.第二个质量分析器(ms2)的检测范围为50～1000质荷比(m/z)，离子积累时间设定为20ms；
22.第三个质量分析器(ms3)的检测范围为50～1000质荷比(m/z)，离子积累时间设定为20ms。
23.进一步的，所述高效液相色谱串联质谱联用检测中的质谱条件为：
24.氮气流速为1.5l/min，cdl温度为200℃，加热块的温度为200℃，pda的检测波长范围为200～700nm，正离子模式下接口电压为4.5kv，负离子模式下接口电压为3.5kv，正负离子模式的检测电压均为1.7kv，干燥气压设定为100kpa，最大压力设定为30mpa。
25.进一步的，所述高效液相色谱串联质谱联用检测中色谱条件为：
26.流动相a为采用浓度为0.04mol/l的乙酸铵溶液，流动相b为乙腈；
27.采用梯度洗脱，所述梯度洗脱的洗脱条件为：
28.0～10min，95％
→
90％流动相a，5％
→
10％流动相b；
29.10～25min，90％
→
60％流动相a，10％
→
40％流动相b；
30.25～50min，60％
→
48％流动相a，40％
→
52％流动相b；
31.50～80min，48％
→
30％流动相a，52％
→
70％流动相b；
32.80～95min，30％
→
15％流动相a，70％
→
85％流动相b。
33.进一步的，所述乙酸铵溶液为氨水调节ph值至10的乙酸铵溶液。
34.进一步的，所述高效液相色谱串联质谱联用检测中色谱柱为phenomenex gemini c18色谱柱、柱温为40℃、流速为1.0ml/min、检测波长为240nm。
35.本发明的附子的多元统计分析方法的有益效果为：
36.本发明通过对附子制品及煻灰炮附子，采用hplc-it-tof-msn结合生物碱含量测
定及多元统计分析方法(主成分分析、偏最小二乘-判别分析、聚类分析)对不同附子制品的化学成分进行了研究和比较，为煻灰炮附子在经典名方中的应用提供理论依据和数据支持；
37.本发明通过多元统计分析方法对附子制品进行快速检测评判，准确将不同品质的附子制品区分处理，能够快速寻找出能够代替古法煻灰炮附子的附子制品；
38.本发明的多元统计分析方法具有快速、简单、灵敏、高效等优点，且不受附子制品种类的影响，可以对不同种类的附子制品进行检测分类；
39.本发明的多元统计分析的方法可用于附子制品的质量控制及质量标准的提高，为附子制品不同品质批次间的质量评价研究提供借鉴；
40.本发明通过调整色谱条件、质谱条件及附子的前处理过程，再采用高效液相色谱串联质谱联用检测方法对附子进行检测，能够有效将活性成分分离并检出；
41.本发明采用高效液相色谱串联质谱联用检测方法对煻灰炮附子及附子制品(生附片、炒附片、蒸附片和烘附子等)中化学成分进行鉴定，结合对照品及质谱特征共指认38种成分，其中煻灰炮附子18种，生附片24种(16种与煻灰炮附子相同)，炒附片22种(15种相同)、蒸附片17种(10种相同)，烘附子20种(15种相同)；
42.本发明采用高效液相色谱法对乌头碱类生物碱的含量进行测定，表明煻灰炮附子(0.0151％～0.0654％)和炒附片(0.0263％～0.0582％)中次乌头碱含量最高，生附片中新乌头碱含量最高(0.0638％～0.1207％)，而蒸附片中苯甲酰次乌头原碱含量最高(0.0106％～0.0213％)；
43.本发明通过对不同附子制品中6种乌头碱类生物碱含量及14个主要化学成分色谱峰面积进行了聚类分析，发现炒附片与煻灰炮附子归为一类的分类距离最小(d＝2.5，d＝6)；对不同附子制品中14个主要化学成分色谱峰面积进行了偏最小二乘法-判别分析，发现新乌头碱、次乌头碱、苯甲酰新乌头原碱的含量对区分不同附子制品及煻灰炮附子的贡献最大(vip≥1)；
44.本发明通过对不同附子制品及煻灰炮附子间的化学成分差异进行多元分析对比，发现炒附片与煻灰炮附子具有15种相同成分，且聚类分析中距离最小，同时两者中次乌头碱的含量均最高，说明炒附片具有与煻灰炮附子最为相似的化学成分物质基础，可以以炒附子作为煻灰炮附子的替代品；
45.本发明通过对5种附子制品(煻灰炮附子、生附子、炒附子、蒸附子、烘附子)进行高效液相色谱串联质谱联用检测，从特征碎片离子信息中鉴定出38种化学成分(均为生物碱)，并指认了33种化学成分；其中，煻灰炮附子中共鉴定出18种化学成分；生附片中共鉴定出24个化学成分，其鉴定出的14-乙酰尼奥灵或其异构体、乌头卡查敏碱b、翠雀碱或异翠雀碱未在其它加工品中检出；炒附片中共鉴定22个化学成分，其鉴定出的10-羟基苯甲酰新乌头原碱、苯甲酰-3，13-双去氧新乌头原碱或其异构体、脱水苯甲酰新乌头原碱、脱水苯甲酰次乌头原碱、唐古替星a或其异构体未在其它炮制方法制备的饮片中检出；蒸附片中共鉴定17个化学成分，其鉴定出的1-乙酰光泽乌头碱或12-乙酰光泽乌头碱、苯甲酰去氧乌头原碱未在其它炮制方法制备的饮片中检出；烘附子中共鉴定20个化学成分，其鉴定出的尼奥灵或10-羟基塔拉乌头胺、森布星a或森布星b、14-乙酰尼奥灵或其异构体、卡米车灵a或其异构体未在其它炮制方法制备的饮片中检出；
46.本发明的高效液相色谱串联质谱联用检测方法具有超高的灵敏度和优异的定性能力；
47.本发明中38种化学成分存在细微结构差异，产生的质谱特征片段存在较大差异，可准确确证活性成分；
48.本发明的多元统计分析方法准确度高、方法简单、分析周期短、稳定性好、重复性良好、活性成分分离效果好，可作为附子制品药理活性和安全性检查的重要方法。
附图说明
49.图1是本发明实施例1中不同附子制品中6种乌头碱类生物碱含量对比图；
50.图2是本发明实施例1中不同附子制品的高效液相色谱对比图；
51.图3是本发明实施例1中蒸附子的高效液相色谱图；
52.图4是本发明实施例1中炒附子的高效液相色谱图；
53.图5是本发明实施例1中生附子的高效液相色谱图；
54.图6是本发明实施例1中煻灰炮附子的高效液相色谱图；
55.图7是本发明实施例1中煻灰炮附子的高效液相色谱图；
56.图8是本发明实施例1中不同附子制品及煻灰炮附子制品的分析色谱峰pca二维空间图；
57.图9是本发明实施例1中不同附子制品及煻灰炮附子制品的分析色谱峰pls-da二维空间图；
58.图10是本发明实施例1中不同附子制品及煻灰炮附子制品的vip值排序图；
59.图11是本发明实施例1中不同附子制品及煻灰炮附子制品的14种分析色谱峰峰面积聚类分析结果图；
60.图12是本发明实施例1中不同附子制品及煻灰炮附子制品的6种乌头碱类生物碱含量聚类分析结果图。
具体实施方式
61.下面对本发明实施例中的技术方案进行清楚、完整地描述。在下面的描述中阐述了很多具体细节以便于充分理解本发明，但是本发明还可以采用其他不同于在此描述的其它方式来实施，本领域技术人员可以在不违背本发明内涵的情况下做类似推广，因此本发明不受下面公开的具体实施例的限制。
62.实施例1一种附子的多元统计分析方法
63.本实施例分别对煻灰炮附子、生附子、炒附子、蒸附子、烘附子这五种不同的饮片进行多元统计分析，具体方法如下：
64.一、6种乌头碱类生物碱含量检测
65.分别精密称定煻灰炮附子、生附子、炒附子、蒸附子的粉末(过3号筛)各2g，置具塞锥形瓶中，加入3ml氨试液(取25％浓氨水400ml，加超纯水使成1000ml，即得)，精密加入50ml体积比为1∶1的异丙醇-乙酸乙酯混合溶液，称定重量，经超声处理(功率300w，频率40khz，水温在25℃以下)30分钟，放至室温，再次称定重量，并用体积比为1∶1的异丙醇-乙酸乙酯混合溶液补足减失的重量，摇匀，滤过，所得续滤液再次精密量取25ml，于40℃以下
减压回收溶剂至干，残余物精密加入3ml体积比为1∶1的异丙醇-二氯甲烷混合溶液溶解，滤过，取续滤液，相应的附子待测液，即煻灰炮附子待测液(制备8次，标记为pfz1)、生附子待测液(制备3次，标记为sfz1)、炒附子待测液(制备3次，标记为cfz1)和蒸附子待测液(制备3次，标记为zfz1)。其中，不同附子待测液中附子浓度均为(2/3)g/ml。
66.取苯甲酰乌头原碱、苯甲酰次乌头原碱、苯甲酰新乌头原碱对照品适量，分别精密称定后置于同一容量瓶中，加入体积比为1∶1的异丙醇-二氯甲烷混合溶液制成每1ml各含10μg苯甲酰乌头原碱、苯甲酰次乌头原碱、苯甲酰新乌头原碱的混合溶液，即得对照品溶液a。其中，苯甲酰乌头原碱、苯甲酰次乌头原碱、苯甲酰新乌头原碱的浓度均为10μg/ml。
67.取乌头碱、次乌头碱、新乌头碱对照品适量，分别精密称定后置于同一容量瓶中，加入体积比为1∶1的异丙醇-二氯甲烷混合溶液制成每1ml各含5μg乌头碱、次乌头碱、新乌头碱的混合溶液，即得对照品溶液b。其中，乌头碱、次乌头碱、新乌头碱的浓度均为5μg/ml。
68.分别煻灰炮附子待测液、生附子待测液、炒附子待测液、蒸附子待测液、对照品溶液a和对照品溶液b进行高效液相色谱检测，即得相应的附子和相应对照品的保留时间及峰面积，然后依据对应的保留时间和峰面积分别计算相应附子中苯甲酰新乌头原碱、苯甲酰乌头原碱、苯甲酰次乌头原碱、新乌头碱、次乌头碱和乌头碱这6种乌头碱类生物碱的含量，具体结果见表1和图1。
69.其中，高效液相色谱检测的检测条件为：
70.色谱柱为bdshypersil c18色谱柱(4.6mm
×
250mm，5μm)；
71.检测波长为235nm；
72.柱温为25℃；
73.流速为1.0ml/min；
74.进样量为10μl；
75.流动相a是0.1mol/l乙酸铵(每1000ml加冰醋酸0.5ml)，流动相b为体积比为25：15的乙腈-四氢呋喃混合溶液；
76.洗脱方式为梯度洗脱，具体洗脱程序为：
77.0～48min，85％
→
74％流动相a，15％
→
26％流动相b；
78.48～49min，74％
→
65％流动相a，26％
→
35％流动相b；
79.49～58min，65％流动相a，35％流动相b；
80.58～65min，65％
→
85％流动相a，35％
→
15％流动相b。
81.含量计算方法具体如下：
82.乌头碱类生物碱含量计算公式(％)＝(附子待测液对应的乌头碱类生物碱峰面积
×
对照品溶液中乌头碱类生物碱的浓度)/(乌头碱类生物碱对照品对应的峰面积
×
附子待测液中附子浓度)
×
100％；
83.比如：煻灰炮附子中苯甲酰新乌头原碱含量(％)＝(煻灰炮附子待测液对应的苯甲酰新乌头原碱峰面积
×
对照品溶液a中苯甲酰新乌头原碱的浓度)/(对照品溶液a中苯甲酰新乌头原碱对应的峰面积
×
煻灰炮附子待测液中附子浓度)
×
100％＝(煻灰炮附子待测液对应的苯甲酰新乌头原碱峰面积
×
1.5
×
10-5
/对照品溶液a中苯甲酰新乌头原碱对应的峰面积)
×
100％。
84.表1附子不同加工品的含量测定结果(％)
[0085][0086][0087]
bma＝苯甲酰新乌头原碱；bac＝苯甲酰乌头原碱；bha＝苯甲酰次乌头原碱；ma＝新乌头碱；ha＝次乌头碱；ac＝乌头碱
[0088]1表示形式为平均值
±
标准差；2n.d.表示实验条件下未检出
[0089]
由于，研究中烘附子样品不足，作为定量研究(含量测定)和统计分析研究(pca、pls-da、ca)的代表性不足，本实施例未对烘附子进行多元分析，但在定性研究中(hplc-it-tof-msn)可以作为不同附子加工品中化学成分差异的参考。
[0090]
二、高效液相色谱串联质谱联用检测(hplc-it-tof-msn)
[0091]
精密称取8批煻灰炮附子粉末(过3号筛)各0.5g，混合均匀，混合粉末共4g，置于具塞锥形瓶中，加3ml氨试液(取25％浓氨水400ml，加超纯水使之成1000ml)，精密加入50ml体积比为1∶1的异丙醇-乙酸乙酯混合溶液，称定重量，再经超声处理(功率300w，频率40khz，水温在25℃以下)30分钟，放至室温(20～25℃)，再次称定重量，用体积比为1∶1的异丙醇-乙酸乙酯混合溶液补足减失的重量，摇匀，滤过，所得续滤液再次精密量取40ml，于40℃以下减压回收溶剂至干，残渣精密加入体积比为1∶1的异丙醇-二氯甲烷混合溶液1ml溶解，滤过，所得滤液即为煻灰炮附子供试品溶液，标记为pfz。
[0092]
取3批生附片粉末(过3号筛)各1.4g，混合均匀，再次取所得混合粉末4g，依照煻灰炮附子供试品溶液的方法制备生附片供试品溶液，标记为sfz。
[0093]
依照生附子供试品溶液的制备方法，分别制备炒附片供试品溶液(标记为cfz)和蒸附片供试品溶液(标记为zfz)。
[0094]
取4g烘附子粉末(过3号筛)，依照煻灰炮附子供试品溶液的方法制备烘附子供试品溶液，标记为hfz。
[0095]
煻灰炮附子供试品溶液pfz、生附片供试品溶液sfz、炒附片供试品溶液cfz、蒸附片供试品溶液zfz、烘附子供试品溶液hfz统称为不同的附子供试品溶液。
[0096]
分别取烘附子供试品溶液hfz、煻灰炮附子供试品溶液pfz、生附片供试品溶液sfz、炒附片供试品溶液cfz、蒸附片供试品溶液zfz、烘附子供试品溶液hfz及上述对照品溶液a和对照品溶液b进行高效液相色谱串联质谱联用检测，即得相应的色谱图和质谱鉴定结果。所得相应的附子色谱图见图2～6，图2为不同附子制品的色谱对比图，其中a测定的是蒸附子供试品溶液，b测定的是炒附子供试品溶液，c测定的是生附子供试品溶液，d测定的是煻灰炮附子供试品溶液；图3为蒸附子的色谱图，图4为炒附子的色谱图，图5为生附子的色谱图，图6为煻灰炮附子的色谱图；图7为烘附子的色谱图。
[0097]
其中，高效液相色谱检测的检测条件为：
[0098]
色谱柱为phenomenex gemini c18色谱柱(4.6mm
×
250mm，5μm)；
[0099]
检测波长为240nm；
[0100]
柱温为40℃；
[0101]
流速为1.0ml/min；
[0102]
进样量为5μl；
[0103]
流动相a是以氨水调节ph值至10、浓度为0.04mol/l的乙酸铵溶液，流动相b乙腈；
[0104]
洗脱方式为梯度洗脱，具体洗脱程序为：
[0105]
0～10min，95％
→
90％流动相a，5％
→
10％流动相b；
[0106]
10～25min，90％
→
60％流动相a，10％
→
40％流动相b；
[0107]
25～50min，60％
→
48％流动相a，40％
→
52％流动相b；
[0108]
50～80min，48％
→
30％流动相a，52％
→
70％流动相b；
[0109]
80～95min，30％
→
15％流动相a，70％
→
85％流动相b。
[0110]
质谱检测条件为：
[0111]
采用esi离子源进行质谱分析，并在正负离子模式下进行检测；
[0112]
第一个质量分析器(ms1)的检测范围为100～1000质荷比(m/z)，离子积累时间设定为30ms；
[0113]
第二个质量分析器(ms2)的检测范围为50～1000质荷比(m/z)，离子积累时间设定为20ms；
[0114]
第三个质量分析器(ms3)的检测范围为50～1000质荷比(m/z)，离子积累时间设定为20ms；
[0115]
氮气流速设定为1.5l/min；
[0116]
cdl温度为200℃；
[0117]
加热块温度为200℃；
[0118]
pda检测波长范围为200～700nm；
[0119]
正离子模式下接口电压为4.5kv，负离子模式下接口电压为3.5kv，正负离子模式的检测电压均为1.7kv；
[0120]
干燥气压力设定为100kpa，最大压力为30mpa。
[0121]
针对不同附子制品所测的色谱图中不同化学特征色谱峰进行对比、统计，并与对照品的化学特征色谱峰进行对比，获得不同附子制品的36个hplc化学特征色谱峰，具体见表2。
[0122]
表2不同附子制品的36个hplc化学特征色谱峰
[0123][0124]
表2中，3号色谱峰对应的化学成分为1-乙酰光泽乌头碱或12-乙酰光泽乌头碱，5号色谱峰对应的化学成分为附子灵、10-羟基苯甲酰新乌头原碱，6号色谱峰对应的化学成分为苯甲酰新乌头原碱、尼奥灵(或其异构体10-羟基塔拉乌头胺)，8号色谱峰对应的化学成分为苯甲酰乌头原碱，9号色谱峰对应的化学成分为14-乙酰尼奥灵或其异构体，10号色谱峰对应的化学成分为苯甲酰次乌头原碱，11号色谱峰对应的化学成分为14-乙酰尼奥灵或其异构体，16号色谱峰对应的化学成分为卡米车灵a或其异构体，17号色谱峰对应的化学成分为苯甲酰-3，13-双去氧新乌头原碱、8-乙酰-14-苯甲酰尼奥灵(或其异构体)、trifoliolasine e，18号色谱峰对应的化学成分为脱水苯甲酰新乌头原碱，20号色谱峰对
应的化学成分为新乌头碱，21号色谱峰对应的化学成分为乌头芬碱，22号色谱峰对应的化学成分为翠雀碱或异翠雀碱，23号色谱峰对应的化学成分为乌头碱，24号色谱峰对应的化学成分为乌头卡查敏碱b，25号色谱峰对应的化学成分为脱水苯甲酰次乌头原碱，27号色谱峰对应的化学成分为次乌头碱，29号色谱峰对应的化学成分为3-去氧乌头碱或印乌头碱，30号色谱峰对应的化学成分为唐古替星a或其异构体，31号色谱峰对应的化学成分为展花乌头碱或3，13-双脱氧乌头碱，34号色谱峰对应的化学成分为8-乙酰-14-苯甲酰尼奥灵或其异构体，36号色谱峰对应的化学成分为1-羰基翠雀素或6，14-二乙酰翠雀拉亭。
[0125]
再次针对不同附子制品所测的质谱鉴定结果进行对比、统计，并与对照品的质谱鉴定结果进行对比，获得不同附子制品鉴定的38种化学成分，且这些化学成分均为生物碱，具体见表3。
[0126]
[0127][0128]
从表3中可以看出，煻灰炮附子中共鉴定18个化学成分；生附片中共鉴定24个化学成分，其中14-乙酰尼奥灵或其异构体、乌头卡查敏碱b、翠雀碱或异翠雀碱未在其他加工品
中检出；炒附片中共鉴定22个化学成分，其中10-羟基苯甲酰新乌头原碱、苯甲酰-3，13-双去氧新乌头原碱或其异构体、脱水苯甲酰新乌头原碱、脱水苯甲酰次乌头原碱、唐古替星a或其异构体未在其他加工品中检出；蒸附片中共鉴定17个化学成分，其中1-乙酰光泽乌头碱或12-乙酰光泽乌头碱、苯甲酰去氧乌头原碱未在其他加工品中检出；烘附子中共鉴定20个化学成分，其中尼奥灵或10-羟基塔拉乌头胺、森布星a或森布星b、14-乙酰尼奥灵或其异构体、卡米车灵a或其异构体未在其他加工品中检出。
[0129]
三、多元统计分析
[0130]
在多元分析过程中，舍弃样本量不足的烘附子。
[0131]
在hplc-it-tof-msn研究的基础上，从煻灰炮附子所具有的18个化学特征色谱峰(图6)中选择14个峰面积较大的主要色谱峰(具体见表4)作为分析色谱峰。其中，33号色谱峰虽然具有最高的峰面积，但在hplc-it-tof-msn液质分析中未得到任何母离子和碎片离子信息，故不作为多元统计分析的分析色谱峰。具体分析色谱峰如下：
[0132]
表4 14种煻灰炮附子的分析色谱峰一览
[0133][0134]
表4中色谱峰序号列中数字与表2中色谱峰序号对应
[0135]
分别将不同附子制品(生附片、炒附片、蒸附片)及煻灰炮附子的对应14种分析色谱峰峰面积及6种乌头碱类生物碱含量导入excel表格，添加“label”行并为相同加工品种的不同批次标上同一标签，保存为csv格式的csv文件。将csv文件导入在线分析软件metaboanalyst 4.0(https：//www.metaboanalyst.ca/)的单因素统计分析工具集中，导入时数据类型选择“峰面积”或“浓度”，数据格式选择“未配对列向排布(sample in column)”。数据过滤项选择“none”，标准化处理时，sample normalization、data transformation及data scaling三项均选择“none”，最后选择相应的分析模块分别对上述数据进行主成分分析(pca)、偏最小二乘-判别分析(pls-da)和聚类分析(ca)，聚类分析方法为层次聚类(hierarchical clustering)，距离标准为pearson，聚类算法为ward，分析结束后将分析结果导出为树状图(dendrogram)。
[0136]
1)主成分分析
[0137]
主成分分析(pca)是一种对数据进行旋转变换的统计学方法，其本质是在线性空间中进行一个基变换，使得变换后的数据投影在一组新的“坐标轴”上的方差最大化，随后，裁剪掉变换后方差很小的“坐标轴”，剩下的新“坐标轴”即被称为主成分(principal component)，它们可以在一个较低维度的子空间中尽可能地表示原有数据的性质。
[0138]
通过pca，我们可以将表征附子制品及煻灰炮附子选取的化学物质基础的n个成分的含量(n个维度)，通过适当的函数变换，降维成2个“主成分”(2个维度)，同时尽可能地保留了原有数据的性质。在新的2维空间坐标中，2个样品点的空间距离越近，表明该样品原有
数据(n个成分的含量)的相近成度越高，反之相近程度越低。
[0139]
主成分分析(pca)结果见图8。图8中，绿色数据点表示煻灰炮附子中14种分析色谱峰峰面积数据经处理后在二维坐标空间中的投影，绿色空间区域表示煻灰炮附子中14个化学特征色谱峰峰面积数据经处理后在二维坐标空间中投影的分布范围情况；蓝色数据点表示生附片样品中14种分析色谱峰峰面积数据经处理后在二维坐标空间中的投影，蓝色空间区域表示生附片样品中14种分析色谱峰峰面积数据经处理后在二维坐标空间中投影的分布范围情况；红色数据点表示炒附片样品中14种分析色谱峰峰面积数据经处理后在二维坐标空间中的投影，红色空间区域表示炒附片样品中14种分析色谱峰峰面积数据经处理后在二维坐标空间中投影的分布范围情况；青色数据点表示蒸附片样品中14种分析色谱峰峰面积数据经处理后在二维坐标空间中的投影，青色空间区域表示蒸附片样品中14种分析色谱峰峰面积数据经处理后在二维坐标空间中投影的分布范围情况。
[0140]
在图8中，若两个样品的数据投影点的空间分布越相近，则表示这两个样品在化学成分上的相似性越高(体现在14种分析色谱峰峰面积的相似性上)、差异性越小；反之，若两个样品的数据投影点的空间分布越相离，则表示这两个样品在化学成分上的相似性越低、差异性越大。同理，在图8中，若附子制品与煻灰炮附子的数据投影空间区域越相近、重合度越大，则表示这种附子制品与煻灰炮附子在化学成分上的相似性越高、差异性越小；反之，若这种附子制品与煻灰炮附子的数据投影空间区域越相离、重合度越小，则表示这种附子制品与煻灰炮附子在化学成分上的相似性越低、差异性越大。
[0141]
根据pca结果可以看出，在图8中，炒附片样品的数据投影空间区域与煻灰炮附子的数据投影空间区域最为相近，它们的数据投影空间区域几乎完全重合，表明在3种现代附子制品中，炒附片的化学成分与煻灰炮附子的相似程度最高；生附片样品的数据投影空间区域与煻灰炮附子的数据投影空间区域分布较为相离，它们的数据投影空间区域仅有小部分重合，表明生附片的化学成分与煻灰炮附子有较大差异，相似程度不如炒附片；蒸附片样品的数据投影空间区域与煻灰炮附子的数据投影空间区域最为相离，它们的数据投影空间区域完全没有发生重合，表明在3种现代附子制品中，蒸附片的化学成分与煻灰炮附子的相似程度最低。
[0142]
2)偏最小二乘-判别分析
[0143]
pls-da(partial least squares discriminant analysis)，即偏最小二乘法判别分析，是多变量数据分析技术中的判别分析法，经常用来处理分类和判别问题。通过对主成分适当的旋转，pls-da可以有效的对组间观察值进行区分，并且能够找到导致组间区别的影响变量。该方法与pca类似。
[0144]
同理，根据pls-da结果，参见图9，可以看出，在图9中，炒附片样品的数据投影空间区域与煻灰炮附子的数据投影空间区域最为相近，它们的数据投影空间区域有很大的重合度，表明在3种现代附子制品中，炒附片的化学成分与煻灰炮附子的相似程度最高；在图9中，生附片样品的数据投影空间区域与煻灰炮附子的数据投影空间区域分布较为相离，它们的数据投影空间区域仅有小部分重合，表明生附片的化学成分与煻灰炮附子有较大差异，相似程度不如炒附片；在图9中，蒸附片样品的数据投影空间区域与“煻灰炮”附子的数据投影空间区域最为相离，它们的数据投影空间区域完全没有发生重合，表明在3种现代附子制品中，蒸附片的化学成分与煻灰炮附子的相似程度最低。
[0145]
根据pls-da的结果，我们发现在3种附子现代加工品中，炒附片在化学成分上与煻灰炮附子具有最高的相似程度，此结果与主成分分析(pca)结果一致。
[0146]
图8和图9综合来看，可以看出炒附片的空间分布区域(红)与煻灰炮附子(绿)具有最高的重合度，而蒸附片的空间分布(青)与煻灰炮附子(绿)无重合。
[0147]
vip值是pls-da模型变量的变量权重值，可用于衡量各化学成分含量差异对各组样本分类判别的影响强度和解释能力，vip≥1为常见的差异化学成分筛选标准。对偏最小二乘-判别分析pls-da获得的结果进行vip值排序，排序结果见图10，发现vip值≥1的化学特征色谱峰有3个，降序排列依次为峰20(tr＝41.224min)、峰27(tr＝59.179min)、峰6(tr＝28.279min)。根据hplc-it-tof液质指认结果，峰20为新乌头碱，峰27为次乌头碱，峰6主要为苯甲酰新乌头原碱。上述结果表明，4种附子不同加工品(“煻灰炮”附子、生附片、炒附片、蒸附片)中新乌头碱、次乌头碱、苯甲酰新乌头原碱的含量差异是导致它们具有不同化学成分特征的重要因素，也是区分“煻灰炮”附子、生附片、炒附片、蒸附片的关键成分。其中，炒附子的vip值与煻灰炮附子最为接近。
[0148]
3)聚类分析
[0149]
聚类分析(clusteranalysis)又称群分析，是根据“物以类聚”的道理，对样品或指标进行分类的一种多元统计分析方法，它们讨论的对象是大量的样品，要求能合理地按各自的特性来进行合理的分类，没有任何模式可供参考或依循，即是在没有先验知识的情况下进行的。聚类是将数据分类到不同的类或者簇这样的一个过程，所以同一个簇中的对象有很大的相似性，而不同簇间的对象有很大的相异性。
[0150]
在聚类分析中，衡量样本间相似性和差异性的方法就是计算两个样本之间的分类距离(d)，分类距离越小，样本间的相似性越大，反之相似性越小。
[0151]
聚类分析的分析结果见图11和图12，其中，pfz表示煻灰炮附子，sfz表示生附片，cfz表示炒附片，zfz表示蒸附片，序号表示不同的批次，即pfz-7表示第7批次的煻灰炮附子。
[0152]
图11为14种分析色谱峰峰面积的聚类分析结果，图12为6种乌头碱类生物碱含量的聚类分析结果，可以看出，以14个主要色谱峰hplc峰面积为变量，发现分类距离d＝5时生附片、炒附片与煻灰炮附子可聚为一类，而d＝20时蒸附片可与煻灰炮附子聚为一类；以6种生物碱含量为变量，发现d＝2.5时炒附片可与煻灰炮附子聚为一类，d＝5时生附片可与煻灰炮附子聚为一类，而d＝7.5时蒸附片可与煻灰炮附子聚为一类。因此，炒附片与煻灰炮附子归为一类的分类距离最小。
[0153]
当主成分分析(pca)、偏最小二乘-判别分析(pls-da)和聚类分析(ca)的分析结果一致，且判定结果均为附子制品与煻灰炮附子相似度较高时，即可判定该附子制品(如：炒附子)可代替煻灰炮附子使用；当分析结果不一致，出现附子制品与煻灰炮附子相似度不高时，则判定该附子制品(如：蒸附子、生附子)无法代替煻灰炮附子使用。
[0154]
显然，所描述的实施例仅仅是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。