展浩 言方榮 陸濤

摘要系統(tǒng)闡述了中藥指紋圖譜數(shù)據(jù)分析處理方法（包括指紋圖譜數(shù)據(jù)前期處理方法、相似度評價(jià)方法以及機(jī)器學(xué)習(xí)方法）在指紋圖譜數(shù)據(jù)處理中的應(yīng)用，以期為中藥指紋圖譜數(shù)據(jù)分析提供參考。

關(guān)鍵詞指紋圖譜；機(jī)器學(xué)習(xí)；數(shù)據(jù)分析方法

中圖分類號S126文獻(xiàn)標(biāo)識碼

A文章編號0517-6611（2017）16-0203-03

Research Progress on Data Analysis Methods in Traditional Chinese Medicine Fingerprints

ZHAN Hao，YAN Fangrong，LU Tao*（School of Sciences，China Pharmaceutical University，Nanjing，Jiangsu 211198）

AbstractThe application of the methods of data analysis and processing of traditional Chinese medicine fingerprints（including the fingerprints data pretreatment methods，similarity evaluation method as well as machine learning methods） in the data processing of fingerprints were systematically expounded，so as to provide references for data analysis of traditional Chinese medicine fingerprints.

Key wordsFingerprints； Machine learning； Data analysis methods

基金項(xiàng)目國家社會科學(xué)基金一般項(xiàng)目（16BTJ021）。

作者簡介展浩（1993—），男，安徽宿州人，碩士研究生，研究方向：生物統(tǒng)計(jì)與計(jì)算藥學(xué)。*通訊作者，教授，博士，博士生導(dǎo)師，從事藥物化學(xué)和藥學(xué)信息學(xué)研究。

收稿日期2017-04-12

近年來，中藥及中藥制劑在世界范圍內(nèi)的需求增長迅速，因此中藥質(zhì)量的評價(jià)和控制顯得尤為重要。然而，由于中藥具有組分復(fù)雜、靶點(diǎn)較多、整體協(xié)同作用等特點(diǎn)，使得中藥質(zhì)量控制與評價(jià)成為制約中藥發(fā)展的一大難題。傳統(tǒng)的針對中藥單一化學(xué)成分或其中幾個(gè)成分的分析已無法準(zhǔn)確對中藥質(zhì)量進(jìn)行整體評價(jià)[1]。

隨著生命科學(xué)以及分析化學(xué)技術(shù)的進(jìn)步和發(fā)展，中藥質(zhì)量控制領(lǐng)域的研究手段越來越豐富。指紋圖譜技術(shù)在中藥優(yōu)劣評價(jià)、真?zhèn)舞b別、活性成分識別以及一致性和穩(wěn)定性評價(jià)方面已成為國際公認(rèn)有效方法。中藥指紋圖譜是指利用質(zhì)量分析及信息采集等手段對中藥材或中成藥進(jìn)行處理而得到的能夠顯示其性質(zhì)的色譜、光譜等數(shù)據(jù)。常用的分析手段有紅外光譜法（IR）、紫外光譜法（UV）、核磁共振波譜法（NMR）、薄層色譜法（TLC）、高效液相色譜法（HPLC）、氣相色譜法（GC）等[2]。然而，中藥指紋圖譜技術(shù)的出現(xiàn)隨之而來的是大量的化學(xué)測量數(shù)據(jù)，如何通過有效的數(shù)據(jù)分析手段對指紋圖譜數(shù)據(jù)進(jìn)行分析，進(jìn)而挖掘出其中有價(jià)值信息已成為中藥指紋圖譜應(yīng)用于中藥質(zhì)量控制中所面臨的新難題。在中藥指紋圖譜數(shù)據(jù)處理過程中，由于圖譜十分復(fù)雜，常見的情況是將指紋圖譜數(shù)據(jù)處理與計(jì)算機(jī)圖譜解析、統(tǒng)計(jì)學(xué)、機(jī)器學(xué)習(xí)等技術(shù)相結(jié)合[3]。筆者對現(xiàn)有的指紋圖譜數(shù)據(jù)分析處理方法進(jìn)行了系統(tǒng)闡述，旨在為研究者在相關(guān)領(lǐng)域的研究提供參考，并選擇適合的研究工具。

