CN112513619A - 光谱校正装置及光谱校正方法 - Google Patents

光谱校正装置及光谱校正方法 Download PDF

Info

Publication number
CN112513619A
CN112513619A CN201980050349.7A CN201980050349A CN112513619A CN 112513619 A CN112513619 A CN 112513619A CN 201980050349 A CN201980050349 A CN 201980050349A CN 112513619 A CN112513619 A CN 112513619A
Authority
CN
China
Prior art keywords
color conversion
signal
calculated
conversion matrix
candidate
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
Application number
CN201980050349.7A
Other languages
English (en)
Other versions
CN112513619B (zh
Inventor
川口洋平
乔治·查尔基迪斯
安田知弘
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi High Tech Corp
Original Assignee
Hitachi High Technologies Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hitachi High Technologies Corp filed Critical Hitachi High Technologies Corp
Publication of CN112513619A publication Critical patent/CN112513619A/zh
Application granted granted Critical
Publication of CN112513619B publication Critical patent/CN112513619B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/27—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands using photo-electric detection ; circuits for computing concentration
    • G01N21/274—Calibration, base line adjustment, drift correction
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/64—Fluorescence; Phosphorescence
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/64—Fluorescence; Phosphorescence
    • G01N21/6408—Fluorescence; Phosphorescence with measurement of decay time, time resolved fluorescence
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/64—Fluorescence; Phosphorescence
    • G01N21/6428—Measuring fluorescence of fluorescent products of reactions or of fluorochrome labelled reactive substances, e.g. measuring quenching effects, using measuring "optrodes"
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/64—Fluorescence; Phosphorescence
    • G01N21/6428—Measuring fluorescence of fluorescent products of reactions or of fluorochrome labelled reactive substances, e.g. measuring quenching effects, using measuring "optrodes"
    • G01N2021/6439—Measuring fluorescence of fluorescent products of reactions or of fluorochrome labelled reactive substances, e.g. measuring quenching effects, using measuring "optrodes" with indicators, stains, dyes, tags, labels, marks
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/64—Fluorescence; Phosphorescence
    • G01N21/6428—Measuring fluorescence of fluorescent products of reactions or of fluorochrome labelled reactive substances, e.g. measuring quenching effects, using measuring "optrodes"
    • G01N2021/6439—Measuring fluorescence of fluorescent products of reactions or of fluorochrome labelled reactive substances, e.g. measuring quenching effects, using measuring "optrodes" with indicators, stains, dyes, tags, labels, marks
    • G01N2021/6441—Measuring fluorescence of fluorescent products of reactions or of fluorochrome labelled reactive substances, e.g. measuring quenching effects, using measuring "optrodes" with indicators, stains, dyes, tags, labels, marks with two or more labels
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2201/00—Features of devices classified in G01N21/00
    • G01N2201/12—Circuits of general importance; Signal processing
    • G01N2201/127—Calibration; base line adjustment; drift compensation

Landscapes

  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Chemical & Material Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Pathology (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Optics & Photonics (AREA)
  • Mathematical Physics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
  • Apparatus Associated With Microorganisms And Enzymes (AREA)

Abstract

本发明提供一种光谱校正装置及光谱校正方法，即使在荧光色素间同时出现峰的情况下也能够高精度地进行光谱校正。光谱校正装置算出用于颜色变换处理的颜色变换矩阵，包括：分光信号取得部，其取得随时间而检测到的荧光的分光信号；候选算出部，其基于所述分光信号，针对依赖于荧光色素的荧光的峰同时出现的频度的参数的每个值，算出颜色变换矩阵的候选；以及选择部，其基于针对各候选算出的评价值来选择颜色变换矩阵。

Description

光谱校正装置及光谱校正方法
技术领域
本发明涉及毛细管电泳型测序仪(Sequencer)的光谱校正装置及光谱校正方法。
背景技术
毛细管电泳型测序仪是一边使核酸试样在被称为毛细管的细管中电泳，一边对利用激发光的照射而产生的荧光进行分光并随时间来检测，使用检测到的信号来分析核酸的装置。具体而言，进行确定构成核酸的碱基序列的测序处理和比较DNA片段的链长的片段解析处理。在测序处理和片段解析处理的过程中实施颜色变换处理。颜色变换处理是将分光后的信号的时间序列数据变换为标识核酸试样的荧光色素的荧光的时间序列数据的处理，将与各荧光色素的光谱响应特性对应的颜色变换矩阵的伪逆矩阵乘以检测到的信号的时间序列数据来进行。
取得颜色变换矩阵的处理被称为光谱校正，在专利文献1中公开了使用由能够识别光谱的多个分子化学种类生成的电泳图的光谱校正方法。即，通过对由电泳图生成的数据矩阵的部分进行秩分析来识别具有光谱响应特性的区域，将根据识别出的区域判定的纯成分光谱响应分成簇，将各簇的代表性的光谱响应与分子化学种类相关联。
现有技术文献
专利文献
专利文献1：国际公开WO01/055699号公报
发明内容
发明要解决的课题
然而，在专利文献1中，假定在电泳图的时间区域中，表示荧光色素的各峰充分地分离，在荧光色素间同时出现峰的情况下，光谱校正的精度降低。在毛细管电泳型测序仪中，随着装置的小型化，毛细管变短，因此在毛细管内难以分离碱基，被检测的信号的峰的重叠增加。
因此，本发明的目的在于提供一种即使在荧光色素间同时出现峰的情况下，也能够高精度地进行光谱校正的光谱校正装置以及光谱校正方法。
用于解决课题的手段
为了实现上述目的，本发明是一种光谱校正装置，其算出用于颜色变换处理的颜色变换矩阵，其中，包括：分光信号取得部，其取得随时间而检测到的荧光的分光信号；候选算出部，其基于所述分光信号，针对依赖于荧光色素的荧光的峰同时出现的频度的参数的每个值，算出颜色变换矩阵的候选；以及选择部，其基于针对各候选算出的评价值来选择颜色变换矩阵。
另外，本发明是一种光谱校正方法，其算出用于颜色变换处理的颜色变换矩阵，包括：取得步骤，取得随时间而检测到的荧光的分光信号；候选算出步骤，基于所述分光信号，针对依赖于荧光色素的荧光的峰同时出现的频度的参数的每个值，算出颜色变换矩阵的候选；以及选择步骤，基于针对各候选算出的评价值来选择颜色变换矩阵。
发明效果
根据本发明，能够提供即使在荧光色素间同时出现峰的情况下，也能够高精度地进行光谱校正的光谱校正装置以及光谱校正方法。
附图说明
图1是光谱校正装置的硬件结构图。
图2是说明检测到的分光信号的例子、分光信号与颜色变换矩阵、荧光的时间强度的关系的图。
图3是实施例1的功能框图。
图4是表示实施例1的处理的流程的一例的图。
图5是表示实施例1的初始值的计算处理的流程的一例的图。
图6是表示实施例1的颜色变换矩阵的候选的计算处理的流程的一例的图。
图7是表示实施例1的候选的选择处理的流程的一例的图。
图8是表示实施例2的候选的选择处理的流程的一例的图。
图9是表示实施例3的候选的选择处理的流程的一例的图。
图10是说明迁移率修正处理的图。
图11是说明对包含上拉或下拉的颜色变换信号实施迁移率修正处理时的例子的图。
具体实施方式
以下，根据附图对本发明的光谱校正装置以及光谱校正方法的优选实施例进行说明。此外，在以下的说明以及附图中，对具有相同功能结构的构成要素标注相同的附图标记，从而省略重复说明。
实施例1
图1是表示光谱校正装置1的硬件结构的图。光谱校正装置1构成为，CPU(CentralProcessing Unit：中央处理器)2、主存储器3、存储装置4、显示存储器5、显示装置6、控制器7、输入装置8、网络适配器10通过系统总线11以能够收发信号的方式而连接。光谱校正装置1经由网络12以能够收发信号的方式与毛细管电泳型测序仪13和分光信号数据库14连接。在此，“能够收发信号”是指以电、光学方式不论有线、无线，能够相互或者从一方向另一方收发信号的状态。
CPU2是控制各构成要素的动作的装置。CPU2将存储在存储装置4中的程序、程序执行所需的数据加载到主存储器3中并执行。存储装置4是存储CPU2执行的程序、程序执行所需的数据的装置，具体而言是硬盘、SSD(Solid State Drive：固态驱动器)等记录装置或者对IC卡、SD卡、DVD等记录介质进行读写的装置。各种数据经由LAN(Local Area Network：局域网)等网络12进行收发。主存储器3是存储CPU2执行的程序、运算处理的中途经过的装置。
显示存储器5是暂时存储用于在液晶显示器等显示装置6上显示的显示数据的装置。输入装置8是操作者对光谱校正装置1进行操作指示的操作设备，具体而言是键盘、鼠标等。鼠标也可以是触控板、跟踪球等其他指点设备。控制器7是检测键盘、鼠标的状态，并将检测到的信息等输出到CPU2的装置。另外，在显示装置6为触摸面板的情况下，触摸面板也作为输入装置8发挥功能。网络适配器10是用于将光谱校正装置1与LAN、电话线路、因特网等网络12连接的装置。
毛细管电泳型测序仪13是一边使核酸试样在毛细管中电泳，一边对利用激发光的照射而产生的荧光进行分光并随时间来检测，使用作为检测到的信号的分光信号来分析核酸的装置。在核酸的分析中，通过对分光信号乘以颜色变换矩阵的伪逆矩阵，实施变换为标识核酸试样的荧光色素的荧光的时间序列数据的颜色变换处理。分光信号数据库14是保管由毛细管电泳型测序仪13取得的分光信号的数据库系统。
使用图2，对分光信号与颜色变换矩阵、荧光的时间强度的关系进行说明。图2的(a)是检测到的分光信号的一例，纵轴是信号强度，横轴是时间，示出了与被分光的波长对应的每个检测通道的波形。
图2的(b)是表示分光信号F与颜色变换矩阵S、荧光的时间强度C的关系的图。在毛细管中通过电泳而移动的核酸试样中，含有四种碱基A(腺嘌呤)、G(鸟嘌呤)、C(胞嘧啶)、T(胸腺嘧啶)，通过知晓各碱基的移动顺序，能够掌握核酸试样中的碱基序列。由于各碱基与不同的荧光色素对应，因此能够根据各荧光色素的荧光的时间强度C知晓各碱基的移动顺序，从而决定碱基序列。
由毛细管电泳型测序仪13检测的分光信号F是对不按照每个荧光色素区分荧光的时间强度C而集中检测到的信号进行分光而得的结果，所以如果使用颜色变换矩阵S，则F＝S×C。即，为了由分光信号F得到荧光的时间强度C，只要将S的伪逆矩阵即S+乘以两边来算出C＝S+×F即可，因此需要取得颜色变换矩阵S的处理即光谱校正。此外，F是波长通道编号n的行和时间编号t的列的矩阵，S是波长通道编号n的行和荧光色素编号k的列的矩阵，C是荧光色素编号k的行和时间编号t的列的矩阵。
为了取得颜色变换矩阵S，虽然有对分光信号F进行公知的独立成分分析的方法，但所得到的颜色变换矩阵S中精度不够。然而，由于在各碱基中标识不同的荧光色素，所以各荧光色素的荧光的峰同时出现的频度低。因此在本实施例中，算出多个荧光的时间强度C中各荧光色素的荧光的峰同时出现的频度变低的颜色变换矩阵S的候选，基于针对各候选计算出的评价值，从多个候选中选择颜色变换矩阵S。
使用图3对本实施例的主要部分进行说明。需要说明的是，这些主要部分既可以由专用的硬件构成，也可以由在CPU2上运行的软件构成。在以下的说明中，对本实施例的主要部分由软件构成的情况进行说明。本实施例具备分光信号取得部31、初始值算出部32、候选算出部33和选择部34。
以下，对各部分进行说明。
分光信号取得部31取得用于光谱校正的分光信号30。分光信号30可以是存储在存储装置4中的信号，也可以是经由网络适配器10从毛细管电泳型测序仪13或分光信号数据库14取得的信号。
初始值算出部32算出颜色变换矩阵S和荧光的时间强度C的初始值。也可以在初始值的计算中使用分光信号30，在计算初始值之前，进行几个预处理。作为预处理，例如进行尖峰的去除、基线的去除、波长通道间的归一化等，由此能够提高以后的处理中的精度。在颜色变换矩阵S的初始值的计算中，使用公知的独立成分分析。使用颜色变换矩阵S的初始值和分光信号F来算出荧光的时间强度C的初始值。
候选算出部33算出颜色变换矩阵S的候选。在颜色变换矩阵S的候选的计算中，使用依赖于在荧光的时间强度C中荧光色素的荧光的峰同时出现的频度即同时出现频度的参数λ。参数λ相对于同时出现频度处于负的相关关系。即，在参数λ小时，同时出现频度变多，当参数λ变大时，同时出现频度变少。针对每个参数λ的值算出颜色变换矩阵S的候选。
选择部34从针对每个参数λ的值算出的候选中，基于针对各候选算出的评价值来选择颜色变换矩阵。针对各候选算出的评价值，例如基于使用颜色变换矩阵的各候选算出的荧光的时间强度C的元素而被计算出。由于在物理上不会发生荧光的时间强度C的元素为负值的情况，所以这样的颜色变换矩阵的候选使用成为低评价值的代价函数。
使用图4，对具备以上各部的光谱校正装置1所执行的处理的流程的一例进行说明。以下，对各步骤进行说明。
(S401)
分光信号取得部31取得分光信号30。分光信号30既可以从存储装置4取得，也可以从毛细管电泳型测序仪13或分光信号数据库14取得。此外，分光信号30是由波长通道编号n的行和时间编号t的列构成的矩阵，将各元素设为x(t，n)。另外，1≤t≤T，T是表示检测时刻的点数的整数，例如T＝1000。并且，1≤n≤N，N是表示波长通道数的整数，例如N＝20。
(S402)
初始值算出部32根据分光信号30算出颜色变换矩阵S和荧光的时间强度C的初始值。此外，颜色变换矩阵S是由波长通道编号n的行和荧光色素编号k的列构成的矩阵，将各元素设为s(n，k)。另外，1≤k≤K，K是表示荧光色素的个数的整数，例如K＝4。并且，荧光的时间强度C是由荧光色素编号k的行和时间编号t的列构成的矩阵。
在图5中示出算出颜色变换矩阵S和荧光的时间强度C的初始值的处理的流程的例子，关于图5的各步骤在以下进行说明。
(S501)
初始值算出部32从分光信号30中去除尖峰信号。尖峰信号的去除例如通过以下的步骤来执行。
(1)首先，通过[数式1]，计算x(t，n)的时间方向的差分信号v(t，n)。
[数式1]
v(t，n)＝x(t，n)-x(t-1，n)
(2)设为Δ＝Δ_max。此外，Δ_max是时间编号的最大值，Δ_min是时间编号的最小值，Δ_tol是时间编号的允许值，θ1是分光信号x(t，n)的阈值，θ2是差分信号v(t，n)的阈值。
(3)搜索满足[数式2]的p，如果有这样的p，则将t∈[p-Δ_tol，p+Δ+Δ_tol]的范围的x(t，n)置换为0。
[数式2]
Figure BDA0002920604540000061
且
Figure BDA0002920604540000062
且
Figure BDA0002920604540000063
(4)从Δ减去1，如果Δ＜Δ_min则结束，否则返回(3)。
通过以上的步骤，取得去除了尖峰信号的分光信号x(t，n)。在分光信号中包含尖峰信号的情况下，由于使以后的处理的精度降低，所以优选去除尖峰信号。
(S502)
初始值算出部32从S501的计算结果x(t，n)去除基线。基线的去除例如使用从x(t，n)减去各时刻t的附近区间的最小值等公知的方法。将去除了基线的信号设为z(t，n)。
(S503)
初始值算出部32对S502的计算结果z(t，n)进行波长通道间强度的归一化处理。在归一化处理中使用[数式3]和[数式4]，取得归一化信号f(t，n)。
[数式3]
Figure BDA0002920604540000071
[数式4]
Figure BDA0002920604540000072
在波长通道间的强度存在较大差异时，变成对强度大的波长通道进行过度评价，对强度小的波长通道进行过小评价的情况，使以后的处理的精度降低。因此，优选将波长通道间的强度归一化。
此外，S501～S503是用于以后的处理的预处理，也可以不一定执行。但是，为了提高以后的处理的计算精度，优选执行。
(S504)
初始值算出部32通过独立成分分析算出颜色变换矩阵S的初始值。在独立成分分析的输入中，优选使用S503的计算结果f(t，n)。另外，在独立成分分析中，使用Fast ICA、Auxiliary-function-based ICA、Natural-gradient-based ICA等公知的方法。
此外，在独立成分分析之前，也可以根据归一化信号f(t，n)算出时间方向的差分信号g(t，n)＝f(t，n)-f(t-1，n)，将差分信号g(t，n)用于独立成分分析。独立成分分析假定输入信号遵循零均值的分布。由于差分信号g(t，n)遵循零均值的分布，所以能够使独立成分分析有效地发挥功能。基于相同的理由，也可以将从归一化信号f(t，n)减去归一化信号f(t，n)的平均值而得到的信号用于独立成分分析。
另外，也可以将对通过独立成分分析得到的颜色变换矩阵S的各元素a(n，k)实施非负化处理而得到的各元素b(n，k)作为颜色变换矩阵S的初始值。在非负化处理中例如使用[数式5]。
[数式5]
b(n，k)＝|a(n，k)|
另外，作为其他非负化处理，也可以实施从各a(n，k)减去a(n，k)的最小值的处理、将a(n，k)中的负值置换为零的处理，从而算出颜色变换矩阵S的初始值。计算出的颜色变换矩阵S(n，k)的初始值被保管在存储装置4等中。
(S505)
初始值算出部32使用S504的计算结果算出荧光的时间强度C的初始值。在C的初始值的计算中使用[数式6]。
[数式6]
C＝S+F
此外，S+是S的伪逆矩阵，F是以归一化信号f(t，n)为元素的矩阵。计算出的荧光的时间强度C(k，t)的初始值被保管在存储装置4等中。
返回图4的说明。
(S403)
候选算出部33算出颜色变换矩阵的候选。颜色变换矩阵的候选按照依存于荧光的时间强度C中荧光色素的荧光的峰同时出现的频度即同时出现频度的参数λ的每个值而被计算出。另外，参数λ例如设定为1.01、1.02、…、1.10。
在图6中示出算出颜色变换矩阵的候选的处理的流程的例子，关于图6的各步骤在以下进行说明。
(S601)
候选算出部33设定依赖于荧光色素的荧光的峰同时出现的频度的参数λ。例如设定λ＝1.01。
(S602)
候选算出部33从存储装置4等读出颜色变换矩阵S(n，k)和荧光的时间强度C(k，t)的初始值。
(S603)
候选算出部33将颜色变换矩阵S(n，k)与波长通道编号n关联地进行归一化。在归一化中例如使用[数式7]，作为归一化的结果计算出P(n|k)。
[数式7]
Figure BDA0002920604540000091
另外，也可以使用[数式8]来代替[数式7]。
[数式8]
Figure BDA0002920604540000092
其中，ξ是预定的正数。
在使用[数式8]的情况下，能够通过ξ来调整相当于颜色变换矩阵S(n，k)的各列的各荧光色素的波长响应特性相互的类似度。
(S604)
候选算出部33将包含参数λ的荧光的时间强度C(k，t)与荧光色素编号k关联地进行归一化。在归一化中例如使用[数式9]，作为归一化的结果计算出P(k|t)。
[数式9]
Figure BDA0002920604540000093
通过给出大于1.0的参数λ并使用[数式9]，根据参数λ的大小来抑制同时出现频度，在下一步骤中更新颜色变换矩阵S(n，k)和荧光的时间强度C(k，t)。
另外，也可以使用[数式10]或[数式11]来代替[数式9]。
[数式10]
Figure BDA0002920604540000101
其中、
Figure BDA0002920604540000102
[数式11]
Figure BDA0002920604540000103
其中、
Figure BDA0002920604540000104
ε是充分小的正的常数。
在使用[数式10]或[数式11]的情况下，由于容易更新为强度小的峰少的荧光的时间强度C(k，t)，所以在噪声大的情况下能够得到精度高的结果。
(S605)
候选算出部33更新颜色变换矩阵S(n，k)和荧光的时间强度C(k，t)。在颜色变换矩阵S(n，k)的更新中使用[数式12]和[数式13]，在荧光的时间强度C(k，t)的更新中使用[数式12]和[数式14]。
[数式12]
Figure BDA0002920604540000105
[数式13]
Figure BDA0002920604540000106
[数式14]
Figure BDA0002920604540000107
(S606)
候选算出部33判定是否达到了预定次数，在达到的情况下使处理进入S607，在未达到的情况下使处理返回到S603。
(S607)
候选算出部33算出与在S601中设定的参数λ对应的颜色变换矩阵R_λ(n，k)＝σ(n)S(n，k)。通过对S(n，k)乘以在波长通道间强度的归一化处理中使用的σ(n)，能够得到与归一化处理前的信号z(t，n)对应的颜色变换矩阵R_λ(n，k)。
(S608)
候选算出部33判定是否更新参数λ，在更新的情况下使处理返回S601来重新设定参数λ，例如设为λ＝1.02，在不更新的情况下结束S403的处理。参数λ例如按照1.01、1.02、…、1.10的顺序进行设定。
此外，在S403中，候选算出部33也可以使用贝叶斯估计来算出颜色变换矩阵的候选。即，通过使用Gibbs采样或Metropolis-Hastings采样等公知的马尔可夫链蒙特卡罗方法，求出将用[数式15]定义的后验概率最大化的颜色变换矩阵S(n，k)和荧光的时间强度C(k，t)。
[数式15]
Figure BDA0002920604540000111
其中、Θ＝(σnoise 2，αk，βk，λ)
这里，αk和βk是伽马分布的参数。在使用贝叶斯估计的情况下，能够解决初始值依赖性的问题，因此能够得到精度高的颜色变换矩阵。
返回图4的说明。
(S404)
选择部34从在S403中计算出的颜色变换矩阵的多个候选R_λ(n，k)中，基于针对各候选算出的评价值来选择颜色变换矩阵。在图7中示出从候选中选择颜色变换矩阵的处理的流程的例子，关于图7的各步骤在以下进行说明。
(S701)
选择部34针对每个参数λ算出颜色变换信号y_λ(t，k)。颜色变换信号y_λ(t，k)使用颜色变换矩阵的候选R_λ(n，k)和归一化处理前的信号z(t，n)，通过[数式16]而被算出。
[数式16]
Yλ＝Rλ +Z
此外，Yλ是以颜色变换信号y_λ(t，k)为元素的矩阵，Rλ +是颜色变换矩阵的候选R_λ(n，k)的伪逆矩阵，Z是以归一化处理前的信号z(t，n)为元素的矩阵。
(S702)
选择部34提取在S701中计算出的颜色变换信号y_λ(t，k)中的下拉。当某个荧光色素的颜色变换信号y_λ(t，k)在某个时刻具有充分大的峰时，将其他荧光色素的颜色变换信号y_λ(t，k)在该时刻具有的负值的峰的情况称为“下拉”。下拉在物理上不会发生，起因于在颜色变换信号y_λ(t，k)的计算中使用的颜色变换矩阵的候选R_λ(n，k)所包含的误差。因此，在本步骤中，提取颜色变换信号y_λ(t，k)中包含的下拉，基于提取出的下拉，在下一步骤以后评价颜色变换矩阵的各候选R_λ(n，k)。此外，在颜色变换矩阵的候选R_λ(n，k)中包含误差的情况下，产生下拉的同时也产生正值的峰即上拉，但在本实施例中不处理。
颜色变换信号y_λ(t，k)中的下拉通过搜索各荧光色素的颜色变换信号y_λ(t，k)为负值、且y_λ(t，k)的时间方向的差分值y_λ(t，k)-y_λ(t-1，k)从负变为正的时刻来提取。
(S703)
选择部34基于在步骤S702中提取的下拉算出代价函数c_λ。代价函数c_λ的计算例如使用[数式17]。
[数式17]
Figure BDA0002920604540000121
在此，m是提取出的下拉的编号，dm是各下拉的强度的绝对值，em是与下拉相同时刻的其他荧光色素的最大强度，D是预先决定的下拉的允许值。根据[数式17]，通过基于最大强度的总和的除法运算，能够避免零除法所导致的运算错误。
此外，基于下拉的代价函数c_λ并不限定于[数式17]，也可以是能够评价下拉的面积的总和的函数。
(S704)
选择部34根据在S703中计算出的代价函数c_λ，从候选中选择颜色变换矩阵R_λ(n，k)。在代价函数c_λ为[数式17]的情况下，选择代价函数c_λ为最小的颜色变换矩阵R_λ(n，k)。
通过以上说明的处理的流程，即使在荧光色素间同时出现峰的情况下，也能够求出适当的颜色变换矩阵R_λ(n，k)。即，即使在分光信号的时间分辨率低的情况下或荧光的峰在不同的荧光色素间重复的频度高的情况下，也能够高精度地进行光谱校正。特别是在毛细管短的小型的毛细管电泳型测序仪或测定时间短的情况下，本实施例有效。
实施例2
在实施例1中，说明了使用基于颜色变换信号y_λ(t，k)中的下拉的代价函数c_λ来选择颜色变换矩阵R_λ(n，k)的情况。在本实施例中，对使用基于对颜色变换信号y_λ(t，k)实施反卷积处理而得到的反卷积信号的代价函数c_λ的情况进行说明。此外，本实施例的整体结构与实施例1相同，因此省略说明。
使用图8对本实施例的候选的选择处理进行说明。此外，图8的处理的流程在图4的S403中执行。
(S801)
与步骤S701相同，选择部34针对每个参数λ算出颜色变换信号y_λ(t，k)。
(S802)
选择部34对在S801中计算出的颜色变换信号y_λ(t，k)进行反卷积处理。在反卷积处理中，使用公知的盲反卷积处理等。通过反卷积处理，针对每个参数λ算出反卷积信号。
(S803)
选择部34基于在步骤S802中计算出的反卷积信号来算出代价函数c_λ。
(S804)
选择部34根据在S803中计算出的代价函数c_λ，从候选中选择颜色变换矩阵R_λ(n，k)。
通过以上说明的处理流程，即使在荧光色素间同时出现峰的情况下，也能够求出适当的颜色变换矩阵R_λ(n，k)。此外，由于在反卷积信号中与颜色变换信号相比分离度变高，所以能够更容易地检测下拉。
实施例3
在实施例1中，说明了使用基于颜色变换信号中的下拉的代价函数来选择颜色变换矩阵的情况。在本实施例中，对使用基于对颜色变换信号实施迁移率修正处理而得到的迁移率修正处理信号的代价函数的情况进行说明。此外，本实施例的整体结构与实施例1相同，因此省略说明。
使用图9对本实施例的候选的选择处理进行说明。此外，图9的处理流程在图4的S403中执行。
(S901)
与步骤S701相同，选择部34针对每个参数λ算出颜色变换信号y_λ(t，k)。
(S902)
选择部34对在S901中计算出的颜色变换信号y_λ(t，k)进行迁移率修正处理。迁移率修正处理是以荧光色素的荧光的峰成为等间隔的方式使各荧光色素的荧光光谱在时间方向上移动的处理。图10中示出针对四个荧光色素A、G、C、T的荧光光谱的迁移率修正处理的一例。图10的(a)是迁移率修正处理前的荧光光谱的例子，荧光的峰不等间隔。图10的(b)是迁移率修正处理后的荧光光谱的例子，使荧光色素A和荧光色素C向时间的正方向移动，使荧光色素G保持原样，使荧光色素T向时间的负方向移动，由此荧光的峰成为等间隔。
但是，在颜色变换矩阵的候选R_λ(n，k)中包含误差的情况下，在颜色变换信号y_λ(t，k)中包含图11的(a)所示的上拉1101或下拉1102。当对包含上拉1101或下拉1102的颜色变换信号y_λ(t，k)实施迁移率修正处理时，如图11的(b)所示，峰间隔紊乱。在本实施例中，为了利用峰间隔的紊乱，在本步骤中进行迁移率修正处理。在迁移率修正处理中使用公知的方法，通过迁移率修正处理，针对每个参数λ算出迁移率修正信号h_λ(t，k)。
(S903)
选择部34算出基于在S902中计算出的迁移率修正信号h_λ(t，k)的代价函数c_λ。说明代价函数c_λ的算出步骤的一例。
(1)在k方向上累计h_λ(t，k)，得到时间序列信号u_λ(t)。即，不按荧光色素区分地处理迁移率修正信号h_λ(t，k)。
(2)对u_λ(t)进行离散傅里叶变换，得到U_λ(ω)。其中，ω是离散化的频率索引。
(3)计算功率谱V_λ(ω)＝|U_λ(ω)|2。
(4)针对规定的范围[ωbegin，ωend]内的ω，求出V_λ(ω)的最大值W_λ。W_λ表示与没有紊乱时的峰间隔相当的频率分量的大小。
(5)在规定的范围[ωbegin，ωend]中，计算V_λ(ω)的在ω方向上的累计值E_λ＝Σ(V_λ(ω))。E_λ是包含迁移率修正信号h_λ(t，k)的峰间隔的全部频率分量的功率，峰间隔越乱，则E_λ越大。
(6)作为代价函数，计算c_λ＝E_λ/W_λ。
(S904)
选择部34根据在S903中计算出的代价函数c_λ，从候选中选择颜色变换矩阵R_λ(n，k)。本实施例的代价函数c_λ是用与没有紊乱时的峰间隔相当的频率分量的大小除以包含全部频率分量的功率而得到的值，因此峰间隔越乱，则代价函数c_λ越大。因此，从候选中选择代价函数c_λ最小的颜色变换矩阵R_λ(n，k)。
通过以上说明的处理的流程，即使在荧光色素间同时出现峰的情况下，也能够求出适当的颜色变换矩阵R_λ(n，k)。此外，在本实施例的代价函数c_λ中不需要下拉或上拉的检测，因此即使在噪声多而难以进行下拉或上拉的检测的情况下，也能够选择适当的颜色变换矩阵R_λ(n，k)。
以上，对本发明的多个实施例进行了说明。本发明并不限定于这些实施例，包括各种变形例。例如，上述实施例是为了容易理解地说明本发明而详细说明的实施例，并不限定于具备所说明的所有结构。另外，能够将某实施例的结构的一部分置换为其他实施例的结构。并且，也能够在某实施例的结构中添加其他实施例的结构。另外，对于各实施例的结构的一部分，能够进行其他结构的追加、删除、置换。
附图标记说明
1：光谱校正装置、2：CPU、3：主存储器、4：存储装置、5：显示存储器、6：显示装置、7：控制器、8：输入装置、10：网络适配器、11：系统总线、12：网络、13：毛细管电泳型测序仪、14：分光信号数据库、30：分光信号、31：分光信号取得部、32：初始值算出部、33：候选算出部、34：选择部、35：颜色变换矩阵、1101：上拉、1102：下拉。

Claims (8)

1.一种光谱校正装置，其算出用于颜色变换处理的颜色变换矩阵，其特征在于，包括：
分光信号取得部，其取得随时间而检测到的荧光的分光信号；
候选算出部，其基于所述分光信号，针对依赖于荧光色素的荧光的峰同时出现的频度的参数的每个值，算出颜色变换矩阵的候选；以及
选择部，其基于针对各候选算出的评价值来选择颜色变换矩阵。
2.根据权利要求1所述的光谱校正装置，其特征在于，
所述选择部从颜色变换信号中提取下拉，使用基于所述下拉的代价函数来算出所述评价值，其中，使用针对所述参数的每个值算出的颜色变换矩阵的各候选来算出所述颜色变换信号。
3.根据权利要求1所述的光谱校正装置，其特征在于，
所述选择部对颜色变换信号执行反卷积处理来取得反卷积信号，使用基于所述反卷积信号的代价函数来算出所述评价值，其中，使用针对所述参数的每个值算出的颜色变换矩阵的各候选来算出所述颜色变换信号。
4.根据权利要求1所述的光谱校正装置，其特征在于，
所述选择部对颜色变换信号执行迁移率修正处理来取得迁移率修正信号，使用基于所述迁移率修正信号的代价函数来算出所述评价值，其中，使用针对所述参数的每个值算出的颜色变换矩阵的各候选来算出所述颜色变换信号。
5.一种光谱校正方法，其算出用于颜色变换处理的颜色变换矩阵，其特征在于，包括：
取得步骤，取得随时间而检测到的荧光的分光信号；
候选算出步骤，基于所述分光信号，针对依赖于荧光色素的荧光的峰同时出现的频度的参数的每个值，算出颜色变换矩阵的候选；以及
选择步骤，基于针对各候选算出的评价值来选择颜色变换矩阵。
6.根据权利要求5所述的光谱校正方法，其特征在于，
在所述选择步骤中，从颜色变换信号中提取下拉，使用基于所述下拉的代价函数来算出所述评价值，其中，使用针对所述参数的每个值算出的颜色变换矩阵的各候选来算出所述颜色变换信号。
7.根据权利要求5所述的光谱校正方法，其特征在于，
在所述选择步骤中，对颜色变换信号执行反卷积处理来取得反卷积信号，使用基于所述反卷积信号的代价函数来算出所述评价值，其中，使用针对所述参数的每个值算出的颜色变换矩阵的各候选来算出所述颜色变换信号。
8.根据权利要求5所述的光谱校正方法，其特征在于，
在所述选择步骤中，对颜色变换信号执行迁移率修正处理来取得迁移率修正信号，使用基于所述迁移率修正信号的代价函数来算出所述评价值，其中，使用针对所述参数的每个值算出的颜色变换矩阵的各候选来算出所述颜色变换信号。
CN201980050349.7A 2018-09-10 2019-08-09 光谱校正装置及光谱校正方法 Active CN112513619B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2018168596A JP7022670B2 (ja) 2018-09-10 2018-09-10 スペクトル校正装置及びスペクトル校正方法
JP2018-168596 2018-09-10
PCT/JP2019/031581 WO2020054292A1 (ja) 2018-09-10 2019-08-09 スペクトル校正装置及びスペクトル校正方法

Publications (2)

Publication Number Publication Date
CN112513619A true CN112513619A (zh) 2021-03-16
CN112513619B CN112513619B (zh) 2024-02-02

Family

ID=69777976

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201980050349.7A Active CN112513619B (zh) 2018-09-10 2019-08-09 光谱校正装置及光谱校正方法

Country Status (6)

Country Link
US (1) US11692937B2 (zh)
JP (1) JP7022670B2 (zh)
CN (1) CN112513619B (zh)
DE (1) DE112019003466B4 (zh)
GB (1) GB2590833B (zh)
WO (1) WO2020054292A1 (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7022670B2 (ja) 2018-09-10 2022-02-18 株式会社日立ハイテク スペクトル校正装置及びスペクトル校正方法
WO2023195077A1 (ja) * 2022-04-05 2023-10-12 株式会社日立ハイテク 塩基配列の解析方法及び遺伝子解析装置
CN115728276B (zh) * 2022-11-14 2024-01-23 中船重工安谱（湖北）仪器有限公司 爆炸物检测方法和检测系统

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS63313035A (ja) * 1987-06-09 1988-12-21 アプライド　バイオシステムズ　インコ−ポレイテツド Ｄｎａ配列決定のための実時間走査電気泳動装置
CN1358868A (zh) * 2000-11-29 2002-07-17 株式会社岛津制作所 核酸序列测定方法
US6863791B1 (en) * 2000-09-11 2005-03-08 Spectrumedix Llc Method for in-situ calibration of electrophoretic analysis systems
CN102411024A (zh) * 2011-06-15 2012-04-11 公安部第一研究所 一种基于空间校正及光谱校正的毛细管阵列电泳检测方法
CN202196015U (zh) * 2011-06-15 2012-04-18 公安部第一研究所 一种基于光谱校正的毛细管阵列电泳检测装置
CN102539397A (zh) * 2010-11-11 2012-07-04 索尼公司 荧光光谱校正方法和荧光光谱测量装置
US20130026391A1 (en) * 2011-07-25 2013-01-31 Sony Corporation Information processing apparatus, information processing method, program, and method of correcting intensity of fluorescence spectrum
CN103529568A (zh) * 2013-10-21 2014-01-22 上海交通大学 基于受激布里渊散射效应的矩形光滤波器系统
CN104870980A (zh) * 2012-12-17 2015-08-26 株式会社日立高新技术 基因型解析装置以及基因型解析方法
CN107003244A (zh) * 2015-01-27 2017-08-01 株式会社日立高新技术 多色荧光分析装置
CN107976478A (zh) * 2017-11-21 2018-05-01 南京溯远基因科技有限公司 一种基于毛细管电泳的多染料集核酸分析方法及其应用
CN108473925A (zh) * 2016-01-28 2018-08-31 株式会社日立高新技术 碱基序列确定装置、毛细管阵列电泳装置和方法

Family Cites Families (37)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5871628A (en) * 1996-08-22 1999-02-16 The University Of Texas System Automatic sequencer/genotyper having extended spectral response
US6738502B1 (en) * 1999-06-04 2004-05-18 Kairos Scientific, Inc. Multispectral taxonomic identification
US6333501B1 (en) 2000-01-27 2001-12-25 Perkin-Elmer Corporation Methods, apparatus, and articles of manufacture for performing spectral calibration
US6807490B1 (en) * 2000-02-15 2004-10-19 Mark W. Perlin Method for DNA mixture analysis
JP2001255301A (ja) * 2000-03-08 2001-09-21 Shimadzu Corp 電気泳動装置
US6734420B2 (en) * 2000-04-06 2004-05-11 Quantum Dot Corporation Differentiable spectral bar code methods and systems
US6618140B2 (en) * 2001-06-18 2003-09-09 Amnis Corporation Spectral deconvolution of fluorescent markers
DE10250100A1 (de) * 2002-10-28 2004-05-13 Leica Microsystems Heidelberg Gmbh Mikroskopsystem und Verfahren zur Analyse und Auswertung von Mehrfachfärbungen eines mikroskopischen Objekts
WO2004113557A2 (en) 2003-06-18 2004-12-29 Applera Corporation Methods and systems for the analysis of biological sequence data
DE10339312A1 (de) * 2003-08-27 2005-03-31 Leica Microsystems Heidelberg Gmbh Verfahren zur Trennung von Fluoreszenzspektren von in einer Probe vorhandenen Farbstoffen
US7417726B2 (en) * 2003-09-19 2008-08-26 Applied Biosystems Inc. Normalization of data using controls
US20090001262A1 (en) 2003-10-22 2009-01-01 Erik Visser System and Method for Spectral Analysis
US7295316B2 (en) * 2004-01-14 2007-11-13 Applera Corporation Fluorescent detector with automatic changing filters
US7383261B2 (en) * 2004-01-16 2008-06-03 Xerox Corporation Reference database and method for determining spectra using measurements from an LED color sensor, and method of generating a reference database
US7471385B2 (en) * 2005-01-13 2008-12-30 Xerox Corporation Systems and methods for selecting a reference database for determining a spectrum of an object based on fluorescence of the object
US8743137B2 (en) * 2006-04-10 2014-06-03 Edgenet, Inc. Method for electronic color matching
EP2044432B1 (en) * 2006-06-28 2023-03-29 Life Technologies Corporation Minimizing effects of dye crosstalk
US7519514B2 (en) * 2006-07-14 2009-04-14 Agilent Technologies, Inc. Systems and methods for removing noise from spectral data
JP4649635B2 (ja) * 2006-08-02 2011-03-16 独立行政法人科学技術振興機構 画像特徴抽出方法および画像圧縮方法
US8306756B2 (en) * 2006-10-26 2012-11-06 Shimadzu Corporation Method of determining base sequence of nucleic acid
US8634640B2 (en) * 2010-10-21 2014-01-21 Hewlett-Packard Development Company, L.P. Method and apparatus for selecting a color palette
EP2783193A4 (en) * 2011-11-03 2015-08-26 Verifood Ltd COST-EFFECTIVE SPECTROMETRIC SYSTEM FOR USER-EATING FOOD ANALYSIS
CA2935690A1 (en) * 2013-12-31 2015-07-09 Memorial Sloan Kettering Cancer Center Systems, methods, and apparatus for multichannel imaging of fluorescent sources in real-time
WO2015101992A2 (en) * 2014-01-03 2015-07-09 Verifood, Ltd. Spectrometry systems, methods, and applications
US9459201B2 (en) * 2014-09-29 2016-10-04 Zyomed Corp. Systems and methods for noninvasive blood glucose and other analyte detection and measurement using collision computing
WO2016125165A2 (en) * 2015-02-05 2016-08-11 Verifood, Ltd. Spectrometry system with visible aiming beam
WO2016125164A2 (en) * 2015-02-05 2016-08-11 Verifood, Ltd. Spectrometry system applications
WO2016162865A1 (en) * 2015-04-07 2016-10-13 Verifood, Ltd. Detector for spectrometry system
US10203246B2 (en) * 2015-11-20 2019-02-12 Verifood, Ltd. Systems and methods for calibration of a handheld spectrometer
MY187707A (en) * 2015-12-17 2021-10-13 Sicpa Holding Sa Security element, security arrangement, method for its production and authentication method using the same
US10254215B2 (en) * 2016-04-07 2019-04-09 Verifood, Ltd. Spectrometry system applications
US10791933B2 (en) * 2016-07-27 2020-10-06 Verifood, Ltd. Spectrometry systems, methods, and applications
WO2018056976A1 (en) * 2016-09-22 2018-03-29 Halliburton Energy Services, Inc. Methods and systems for obtaining high-resolution spectral data of formation fluids from optical computing device measurements
EP3551764A4 (en) * 2016-12-09 2020-08-05 IntegenX Inc. ELECTROPHEROGRAM ANALYSIS
WO2018151843A2 (en) 2017-02-17 2018-08-23 Life Technologies Corporation Automated quality control and spectral error correction for sample analysis instruments
JP7272349B2 (ja) * 2018-03-22 2023-05-12 凸版印刷株式会社 色対応情報生成システム、プログラム及び色対応情報生成方法
JP7022670B2 (ja) 2018-09-10 2022-02-18 株式会社日立ハイテク スペクトル校正装置及びスペクトル校正方法

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS63313035A (ja) * 1987-06-09 1988-12-21 アプライド　バイオシステムズ　インコ−ポレイテツド Ｄｎａ配列決定のための実時間走査電気泳動装置
US6863791B1 (en) * 2000-09-11 2005-03-08 Spectrumedix Llc Method for in-situ calibration of electrophoretic analysis systems
CN1358868A (zh) * 2000-11-29 2002-07-17 株式会社岛津制作所 核酸序列测定方法
CN102539397A (zh) * 2010-11-11 2012-07-04 索尼公司 荧光光谱校正方法和荧光光谱测量装置
CN102411024A (zh) * 2011-06-15 2012-04-11 公安部第一研究所 一种基于空间校正及光谱校正的毛细管阵列电泳检测方法
CN202196015U (zh) * 2011-06-15 2012-04-18 公安部第一研究所 一种基于光谱校正的毛细管阵列电泳检测装置
US20130026391A1 (en) * 2011-07-25 2013-01-31 Sony Corporation Information processing apparatus, information processing method, program, and method of correcting intensity of fluorescence spectrum
CN104870980A (zh) * 2012-12-17 2015-08-26 株式会社日立高新技术 基因型解析装置以及基因型解析方法
US20150337360A1 (en) * 2012-12-17 2015-11-26 Hitachi High-Technologies Corporation Device for genotypic analysis and method for genotypic analysis
CN103529568A (zh) * 2013-10-21 2014-01-22 上海交通大学 基于受激布里渊散射效应的矩形光滤波器系统
CN107003244A (zh) * 2015-01-27 2017-08-01 株式会社日立高新技术 多色荧光分析装置
CN108473925A (zh) * 2016-01-28 2018-08-31 株式会社日立高新技术 碱基序列确定装置、毛细管阵列电泳装置和方法
CN107976478A (zh) * 2017-11-21 2018-05-01 南京溯远基因科技有限公司 一种基于毛细管电泳的多染料集核酸分析方法及其应用

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
LI, B ETAL: "STR Fluorescence Spectral Analysis Based on Matrix Analysis", SPECTROSCOPY AND SPECTRAL ANALYSIS, vol. 35, no. 2, pages 429 - 432 *
MUNOZ-BARRUTIA ETAL: "Blind spectral unmixing of M-FISH images by non-negative matrix factorization", PROCEEDINGS OF THE 29TH ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE EMBS, pages 6247 - 6250 *

Also Published As

Publication number Publication date
GB202102022D0 (en) 2021-03-31
WO2020054292A1 (ja) 2020-03-19
JP7022670B2 (ja) 2022-02-18
US20210310951A1 (en) 2021-10-07
JP2020041876A (ja) 2020-03-19
GB2590833B (en) 2022-06-22
GB2590833A (en) 2021-07-07
CN112513619B (zh) 2024-02-02
DE112019003466T5 (de) 2021-04-22
DE112019003466B4 (de) 2023-03-30
US11692937B2 (en) 2023-07-04

Similar Documents

Publication Publication Date Title
US8392126B2 (en) Method and system for determining the accuracy of DNA base identifications
Karabiber et al. QuShape: rapid, accurate, and best-practices quantification of nucleic acid probing information, resolved by capillary electrophoresis
CN112513619B (zh) 光谱校正装置及光谱校正方法
US6681186B1 (en) System and method for improving the accuracy of DNA sequencing and error probability estimation through application of a mathematical model to the analysis of electropherograms
US6950755B2 (en) Genotype pattern recognition and classification
US20030055615A1 (en) System and methods for processing biological expression data
US6334099B1 (en) Methods for normalization of experimental data
Urfer et al. Statistics for proteomics: a review of tools for analyzing experimental data
Jin et al. Robust discriminant analysis and its application to identify protein coding regions of rice genes
US20180300451A1 (en) Techniques for fractional component fragment-size weighted correction of count and bias for massively parallel DNA sequencing
WO2017079398A1 (en) A system and method for compensating noise in sequence data for improved accuracy and sensitivity of dna testing
Maji et al. Multimodal Omics Data Integration Using Max Relevance--Max Significance Criterion
Touati et al. New methodology for repetitive sequences identification in human X and Y chromosomes
CN107077535B (zh) 用于在基因材料样本中检测微小变异体的方法和系统
Vishwakarma et al. A weight function method for selection of proteins to predict an outcome using protein expression data
JP6514369B2 (ja) 塩基配列決定装置、キャピラリアレイ電気泳動装置及び方法
US20200105374A1 (en) Mixture model for targeted sequencing
Sharma et al. Integer period discrete Fourier transform-based algorithm for the identification of tandem repeats in the DNA sequences
Lauria Rank-based miRNA signatures for early cancer detection
JP4894860B2 (ja) 核酸塩基配列信頼度の算定方法
CN116434830B (zh) 基于ctDNA多位点甲基化的肿瘤病灶位置识别方法
KR20190061771A (ko) Geo에 공개된 차세대 염기서열 데이터를 이용한 유전체 분석 방법
CN116134526A (zh) 具有合成等位基因阶梯库的dna分析仪
EP1190365A2 (en) Methods and system for amplitude normalization and selection of data peaks
Myllymäki et al. Information Theoretic Methods for Bioinformatics

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