CN106990214A - 一种评价中药材质量的方法 - Google Patents
一种评价中药材质量的方法 Download PDFInfo
- Publication number
- CN106990214A CN106990214A CN201710316615.9A CN201710316615A CN106990214A CN 106990214 A CN106990214 A CN 106990214A CN 201710316615 A CN201710316615 A CN 201710316615A CN 106990214 A CN106990214 A CN 106990214A
- Authority
- CN
- China
- Prior art keywords
- chinese medicine
- traditional chinese
- fingerprint
- quality
- similarity
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000003814 drug Substances 0.000 title claims abstract description 104
- 238000000034 method Methods 0.000 title claims abstract description 55
- 239000000463 material Substances 0.000 claims abstract description 36
- 239000003795 chemical substances by application Substances 0.000 claims abstract description 31
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 12
- 239000000126 substance Substances 0.000 claims description 29
- 238000001228 spectrum Methods 0.000 claims description 25
- 239000011159 matrix material Substances 0.000 claims description 20
- 238000005516 engineering process Methods 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 14
- 230000004927 fusion Effects 0.000 claims description 13
- 239000013598 vector Substances 0.000 claims description 12
- 239000003550 marker Substances 0.000 claims description 8
- 238000004809 thin layer chromatography Methods 0.000 claims description 8
- 238000004817 gas chromatography Methods 0.000 claims description 7
- 239000002028 Biomass Substances 0.000 claims description 6
- 238000004566 IR spectroscopy Methods 0.000 claims description 6
- 238000012549 training Methods 0.000 claims description 6
- 238000000870 ultraviolet spectroscopy Methods 0.000 claims description 6
- 238000003745 diagnosis Methods 0.000 claims description 3
- 230000002194 synthesizing effect Effects 0.000 claims description 2
- 229940079593 drug Drugs 0.000 abstract description 6
- 238000001514 detection method Methods 0.000 abstract description 5
- 238000011156 evaluation Methods 0.000 abstract description 5
- 238000013441 quality evaluation Methods 0.000 abstract description 5
- 238000002156 mixing Methods 0.000 abstract description 2
- 230000001737 promoting effect Effects 0.000 abstract description 2
- 206010037660 Pyrexia Diseases 0.000 abstract 2
- 230000015572 biosynthetic process Effects 0.000 abstract 1
- 238000011524 similarity measure Methods 0.000 abstract 1
- 238000003786 synthesis reaction Methods 0.000 abstract 1
- 238000004458 analytical method Methods 0.000 description 13
- 239000000523 sample Substances 0.000 description 10
- OKKJLVBELUTLKV-UHFFFAOYSA-N Methanol Chemical compound OC OKKJLVBELUTLKV-UHFFFAOYSA-N 0.000 description 9
- 229940126680 traditional chinese medicines Drugs 0.000 description 8
- 230000007547 defect Effects 0.000 description 7
- WEVYAHXRMPXWCK-UHFFFAOYSA-N Acetonitrile Chemical compound CC#N WEVYAHXRMPXWCK-UHFFFAOYSA-N 0.000 description 6
- 230000006399 behavior Effects 0.000 description 6
- 230000008859 change Effects 0.000 description 6
- 230000000694 effects Effects 0.000 description 6
- 230000006870 function Effects 0.000 description 6
- 238000003909 pattern recognition Methods 0.000 description 6
- 239000000047 product Substances 0.000 description 6
- 238000013528 artificial neural network Methods 0.000 description 5
- 235000008216 herbs Nutrition 0.000 description 5
- 230000014759 maintenance of location Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000002360 preparation method Methods 0.000 description 4
- 238000000513 principal component analysis Methods 0.000 description 4
- 238000003908 quality control method Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 230000036541 health Effects 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 241000132012 Atractylodes Species 0.000 description 2
- NBIIXXVUZAFLBC-UHFFFAOYSA-N Phosphoric acid Chemical compound OP(O)(O)=O NBIIXXVUZAFLBC-UHFFFAOYSA-N 0.000 description 2
- 230000009471 action Effects 0.000 description 2
- 239000004480 active ingredient Substances 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000004587 chromatography analysis Methods 0.000 description 2
- 238000007621 cluster analysis Methods 0.000 description 2
- 150000001875 compounds Chemical class 0.000 description 2
- 238000003066 decision tree Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000000857 drug effect Effects 0.000 description 2
- 230000005518 electrochemistry Effects 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 238000004128 high performance liquid chromatography Methods 0.000 description 2
- 238000002329 infrared spectrum Methods 0.000 description 2
- 238000004811 liquid chromatography Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 108090000623 proteins and genes Proteins 0.000 description 2
- 230000024053 secondary metabolic process Effects 0.000 description 2
- 238000002211 ultraviolet spectrum Methods 0.000 description 2
- 229940126673 western medicines Drugs 0.000 description 2
- 244000061520 Angelica archangelica Species 0.000 description 1
- 241000196324 Embryophyta Species 0.000 description 1
- 235000011201 Ginkgo Nutrition 0.000 description 1
- 244000194101 Ginkgo biloba Species 0.000 description 1
- 235000008100 Ginkgo biloba Nutrition 0.000 description 1
- 239000009429 Ginkgo biloba extract Substances 0.000 description 1
- 235000001287 Guettarda speciosa Nutrition 0.000 description 1
- 241000241413 Propolis Species 0.000 description 1
- GAMYVSCDDLXAQW-AOIWZFSPSA-N Thermopsosid Natural products O(C)c1c(O)ccc(C=2Oc3c(c(O)cc(O[C@H]4[C@H](O)[C@@H](O)[C@H](O)[C@H](CO)O4)c3)C(=O)C=2)c1 GAMYVSCDDLXAQW-AOIWZFSPSA-N 0.000 description 1
- 239000013543 active substance Substances 0.000 description 1
- 229910000147 aluminium phosphate Inorganic materials 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000003153 chemical reaction reagent Substances 0.000 description 1
- 238000013375 chromatographic separation Methods 0.000 description 1
- 238000007635 classification algorithm Methods 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 238000007418 data mining Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000004438 eyesight Effects 0.000 description 1
- 239000012467 final product Substances 0.000 description 1
- 229930003944 flavone Natural products 0.000 description 1
- 150000002212 flavone derivatives Chemical class 0.000 description 1
- 235000011949 flavones Nutrition 0.000 description 1
- 229930003935 flavonoid Natural products 0.000 description 1
- 150000002215 flavonoids Chemical class 0.000 description 1
- 235000017173 flavonoids Nutrition 0.000 description 1
- 238000007667 floating Methods 0.000 description 1
- 238000002290 gas chromatography-mass spectrometry Methods 0.000 description 1
- 238000007429 general method Methods 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 229940068052 ginkgo biloba extract Drugs 0.000 description 1
- 235000020686 ginkgo biloba extract Nutrition 0.000 description 1
- 229930184727 ginkgolide Natural products 0.000 description 1
- 150000004676 glycans Chemical class 0.000 description 1
- 238000003306 harvesting Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 150000002596 lactones Chemical class 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000007477 logistic regression Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 238000012567 pattern recognition method Methods 0.000 description 1
- 230000000144 pharmacologic effect Effects 0.000 description 1
- 239000006187 pill Substances 0.000 description 1
- 229920001282 polysaccharide Polymers 0.000 description 1
- 239000005017 polysaccharide Substances 0.000 description 1
- 238000012847 principal component analysis method Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 229940069949 propolis Drugs 0.000 description 1
- 102000004169 proteins and genes Human genes 0.000 description 1
- 239000008213 purified water Substances 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 239000012488 sample solution Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000002791 soaking Methods 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 238000013112 stability test Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 239000012085 test solution Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
- VHBFFQKBGNRLFZ-UHFFFAOYSA-N vitamin p Natural products O1C2=CC=CC=C2C(=O)C=C1C1=CC=CC=C1 VHBFFQKBGNRLFZ-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/15—Medicinal preparations ; Physical properties thereof, e.g. dissolubility
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Food Science & Technology (AREA)
- Analytical Chemistry (AREA)
- Pharmacology & Pharmacy (AREA)
- Biophysics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Medicinal Chemistry (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Medicines Containing Plant Substances (AREA)
Abstract
本发明公开了一种评价中药材质量的方法,包括以下步骤:对中药获取的指纹图谱釆用偏最小二乘判别模型进行处理,获取相应的中药寒热特征标记,针对这些中药寒热特征标记,采用Multi‑Agent融合算法进行分类,采用改进的灰色综合关联度进行指纹图谱的相似度计算。本发明能够定性和定量评定图谱间的相似度,进而评价中药材质量,并比较其稳定性、一致性对实验结果进行鉴定、评价,是一种评价中药指纹图谱的有效手段,为云南道地药材的质量评价检测提供相应依据。参照评价结果能够科学、准确地判断相关药材的真伪、道地性、质量的优劣。这对于民族药质量标准化、走向国际市场起到推动作用。
Description
技术领域
本发明涉及医药领域,具体是一种评价中药材质量的方法。
背景技术
中药是中国医药学宝库中的一颗明珠,是中华民族的瑰宝,为中华民族的健康和繁衍做出了巨大贡献,在几千年的医疗实践中逐渐形成了自己独特的理论体系,其卓越的临床疗效众口皆碑,并远传东南亚乃至全世界。随着21世纪生命科学的蓬勃发展与人们崇尚回归自然的热潮,同时,西药的研制难度越来越大,周期越来越长,投人越来越多,迫使许多国家将目光转移至中药,这给中药走出国门带来了难得的机遇。
虽然我国对中药的应用已有很久的历史,但是我国中药的质量标准体系还不够完善,质量检测方法及控制技术比较落后。在中药材的鉴别方面,目前从以人为经验鉴别为主,发展到基原鉴别、性状鉴别、理化鉴别等鉴别方法。现在随着药材鉴别科学化,我们采用更科学的方法对药材质量进行鉴别,对其有效成分进行分析。利用紫外光谱(UV)、红外光谱(IR)、薄层色谱(TLC),气相色谱(GC),电化学等手段获取药材成分的指纹谱图,对谱图预处理获取特征值信息,通过特征值信息采用各种分类器分类方法组和,结合总体方法分析以便于对与对药材产地、质量、地道性进行控制。
以指纹图谱作为中药(天然药物)提取物及其制剂的质量控制方法,已成为目前国际共识,各种符合中药(天然药物)特色的指纹图谱控制技术体系正在研究和建立。美国食品药品管理局(FDA)允许草药保健品申报资料中提供色谱指纹图谱;世界卫生组织(WHO)在1996年草药评价指导原则中也规定,如果草药的活性成分不明确,可以提供色谱指纹图谱以证明产品质量的一致;欧共体在草药质量指南中亦称,单靠测定某种有效成分考查质量的稳定性是不够的,因为草药及其制剂是以整体为活性物质。色谱指纹图谱尤其是薄层色谱的鲜明的指纹图谱是很有用的。国外指纹图谱的应用,目的在于解决成分复杂,有效成分不明确的植物药质量检测和产品批次间质量差异的问题。其中德国研制的银杏叶提取物制剂是一个突出的例子。他们应用指纹图谱制定了相应的标准,该图谱体现了制剂所含的33个化学成分(主要为黄酮类和内酯类)和各自的含量。经化学成分和药效相关性研究,发现约24%银杏黄酮和约6%银杏内酯组成的提取物具有最佳疗效。此外,采用“混批勾兑”法,可使最终产品质量稳定,指纹图谱重现性良好,含量浮动范围为5%。
在国内一些专家提出了一些应用指纹图谱对中药材质量控制的方法,主要有:
(1)相似度评价法
相似度可以量化描述指纹图谱相似性,即客观地数字化描述。相似度分析常采用的方法是相关系数和夹角余弦。相关系数是测度样品间在特征变量的变化模式上相似形状的相似度,它忽略了变量值大小之间的差异,可鉴别样品真伪,提供定性信息的相似度。夹角余弦法是在多维空间向量夹角的基础上,计算样品指纹图谱特征向量与共有模式向量之间的夹角余弦相似度,可以提供样品真伪相似性的信息。王龙星等用向量夹角法,利用计算机对11个不同产地及炮制方法的吴茱萸样品的液相色谱指纹图谱分析,以往对中药“有效成分”的定量分析并不能完全表现出中药的内在质量。随着学科的发展,过去被认为是无效的成分如多糖、蛋白等,现在被发现其新的药用活性,所以对中药质量的评价,最好建立在缜密、系统、全面的基础上,结合药理活性等指标,做出综合评价。
(2)模糊聚类分析
中药的质量具有模糊性,图谱的相似性只有相对的意义。田兰等采用聚类分析和逐步判别分析对白术高效液相色谱进行化学模式识别研究,根据识别结果,初步建立了评价白术真伪优劣的新方法。
(3)人工神经网络法
反向传播模型(BP)是目前应用最广泛的人工神经网络模型,具有判别分类可靠性大,速度快等优点。汤丹等建立了广藿香GC-MS指纹图谱解析的人工神经网络模型,并通过BP网络对不同产地广藿香的识别,证明其有较好的识别功能。
(4)主成分分析法
郝燕等运用统计软件SPSS对10批次不同产地的当归指纹图谱的37个多维多息特征参数进行主成分分析。从对各主成分的命名中发现了4个主成分能够反映中药色谱指纹图谱信息的规律,证实了主成分分析能达到降维目的,使繁多的求解目标简化,可用于中药指纹图谱的数据挖掘。余杰等通过主成分分析和空间投影变换方法获取中药材质量的特征指纹图谱,对野生、栽培桑白皮及其伪品共42个样品进行了测定。结果表明,样品鉴别的准确率达到90.5%。
中药指纹图谱技术及模式识别技术是解决中药质评价和分类鉴别的有效技术手段。模式识别技术是近年来发展较快的新技术,广泛的应用于多种领域。有助于中药质量控制。运用模式识别根据物质所含的化学成分,用计算机对其进行分类或描述。可利用TLC,UV,IR,HPLC,GC,MC等方法获得其化学数据,再根据一定的方法进行识别。中药质量的模式识别,既体现了中药成分复杂,多成分多靶点综和作用的机理。又具有较强的科学性和实践性。正如每个人都有指纹而且指纹各不相同一样,每一味中药的特性和有效成分也千差万别。借助计算机和现代分析技术,可以将中药的特性和有效成分采用图谱的形式描绘出来,使每味中药都拥有如人的指纹一样的标准图谱,指纹图谱的计算机模式识别系统是实现将数字化的中药推向世界的最佳手段,可以避免人工比较所造成的工作繁重、准确性差等弊端。计算机模式识别系统可识别不同检测方法的指纹图谱以及各种基因指纹图谱,并能评价指纹图谱的相关性,以相关性判断提供试品的合格与否,方法简便、客观、实用。它能够成为中药质量评价的一种科学的、全面的、准确的方法。
目前传统中药在有效成分、药效、质控及作用机理等方面还缺乏健全的科学量化指标,在中药材化学成分的定性定量与药效,以及不同药材配伍的组效间的相互关系等方面还缺乏足够的理论深度,在安全、有效、均匀、稳定等质量方面还缺乏标准化的科学质量评价体系,难以取得国际市场的认同。有关传统中药定性定量的安全性和有效性数据还远远达不到可使之成为国际通用药物的标准。药材的道地性即药材地理区域性,是指许多药材质量特征和药效完全或主要取决于特殊地域的地理环境和自然条件。近年来,我国不断加大了对道地性药材的种植与研究开发的力度,发展了很多道地性药材基地。对于中药,其单味药的化学组成都已非常复杂并且大多为未知,少则几十种,多则成百上千种,传统的物理或化学工作量极大对于中药的信息获取可能也不充分,而且结果可能还会存在大的偏差,模式识别在中药上的运用则可能获取中药这一复杂化学体系的整体特征信息的差异,从不同的中药材或复方中药样本中寻找出决定样本差异的隐变量,从而达到对不同属性的中药材或复方中药样本进行判别的目的,可能更为简便且可靠。
发明内容
本发明的目的在于提供一种评价中药材质量的方法,以解决上述背景技术中提出的问题。
为实现上述目的,本发明提供如下技术方案:
一种评价中药材质量的方法,包括以下步骤:
一、中药指纹图谱的中药寒热特征标记
1)以单一指纹图谱技术的中药图谱数据为单元,釆用偏最小二乘判别模型识别中药寒热特征标记,进而挑选出不同指纹图谱技术下的中药寒热特征标记;
2)将挑选出的中药寒热特征标记组成中药寒热特征标记物质组群向量,以此向量分别在无机元素、初生物质和次生物质三个层面上再建立偏最小二乘判别模型,以识别各物质层面上的中药寒热特征标记;
3)综合三个物质层面上的所有中药寒热特征标记建立中药寒热整体药性的偏最小二乘模型,以显示中药整体寒热特征标记物质组群及其分布特征;
二、对中药指纹图谱进行分类
对中药指纹图谱的中药寒热特征标记采用Multi-Agent融合算法进行分类;
三、指纹图谱相似度计算方法
采用改进的灰色综合关联度进行指纹图谱的相似度计算;具体计算步骤如下:
设Xo={x0(1),x0(2),···,x0(n)}是数据行为参考序列,
X1={x1(1),x1(2),···,x1(n)}
······
Xi={xi(1),xi(2),···,xi(n)}
······
Xm={xm(1),xm(2),···,xm(n)}
是数据行为比较序列;
设X0,Xi的长度相同,且初值都不是零值,ε0i和r0i分别是灰色绝对关联度和灰色相对关联度,θ∈[0,1],那么称ρ0i=θε0i+(1-θ)r0i是X0,Xi的灰色综合关联度;
令
其中,
作为本发明进一步的方案:单一指纹图谱技术包括但不限于红外光谱、紫外光谱、薄层色谱与气相色谱。
作为本发明进一步的方案:每个智能体包含一个单分类器。
作为本发明进一步的方案:Multi-Agent融合算法的实现过程步骤如下:
Step 1:通过训练集,统计单分类器对训练样本分类所产生的结果信息,得到决策相关性矩阵D和可信度矩阵S的初始值;
Step 2:用v表示各Agent对样本x裁决置信度,定义为可信度矩阵的元素最大值;l表示各智能体Agent最终决策结果;
Step 3:若v大于一个阈值T,则表示各Agent基本达成共识,则进行Step 5;否则用可信度矩阵S调节溯源矩阵,转Step 4;
Step 4:将可信度矩阵S的每一行进行归一化,保证每行的和为1,重新计算v和l的值,转Step 3;
Step 5:l的值为各Agent的最终达成的一致意见,即为Multi-Agent融合算法的最终诊断结果。
与现有技术相比,本发明的有益效果是:
本发明是在利用紫外光谱(UV)、红外光谱(IR)、薄层色谱(TLC),气相色谱(GC),电化学等手段获取药材指纹图谱的基础上进行研究。首先,对本地民族中药材获取的指纹图谱进行预处理,获取相应的特征值,针对这些特征值(大量数据),为了克服多维海量数据对分类的影响其特点。拟采用模式识别中的多分类器融合技术,及指纹图谱总体相似度判定,然后与传统的专家经验相结合验证分类效果。从而为云南道地药材的质量评价检测提供相应依据。参照评价结果能够科学、准确地判断相关药材的真伪、道地性、质量的优劣。这对于民族药质量标准化、走向国际市场起到推动作用。
总之,本发明能够定性和定量评定图谱间的相似度,进而评价中药材质量,并比较其稳定性、一致性对实验结果进行鉴定、评价。是一种评价中药指纹图谱的有效手段。
附图说明
图1是本发明方法示意图。
具体实施方式
下面将结合本发明实施例,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明要解决的问题:
1、确定有效成分的提取工艺并进行优化,同时最适宜的色谱分离条件,获取民族中药的指纹图谱的基础上对指纹图谱的数据进行科学合理的筛选与比较,进而鉴定判别产品的质量,并比较其稳定性、一致性对实验结果进行鉴定、评价。
2、研究了多类分类问题中的分类器组合架构,为克服已有架构的不足,分析集成的整体性能与成员分类器性能和多样性之间关系,提高原系统的多样性及分类性。
3、定性和定量评定图谱间的相似度:目前为辨别中药材的真伪和质量优劣,主要参考的指标是中药指纹图谱的相似度。对于从整体上表现中药材复杂体系所含的化学成分的变化情况,指纹图谱的相似度有着关键的作用。相似度能够全面,客观地反映中药指纹图谱间的相似程度,进而评价中药材质量,是一种评价中药指纹图谱的有效手段。
为了解决上述问题,本发明的技术方案如下:
(1)中药指纹图谱的中药寒热特征标记的处理方法
中药指纹图谱形象地反映了物种具有遗传特性的次生代谢的“整体确定性(共有特征)”,而地域、生长环境、采收等多种不定因素对中药材的次生代谢存在重要影响,具有统计学中多元随机分布的“模糊性(非共有特征)”,中药多维多息指纹图谱数据具有明显的“高维灾难性”特征,即数据矩阵的行(即中药样本)往往很小,而数据矩阵的列(即图谱)往往很大。此类维数灾难性数据对中药药性特征标记的识别提出了挑战。由于维数灾难及其高度共线性的影响,传统的统计模式识别方法(例如Hsher判别,Logistic回归判别等)往往难以处理此类数据。
本发明釆用偏最小二乘判别模型,基于多种种图谱数据对中药特征标记进行综合研究,以从整体上显示中药特征标记物质组群。首先以单一指纹图谱技术(包括红外光谱、紫外光谱等)的中药指纹图谱数据为单元釆用偏最小二乘判别模型识别中药寒热特征标记,进而挑选出不同图谱技术下的特征标记。然后将挑选出的特征标记组成中药寒热特征标记物质组群向量,以此向量分别在无机元素、初生物质和次生物质三个层面上再建立偏最小二乘判别模型,以识别各物质层面上的中药寒热特征标记;最后综合三个物质层面上的所有特征标记建立中药寒热整体药性的偏最小二乘模型,以显示中药整体寒热特征标记物质组群及其分布特征。由于该方法是一种稳健的判别分析统计方法,它将主成分分析,典型相关分析和线性判别分析有机结合,既克服了线性判别分析难以处理高度共线性数据的问题,又克服了主成分分析的无监督缺陷,还通过典型相关,充分利用了协变量与分类标示变量之间的整体相关性。因此,它是处理高维灾难且高度共线性多维多息指纹图谱数据的有效方法。
(2)基于指纹图谱分类的方案
目前研究中药材分类算法的研究已经很多,主要有神经网络、决策树、k-最近邻、模糊聚类法及贝叶斯等,但每种方法都有一定的偏好和局限性。特定的方法往往在一个问题或数据集上取得好的效果,而在另一个问题或数据集上不尽如人意。本发明是针对单分类器存在的问题而提出的,具体方案为,采用多智能体(Multi-Agent)组融合技术,每个智能体包含一个单分类器。其多个智能体之间为共同完成一个任务,具有信息相互协商,相互服务的行为。通过这种信息的交互作用使每个智能体的行为得以更好的表现,同时采用Multi-Agent多分类器融合算法克服经典的MA多分类器融合算法在多个Agent(智能体)进行相互协商时,没有考虑各单分类器的分类效果、识别性能有所差别的情况,在决策的时候没有考虑到迭代次数不收敛,拒绝决策的情况。
Multi-Agent融合算法的主要实现过程步骤如下:
Step 1:通过训练集,统计单分类器对训练样本分类所产生的结果信息,得到决策相关性矩阵D和可信度矩阵S的初始值;
Step 2:用v表示各Agent对样本x裁决置信度,定义为可信度矩阵的元素最大值;l表示各智能体Agent最终决策结果。
Step 3:若v大于一个阈值T,则表示各Agent基本达成共识,则进行Step 5;否则用可信度矩阵S调节溯源矩阵,转Step 4。
Step 4:将可信度矩阵S的每一行进行归一化,保证每行的和为1,重新计算v和l的值,转Step 3。
Step 5:l的值为各Agent的最终达成的一致意见,即为Multi-Agent融合算法的最终诊断结果。
(3)指纹图谱相似度计算方法
相似度是复杂的中药材进行质量控制,判别真伪的一个重要指标,是中药指纹图谱技术的一个有效手段,符合中药指纹图谱的模糊性和整体性的属性。因此,相似度的计算方法非常重要。
目前,有多种相似度计算方法。可以分为两大类,定性相似度和定量相似度。定性相似度能够反映中药化学成分的分布比例的相似情况,但是不具有定量评价功能,即不能够反映化学成分含量之间的差异。定性相似度主要有,相关系数法,夹角余弦,Nei系数法,改进Nei系数法等。定量相似度能够从整体上反映化学成分含量之间的差异情况,主要有程度相似度,投影含量相似度,含量相似度等。
以往的相似度计算方法只有定性相似度和定量相似度,没有综合考虑定性和定量的相似度计算方法,针对这一缺陷,本发明将灰色综合关联度进行了改进,应用于相似度计算中。
设Xo={x0(1),x0(2),···,x0(n)}是数据行为参考序列,
X1={x1(1),x1(2),···,x1(n)}
······
Xi={xi(1),xi(2),···,xi(n)}
······
Xm={xm(1),xm(2),···,xm(n)}
是数据行为比较序列。
设X0,Xi的长度相同,且初值都不是零值,ε0i和r0i分别是的灰色绝对关联度和灰色相对关联度,θ∈[0,1],那么称
ρ0i=θε0i+(1-θ)r0i
是X0,Xi的灰色综合关联度。
综合关联度既反映了序列折线之间的相似程度,又体现了折线相对于始点的变化速率的接近程度,是比较全面的反映序列之间相似程度的一个数量指标。
在灰色综合关联度公式中,θ的一般取值为定值。为了使灰色综合关联度更加准确的反映序列之间的相似性,对θ的取值进行了改进,如下所示:
其中,
以往的相似度计算方法中,没有综合考虑定性与定量功能的相似度,针对这一缺陷,故将改进的灰色综合相似度应用于相似度的计算中。将对照指纹图谱看作数据行为参考序列,将作比较的样品指纹图谱看作数据行为比较列,ε0i反映了定量功能,r0i反映了定性功能,综合考虑了定性与定量。
中药指纹图谱可采用以下仪器及条件进行获得。
(1)仪器、试剂与药材
LC-20AD高效液相色谱仪;岛津Solution色谱工作站;紫外检测器;KH5200型超声波清洗器;电子分析天平;旋转蒸发仪;过滤器。甲醇(分析纯);甲醇(色谱纯);乙腈(色谱纯);水为纯净水。实验采用云南省药物研究所提供的本地民族中药材为样品并提供鉴定。
(2)色谱条件
采用冷浸超声法提取云南蜂胶样品,液相色谱条件为:色谱柱:AgilentXDB C18柱(4.6×150mm,5μm);流动相:乙腈、甲醇和0.001%磷酸水溶液,流速为0.8mL/min,流动相按一定比例进行梯度洗脱;柱温:30℃;扫描波长:320nm。
(3)方法考察
以样品的主要色谱峰的相对保留时间及相对保留峰面积来考察分析方法的稳定性,仪器精密度及实验方法的重现性。
①精密度试验
取样品的同一份供试品溶液连续进样6次,记录指纹图谱,结果表明各主要色谱峰相对保留时间和相对峰面积比值无明显变化。
②稳定性试验
在0,4,8,12,36,48h进样20L,记录指纹图谱,结果表明各主要色谱峰相对保留时间和相对峰面积比值无明显变化。
③重复性试验
取同一批样品按上述方法制备6份供试液,按上述色谱条件测定指纹图谱,结果表明各主要色谱峰相对保留时间和相对峰面积比值无明显变化,表明重复性良好,符合指纹图谱要求。
通过指纹图谱,经过数据处理,拟采用其出峰时间、峰高、峰面积等数据作为特征向量,对其进行辨别。
本发明的特点:
(1)由于中药材在不同的处方中显示不同的功效,有效部位也不是一个,一张色谱指纹图谱很难代表中药材的整体质量,中药材的“标准指纹图谱”在各实验室中量值的传递也会因仪器类型的不同而存在一定的问题。因此,色谱指纹图谱的相似性比较,并不完全适用于中药材质量标准中。本发明提出用灰色综合关联度进行了改进,应用于相似度计算中,克服以往相似度计算方法只有定性相似度和定量相似度,没有综合考虑定性和定量的缺点。为云南民族药材的鉴别提供有效的技术支持,为其他药材的分析提供参考。同时模式识别分析法对不同产地的民族药材进行较好的聚类分析及质量评价。
(2)采用Multi-Agent融合算法进行分类,继承神经网络、决策树、模糊聚类法及贝叶斯等分类器的优点,克服它们对一定数据和特征的偏好和局限性。本发明构建一种基于相关成分分析的多分类器,并将其用于中药材质量模式分类。通过从色谱分析所得到的高维数据集中分步提取分类相关成分,获取化学模式特征向量。
对于本领域技术人员而言,显然本发明不限于上述示范性实施例的细节,而且在不背离本发明的精神或基本特征的情况下,能够以其他的具体形式实现本发明。因此,无论从哪一点来看,均应将实施例看作是示范性的,而且是非限制性的,本发明的范围由所附权利要求而不是上述说明限定,因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。
此外,应当理解,虽然本说明书按照实施方式加以描述,但并非每个实施方式仅包含一个独立的技术方案,说明书的这种叙述方式仅仅是为清楚起见,本领域技术人员应当将说明书作为一个整体,各实施例中的技术方案也可以经适当组合,形成本领域技术人员可以理解的其他实施方式。
Claims (4)
1.一种评价中药材质量的方法,其特征在于,包括以下步骤:
一、中药指纹图谱的中药寒热特征标记
1)以单一指纹图谱技术的中药图谱数据为单元,釆用偏最小二乘判别模型识别中药寒热特征标记,进而挑选出不同指纹图谱技术下的中药寒热特征标记;
2)将挑选出的中药寒热特征标记组成中药寒热特征标记物质组群向量,以此向量分别在无机元素、初生物质和次生物质三个层面上再建立偏最小二乘判别模型,以识别各物质层面上的中药寒热特征标记;
3)综合三个物质层面上的所有中药寒热特征标记建立中药寒热整体药性的偏最小二乘模型,以显示中药整体寒热特征标记物质组群及其分布特征;
二、对中药指纹图谱进行分类
对中药指纹图谱的中药寒热特征标记采用Multi-Agent融合算法进行分类;
三、指纹图谱相似度计算方法
采用改进的灰色综合关联度进行指纹图谱的相似度计算;具体计算步骤如下:
设Xo={x0(1),x0(2),···,x0(n)}是数据行为参考序列,
X1={x1(1),x1(2),···,x1(n)}
······
Xi={xi(1),xi(2),···,xi(n)}
······
Xm={xm(1),xm(2),···,xm(n)}
是数据行为比较序列;
设X0,Xi的长度相同,且初值都不是零值,ε0i和r0i分别是灰色绝对关联度和灰色相对关联度,θ∈[0,1],那么称ρ0i=θε0i+(1-θ)r0i是X0,Xi的灰色综合关联度;
令
其中,
2.根据权利要求1所述的评价中药材质量的方法,其特征在于,单一指纹图谱技术包括但不限于红外光谱、紫外光谱、薄层色谱与气相色谱。
3.根据权利要求1所述的评价中药材质量的方法,其特征在于,每个智能体包含一个单分类器。
4.根据权利要求1所述的评价中药材质量的方法,其特征在于,Multi-Agent融合算法的实现过程步骤如下:
Step 1:通过训练集,统计单分类器对训练样本分类所产生的结果信息,得到决策相关性矩阵D和可信度矩阵S的初始值;
Step 2:用v表示各Agent对样本x裁决置信度,定义为可信度矩阵的元素最大值;l表示各智能体Agent最终决策结果;
Step 3:若v大于一个阈值T,则表示各Agent基本达成共识,则进行Step5;否则用可信度矩阵S调节溯源矩阵,转Step 4;
Step 4:将可信度矩阵S的每一行进行归一化,保证每行的和为1,重新计算v和l的值,转Step 3;
Step 5:l的值为各Agent的最终达成的一致意见,即为Multi-Agent融合算法的最终诊断结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710316615.9A CN106990214B (zh) | 2017-05-08 | 2017-05-08 | 一种评价中药材质量的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710316615.9A CN106990214B (zh) | 2017-05-08 | 2017-05-08 | 一种评价中药材质量的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106990214A true CN106990214A (zh) | 2017-07-28 |
CN106990214B CN106990214B (zh) | 2019-11-15 |
Family
ID=59417584
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710316615.9A Expired - Fee Related CN106990214B (zh) | 2017-05-08 | 2017-05-08 | 一种评价中药材质量的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106990214B (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107449849A (zh) * | 2017-10-11 | 2017-12-08 | 天津中医药大学 | 一种中药鉴别方法 |
CN108226347A (zh) * | 2018-02-09 | 2018-06-29 | 中南民族大学 | 七叶皂苷钠化学成分组分配比及其筛选方法 |
CN109376957A (zh) * | 2018-12-03 | 2019-02-22 | 渤海大学 | 一种火电厂负荷的预测方法 |
CN109668850A (zh) * | 2019-02-28 | 2019-04-23 | 山东中医药大学 | 基于紫外指纹图谱的中药药性识别方法及系统 |
CN110031595A (zh) * | 2019-03-06 | 2019-07-19 | 湖南中医药大学 | 一种动态性与传递性综合控制的中药物质基准制备方法及其产品 |
CN110376142A (zh) * | 2019-08-20 | 2019-10-25 | 陕西中医药大学 | 中药红花质量等级的检测方法 |
CN106990214B (zh) * | 2017-05-08 | 2019-11-15 | 云南民族大学 | 一种评价中药材质量的方法 |
CN110838343A (zh) * | 2019-11-15 | 2020-02-25 | 山东中医药大学 | 一种基于多模态指纹图谱的中药药性识别方法及系统 |
CN112161965A (zh) * | 2020-09-22 | 2021-01-01 | 复旦大学附属妇产科医院 | 检测中药药性的方法、装置、计算机设备和存储介质 |
CN113567587A (zh) * | 2021-07-29 | 2021-10-29 | 山东中医药大学 | 基于挥发油成分特征的中药寒热药性识别方法 |
CN113567639A (zh) * | 2021-07-13 | 2021-10-29 | 中国食品药品检定研究院 | 中药材质量综合评价方法 |
CN115407739A (zh) * | 2022-10-31 | 2022-11-29 | 天津有容蒂康通讯技术有限公司 | 一种电缆制造的生产设备控制方法及系统 |
CN117330678B (zh) * | 2023-12-01 | 2024-02-09 | 内蒙古亿利制药有限公司 | 蓉蛾益肾口服液的质量检测方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101271092A (zh) * | 2007-03-19 | 2008-09-24 | 沈阳药科大学 | 中药色谱指纹图谱特征数字化和全定性全定量质量控制方法 |
CN101532954A (zh) * | 2008-03-13 | 2009-09-16 | 天津天士力现代中药资源有限公司 | 一种用红外光谱结合聚类分析鉴定中药材的方法 |
CN101676717A (zh) * | 2008-09-19 | 2010-03-24 | 天津天士力制药股份有限公司 | 一种中药制品的质量评价方法 |
CN103411912A (zh) * | 2013-05-22 | 2013-11-27 | 首都师范大学 | 一种利用THz-TDS结合模糊规则专家系统鉴定中草药的方法 |
CN104345045A (zh) * | 2014-11-04 | 2015-02-11 | 天津工业大学 | 一种基于化学模式识别和近红外光谱的相似药材鉴别方法 |
CN104372075A (zh) * | 2014-08-22 | 2015-02-25 | 贵州茅台酒股份有限公司 | 一种构建判别模型鉴别大曲质量的方法 |
CN105138861A (zh) * | 2015-05-31 | 2015-12-09 | 青岛市食品药品检验研究院 | 一种大黄药材真伪及基源种属类别预测模型的构建方法 |
CN105608460A (zh) * | 2014-11-24 | 2016-05-25 | 中国电信股份有限公司 | 多分类器融合方法和系统 |
CN106560695A (zh) * | 2016-10-20 | 2017-04-12 | 中国计量大学 | 联合三种检测方法的武夷岩茶产地鉴别方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106990214B (zh) * | 2017-05-08 | 2019-11-15 | 云南民族大学 | 一种评价中药材质量的方法 |
-
2017
- 2017-05-08 CN CN201710316615.9A patent/CN106990214B/zh not_active Expired - Fee Related
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101271092A (zh) * | 2007-03-19 | 2008-09-24 | 沈阳药科大学 | 中药色谱指纹图谱特征数字化和全定性全定量质量控制方法 |
CN101532954A (zh) * | 2008-03-13 | 2009-09-16 | 天津天士力现代中药资源有限公司 | 一种用红外光谱结合聚类分析鉴定中药材的方法 |
CN101676717A (zh) * | 2008-09-19 | 2010-03-24 | 天津天士力制药股份有限公司 | 一种中药制品的质量评价方法 |
CN103411912A (zh) * | 2013-05-22 | 2013-11-27 | 首都师范大学 | 一种利用THz-TDS结合模糊规则专家系统鉴定中草药的方法 |
CN104372075A (zh) * | 2014-08-22 | 2015-02-25 | 贵州茅台酒股份有限公司 | 一种构建判别模型鉴别大曲质量的方法 |
CN104345045A (zh) * | 2014-11-04 | 2015-02-11 | 天津工业大学 | 一种基于化学模式识别和近红外光谱的相似药材鉴别方法 |
CN105608460A (zh) * | 2014-11-24 | 2016-05-25 | 中国电信股份有限公司 | 多分类器融合方法和系统 |
CN105138861A (zh) * | 2015-05-31 | 2015-12-09 | 青岛市食品药品检验研究院 | 一种大黄药材真伪及基源种属类别预测模型的构建方法 |
CN106560695A (zh) * | 2016-10-20 | 2017-04-12 | 中国计量大学 | 联合三种检测方法的武夷岩茶产地鉴别方法 |
Non-Patent Citations (2)
Title |
---|
张新新: "基于多维多息指纹图谱的中药药性特征标记模式识别研究", 《中国优秀硕士学位论文全文数据库医药卫生科技辑》 * |
耿洁: "中药指纹图谱的相似度评价方法", 《中国优秀硕士学位论文全文数据库医药卫生科技辑》 * |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106990214B (zh) * | 2017-05-08 | 2019-11-15 | 云南民族大学 | 一种评价中药材质量的方法 |
CN107449849B (zh) * | 2017-10-11 | 2020-07-03 | 天津中医药大学 | 一种中药鉴别方法 |
CN107449849A (zh) * | 2017-10-11 | 2017-12-08 | 天津中医药大学 | 一种中药鉴别方法 |
CN108226347A (zh) * | 2018-02-09 | 2018-06-29 | 中南民族大学 | 七叶皂苷钠化学成分组分配比及其筛选方法 |
CN109376957A (zh) * | 2018-12-03 | 2019-02-22 | 渤海大学 | 一种火电厂负荷的预测方法 |
CN109376957B (zh) * | 2018-12-03 | 2020-11-20 | 渤海大学 | 一种火电厂负荷的预测方法 |
CN109668850A (zh) * | 2019-02-28 | 2019-04-23 | 山东中医药大学 | 基于紫外指纹图谱的中药药性识别方法及系统 |
CN110031595A (zh) * | 2019-03-06 | 2019-07-19 | 湖南中医药大学 | 一种动态性与传递性综合控制的中药物质基准制备方法及其产品 |
CN110376142A (zh) * | 2019-08-20 | 2019-10-25 | 陕西中医药大学 | 中药红花质量等级的检测方法 |
CN110376142B (zh) * | 2019-08-20 | 2022-07-01 | 陕西中医药大学 | 中药红花质量等级的检测方法 |
CN110838343A (zh) * | 2019-11-15 | 2020-02-25 | 山东中医药大学 | 一种基于多模态指纹图谱的中药药性识别方法及系统 |
CN110838343B (zh) * | 2019-11-15 | 2022-03-01 | 山东中医药大学 | 一种基于多模态指纹图谱的中药药性识别方法及系统 |
CN112161965A (zh) * | 2020-09-22 | 2021-01-01 | 复旦大学附属妇产科医院 | 检测中药药性的方法、装置、计算机设备和存储介质 |
CN112161965B (zh) * | 2020-09-22 | 2022-08-05 | 复旦大学附属妇产科医院 | 检测中药药性的方法、装置、计算机设备和存储介质 |
CN113567639A (zh) * | 2021-07-13 | 2021-10-29 | 中国食品药品检定研究院 | 中药材质量综合评价方法 |
CN113567587A (zh) * | 2021-07-29 | 2021-10-29 | 山东中医药大学 | 基于挥发油成分特征的中药寒热药性识别方法 |
CN115407739A (zh) * | 2022-10-31 | 2022-11-29 | 天津有容蒂康通讯技术有限公司 | 一种电缆制造的生产设备控制方法及系统 |
CN117330678B (zh) * | 2023-12-01 | 2024-02-09 | 内蒙古亿利制药有限公司 | 蓉蛾益肾口服液的质量检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106990214B (zh) | 2019-11-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106990214B (zh) | 一种评价中药材质量的方法 | |
Chen et al. | Quality control and original discrimination of Ganoderma lucidum based on high-performance liquid chromatographic fingerprints and combined chemometrics methods | |
Goodarzi et al. | Similarity analyses of chromatographic herbal fingerprints: a review | |
Höft et al. | Quantitative ethnobotany: Application of multivariate and statistical analysis in ethnobotany | |
Clarke | Extracting process-related information from pharmaceutical dosage forms using near infrared microscopy | |
Cao et al. | Classification and identification of Rhodobryum roseum Limpr. and its adulterants based on fourier-transform infrared spectroscopy (FTIR) and chemometrics | |
Peng et al. | Fingerprint profile of active components for Artemisia selengensis Turcz by HPLC–PAD combined with chemometrics | |
CN101676717B (zh) | 一种中药制品的检测方法 | |
CN102087212B (zh) | 基于主成份分析的葛粉掺假鉴别方法 | |
CN101861157A (zh) | 用于评价传统药物的系统和方法 | |
CN105181916A (zh) | 一种利用电子鼻、电子舌传感器技术检测量化中药气味及其对应证的方法 | |
Dong et al. | Deep learning for geographical discrimination of Panax notoginseng with directly near-infrared spectra image | |
CN114595956A (zh) | 一种基于灰度关联法模糊聚类算法的桉树土壤肥力分析方法 | |
Cheng et al. | Infrared spectroscopic studies of Chinese medicines | |
Jiang et al. | Near-infrared spectroscopy combined with machine learning for rapid identification of Atractylodis rhizoma decoction pieces | |
Dahl et al. | Otolith morphometry and Fourier transform near-infrared (FT-NIR) spectroscopy as tools to discriminate archived otoliths of newly detected cryptic species, Etelis carbunculus and Etelis boweni | |
CN116665848A (zh) | 一种中药配比设备及其中药配比方法 | |
CN109870533A (zh) | 一种基于分子地图的生物样本快速智能识别方法 | |
AU2018101606A4 (en) | A method for identifying meconopsis quintuplinervia regel from different geographical origins | |
CN111638294A (zh) | 一种小叶金花草的鉴别方法及其应用 | |
Hu et al. | A review of the application of vibrational spectroscopy combined with chemometrics in genus Panax | |
Liu et al. | Quality Consistency Evaluation of Traditional Chinese Medicines: Current Status and Future Perspectives | |
San Nicolas et al. | NIR-hyperspectral imaging and machine learning for non-invasive chemotype classification in Cannabis sativa L | |
Lu et al. | Rapid analysis of Radix Astragali using a portable Raman spectrometer with 1064-nm laser excitation and data fusion with PLS-DA | |
Chen et al. | Pipeline for high-throughput modeling of marijuana and hemp extracts |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20191115 Termination date: 20210508 |