1中藥指紋圖譜前期校正方法

指紋圖譜相關(guān)研究通?？煞譃?個(gè)部分：指紋圖譜生成以及指紋圖譜的分析計(jì)算。在指紋圖譜的生成過程中，由于試驗(yàn)條件的限制以及儀器的差別，經(jīng)常導(dǎo)致所獲得的指紋圖譜都會有一定程度的噪聲、基線漂移、保留時(shí)間漂移等現(xiàn)象，因此指紋圖譜在進(jìn)行分析計(jì)算前，必須對所獲得的指紋圖譜數(shù)據(jù)進(jìn)行前期處理，從而解決上述問題，使后續(xù)的分析計(jì)算更加準(zhǔn)確有效[4]。

1.1基于小波變換的基線校正和噪聲去除

小波變換（Wavelet transform ，WT）是一種基于時(shí)頻域的信號處理方法，該方法具有時(shí)頻局部性良好、選基靈活等優(yōu)點(diǎn)[5]。在中藥指紋圖譜數(shù)據(jù)處理過程中，通常使用離散小波變換方法， Mallat快速離散小波變換算法如下：

f（t）=kzcJ，kφJ(rèn)，k（t）+Jj=1

kzdj，kΨJ，k（t） （1）

式中，cJ，k表示f（t）在第J+1個(gè)頻段上第k個(gè)時(shí)段上的分量，也就是圖譜的低頻部分，dj，k表示f（t）在第j個(gè)頻段上的第k個(gè)時(shí)段的分量（1≤j≤J），代表圖譜的高頻部分。

圖譜數(shù)據(jù)的信號噪聲一般集中于高頻段的小波系數(shù)中，基線漂移一般表現(xiàn)在低頻信號中，所以在對圖譜基線進(jìn)行校正時(shí)一般是對低頻段處理，信號噪聲的去除一般是針對高頻段的處理，二者之間相互獨(dú)立。通過將圖譜數(shù)據(jù)在低頻段的小波系數(shù)置為0，可解決基線漂移問題，將較高頻段的小波系數(shù)進(jìn)行閾值處理可實(shí)現(xiàn)圖譜噪聲去除[6]。

在處理指紋圖譜噪聲和基線漂移問題的研究中，邵利民等[7]使用Haar小波對乳酸-稀土絡(luò)合物圖譜進(jìn)行處理，成功去除了其信號中的噪聲，獲得了光滑的色譜曲線。袁海龍等[8]運(yùn)用Daubechies8小波方法處理茵陳注射液HPLC數(shù)據(jù)，成功解決了色譜的基線漂移，得到了平直的基線。

1.2保留時(shí)間漂移校正

從理論上講，多種樣本經(jīng)同一色譜儀分析，在得到的色譜圖中相同物質(zhì)應(yīng)該有相同的保留時(shí)間。然而，因?yàn)樵囼?yàn)條件等因素的限制，不同樣本使用同一色譜儀所測得的色譜圖中，相同物質(zhì)的保留時(shí)間往往存在漂移現(xiàn)象。目前已提出的保留時(shí)間漂移校正算法主要包括相關(guān)最優(yōu)化調(diào)整（Correlation optimized warping，COW）[9]、基于遺傳算法的校準(zhǔn)方法（Peak alignment by a genetic algorithm，PAGA）[10]、動(dòng)態(tài)時(shí)間調(diào)整（Dynamic time warping）[11]、參數(shù)化時(shí)間調(diào)整（Parametric time warping，PTW）[12]、半?yún)?shù)化時(shí)間調(diào)整（Semiparametric time warping，STW）[13]、模糊規(guī)整（Fuzzy alignment，F(xiàn)W）[14]等，其中相關(guān)最優(yōu)化調(diào)整、動(dòng)態(tài)時(shí)間調(diào)整以及基于遺傳算法的校準(zhǔn)方法在指紋圖譜保留時(shí)間漂移校正方面的應(yīng)用尤為廣泛。LucioGutiérrez等[15]在對38批特納草HPLC指紋圖譜數(shù)據(jù)使用偏最小二乘回歸模型建模之前，采用相關(guān)最優(yōu)化調(diào)整方法對其進(jìn)行保留時(shí)間漂移校正，偏峰問題明顯得到解決。

