CN105527361A - 一种高效微流电色谱相似度分析以及计算影响度的方法 - Google Patents
一种高效微流电色谱相似度分析以及计算影响度的方法 Download PDFInfo
- Publication number
- CN105527361A CN105527361A CN201410514223.XA CN201410514223A CN105527361A CN 105527361 A CN105527361 A CN 105527361A CN 201410514223 A CN201410514223 A CN 201410514223A CN 105527361 A CN105527361 A CN 105527361A
- Authority
- CN
- China
- Prior art keywords
- peak
- similarity
- sample
- time
- alignment
- 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
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明为一种高效微流电色谱相似度分析以及计算影响度的方法,该方法包括:通过对齐步骤将不同图谱上的波峰进行对齐,之后再通过相似度公式的计算,得到相关色谱图中的各检测物质的相似度数值,根据得到的数值来判断被检测物质与标准样品是否属于同一种物质。最后,通过影响程度公式来计算出各个峰在相似度计算时的影响程度。本发明可以有效的比对被检测物质与标准样品是否属于同一种物质,同时还能获取同一个被检测物质中的各个峰对相似度计算的影响程度,从而提高了高效微流电色谱的分析、辨别样品的能力。
Description
技术领域
本发明涉及高效微流电色谱样品分析领域,具体说涉及一种高效微流电色谱相似度分析以及计算影响度的方法。
背景技术
高效微流电色谱是一种非常重要的分离分析仪器,它的最大作用就是确定待测样品的名称,以及分析待测样品的成分组成。
为了能够让高效微流电色谱准确地分析出待测样品的名称,目前色谱分析领域常常在计算模块引入相似度的概念,如中国专利CN104007212。由于仪器、待测样品、色谱柱以及溶剂的不确定性,在进行相似度计算之前,往往先要把所有待测样品得到的色谱数据进行一个对齐操作,保证所有样品在同一个条件下进行相似度计算,如中国专利CN101929989A。
同时,为了能够了解到同一个样品中各个峰在相似度计算时的影响程度,需要引入影响度计算概念。影响度计算方法的引入,可以直观地让用户了解待测物质与标样之间存在差异的主要原因,帮助用户提升物质检测的效率和速度。
然而,目前现有技术还没有将对齐步骤,相似度计算,影响度计算等结合起来组成完整的技术方案,而影响度的计算法方式也比较简单,缺少高信度的计算方法。
发明内容
本发明的目的是为了提高高效微流电色谱在比对样品时的准确率以及帮助用户了解各个峰对相似度计算的影响程度,提高高效微流电色谱物质检测的准确率,提升高效微流电色谱物质检测的效率。
为了实现本发明的目的而采用的技术方案是:一种高效微流电色谱相似度分析以及计算影响度的方法,包括以下步骤:
S1,数据获取步骤,采集标准样品与未知样品的色谱信号数据;
S2,数据对齐步骤,通过对齐步骤将标准样品与未知样品的数据进行对齐;
S3,相似度计算步骤,根据设定的相似度计算公式,得到未知样品的相似度数值;
S4,影响度计算步骤,通过设定影响度计算公式,得到各个色谱峰对相似度计算的影响程度。
优选地,所述数据对齐步骤包括:
S21,从数据源选择数据列表;
S22,选择要对齐的峰,初始值为最后的峰,每一行只能选一个峰;
S23,以第一个数据为基准,进行对齐运算:
Ty=Tx*T1/T2;
T1:第一行对齐峰的响应时间;
T2:当前行对齐峰的响应时间;
Tx:对齐前峰的响应时间;
Ty:对齐后峰的响应时间;
S24,显示对齐后的数据列表。
优选地,所述相似度计算步骤包括:
S31.导入需要分析的所有样品数据,需要的项目是响应时间、峰名称、峰高/峰面积;
S32.根据时间窗/时间带的设置,遍历所有样品,计算出所有的响应时间,若某一样品在某个响应时间内没有峰的话,则峰高/峰面积设为0;
S33.以标准样品为基础,计算所有样品的相似度,预设相似度公式如下:
其中xi是第一个样品的峰高/峰面积,yi是另一个样品的峰高/峰面积,n是响应时间的个数。当其中某一项为0的时候xi/yi设为0,当两项都为0的时候xi/yi设为1。当xi>2yi时xi/yi设为0,反之当yi>2xi时xi/yi也设为0。
优选地,所述影响度计算步骤包括:
S41.取所有样品的相似度的平均值,大于该平均值的分为第一类,小于该值的分为另一类;
S42.分别计算出这两类的各个响应时间的峰高/峰面积的平均值;
S43.计算出各个峰在相似度计算时的影响程度,预设影响度公式如下:
其中xi是第一类样品的峰高/峰面积的平均值,yi是另一类样品的峰高/峰面积的平均值。
附图说明
图1是标准样品与被测样品进行相似度计算和影响度计算的步骤流程图。
具体实施方式
实施例:
本实施例的实验环境为:流动相:0.02mol/L磷酸二氢钠二乙胺缓冲溶液+甲醇(70:30);流速:1mL/min;检测波长:220nm;进样量:50μL;柱温:30℃;色谱柱:C18(l50mm×4.6mm,粒径5μm)。
本实施例具体步骤为:
1.从数据源选择数据列表;给出数据源和具体列表。同理,下同的每一步都要给出具体的数值。
2.选择要对齐的峰,初始值为最后的峰,每一行只能选一个峰;
3.以第一个数据为基准,进行对齐运算:
T1:第一行对齐峰的响应时间;
T2:当前行对齐峰的响应时间;
Tx:对齐前峰的响应时间;
Ty:对齐后峰的响应时间;
Ty=Tx*T1/T2;
4.显示对齐后的数据列表。
完成了标准样品与待测样品相关数据的对齐之后,即进入到相似度计算的步骤。相似度计算的步骤为:
1.导入需要分析的所有样品数据,需要的项目是响应时间、峰名称、峰高/峰面积;
2.根据时间窗/时间带的设置,遍历所有样品,计算出所有的响应时间,若某一样品在某个响应时间内没有峰的话,则峰高/峰面积设为0;
3.以标准样品为基础,计算所有样品的相似度,公式如下:
其中xi是第一个样品的峰高/峰面积,yi是另一个样品的峰高/峰面积,n是响应时间的个数。当其中某一项为0的时候xi/yi设为0,当两项都为0的时候xi/yi设为1。
当xi>2yi时xi/yi设为0,反之当yi>2xi时xi/yi也设为0。
完成相似度计算之后,如果相似度的数值越接近于1,说明被测样品与标准样品越有可能是同一种物质。反之,如果相似度的数值与1的差距越大,说明被测样品与标准样品越不可能是同一种物质。
完成以上这些步骤之后,即可开始进行影响度的计算,方法如下:
1.取所有样品的相似度的平均值,大于该平均值的分为第一类,小于该值的分为另一类;
2.分别计算出这两类的各个响应时间的峰高/峰面积的平均值;
3.计算出各个峰在相似度计算时的影响程度,公式如下:
其中xi是第一类样品的峰高/峰面积的平均值,yi是另一类样品的峰高/峰面积的平均值。
得到影响度计算的结果之后,如果结果越大,说明这个色谱峰对相似度计算的结果影响越大。反之,如果结果越小,说明这个色谱峰对相似度计算的结果影响越小。通过影响度的计算,用户可以很清楚地知道被测样品中哪一种物质对样品的测定起到了最大的影响。
通过对齐步骤、相似度计算步骤以及影响度计算步骤,用户能够很直观地了解标准样品与被测样品之间所存在的联系,以及是何种物质对样品的测定造成了最大的影响。从而提高了高效微流电色谱在比对样品时的准确率以及帮助用户了解了各个峰对相似度计算的影响程度。
以上所述仅为本发明的较佳实施方式,本发明的保护范围并不以上述实施方式为限,但凡本领域普通技术人员根据本发明所揭示内容所作的等效修饰或变化,皆应纳入权利要求书中记载的保护范围内。
Claims (4)
1.一种高效微流电色谱相似度分析以及计算影响度的方法,其特征在于包括以下步骤:
S1,数据获取步骤,采集标准样品与未知样品的色谱信号数据;
S2,数据对齐步骤,通过对齐步骤将标准样品与未知样品的数据进行对齐;
S3,相似度计算步骤,根据设定的相似度计算公式,得到未知样品的相似度数值;相似度公式如下:
其中xi是第一个样品的峰高/峰面积,yi是另一个样品的峰高/峰面积,n是响应时间的个数。当其中某一项为0的时候xi/yi设为0,当两项都为0的时候xi/yi设为1。当xi>2yi时xi/yi设为0,反之当yi>2xi时xi/yi也设为0;
S4,影响度计算步骤,通过设定影响度计算公式,得到各个色谱峰对相似度计算的影响程度;影响度公式如下:
其中xi是第一类样品的峰高/峰面积的平均值,yi是另一类样品的峰高/峰面积的平均值。
2.根据权利要求1所述的高效微流电色谱相似度分析以及计算影响度的方法,其特征在于:所述数据对齐步骤包括:
S21,从数据源选择数据列表;
S22,选择要对齐的峰,初始值为最后的峰,每一行只能选一个峰;
S23,以第一个数据为基准,进行对齐运算:
Ty=Tx*T1/T2;
T1:第一行对齐峰的响应时间;
T2:当前行对齐峰的响应时间;
Tx:对齐前峰的响应时间;
Ty:对齐后峰的响应时间;
S24,显示对齐后的数据列表。
3.根据权利要求1所述的高效微流电色谱相似度分析以及计算影响度的方法,其特征在于:所述相似度计算步骤包括:
S31.导入需要分析的所有样品数据,需要的项目包括响应时间、峰名称、峰高/峰面积;
S32.根据时间窗/时间带的设置,遍历所有样品,计算出所有的响应时间,若某一样品在某个响应时间内没有峰的话,则峰高/峰面积设为0;
S33.以标准样品为基础,计算所有样品的相似度,公式如下:
其中xi是第一个样品的峰高/峰面积,yi是另一个样品的峰高/峰面积,n是响应时间的个数。当其中某一项为0的时候xi/yi设为0,当两项都为0的时候xi/yi设为1。当xi>2yi时xi/yi设为0,反之当yi>2xi时xi/yi也设为0。
4.根据权利要求1所述的高效微流电色谱相似度分析以及计算影响度的方法,其特征在于:所述影响度计算步骤包括:
S41.取所有样品的相似度的平均值,大于该平均值的分为第一类,小于该值的分为另一类;
S42.分别计算出这两类的各个响应时间的峰高/峰面积的平均值;
S43.计算出各个峰在相似度计算时的影响程度,公式如下:
其中xi是第一类样品的峰高/峰面积的平均值,yi是另一类样品的峰高/峰面积的平均值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410514223.XA CN105527361B (zh) | 2014-09-29 | 2014-09-29 | 一种高效微流电色谱相似度分析以及计算影响度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410514223.XA CN105527361B (zh) | 2014-09-29 | 2014-09-29 | 一种高效微流电色谱相似度分析以及计算影响度的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105527361A true CN105527361A (zh) | 2016-04-27 |
CN105527361B CN105527361B (zh) | 2017-06-20 |
Family
ID=55769715
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410514223.XA Active CN105527361B (zh) | 2014-09-29 | 2014-09-29 | 一种高效微流电色谱相似度分析以及计算影响度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105527361B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110161161A (zh) * | 2019-07-01 | 2019-08-23 | 汕头出入境检验检疫局检验检疫技术中心 | 基于高效液相指纹图谱和聚类分析的红木鉴定方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07218491A (ja) * | 1994-01-31 | 1995-08-18 | Shimadzu Corp | クロマトグラフ用検出装置 |
US20060292631A1 (en) * | 2003-09-25 | 2006-12-28 | Per Broberg | Elastin peptide fingerprints and analysis methods for mmp12 related to copd |
US20090306932A1 (en) * | 2008-06-10 | 2009-12-10 | National University Of Ireland, Galway | Similarity index: a rapid classification method for multivariate data arrays |
CN101603959A (zh) * | 2009-06-29 | 2009-12-16 | 上海应用技术学院 | 一种快速测定猪肉香精相似度的方法 |
CN101929989A (zh) * | 2009-06-24 | 2010-12-29 | 陈喆 | 一种代谢谱峰位置分辨与对齐的方法 |
CN104007212A (zh) * | 2014-06-23 | 2014-08-27 | 华中科技大学 | 一种基于局部相似性匹配的色谱特征峰波段提取与分析方法 |
-
2014
- 2014-09-29 CN CN201410514223.XA patent/CN105527361B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07218491A (ja) * | 1994-01-31 | 1995-08-18 | Shimadzu Corp | クロマトグラフ用検出装置 |
US20060292631A1 (en) * | 2003-09-25 | 2006-12-28 | Per Broberg | Elastin peptide fingerprints and analysis methods for mmp12 related to copd |
US20090306932A1 (en) * | 2008-06-10 | 2009-12-10 | National University Of Ireland, Galway | Similarity index: a rapid classification method for multivariate data arrays |
CN101929989A (zh) * | 2009-06-24 | 2010-12-29 | 陈喆 | 一种代谢谱峰位置分辨与对齐的方法 |
CN101603959A (zh) * | 2009-06-29 | 2009-12-16 | 上海应用技术学院 | 一种快速测定猪肉香精相似度的方法 |
CN104007212A (zh) * | 2014-06-23 | 2014-08-27 | 华中科技大学 | 一种基于局部相似性匹配的色谱特征峰波段提取与分析方法 |
Non-Patent Citations (3)
Title |
---|
乔善磊等: "中药指纹图谱相似度计算的规范化研究", 《第二军医大学学报》 * |
刘永锁: "相似系统理论用于中药色谱指纹图谱的相似度评价", 《色谱》 * |
程翼宇: "化学指纹图谱的相似性测度及其评价方法", 《化学学报》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110161161A (zh) * | 2019-07-01 | 2019-08-23 | 汕头出入境检验检疫局检验检疫技术中心 | 基于高效液相指纹图谱和聚类分析的红木鉴定方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105527361B (zh) | 2017-06-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2020106686A8 (en) | Method and system for enhancing a user interface for a web application | |
TWI445958B (zh) | 氣體檢測系統、裝置及方法 | |
CN105223264A (zh) | 一种质谱定量分析的模拟内标方法、装置及应用 | |
NZ629703A (en) | Methods and compositions for detecting multiple analytes with a single signal | |
NZ701972A (en) | Test strip and detecting device | |
WO2015096788A1 (zh) | 用于检测容装体中的样品的拉曼光谱检测方法 | |
CN103558257B (zh) | 一种基于阵列式传感器的农药多残留检测仪 | |
CN203881678U (zh) | 一种煤矿井下便携式sf6检测装置 | |
EP2698727A3 (en) | Terminal and method for determining type of input method editor | |
MX2018001134A (es) | Aparato y metodo de deteccion de gas. | |
CN102736055B (zh) | 三相智能电表计量准确度的在线检测方法 | |
WO2016182885A3 (en) | Apparatus for determining nitrate levels, and method for measuring ion concentration with a standard deviation correction | |
CN103743844B (zh) | 甲醇中乙醇含量的测定方法 | |
CN104569886B (zh) | 基于时频参数标准仪器方式的信号检测设备校准方法 | |
CN103105362A (zh) | 一种多量程水质在线检测仪 | |
CN103308634A (zh) | 含气量测量装置及方法 | |
CN1979172A (zh) | 一种气体分析方法及其装置 | |
CN105527361A (zh) | 一种高效微流电色谱相似度分析以及计算影响度的方法 | |
CN103954691A (zh) | 一种材料成分分数无损检测方法 | |
CN205719955U (zh) | 非甲烷总烃的连续监测装置 | |
CN103376397B (zh) | 一种异步电路的检测系统 | |
Laursen et al. | Enhanced monitoring of biopharmaceutical product purity using liquid chromatography–mass spectrometry | |
Li et al. | Spectral correlation of high-performance liquid chromatography-diode array detection data from two independent chromatographic runs: Peak tracking in pharmaceutical impurity profiling | |
CN109932278B (zh) | 一种固定床反应器气相停留时间分布测定系统及方法 | |
BR112018010634A2 (pt) | método analítico |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |