CN103278591A - 一种色谱指纹图谱相似度评价方法 - Google Patents
一种色谱指纹图谱相似度评价方法 Download PDFInfo
- Publication number
- CN103278591A CN103278591A CN2013101831648A CN201310183164A CN103278591A CN 103278591 A CN103278591 A CN 103278591A CN 2013101831648 A CN2013101831648 A CN 2013101831648A CN 201310183164 A CN201310183164 A CN 201310183164A CN 103278591 A CN103278591 A CN 103278591A
- Authority
- CN
- China
- Prior art keywords
- peak
- sample
- similarity
- finger
- 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
- 238000011156 evaluation Methods 0.000 title abstract description 13
- 239000013598 vector Substances 0.000 claims abstract description 78
- 238000000034 method Methods 0.000 claims abstract description 39
- 239000003814 drug Substances 0.000 claims abstract description 18
- 235000013305 food Nutrition 0.000 claims abstract description 11
- 230000000694 effects Effects 0.000 claims abstract description 8
- 230000005477 standard model Effects 0.000 claims description 35
- 239000000203 mixture Substances 0.000 claims description 30
- 230000008569 process Effects 0.000 claims description 5
- 230000008878 coupling Effects 0.000 claims description 4
- 238000010168 coupling process Methods 0.000 claims description 4
- 238000005859 coupling reaction Methods 0.000 claims description 4
- 238000005251 capillar electrophoresis Methods 0.000 claims description 2
- 238000004817 gas chromatography Methods 0.000 claims description 2
- 238000004128 high performance liquid chromatography Methods 0.000 claims description 2
- 229940079593 drug Drugs 0.000 abstract description 3
- 239000000523 sample Substances 0.000 description 100
- 238000011160 research Methods 0.000 description 11
- 238000005516 engineering process Methods 0.000 description 9
- 239000013068 control sample Substances 0.000 description 7
- 230000006872 improvement Effects 0.000 description 7
- 239000008862 fructus schizandrae, radix ginseng, radix ophiopogonis drug combination Substances 0.000 description 6
- 238000002347 injection Methods 0.000 description 6
- 239000007924 injection Substances 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 239000000126 substance Substances 0.000 description 5
- 230000002159 abnormal effect Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 230000000857 drug effect Effects 0.000 description 3
- 238000003908 quality control method Methods 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 241000894007 species Species 0.000 description 3
- 230000000295 complement effect Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 235000002789 Panax ginseng Nutrition 0.000 description 1
- 240000005373 Panax quinquefolius Species 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 239000000599 controlled substance Substances 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000010790 dilution Methods 0.000 description 1
- 239000012895 dilution Substances 0.000 description 1
- 238000004836 empirical method Methods 0.000 description 1
- 239000013604 expression vector Substances 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 239000004615 ingredient Substances 0.000 description 1
- 239000002398 materia medica Substances 0.000 description 1
- 238000000491 multivariate analysis Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000001988 toxicity Effects 0.000 description 1
- 231100000419 toxicity Toxicity 0.000 description 1
- 229940126680 traditional chinese medicines Drugs 0.000 description 1
Images
Landscapes
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
一种色谱指纹图谱相似度评价方法,包括以下步骤:采用标准样品建立对照指纹图谱;采用与建立对照指纹图谱相同的色谱条件建立待评价样品的样品指纹图谱;将样品指纹图谱的每一个峰值与对照指纹图谱中对应峰值相比,得到样品的比值指纹图谱向量和对照比值指纹图谱向量;计算这两个比值向量的相似度S。本发明既能反映待评价样品指纹图谱与标准指纹图谱之间共有峰的相似性,也考虑了非共有峰数目及其大小对产品质量的影响,能整体反应各样品指纹图谱之间相似程度,该相似度的大小能够灵敏地定量表达样品指纹图谱与标准指纹图谱定性和定量差异,可靠地鉴别药品、食品的真伪,更好地从质量稳定性和安全性上控制药品、食品的质量。
Description
技术领域
本发明属于药品、食品领域,确切地说属于药品、食品色谱指纹图谱质量控制技术领域。
背景技术
指纹图谱技术是目前公认的全面控制中药、食品等复杂成分体系质量的最有效最可靠的技术。
指纹图谱质控技术进行质量控制的关键是相似度评价理论及其软件化。目前应用最广泛的评价色谱指纹图谱相似度方法主要有Tanimoto、峰重叠率及改进方法、向量夹角余弦、相关系数、组合相似度、基于相似系统理论的多种程度相似度、比率定性相似度、总量统计矩相似度、距离相似度和加权马氏距离相似度、以及用相对熵的歧异值来衡量色谱指纹图谱的相似度等方法,专利CN200810087195.2也提出了37个特性指标用于中药质量的控制。色谱指纹图谱相似度评价的实质是考察要比较的两个指纹图谱的组成峰数、峰匹配以及相应峰量化关系的一致性:即一方面从产品质量的稳定性和一致性考虑,相似度大小应能灵敏反映多元化学成分分布比例的差异和整体含量的差异;另一方面从安全性考虑,相似度大小应能考察非共有峰数目和强度大小对产品安全性的影响。
现有的指纹图谱相似度算法难于全面满足上述对药品、食品质量稳定性和安全性的要求。
发明内容
本发明的目的是针对现有的指纹图谱相似度算法难于全面满足对药品、食品质量稳定性和安全性要求的问题,提出一种色谱指纹图谱相似度评价方法。
本发明的技术方案是:
一种色谱指纹图谱相似度评价方法,它包括以下步骤:
(2)、采用与建立对照指纹图谱相同的色谱条件建立待评价样品的样品指纹图谱,样品指纹图谱的向量为:
其中,1≤i≤n,n为样品指纹图谱中与对照指纹图谱共有峰的色谱峰数,xi代表待评价样品中i成分即样品指纹图谱中与对照指纹图谱共有峰的峰面积或峰高值,其中,对照指纹图谱中存在而在样品指纹图谱中不存在的成分,样品指纹图谱中该共有峰的取值为xi=0;
1≤j≤m,m为样品指纹图谱中与对照指纹图谱非共有峰的色谱峰数,xaj代表待评价样品中j成分即样品指纹图谱中与对照指纹图谱非共有峰的峰面积或峰高值,;
(3)、将样品指纹图谱的每一个峰值与对照指纹图谱中对应峰值相比,样品的比值指纹图谱表示为n+1维向量: 对照比值指纹图谱表示为n+1维向量:采用下述公式计算这两个比值向量的相似度S:
其中:
δi、δaj分别是各共有峰、非共有峰总面积的权重系数,权重的取值范围δi≥0、δa≥0;
表示样品指纹图谱中样品指纹图谱与对照指纹图谱非共有峰的总面积或峰总高值,作为求非共有峰比值的标准,表示对照指纹图谱峰总面积或峰总高值,nt表示实际的共有峰峰数,nt≤n,当对照指纹图谱存在的成分在样品指纹图谱中不存在时,样品指纹图谱中该共有峰的取值为xi=0;
本发明的步骤(1)中,采用标准样品建立对照指纹图谱包括以下步骤:
(a)、获取标准样品色谱指纹图谱的步骤:取标准样品,采用任一色谱方法按照色谱指纹图谱建立的规范要求建立能表征产品化学组成特征的中药或食品特征的标准样品色谱指纹图谱;
(b)、标准样品色谱指纹图谱预处理的步骤:取若干批次标准样品,采用前述方法获取各批次标准样品的色谱指纹图谱,对数据进行预处理和色谱峰匹配,以平均矢量法或中位数矢量法建立该产品的对照指纹图谱,所述对照指纹图谱的向量为
本发明的色谱方法包括高效液相色谱、气相色谱或高效毛细管电泳色谱。
本发明的步骤(b)中,标准样品的批次总数为k,对于对照指纹图谱的向量按照步骤(3)计算k个批次标准样品的比值指纹图谱向量Yk与此次确定的对照比值指纹图谱向量的相似度采用箱线图法对数据按从小到大排序,并求出该序列数据的统计量第一四分位数Q1、第三四分位数Q3,并求出这两个四分位数差值df=Q3-Q1,将所有数据与Q1-1.5df比较,所有小于Q1-1.5df的数据判断为异常值而将对应批次的标准样品去除,重复前述步骤,直至剩余的不再出现异常值,对剩余批次标准样品指纹图谱进行预处理和色谱峰匹配,以平均矢量法或中位数矢量法建立对照指纹图谱,所述对照指纹图谱的向量为
本发明的还包括合格样品的确定步骤:将最终剩余批次标准样品指纹图谱与最终确定的对照指纹图谱按照权利要求1的步骤(3)计算相似度S,其中S的最小值即为判断阈值S限;待评价样品的指纹图谱与最终确定的对照指纹图谱的相似度按照权利要求1的步骤(3)计算相似度S大于S限即为合格样品。
本发明的步骤(3)中,δi=1。
本发明的步骤(3)中,δi的取值原则是:相应成分对产品有效性安全性影响非常明显时δi>1,该成分对产品有效性安全性影响较小时1>δi≥0。
本发明的有益效果:
本发明是一种综合评价样品指纹图谱与标准指纹图谱相似性的方法,既能反映待评价样品指纹图谱与标准指纹图谱之间共有峰的相似性,也考虑了非共有峰数目及其大小对产品质量的影响,能整体反应各样品指纹图谱之间相似程度,该相似度的大小能够灵敏地定量表达样品指纹图谱与标准指纹图谱定性和定量差异,可靠地鉴别药品、食品的真伪,更好地从质量稳定性和安全性上控制药品、食品的质量。
本发明中,通过向量X与Y间的相似度S大小表征待评价样品与标准样品间质量的稳定性和一致性;相似度越大表示化学成分组成的一致性越好;相似度越小,表示化学成分组成的差异越大,安全性可能得不到保证。应用时,可以通过经验或有效性安全性等试验确定相似度阈值S限,凡待评价样品相似度S大于阈值S限的表示质量稳定,安全性得以保证;根据S的大小也可用于产品真伪的鉴别。
附图说明
图1是本发明实施例中一个样品和对照样品指纹图谱向量图。
图2是本发明的实施例中不同样品和对照样品指纹图谱向量及投影图。
具体实施方式
下面结合附图和实施例对本发明作进一步的说明。
一种色谱指纹图谱相似度评价方法,它包括以下步骤:
(2)、采用与建立对照指纹图谱相同的色谱条件建立待评价样品的样品指纹图谱,样品指纹图谱的向量为:
其中,1≤i≤n,n为样品指纹图谱中与对照指纹图谱共有峰的色谱峰数,xi代表待评价样品中i成分即样品指纹图谱中与对照指纹图谱共有峰的峰面积或峰高值,其中,对照指纹图谱中存在而在样品指纹图谱中不存在的成分,样品指纹图谱中该共有峰的取值为xi=0;
1≤j≤m,m为样品指纹图谱中与对照指纹图谱非共有峰的色谱峰数,xaj代表待评价样品中j成分即样品指纹图谱中与对照指纹图谱非共有峰的峰面积或峰高值,;
其中:
δi、δaj分别是各共有峰、非共有峰总面积的权重系数,权重的取值范围δi≥0、δa≥0;
表示样品指纹图谱中样品指纹图谱与对照指纹图谱非共有峰的总面积或峰总高值,作为求非共有峰比值的标准,表示对照指纹图谱峰总面积或峰总高值,nt表示实际的共有峰峰数,nt≤n,当对照指纹图谱存在的成分在样品指纹图谱中不存在时,样品指纹图谱中该共有峰的取值为xi=0。
相似度算法的提出如下:
本研究构建的相似度算法是基于事先确定的对照指纹图谱为计算参照标准;以中药色谱指纹图谱峰面积值或峰高为相似度评价信息参数;以图谱矢量化空间模式向量间相似性计算为基本出发点。
为说明相似度新算法构建的原理,现假设对照指纹图谱向量为 yi代表标准样品中i成分的峰面积或峰高值,n为色谱峰数;样品指纹图谱向量为xi代表样品中i成分的峰面积或峰高值,n为色谱峰数,如图1所示。用向量表示向量与向量的差向量:Z=(x1-y1,x2-y2,x3-y3,…,xi-yi,…xn-yn),则向量的模长以|Z|表示,其值表示向量与在多维空间的欧几里得距离:向量的模长用|Y|表示:显然|Z|与|Y|的比值|Z|/|Y|的大小在一定程度上反应向量与向量的相似程度,降低了变量单位不同以及检测器响应不同的对相似性大小的影响;由于其取值范围在0~∞之间,数值为零时表示两向量完全相同,数值越大两向量差异越大,考虑相似性定义的习惯,1表示两向量完全相同,初步定义相似度S计算公式(1):
S取值范围在-∞~1,数值越小两向量差异越大,这里相似度为负值被赋予实际意义,表示向量差异非常显著。
相似度算法的优化
当样品指纹图谱色谱峰比对照指纹图谱对应峰都偏大时,样品与对照样品的化学成分分布比例较一致,而如果样品指纹图谱色谱峰比对照指纹图谱对应峰一部分偏大另一部分偏小时,二者的化学成分分布比例相差较大,采用式(1)计算上述两种情况的相似度时不能灵敏反应这种差异,必须对公式(1)进行修正,修正后的计算方法见公式(2)。
式中AD表示与向量差值和的绝对值:MD表示与向量差绝对值的和:显然AD的取值范围为0~MD,AD与MD的比值AD/MD取值范围为0~1,当样品指纹图谱色谱峰比对照指纹图谱对应峰都偏大时AD/MD比值为1;如果样品指纹图谱色谱峰比对照指纹图谱对应峰一部分偏大另一部分偏小时AD/MD比值小于1,特殊情况下AD/MD比值可能为0,因此用项对(1)式修正后能够较准确反映指纹图谱峰值比例的变化。
然而用式(2)实际计算相似度时发现对样品指纹图谱中大峰与小峰有相同的绝对差异仍不能区分,如对2±2与10±2两种情形具有相同的评价结果,二者若从变量本身基数大小来考虑时,前者为基数基础上100%的变化,而后者为20%的变化,二者并不具有完全的等同性。为解决该问题可以这样考虑:既然计算相似度时是以对照指纹图谱为标准,那么可以将样品指纹图谱的每一个峰值与对照指纹图谱中对应峰值相比,可得到称之为样品的比值指纹图谱:显然对照比值指纹图谱为:y′=(1,1,1,...,1,...,1),再用公式(2)计算样品比值指纹图谱x′与对照比值指纹图谱y′的相似度,上述问题得以解决。
上述这种数据处理方法对共有峰模式指纹图谱是可行的,若样品指纹图谱中存在非共有峰(不能与对照指纹图谱相匹配的峰)时就出现问题:如何确定非共有峰比值。为表述方便,用n表示对照指纹图谱色谱峰总个数,m表示样品指纹图谱非共有峰的个数,xi代表样品中共有峰i成分的峰面积,xaj表示样品指纹图谱中非共有峰j的峰面积,表示样品指纹图谱的非共有峰总面积,参照中药注射剂指纹图谱研究的技术要求对指纹图谱共有峰控制每个峰比值的波动范围,而对非共有峰控制总面积的思路,选择对照指纹图谱总面积作为求非共有峰比值的标准,将所有非共有峰合并为一个峰组,其值就是非共有峰总面积那么样品的比值指纹图谱可表示为n+1维向量: 对照比值指纹图谱可表示为n+1维向量:y′=(1,1,1,...,1,...,1,0),以公式(2)来计算这两个比值向量的相似度,具体计算方法见公式(3):
公式(3)理论取值范围为-∞~1,相似度为-∞表示待比较样品指纹图谱与对照指纹图谱没有共有峰,也就是说样品指纹图谱向量与对照指纹图谱向量是正交的,显然特殊情况下当远小于时,用公式(3)计算它们间的相似度几乎为0,这显然是不合理的;公式(3)中的n为对照指纹图谱色谱峰峰数,是假定样品指纹图谱与对照指纹图谱的共有峰峰数为n,实际上当样品指纹图谱中缺少与对照指纹图谱相匹配的峰时,图谱向量化时将样品指纹图谱向量相应元素以0代替,若以nt表示共有峰数,显然nt≤n;若以实际的共有峰峰数nt代替公式(3)中n,就不存在上述不合理的情况了,当没有共有峰时即nt=0,S为-∞。相似度计算公式可由公式(4)表示。
由于中药组成成分复杂,各成分的结构及药效等不尽相同,故对中药指纹图谱相似性评价时,各成分可能存在有不同的比重系数。为此,实际应用时可考虑根据药效、毒性研究等确定中药不同成分相似度评价的药效影响因子,对指纹图谱相似度评价中的特征变量进行加权运算。加权运算公式为公式(5):
式中δi、δa分别是各共有峰、非共有峰总面积的权重系数,权重的取值范围δi≥0或δa≥0,δi、δa取值原则是:①向量的各维变量对产品有效性安全性影响不明确时δi及δa可全取1;②某些变量对产品有效性安全性影响非常明显时δi>1或δa>1,具体大小通过实验确定;③某些变量对产品有效性安全性影响较小时1>δi≥0或1>δa≥0,具体大小通过实验确定。确定δi、δa大小的具体试验方法有待进一步研究。
定性比较新算法与其它相似度算法:
假设有某一产品的对照指纹图谱向量和三个待比较样品指纹图谱向量如图2所示,是三个样品指纹图谱向量与的差向量,是向量在向量的投影,是在向量上的投影,|Z1|、|Z2|、|Z3|分别表示向量与之间的欧几里得距离,|Y|表示向量的模长,|C|、|D|表示向量的模长。从图1可以看出向量与向量的夹角余弦值相同,而|Z1|值大于|Z2|值,显然从本研究定义的相似度公式可以发现与的相似度大于与与的相似度,这表明本研究定义的相似度比夹角余弦法能更灵敏地反应向量之间的差异;图中也可以看出向量明显不同,而它们在向量上的投影向量的模长|C|(包括向量在向量上的投影向量模长|D|)是相等的,表明文献(SUN GX,HOU ZF,ZHANG CL,et al.Comparison between thequalitative similarity and the quantitative similarity of chromatographic fingerprints oftraditional Chinese medicines[J].Acta Pharmaceutica Sinica(药学学报),2007,42(1):75-80)提出的投影含量相似度C%及投影含量相似度误差△C%不能反应的差异,图中可以看出|Z1|、|Z3|是不同的,表明本发明定义的相似度能够反应的差异。
定量比较新算法与其它相似度算法:
表1是采用差异明确的模拟数据来考察4种相似度的计算结果。1~6号样品是对照样品的稀释或浓缩,即样品中的每个峰的大小与对照样品中相应峰大小关系为线性关系:xi=k*yi,k分别为0.5、0.9、1.1、1.5、2.0和3.0;7~12号样品比对照样品多出一个非共有峰,非共有峰峰大小依次由1变为25,共有峰对应峰大小相同;13号样品与3号样品的各峰的相对差异都为10%,但变化方向不同,3号样品所有峰都比对照指纹图谱峰偏大,而13号样品比对照指纹图谱峰一部分偏大,另一部分偏小;14号样品的一个小峰与15号样品的一个大峰相对对照指纹图谱的相对差异相同,都为10%;16号样品的小峰与17号样品的大峰相对对照指纹图谱的绝对差异相同,都为3。
表1.模拟数据的四种相似度结果
表1中1~6号样品的相似度数据表明当峰值大小的总差异不超过100%时,本研究定义的相似度S、改良程度相似度q、新改良程度相似度q′数值大小与样品和参照样品的相对差异大小一致,直接给出直观的综合定量评价结果;即使相对差异超过100%,数值的大小也能反映样品和参照样品的差异,本研究计算的相似度数值变为负值时,其绝对值越大,样品与对照样品的差异越大;而夹角余弦相似度cosθ全为1,不能反应样品与参照样品的这种差异。7~12号样品与对照样品相似度数据显示当非共有峰峰值由1变为25时,本研究定义相似度S由0.9892变化为0.7314,差值为0.2578;而夹角余弦相似度cosθ由0.9996变化0.8293,差值为0.1703;表明本发明定义相似度S比夹角余弦相似度cosθ能更灵敏地反映非共有峰峰大小的影响,而改良相似度q及新改良程度相似度q′由于只能反应共有峰的变化对相似度的影响,相似度全为1。相对于对照指纹图谱而言,13号样品峰比例关系的一致性比3号样品峰比例关系的一致性要差,按公式(3)计算的相似度S,13号样品的0.8184小于3号样品相似度0.9000,二者相似度△S的差值比夹角余弦值差值大,而改良程度相似度q及新改良程度相似度q′不能反映这种差异。表1中14号、15号样品相似度数据表明按公式(4)计算的相似度无论大峰还是小峰,只要相对差异相同对相似度的影响是一样的,而夹角余弦值对大小峰的灵敏度是不同的。16号、17号样品相似度数据表明按公式(4)计算的相似度能灵敏反应大小峰产生相同的绝对差异时对相似度不同的影响,对16号样品小峰的绝对差异为3时相对差异为60%,而17号样品大峰的绝对差异为3时相对差异为12%,因此17号样品的相似度远大于16号样品的相似度。
新算法评价中药制剂质量稳定性的应用
表2是文献(吴昊,田燕华,郭平平.多元统计学在参麦注射液指纹图谱中的应用[J],中成药,2002,Vol.24No.1:3-6.)中23个参麦注射液样品指纹图谱数据及采用式(4)计算的相似度结果;1~11号是以工艺A制得的参麦注射液样品,作为合格品,12~17号以工艺B制得的参麦注射液样品,2’、3’、4’、5’、9’号是2、3、4、5、9号参麦注射液样品的重复进样获得的指纹图谱,18号是以西洋参代替红参所作的伪品;对照指纹图谱是以中位矢量法用1~11号样品指纹图谱数据确定的。根据本研究确定相似度判断阈值S限的方法,发现7号样品和11号样品为异常样品,剔除7、11号样品后,将1~6号、9~10号作为合格样本,以中位法重新确定对照指纹图谱,依次以公式(4)计算1~6号、9~10号、2’、3’、4’、5’、9’号、7号、11号及12~18号样品相对新对照指纹图谱的相似度,结果按顺序为0.7671,0.8029,0.7978,0.5627,0.9052,0.6964、0.6549,0.7622、0.8118、0.7771、0.8024、0.5783、0.8972、0.7604、0.0247、0.2812、-1.132,-1.700,-1.941,-0.9393,-2.726,-2.263,-10.89;1~6号、9~10号合格样本的相似度未发现异常值,4号样本的相似度0.5627为最小,确定为判断阈值S限;显然12~18号相似度S值均小于0.5627,可判定为不合格品(工艺B与工艺A有显著差别),2’、3’、4’、5’、9’号相似度的均大于0.5627,可判定为合格品,特别是18号伪品其相似度值最低为-10.89,按新的阈值S限判断,7、11号样品依然为异常样品,分析7、11号样品指纹图谱数据可发现,二者的指纹图谱第一个色谱峰值分别为2.72、3.05,明显比其它合格样本对应峰值大得多,说明判断是合理的;即使在剔除7、11号样品前计算的相似度,从表1的数据也可以看出1~11号样本的相似度也远大于12~18号样品的相似度,表明公式(4)算法计算的相似度能够反应两种不同工艺的差别,从另一个角度说明公式(4)算法计算能够反应中药制剂质量的稳定性均匀性。
表2参麦注射液指纹图谱数据及相似度结果
综上所述,本研究构建的相似度算法扩展了相似度的取值范围为-∞~1,当相似度为负值时从原理和现实都能得到合理解释:样品质量与对照样本质量差别越显著相似度数值为负值的绝对值越大;模拟数据计算结果显示本研究构建的相似度算法不仅能够准确灵敏地反应指纹图谱共有峰峰分布比例关系的变化和峰面积总量的波动,也能够定量反应非共有峰总面积大小对相似度的影响,具有较好的综合评价能力;本研究根据箱线图法原理提出的确定相似度判断阈值的经验方法具有较强的可操作性和实用性。
本发明未涉及部分均与现有技术相同或可采用现有技术加以实现。
Claims (7)
1.一种色谱指纹图谱相似度评价方法,其特征是它包括以下步骤:
(2)、采用与建立对照指纹图谱相同的色谱条件建立待评价样品的样品指纹图谱,样品指纹图谱的向量为:
其中,1≤i≤n,n为样品指纹图谱中与对照指纹图谱共有峰的色谱峰数,xi代表待评价样品中i成分即样品指纹图谱中与对照指纹图谱共有峰的峰面积或峰高值,其中,对照指纹图谱中存在而在样品指纹图谱中不存在的成分,样品指纹图谱中该共有峰的取值为xi=0;
1≤j≤m,m为样品指纹图谱中与对照指纹图谱非共有峰的色谱峰数,xaj代表待评价样品中j成分即样品指纹图谱中与对照指纹图谱非共有峰的峰面积或峰高值,;
其中:
δi、δaj分别是各共有峰、非共有峰总面积的权重系数,权重的取值范围δi≥0、δa≥0;
3.根据权利要求2所述的色谱指纹图谱相似度评价方法,其特征是所述的色谱方法包括高效液相色谱、气相色谱或高效毛细管电泳色谱。
4.根据权利要求2所述的色谱指纹图谱相似度评价方法,其特征是所述的步骤(b)中,标准样品的批次总数为k,对于对照指纹图谱的向量按照步骤(3)计算k个批次标准样品的比值指纹图谱向量Yk与此次确定的对照比值指纹图谱向量的相似度采用箱线图法对数据按从小到大排序,并求出该序列数据的统计量第一四分位数Q1、第三四分位数Q3,并求出这两个四分位数差值df=Q3-Q1,将所有数据与Q1-1.5df比较,所有小于Q1-1.5df的数据判断为异常值而将对应批次的标准样品去除,重复前述步骤,直至剩余的不再出现异常值,对剩余批次标准样品指纹图谱进行预处理和色谱峰匹配,以平均矢量法或中位数矢量法建立对照指纹图谱,所述对照指纹图谱的向量为
6.根据权利要求1所述的色谱指纹图谱相似度评价方法,其特征是所述的步骤(3)中,δi=1。
7.根据权利要求1所述的色谱指纹图谱相似度评价方法,其特征是所述的步骤(3)中,δi的取值原则是:相应成分对产品有效性安全性影响非常明显时δi>1,该成分对产品有效性安全性影响较小时1>δi≥0。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310183164.8A CN103278591B (zh) | 2013-05-16 | 2013-05-16 | 一种色谱指纹图谱相似度计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310183164.8A CN103278591B (zh) | 2013-05-16 | 2013-05-16 | 一种色谱指纹图谱相似度计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103278591A true CN103278591A (zh) | 2013-09-04 |
CN103278591B CN103278591B (zh) | 2015-08-26 |
Family
ID=49061171
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310183164.8A Expired - Fee Related CN103278591B (zh) | 2013-05-16 | 2013-05-16 | 一种色谱指纹图谱相似度计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103278591B (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104777143A (zh) * | 2014-01-15 | 2015-07-15 | 中国人民解放军第二军医大学 | 一种用于过期药的、基于拉曼光谱的相似度鉴定方法 |
CN105651875A (zh) * | 2015-12-31 | 2016-06-08 | 河北中烟工业有限责任公司 | 一种指纹图谱的相似度评价算法 |
CN106483187A (zh) * | 2016-10-11 | 2017-03-08 | 王海燕 | 基于离子迁移谱和相似度的中药材与其易混伪品鉴别方法 |
CN106650779A (zh) * | 2016-10-17 | 2017-05-10 | 浙江和谱生物科技有限公司 | 谱图相似度的计算方法 |
CN106788962A (zh) * | 2016-12-13 | 2017-05-31 | 电子科技大学 | 隐私保护下的向量相似性判断方法 |
CN107784192A (zh) * | 2017-09-15 | 2018-03-09 | 山东大学 | 指纹图谱相似度计算方法、装置和样品质量评价系统 |
CN109100477A (zh) * | 2017-06-20 | 2018-12-28 | 香港理工大学 | 食用油分析方法、识别系统、产生库的方法及数据载体 |
CN109668850A (zh) * | 2019-02-28 | 2019-04-23 | 山东中医药大学 | 基于紫外指纹图谱的中药药性识别方法及系统 |
CN110987856A (zh) * | 2019-12-19 | 2020-04-10 | 中国检验检疫科学研究院 | 基于配方体系和指纹图谱的化妆品质量快速鉴定方法 |
CN111426648A (zh) * | 2020-03-19 | 2020-07-17 | 甘肃省交通规划勘察设计院股份有限公司 | 一种红外谱图的相似度确定方法及系统 |
CN116973495A (zh) * | 2023-09-21 | 2023-10-31 | 山东鲁地源天然药物有限公司 | 基于气相色谱仪的中药饮片检测数据分析管理系统 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002286704A (ja) * | 2001-03-28 | 2002-10-03 | Shimadzu Corp | クロマトグラフ用データ処理装置 |
US20040034477A1 (en) * | 2002-08-19 | 2004-02-19 | Mcbrien Michael | Methods for modeling chromatographic variables |
CN101256177A (zh) * | 2007-11-21 | 2008-09-03 | 皖南医学院 | 中药数字化色谱指纹谱相似度评价系统 |
CN101271092A (zh) * | 2007-03-19 | 2008-09-24 | 沈阳药科大学 | 中药色谱指纹图谱特征数字化和全定性全定量质量控制方法 |
CN101676717A (zh) * | 2008-09-19 | 2010-03-24 | 天津天士力制药股份有限公司 | 一种中药制品的质量评价方法 |
WO2011058883A1 (ja) * | 2009-11-13 | 2011-05-19 | 株式会社日立ハイテクノロジーズ | 液体クロマトグラフィー質量分析装置、及び液体クロマトグラフィー質量分析装置を用いた分析方法 |
CN102507815A (zh) * | 2011-10-27 | 2012-06-20 | 浙江工业大学 | 一种色谱指纹谱图的相似度测定方法 |
CN103018382A (zh) * | 2012-12-07 | 2013-04-03 | 南京中医药大学 | 一种指纹图谱相似度的检测方法 |
-
2013
- 2013-05-16 CN CN201310183164.8A patent/CN103278591B/zh not_active Expired - Fee Related
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002286704A (ja) * | 2001-03-28 | 2002-10-03 | Shimadzu Corp | クロマトグラフ用データ処理装置 |
US20040034477A1 (en) * | 2002-08-19 | 2004-02-19 | Mcbrien Michael | Methods for modeling chromatographic variables |
CN101271092A (zh) * | 2007-03-19 | 2008-09-24 | 沈阳药科大学 | 中药色谱指纹图谱特征数字化和全定性全定量质量控制方法 |
CN101256177A (zh) * | 2007-11-21 | 2008-09-03 | 皖南医学院 | 中药数字化色谱指纹谱相似度评价系统 |
CN101676717A (zh) * | 2008-09-19 | 2010-03-24 | 天津天士力制药股份有限公司 | 一种中药制品的质量评价方法 |
WO2011058883A1 (ja) * | 2009-11-13 | 2011-05-19 | 株式会社日立ハイテクノロジーズ | 液体クロマトグラフィー質量分析装置、及び液体クロマトグラフィー質量分析装置を用いた分析方法 |
CN102507815A (zh) * | 2011-10-27 | 2012-06-20 | 浙江工业大学 | 一种色谱指纹谱图的相似度测定方法 |
CN103018382A (zh) * | 2012-12-07 | 2013-04-03 | 南京中医药大学 | 一种指纹图谱相似度的检测方法 |
Non-Patent Citations (10)
Title |
---|
FENG GAN ET AL: "New approach on similarity analysis of chromatographic fingerprint of herbal medicine", 《JOURNAL OF CHROMATOGRAPHY A》, vol. 1104, no. 12, 3 February 2006 (2006-02-03) * |
JUN YE: "Cosine similarity measures for intuitionistic fuzzy sets and their applications", 《MATHEMATICAL AND COMPUTER MODELLING》, vol. 53, no. 12, 31 January 2011 (2011-01-31) * |
KAI-TAI FANG ET AL: "Critical value determination on similarity of fingerprints", 《CHEMOMETRICS AND INTELLIGENT LABORATORY SYSTEMS》, vol. 82, no. 12, 26 May 2006 (2006-05-26) * |
KATTY X. WAN ET AL: "Comparing Similar Spectra: From Similarity Index to Spectral Contrast Angle", 《JOURNAL OF THE AMERICAN SOCIETY FOR MASS SPECTROMETRY》, vol. 13, no. 1, 31 January 2002 (2002-01-31), XP004326927, DOI: doi:10.1016/S1044-0305(01)00327-0 * |
SHI-KAI YAN ET AL: "An approach to develop two-dimensional fingerprint for the quality control of Qingkailing injection by high-performance liquid chromatography with diode array detection", 《JOURNAL OF CHROMATOGRAPHY A》, vol. 1090, no. 12, 7 October 2005 (2005-10-07) * |
刘永锁 等: "相似系统理论用于中药色谱指纹图谱的相似度评价", 《色谱》, vol. 23, no. 2, 30 March 2005 (2005-03-30) * |
孟庆华 等: "色谱指纹图谱相似度的新算法及其应用", 《中成药》, vol. 25, no. 1, 25 January 2003 (2003-01-25) * |
詹雪艳 等: "基于相似系统理论的相似度计算方法的改进", 《分析化学》, vol. 38, no. 2, 15 February 2010 (2010-02-15) * |
詹雪艳 等: "色谱指纹图谱相似度方法的研究进展", 《中国实验方剂学杂志》, vol. 17, no. 2, 20 January 2011 (2011-01-20) * |
詹雪艳 等: "色谱指纹图谱组合相似度的算法", 《色谱》, vol. 28, no. 11, 28 November 2010 (2010-11-28) * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104777143A (zh) * | 2014-01-15 | 2015-07-15 | 中国人民解放军第二军医大学 | 一种用于过期药的、基于拉曼光谱的相似度鉴定方法 |
CN105651875A (zh) * | 2015-12-31 | 2016-06-08 | 河北中烟工业有限责任公司 | 一种指纹图谱的相似度评价算法 |
CN106483187A (zh) * | 2016-10-11 | 2017-03-08 | 王海燕 | 基于离子迁移谱和相似度的中药材与其易混伪品鉴别方法 |
CN106650779B (zh) * | 2016-10-17 | 2019-10-25 | 浙江和谱生物科技有限公司 | 谱图相似度的计算方法 |
CN106650779A (zh) * | 2016-10-17 | 2017-05-10 | 浙江和谱生物科技有限公司 | 谱图相似度的计算方法 |
CN106788962A (zh) * | 2016-12-13 | 2017-05-31 | 电子科技大学 | 隐私保护下的向量相似性判断方法 |
CN106788962B (zh) * | 2016-12-13 | 2020-04-14 | 电子科技大学 | 隐私保护下的向量相似性判断方法 |
CN109100477A (zh) * | 2017-06-20 | 2018-12-28 | 香港理工大学 | 食用油分析方法、识别系统、产生库的方法及数据载体 |
US10859552B2 (en) | 2017-06-20 | 2020-12-08 | The Hong Kong Polytechnic University | Edible oil analysis system and method |
CN107784192A (zh) * | 2017-09-15 | 2018-03-09 | 山东大学 | 指纹图谱相似度计算方法、装置和样品质量评价系统 |
CN109668850A (zh) * | 2019-02-28 | 2019-04-23 | 山东中医药大学 | 基于紫外指纹图谱的中药药性识别方法及系统 |
CN110987856A (zh) * | 2019-12-19 | 2020-04-10 | 中国检验检疫科学研究院 | 基于配方体系和指纹图谱的化妆品质量快速鉴定方法 |
CN111426648A (zh) * | 2020-03-19 | 2020-07-17 | 甘肃省交通规划勘察设计院股份有限公司 | 一种红外谱图的相似度确定方法及系统 |
CN116973495A (zh) * | 2023-09-21 | 2023-10-31 | 山东鲁地源天然药物有限公司 | 基于气相色谱仪的中药饮片检测数据分析管理系统 |
CN116973495B (zh) * | 2023-09-21 | 2023-12-15 | 山东鲁地源天然药物有限公司 | 基于气相色谱仪的中药饮片检测数据分析管理系统 |
Also Published As
Publication number | Publication date |
---|---|
CN103278591B (zh) | 2015-08-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103278591B (zh) | 一种色谱指纹图谱相似度计算方法 | |
Bailey et al. | Exponent of cross‐sectional dependence: Estimation and inference | |
Ikeda et al. | Mode-coupling theory as a mean-field description of the glass transition | |
CN106560697A (zh) | 联合近红外光谱和微量元素的武夷岩茶产地鉴别方法 | |
CN106021685B (zh) | 一种考虑测量误差的退化可靠性分析方法 | |
CN103268450A (zh) | 基于测试的移动智能终端系统安全评估体系模型与评估方法 | |
CN107390259A (zh) | 一种基于svd和svm的核素识别方法 | |
Smarandache | α-discounting method for multi-criteria decision making (α-d MCDM) | |
Yujun et al. | Multiscale multifractal multiproperty analysis of financial time series based on Rényi entropy | |
Chen et al. | Testing and support recovery of correlation structures for matrix-valued observations with an application to stock market data | |
de Carvalho Rocha et al. | A comparison of three procedures for robust PCA of experimental results of the homogeneity test of a new sodium diclofenac candidate certified reference material | |
CN101566569B (zh) | 一种通过特征参量鉴别多种荧光光谱混叠物质的系统和方法 | |
Zhu et al. | Conditional tests for elliptical symmetry | |
CN103186696B (zh) | 面向高维非线性软测量模型的辅助变量约简方法 | |
De Almeida et al. | A gage study through the weighting of latent variables under orthogonal rotation | |
Chhabra et al. | Sequential decision process for tradespace exploration by bounding probabilistic decision criteria using mean-risk analysis | |
CN104573732A (zh) | 一种目标光谱匹配方法 | |
CN105991609B (zh) | 一种风险事件确定方法及装置 | |
Zhao et al. | Distribution of eigenvalues of detrended cross-correlation matrix | |
Lingke et al. | A new method to solve robust data reconciliation in nonlinear process | |
CN102282560A (zh) | 用于产生马库什结构专利权利要求内的虚拟化合物链接库的方法 | |
Caldarelli et al. | Weighted networks as randomly reinforced urn processes | |
Bonamente | Systematic errors in the maximum-likelihood regression of Poisson count data: introducing the overdispersed χ2 distribution | |
US20130197813A1 (en) | Similarity evaluating method, similarity evaluating program, and similarity evaluating device for collective data | |
Ekong et al. | A comparative study of the impact of dummy variables on regression coefficients and canonical correlation indices: An empirical perspective |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150826 |