2基于指紋圖譜數(shù)據(jù)的中藥相似度計(jì)算方法

相似度是指中藥指紋圖譜之間的相似程度，已被認(rèn)定成為評價(jià)中藥指紋圖譜的重要指標(biāo)。常用的中藥相似度計(jì)算方法包括以下4種：夾角余弦法（Vector cosine）、相關(guān)系數(shù)法（Correlation coefficient）、峰重疊率法（Nei系數(shù)法）以及峰重疊率與共有峰強(qiáng)度結(jié)合法（改進(jìn)Nei系數(shù)法）。

2.1夾角余弦法

中藥指紋圖譜數(shù)據(jù)是由各保留時(shí)間對應(yīng)的峰面積或峰高的數(shù)值組成的一個(gè)n維向量X=[x1，x2，…，xn]，其中n為譜峰數(shù)或原始數(shù)據(jù)點(diǎn)數(shù)。評價(jià)2個(gè)中藥的相似度，就可以通過計(jì)算2個(gè)n維向量X1和X2的相似度來實(shí)現(xiàn)。在評價(jià)向量相似度時(shí)，一般使用向量間夾角余弦值來表示，其夾角余弦值越趨近于1意味著向量之間相似度越高，反之相似度越低。向量X1和X2之間夾角余弦值的計(jì)算公式如下：

r（X1，X2）=cos（X1X2）=X′1X2（X′1X2）（X′2X1） （2）

王龍星等[16]通過計(jì)算11個(gè)不同產(chǎn)地及炮制方法的吳茱萸樣品指紋圖譜之間的夾角余弦值來評價(jià)樣品之間的相似度，結(jié)果顯示該方法可反映樣品之間的相似性。

2.2相關(guān)系數(shù)法

相關(guān)系數(shù)法是使用向量之間的相關(guān)系數(shù)來反映向量之間相似程度的方法，其本質(zhì)上是對夾角余弦法的中心化變換，其計(jì)算公式如下：

r（X1，X2）=（X1-1）′（X2-2）（X1-1）′（X1-1）（X2-2）′（X2-2） （3）

聶磊等[17]使用相關(guān)系數(shù)法等4種相似度評價(jià)算法對痛必定粉針指紋圖譜進(jìn)行相似度評價(jià)，結(jié)果顯示相關(guān)系數(shù)法對大峰缺失表現(xiàn)得比較敏感，而對小峰缺失問題不夠敏感。

2.3Nei系數(shù)法與改進(jìn)Nei系數(shù)法

Nei系數(shù)法原是針對個(gè)體間遺傳相似性評價(jià)所提出的方法，現(xiàn)在應(yīng)用于中藥指紋圖譜數(shù)據(jù)中，用于評價(jià)中藥之間的相似度，其計(jì)算公式如下：

r=2n0n1+n2×100%（4）

式中，n0表示兩待分析圖譜之間的共有峰的數(shù)量，n1和n2表示兩待分析圖譜的所有譜峰數(shù)。張聰?shù)萚18]對11批紅參甲醇提取液的指紋圖譜使用Nei系數(shù)法計(jì)算其重疊率，通過對八強(qiáng)峰的分析比較，結(jié)果發(fā)現(xiàn)國產(chǎn)紅參與高麗紅參品質(zhì)相近。但是，Nei系數(shù)法只是考慮到圖譜共有峰的數(shù)量問題，并未考慮共有峰其峰強(qiáng)度的影響。孟慶華等[19]在峰重疊率的基礎(chǔ)上引入峰強(qiáng)度的信息提出了改進(jìn)Nei系數(shù)法，其公式如下：

r=2n0n1+n2-2n1+n2h10-h20

h10+h20（5）

式中，h10、h20分別為2個(gè)待分析圖譜之間共有峰的峰強(qiáng)度。

3機(jī)器學(xué)習(xí)方法在指紋圖譜數(shù)據(jù)中的應(yīng)用

機(jī)器學(xué)習(xí)（Machine learning，ML）是一類可以從數(shù)據(jù)中自動(dòng)學(xué)習(xí)獲得規(guī)律，并且利用學(xué)習(xí)到的規(guī)律對未知數(shù)據(jù)進(jìn)行預(yù)測的算法。根據(jù)樣本數(shù)據(jù)有無標(biāo)記，機(jī)器學(xué)習(xí)算法可分為無監(jiān)督學(xué)習(xí)和監(jiān)督學(xué)習(xí)2大類。根據(jù)輸出變量是否為連續(xù)變量，監(jiān)督學(xué)習(xí)方法又可分為分類算法和回歸算法。

3.1無監(jiān)督學(xué)習(xí)

無監(jiān)督學(xué)習(xí)方法是一類用于無標(biāo)記數(shù)據(jù)的機(jī)器學(xué)習(xí)方法，主要包括主成分分析法（Principal component analysis，PCA）和聚類分析（Cluster analysis，CA）等。

3.1.1主成分分析法。

PCA可用于提取主成分，常用于高維數(shù)據(jù)的降維、變量的選擇以及去除變量間的相關(guān)性。由于中藥指紋圖譜數(shù)據(jù)往往為高維數(shù)據(jù)，因此使用PCA對其降維十分必要。在數(shù)據(jù)處理中PCA常與其他分類或回歸方法聯(lián)用，將其提取的主成分作為其他模型的輸入變量。馮慧萍等[20]對產(chǎn)自浙江麗水的15個(gè)厚樸樣品建立了HPLC-DAD指紋圖譜，利用PCA對其共有峰的相對峰面積進(jìn)行分析，并將前4個(gè)主成分的得分作為人工神經(jīng)網(wǎng)絡(luò)的輸入變量，厚樸樣品清除活性氧指標(biāo)作為輸出變量，建立厚樸指紋圖譜與其清除活性氧能力之間的神經(jīng)網(wǎng)絡(luò)。結(jié)果表明，使用主成分分析+神經(jīng)網(wǎng)絡(luò)技術(shù)構(gòu)建的網(wǎng)絡(luò)可準(zhǔn)確預(yù)測厚樸清除活性氧能力的IC50值。

3.1.2聚類分析。

聚類分析又稱群集分析，它是根據(jù)“物以類聚”的思想把相似的對象通過靜態(tài)分類的方法分成不同的組別或者更多的子集，使在同一個(gè)子集中的樣本都具有一些相似的屬性。該方法通過對中藥指紋圖譜數(shù)據(jù)處理，可解決中藥品種的分類，真?zhèn)舞b別、質(zhì)量評價(jià)、新舊工藝或不同炮制方法比較等問題。田蘭等[21]用各色譜峰對內(nèi)標(biāo)峰的相對保留時(shí)間定性獲得45個(gè)色譜峰，計(jì)算出各色譜峰相對于內(nèi)標(biāo)物的含量得到45個(gè)特征，并對所得特征進(jìn)行標(biāo)準(zhǔn)化后使用聚類分析等方法，發(fā)現(xiàn)白術(shù)優(yōu)質(zhì)品種主要集中在浙江、湖南、四川一帶。

3.2監(jiān)督學(xué)習(xí)

3.2.1分類算法。分類算法適用于樣本標(biāo)記為分類變量的數(shù)據(jù)，常用的分類算法主要包括K近鄰法（knearest neighbor，kNN）、支持向量機(jī)（Support vector machine，SVM）、人工神經(jīng)網(wǎng)絡(luò)（Artificial neural network，ANN）等，其在中藥指紋圖譜數(shù)據(jù)中一般用于解決中藥材的分類及歸屬問題。

3.2.1.1K近鄰法。

K近鄰法是由Cover與Hart于1967年提出的一種分類方法，其算法思想是若1個(gè)樣本的k個(gè)與其最相鄰的樣本大多屬于某一個(gè)類別，則該樣本也屬于這個(gè)類別，并且具有該類別樣本的某些特性[22]。在中藥指紋圖譜數(shù)據(jù)處理中，首先按照指紋圖譜特征參數(shù)相似度最相似的原則，取未知樣本的k個(gè)最近鄰樣本，然后所取的k個(gè)近鄰樣本大多數(shù)的歸類即為該未知樣本的類別。莊花等[23]使用K近鄰法對63個(gè)由安徽、山東、河北、貴州4省所生產(chǎn)的合歡花樣品的指紋圖譜數(shù)據(jù)及來源數(shù)據(jù)建立分類模型，實(shí)現(xiàn)了南北方樣品正確分類。

3.2.1.2支持向量機(jī)。

支持向量機(jī)是由Corinna Cortes和Vapnik等于1995年首先提出的解決小樣本、非線性及高維數(shù)據(jù)等問題的機(jī)器學(xué)習(xí)算法。當(dāng)數(shù)據(jù)線性可分時(shí)，該方法通過樣本數(shù)據(jù)訓(xùn)練學(xué)習(xí)線性支持向量機(jī)，找出間隔最大的的“超平面”，將樣本分類；當(dāng)數(shù)據(jù)線性不可分時(shí)，通過使用核函數(shù)學(xué)習(xí)非線性支持向量機(jī)，相當(dāng)于隱式的在高維特征空間中學(xué)習(xí)線性支持向量機(jī)實(shí)現(xiàn)樣本分類。王曉燕等[24]對10種寒性中藥和10種熱性中藥提取和精制多糖，使用三氟乙酸將得到的多糖水解成單糖，并測定了多糖的單糖組成HPLC指紋圖譜；然后利用支持向量機(jī)對指紋圖譜數(shù)據(jù)及藥性數(shù)據(jù)建立分類模型，實(shí)現(xiàn)了對20種中藥寒、熱藥性的準(zhǔn)確分類，分類正確率達(dá)100%。

3.2.1.3人工神經(jīng)網(wǎng)絡(luò)。

人工神經(jīng)網(wǎng)絡(luò)是一種抽象人腦神經(jīng)的一類機(jī)器學(xué)習(xí)算法，它是由大量的節(jié)點(diǎn)之間相互聯(lián)接構(gòu)成。其中，每個(gè)節(jié)點(diǎn)代表1個(gè)激勵(lì)函數(shù)（Activation function），每2個(gè)節(jié)點(diǎn)間的連接都代表一個(gè)對于通過該連接信號的權(quán)重。應(yīng)用最廣泛的神經(jīng)網(wǎng)絡(luò)為BP神經(jīng)網(wǎng)絡(luò)，其在結(jié)構(gòu)上分為3層：輸入層、隱含層和輸出層。李芳等[25]使用BP神經(jīng)網(wǎng)絡(luò)對采集的12個(gè)主產(chǎn)縣90批黃芪藥材的FTIR指紋圖譜建立了產(chǎn)地的鑒別及預(yù)測模型，模型成功地對測試集中的黃芪藥材產(chǎn)地進(jìn)行了預(yù)測，其預(yù)測準(zhǔn)確率達(dá)83%。

3.2.2回歸算法。

回歸算法適用于樣本標(biāo)記為連續(xù)變量的數(shù)據(jù)，在處理高維數(shù)據(jù)時(shí)，常用算法為偏最小二乘回歸（Partial least squares regression，PLSR）。在處理中藥指紋圖譜數(shù)據(jù)時(shí)，一般結(jié)合相應(yīng)的藥效數(shù)據(jù)，根據(jù)定量組效關(guān)系（Quantitative composition activity relationship，QCAR）來發(fā)現(xiàn)中藥材中相應(yīng)的活性成分信息[26]。

偏最小二乘回歸是一種將主成分分析、回歸分析及典型相關(guān)分析有機(jī)結(jié)合起來的一種多元線性回歸方法。該方法較好地解決樣本個(gè)數(shù)少于變量個(gè)數(shù)時(shí)的建模問題，同時(shí)也克服了變量間多重共線性問題，因此非常適用于對中藥指紋圖譜數(shù)據(jù)進(jìn)行回歸分析。在使用偏最小二乘回歸預(yù)測中藥材活性成分時(shí)，首先根據(jù)指紋圖譜數(shù)據(jù)與相應(yīng)藥效數(shù)據(jù)建立回歸模型，然后使用顯著性檢驗(yàn)計(jì)算回歸系數(shù)顯著性，顯著系數(shù)所對應(yīng)的峰即為潛在的活性成分。Liu等[27]通過偏最小二乘回歸對來自29個(gè)省份的山香圓葉HPLC指紋圖譜數(shù)據(jù)及其抗氧化活性指標(biāo)進(jìn)行分析，構(gòu)建了偏最小二乘回歸模型，通過對回歸系數(shù)進(jìn)行jackknife檢驗(yàn)，成功預(yù)測出山香圓葉中7種抗氧化活性成分。

45卷16期展 浩等中藥指紋圖譜數(shù)據(jù)分析方法研究進(jìn)展

4小結(jié)與展望

目前，中藥指紋圖譜發(fā)展尚處于初級階段，數(shù)據(jù)量少、數(shù)據(jù)處理過程中分析方法的不合理應(yīng)用制約著其在中藥質(zhì)量控制等相關(guān)領(lǐng)域的發(fā)展。隨著分析化學(xué)手段的不斷豐富，計(jì)算機(jī)科學(xué)與數(shù)據(jù)科學(xué)的不斷發(fā)展，將會產(chǎn)生更多的指紋圖譜數(shù)據(jù)以及更加精確的指紋圖譜數(shù)據(jù)分析方法。從長遠(yuǎn)來看，這種趨勢不僅將加速中藥指紋圖譜數(shù)據(jù)在中藥質(zhì)量控制方面的應(yīng)用，而且會促進(jìn)其在中藥藥效、體內(nèi)代謝變化以及臨床用藥療效等研究領(lǐng)域的應(yīng)用與發(fā)展。

參考文獻(xiàn)

[1]

朱曉勤，尹蓮.化學(xué)模式識別在中藥分類和質(zhì)量評價(jià)中的應(yīng)用[J].中華中醫(yī)藥學(xué)刊，2007，25（10）：2111-2114.

[2] 袁琴琴.中藥材指紋圖譜研究進(jìn)展[J].安徽農(nóng)業(yè)科學(xué)，2017，45（4）：132-134.

[3] 屈景輝，廖琪梅，張星.指紋圖譜數(shù)據(jù)庫建立技術(shù)[J].醫(yī)學(xué)信息，2006，19（2）：190-191.

[4] 張鋒.中藥指紋圖譜數(shù)據(jù)預(yù)處理及相似度計(jì)算評價(jià)軟件的開發(fā)[D].廣州：華南理工大學(xué)，2006.

[5] 朱來東，廉小親，江遠(yuǎn)志.小波變換在信號降噪中的應(yīng)用及MATLAB實(shí)現(xiàn)[J].北京工商大學(xué)學(xué)報(bào)（自然科學(xué)版），2009，27（2）：46-49.

[6] 方勇華，孔超，蘭天鴿，等.應(yīng)用小波變換實(shí)現(xiàn)光譜的噪聲去除和基線校正[J].光學(xué)精密工程，2006，14（6）：1088-1092.

[7] 邵利民，唐兵，邵學(xué)廣，等.小波變換用于高效液相色譜的噪聲濾除[J].分析化學(xué)，1997，25（1）：15-18.

[8] 袁海龍，雷長海，肖小河，等.小波變換校正茵陳注射液HPLC指紋圖譜基線的研究[J].中國新醫(yī)藥，2003，2（9）：13-14.

[9] NIELSEN N P V，CARSTENSEN J M，SMEDSGAARD J.Aligning of single and multiple wavelength chromatographic profiles for chemometric data analysis using correlation optimised warping[J].Journal of chromatography A，1998，805（1/2）：17-35.

[10] FORSHED J，SCHUPPEKOISTINEN I，JACOBSSON S P.Peak alignment of NMR signals by means of a genetic algorithm[J].Analytica chimica acta，2003，487（2）：189-199.

[11] KASSIDAS A，MACGREGOR J F，TAYLOR P A.Synchronization of batch trajectories using dynamic time warping[J].AIChE Journal，1998，44（4）：864-875.

[12] EILERS P H C.Parametric time warping[J].Analytical chemistry，2004，76（2）：404-411.

[13] VAN NEDERKASSEL A M，XU C J，LANCELIN P，et al.Chemometric treatment of vanillin fingerprint chromatograms.Effect of different signal alignments on principal component analysis plots[J].Journal of chromatography A，2006，1120（1/2）：291-298.

[14] WALCZAK B，WU W.Fuzzy warping of chromatograms[J].Chemometrics & intelligent laboratory systems，2005，77（1/2）：173-180.

[15] LUCIOGUTIRREZ J R，GARZAJUREZ A，COELLO J，et al.Multiwavelength highperformance liquid chromatographic fingerprints and chemometrics to predict the antioxidant activity of Turnera diffusa，as part of its quality control[J].Journal of chromatography A，2012，1235（8）：68-76.

[16] 王龍星，肖紅斌，梁鑫淼，等.一種評價(jià)中藥色譜指紋譜相似性的新方法：向量夾角法[J].藥學(xué)學(xué)報(bào)，2002，37（9）：713-717.

[17] 聶磊，曹進(jìn)，羅國安，等.中藥指紋圖譜相似度評價(jià)方法的比較[J].中成藥，2005，27（3）：249-252.

[18] 張聰，王智華，金德莊.中國紅參與高麗紅參的指紋譜（HPLC-FPS）比較研究[J].中成藥，2001，23（3）：160-163.

[19] 孟慶華， 劉永鎖， 王健松，等. 色譜指紋圖譜相似度的新算法及其應(yīng)用[J]. 中成藥， 2003， 25（1）：4-8.

[20] 馮慧萍， 楊中林， 尹小英，等. 應(yīng)用PCA-ANN技術(shù)研究中藥厚樸色譜指紋譜與抗氧化活性的關(guān)系[J]. 中國天然藥物， 2005， 3（6）：377-381.

[21] 田蘭，畢開順，孫穩(wěn)健，等.白術(shù)的化學(xué)模式識別[J].中國中藥雜志，2003，28（2）：143-146.

[22] COVER T，HART P.Nearest neighbor pattern classification[J].IEEE Transactions on Information Theory，1967，13（1）：21-27.

[23] 莊花，倪永年.不同產(chǎn)地的合歡花藥材高效液相指紋圖譜[J].南昌大學(xué)學(xué)報(bào)（理科版），2013，37（1）：59-63.

[24] 王曉燕，李峰.20種中藥多糖水解成分HPLC指紋圖譜與寒熱藥性關(guān)系的SVM分析研究[J].山東中醫(yī)藥大學(xué)學(xué)報(bào)，2012，36（5）：439-442.

[25] 李芳，李欽，顧志榮，等.基于FTIR的黃芪藥材產(chǎn)地鑒別及預(yù)測研究[J].天然產(chǎn)物研究與開發(fā)，2016（11）：1712-1717.

[26] 周立東.建議在天然藥物研究中建立“定量組效關(guān)系”（QCAR）概念[J].世界科學(xué)技術(shù)：中醫(yī)藥現(xiàn)代化，1999（2）：33-34.

[27] LIU X，ZHAN H，QIAO Z，et al.Chemometric analysis based on HPLC multiwavelength fingerprints for prediction of antioxidant components in Turpiniae Folium[J].Chemometrics & intelligent laboratory systems，2016，152：54-61.