CN1091830A - 粒子分析装置及方法 - Google Patents
粒子分析装置及方法 Download PDFInfo
- Publication number
- CN1091830A CN1091830A CN93121213A CN93121213A CN1091830A CN 1091830 A CN1091830 A CN 1091830A CN 93121213 A CN93121213 A CN 93121213A CN 93121213 A CN93121213 A CN 93121213A CN 1091830 A CN1091830 A CN 1091830A
- Authority
- CN
- China
- Prior art keywords
- particle
- distributed data
- data
- similar degree
- mathematical expression
- 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
- 239000002245 particle Substances 0.000 title claims abstract description 111
- 238000000034 method Methods 0.000 title claims description 13
- 238000009826 distribution Methods 0.000 claims abstract description 34
- 201000010099 disease Diseases 0.000 claims abstract description 14
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 claims abstract description 14
- 238000006243 chemical reaction Methods 0.000 claims description 14
- 238000004458 analytical method Methods 0.000 claims description 8
- 210000000601 blood cell Anatomy 0.000 claims description 6
- 238000001514 detection method Methods 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 4
- 238000001303 quality assessment method Methods 0.000 claims description 3
- 230000006870 function Effects 0.000 description 8
- 239000000463 material Substances 0.000 description 8
- 238000007405 data analysis Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 230000008859 change Effects 0.000 description 5
- 238000007726 management method Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 230000035945 sensitivity Effects 0.000 description 4
- 241001269238 Data Species 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 238000007689 inspection Methods 0.000 description 3
- 210000000265 leukocyte Anatomy 0.000 description 3
- 238000002203 pretreatment Methods 0.000 description 3
- 210000004369 blood Anatomy 0.000 description 2
- 239000008280 blood Substances 0.000 description 2
- 239000012141 concentrate Substances 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 238000005315 distribution function Methods 0.000 description 2
- 230000003203 everyday effect Effects 0.000 description 2
- 210000003714 granulocyte Anatomy 0.000 description 2
- 238000011835 investigation Methods 0.000 description 2
- 210000004698 lymphocyte Anatomy 0.000 description 2
- 238000013316 zoning Methods 0.000 description 2
- 241000345998 Calamus manan Species 0.000 description 1
- 206010018910 Haemolysis Diseases 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 238000003556 assay Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000008588 hemolysis Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 210000001616 monocyte Anatomy 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 239000011164 primary particle Substances 0.000 description 1
- 235000012950 rattan cane Nutrition 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000011282 treatment Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/10—Investigating individual particles
- G01N15/1031—Investigating individual particles by measuring electrical or magnetic effects
- G01N15/12—Investigating individual particles by measuring electrical or magnetic effects by observing changes in resistance or impedance across apertures when traversed by individual particles, e.g. by using the Coulter principle
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/10—Investigating individual particles
- G01N15/14—Optical investigation techniques, e.g. flow cytometry
- G01N15/1456—Optical investigation techniques, e.g. flow cytometry without spatial resolution of the texture or inner structure of the particle, e.g. processing of pulse signals
- G01N15/1459—Optical investigation techniques, e.g. flow cytometry without spatial resolution of the texture or inner structure of the particle, e.g. processing of pulse signals the analysis being performed on a sample stream
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/01—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials specially adapted for biological cells, e.g. blood cells
- G01N2015/016—White blood cells
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/10—Investigating individual particles
- G01N2015/1019—Associating Coulter-counter and optical flow cytometer [OFC]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/10—Investigating individual particles
- G01N15/1031—Investigating individual particles by measuring electrical or magnetic effects
- G01N15/12—Investigating individual particles by measuring electrical or magnetic effects by observing changes in resistance or impedance across apertures when traversed by individual particles, e.g. by using the Coulter principle
- G01N2015/139—Measuring the ratio of AC/DC impedances
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/10—Investigating individual particles
- G01N15/14—Optical investigation techniques, e.g. flow cytometry
- G01N2015/1477—Multiparameters
Landscapes
- Chemical & Material Sciences (AREA)
- Dispersion Chemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Biological Materials (AREA)
Abstract
公开了一种能定量而且简单容易对患有某种疾
病的可能性进行判断以及对装置的精度管理及性能
进行评价的粒子分析装置。在该装置中,粒子检测器
能对多个粒子中的每一个粒子进行测定,并把对每个
粒子测得的2种数据,通过变换单元22变换成二元
分布数据。将相对于粒子的基准分布数据存储在数
据库存储单元26中,根据二元分布数据和上述基准
分布数据,类似度计算单元28能计算出二元分布数
据相对于基准分布数据的类似度δ,由显示装置20
示出该类似度。
Description
本发明涉及粒子分析装置及方法,尤其涉及可以将作为粒子分析结果获得的数据,与作为基准分布数据的类似程度数值化的装置及方法。
以前,作为粒子分析装置,有例如对多个粒子的每一个粒子分别测定来获得数据,根据这些测定数据的不同将粒子分类的装置。例如,在粒子是血球的情况下,将悬浮着血球的检样顺次供给粒子检测装置,每检测一次粒子,粒子检测装置例如会产生两种数据X1,X2。将这些粒子各自获得数据X1,X2变换成特征参数空间中的分布数据。此处,分布数据,是将数据X1,X2各自作为座标轴的二元点图F(X1,X2),根据该点图,将血球例如分类为粒细胞,淋巴细胞,单核细胞等。
在上述粒子分布装置中,如果观察点图的分布图形,至少可以判断它是正常检样,或是异常检样;如果仔细研究分布图形,还可以推断疾病。因此,以前,对使用者提供记载有各病例中典型点图例的病例集,用户通过比较由粒子分析装置获得的点图和病例集中所示的点图,观察两者的类似性,用作为检查资料。
然而,依赖这种人类的主观性,判断时要花费很多工夫,而且判断结果也不太明确,因此很难说这种状况能有效地利用。因此,希望研究出一种能够通过用数值表示粒子分析装置得到的点图和疾病典型例点图的类似度,定量地而且简单地了解患疾病的可能性。
而且,作为粒子分析装置,有必要对该粒子分析装置作出是否正常动作的判断的管理,即精度管理。以前,该精度管理,是对粒子分析装置的每个测定项目都求出作为精度管理指标的参数,通过检查该参数的变动来进行。然而,如果想提高该精度管理的精度(是否正常动作的检查能力),就必须增加参数的数目,由粒子分析获得的分布数据尽可能多元以致从二元到三元,管理的参数也增多,因此,这种管理既费事又费时。
因此,希望研究出一种精度管理耗时短而且容易进行的粒子分析装置。
而且,制造粒子分析装置时,还有必要评价由该粒子分析装置得到的分布数据是否有再现性和稳定性(同一检样在多次分析时是否可得到同样的分布数据)。也就是说必需进行性能评价,以鉴定制得的粒子分析装置,是否能作为制品在市场上销售。历来在进行这种评价时都是按照与精度管理同样的方法来进行,但是这种评价既费事又费时。
因此,希望研究出一种能在短时间内而且进行性能评价的粒子分析装置。
本发明的目的是提供一种能满足上述各项要求的粒子分析装置及方法。能达到上述目的的粒子分析装置的特征是,在具有通过对多个粒子中的每一个粒子进行测定,对上述每一个粒子获得至少一种数据的粒子检测单元,将这些数据变换成分布数据的单元,并根据上述分布数据将上述各种粒子进行分类的分类单元的粒子分析装置中,具备存储有相对于上述粒子的基准分布数据的数据库存储单元,和根据上述分布数据和上述基准分布数据,计算出上述分布数据相对于上述基准分布数据的分布类似程度δ的类似度计算单元。
而且,上述类似度计算单元可以对上述分布数据F(X1,X2……Xn)和上述基准分布数据G(X1,X2……Xn)的指定区域的数据,根据〔数学式3〕,计算出类似度。此外,还可以设上述粒子为血球,将上述基准分布数据,作为某种疾患中典型的分布数据(血球分布数据)。而且,还可将上述基准分布数据,作为粒子分析装置的精度管理用的基准分布数据,或者作为质量评价用的基准分布数据。
按照本发明的粒子分析方法,它具备通过对多个粒子中的每一个粒子进行测定,对每个粒子获得至少一种数据的过程,和将这些数据变换成分布数据的过程以及根据上述分布数据和基准分布数据计算出上述分布数据相对于上述基准分布数据的类似度δ的过程。还可以通过计算该类似度δ的过程,即相对于上述分布数据F(X1,X2……Xn)和上述基准分布数据G(X1,X2……Xn)的指定区域,根据〔数学式4〕,计算出类似度。
利用本发明的粒子分析装置和方法,一旦获得分布数据,就可以根据该分布数据和基准分布数据,获得数值形式的分布数据相对于基准分布数据的类似度δ。
如果采用某种疾病的分布数据作为基准分布数据,就可以从计算出的类似度δ,获得检样是否患该疾病的可能性,即诊断资料信息。
此外,如果采用粒子分析装置的精度管理用的基准分布数据,或者是质量评价用的基准分布数据,作为基准分布数据,从类似度δ可以知道装置是否正常,或者是偏离正常的程度,从而进行精度管理和对装置进行评价。
实施例
按照本发明的粒子分析装置的一个实施例,如图2所示,具有粒子检测器10。该粒子检测器10,是用于测定那些为了测定而进行过前处理的检样,所说的前处理,例如有;对人的血液进行稀释和溶血处理等,然后再将白血球作为粒子使其悬浮,经过这些前处理的检样,就可用于测定。该粒子检测器10,当粒子每次通过它时,都能发出多种信号。
例如,作为粒子检测器10,如果使用流体观察窗,则每个粒子都产生散乱光和萤光信号。而且,作为粒子检测器10,设有粒子通过的微孔,如果使用中同时供给该微孔以直流电流和高频电流,则每当粒子通过微孔时,就可检测出基于直流阻抗变化的信号,和基于高频阻抗变化的信号。
由粒子检测器10测出的多种信号,经过放大电路12放大后,在波形处理电路14中经过波形处理,再在A/D转换器16中进行A/D变换。
来自该A/D转换器16的数字信号,供给数据分析装置。数字化的粒子检测信号,在该数据分析装置18中被分析,根据其分析结果进行分类,其结果由显示装置20显示出来。
该数据分析装置18,例如由微型计算机构成,通过软件实现图1所示的单元。输到该数据分析装置18的数字信号,包括信号X1和信号X2,其中,用粒子检测器10对每一个粒子测定其直流阻抗,然后把基于这种阻抗变化的信号进行数字化后即为信号X1,而把基于高频阻抗变化的信号进行数字化后即为信号X2。
数据分析装置18,具有数据变换单元22。该数据变换单元22,将每个粒子发出的信号X1,X2转换成特征参数空间的分布数据F(X1,X2)。
也就是说,数据变换单元22具有刻度,该刻度如图3所示有从0至m1合计为m1+1条通道(序列)构成的参数X1,和从0至m2合计为m2+1条通道(序列)构成的参数X2。因此,该刻度被划分成〔(m1+1)*(m2+1)〕的基本要素。
而且,每当X1,X2信号输入进数据变换单元22中,与该输入的X1的大小和X2的大小相对应的基本要素记忆值随之增加。因此,当对应于全部的X1,X2信号都输进数据换单元22中时,例如,如果F(1,1)是2,则参数X1是1,参数X2是1的粒子,表示全部粒子中有2个。
如此这般,通过数据变换单元22获得的分布数据F(X1,X2),供给群集(Cluster ring)单元24,进行聚集成组(分类),例如进行白血球单细胞、粒细胞、淋巴细胞等的分类。
在该数据分析装置18中,还设有数据库存储单元26,在其中存储有作为基准的基准分布数据G(X1,X2)。图4示出该基准分布数据G(X1,X2)的一例,其构造与数据变换单元22的刻度相同,例如G(1,1)表示X1参数为1,X2参数也为1的粒子有10个存在于全部粒子之中。
作为该基准分布数据,例如可以使用某种疾病中典型的二元分布数据,和粒子控制装置是正常状态下分析某种精度管理用物质而得到的二元分布数据(精度管理用二元分布数据)。
根据从数据变换单元22得到的二元分布数据F(X1,X2),和数据库存储单元26的基准分布数据G(X1,X2),由类似度计算单元28算出二元分布数据F(X1,X2)相对于基准分布数据G(X1,X2)的类似度δ。该计算结果由显示装置20显示出来。
具体地说,根据〔数学式5〕,类似度计算单元28算出类似度。
数学式5中,类似度δ的计算区域,对参数X1规定为a1至b1通道,对参数X2规定为a2至b2通道。类似度δ的计算区域,没有必要是特征参数空间的全部(对参数X1来说是从0通道至m1通道,对参数X2来说是从0通道到m2通道),还可以根据目的设定最适宜的区域。例如,在可以预先判断没有数据分布的区域时,也可以仅就该区域之外的部分计算类似度。
具体地说,如果根据图3、图4的分布数据F(X1,X2)和基准分布数据G(X1,X2)来计算类似度δ,则可按下述方法进行。如图3及图4所示,设定a1=1通道,b1=2通道,a2=1通道,b2=4通道,如此进行计算。
数字式5的分子,用〔数学式6)表示。
而分母的前部用〔数学式7〕表示。
该分母的后半部用〔数学式8〕表示。
因此,类似度δ,为265/(316.416)1/2=0.731。
数学式5表示函数F(X1,X2)、G(X1,X2)的相关系数。当分布数据F(X1,X2)、和基准数据相似时,即F(X1,X2)=K·G(X1,X2)时(但K是常数),类似度δ为最大值1。如两者的分布是分离地排列(或者说是其分布呈无类似性排列),则这时δ值变小,接近0。即,类似度δ定量地表示二元分布的类似性。
以下从理论上进行说明,将数学式5相对于分布的变化会显示出什么样的行为进行了分析。该分析装置中使用的数据是作为离散数据的数字信号,而为简便起见将分布数据和基准分布数据都假定为连续分布,则类似度δ按照二元分布,可用〔数学式9〕表示。
为了说明方便,首先考虑为连续的一元分布,类似度δ用〔数学式10〕表示(参见株式会社昭晃堂发行,土井康弘,安藤繁共著、〔画像处理理论〕第11页)。
如果设函数元的分布函数f(X)为正规分布N(m,σ),则f(x)用〔数学式11〕表示。
如果以相对于f(x)将σ变换成aσ,并将m变换成am+b的函数作为g(x),则该g(x)为正规分布N(am+b、a2σ2),用〔数学式12〕表示(参看图5)。
f(x)、g(x)从负无限大到无限大的积分值都是1。
使用数学式11,计算数学式10的分母前半部分时,可用〔数学式13〕表示。
同样,用数学式12计算数学式10的分母后半部分时,可用〔数学式14〕表示。
同样,用数学式11、数学式12计算数学式10的分子时,即成为〔数学式15〕。
式中,如果C的定义如同〔数学式16〕表示,则数学式15可用〔数学式17〕表示。因此,δ可用〔数学式18〕表示。
如果考虑是在移动的情况下,也就是说,相对于函数元的分布数据f(x)来说,g(x)的分散σ和f(x)相同(a=1),其平均值从f(x)只移动到b的m+b的情况下,σ用〔数学式19〕表示。
〔表1〕中示出设m=100,σ=10时移动量b和类似度σ的关系。由该表可看出,δ的移动量b越大,δ就越急剧地减少。
此外,灵敏度变化时,也就是说,相对于函数元的分布数据f(x)来说,g(x)的分散是aσ,其平均值是f(x)的a倍,即am(b=0)的情况下,δ用〔数学式20〕表示。
与上述相同,〔表2〕中示出设m=100,σ=10时灵敏度a和类似度δ的关系。从表2可看出,随着a从1逐渐增大,δ急剧地减少。
以下是考虑二元分布的情况。δ用数学式9表示,设函数元的分布函数为正规分布时,f(x)可用〔数学式21〕表示。但是,X=(x,y)。式中,通过适当的变数变换,即使设K=σ2I(I为单位行列),也不会失去〔数学式21〕的普遍性。因此,数学式9,可用〔数学式22〕表示。
相对于该f(x,y),在x方向将σ变换成axσ,将mx变换成axm+bx,在y方向将σ换成ayσ,将my变换成aymy+by的数据变化函数规定为g(x,y)。该函数用〔数学式23〕表示(参看图6)。当然,不用说就是〔数学式24〕。
当f(x,y)用〔数学式25〕表示时,〔数学式26〕成立,因此,数学式9的分母前半部分,可用〔数学式27〕表示。而且,数学式9的分母后半部分,可用〔数学式28〕表示。因此数学式9的分母,可用〔数学式29〕表示。
根据数学式22、23,数学式9的分子部分可用〔数学式30〕表示。但是,如果设C1、C2为〔数学式31〕中所示,则可被表示成〔数学式32〕。
因此,数学式9可用〔数学式33〕表示,与一元的情况相同。
与一元的情况相同,如果考虑移动和灵敏度变化,则移动的情况(ax=ay=1)下,δ可用〔数学式34〕表示。
此外,灵敏度变化的情况(bx=by=0)下,δ可用〔数学式35〕表示。因此,二元分布f(x,y),g(x,y)距离越远,δ越是急剧地减少。不用说,三元分布也是同样的。
因此,如果引出图7(a)、(b)中所示的对角线,即使是分布位置的差异很难判别的分布,也可以通过求出类似度δ来明确地判断两者的差异。
图8中示出,关于用该粒子分析装置对健康正常的100个人的血液检样进行分析而得到的白血球二元分布数据,调查分布数据相互间的类似度δ(数据数100C2=4950)的结果。健康正常人相互间的类似度δ集中在0.8乃至0.9,正常检样和异常检样之间的类似度δ成为比上述0.8乃至0.9更低的值。但是,低到如何程度,根据检样而不同。
将该粒子分析装置用于帮助诊断时,在数据库存储单元24中,存储有各种疾病的典型分布数据Gi,和对应于该分布数据Gi的病例名。也就是说,典型的分布数据记忆着每一种相应的疾病。各典型的数据Gi(病例分布数据),是由粒子分析装置实际地分析各种疾病名称清楚的检样,并将此时得到的分布数据直接存储在数据库存储单元24中。
根据测得的分布数据,和数据库存储单元24的病例分布数据,就可以逐次求出测得的分布数据相对于各病例分布数据Gi的类似度δ。将与类似度δ显示出较大数值的病例分布数据Gi相对应的病例名作为新的检查资料,和测定结果一起在显示装置20中显示出来。
所谓类似度δ大,是指被测定的分布数据和比较的病例分布数据Gi相类似,被测定的检样,可判断为符合其病例i的可能性高。于是,通过将各种病例的分布数据Gi存储在数据库存储单元24中,就可以由测定的分布数据检索病例。
而且,相对于一个病例i,预先准备了例如Gi1,Gi2,Gi3等许多个病例分布数据,计算出相对于各病例分布数据Gi1,Gi2,Gi3的被测定分布数据的类似度δi1,δi2,δi3,进而求出其平均δi=(δi1+δi2+δi3)/3,就可以获得每一种疾病的高精度的类似度。
此外,用于精度管理时,粒子分析装置在正常状态,将某种精度管理用物质经测定得到的分布数据,作为一个基本分布数据,预先存储在数据库存储单元24中。而且,例如每天,在开始用该粒子分析装置进行分析之前,在该粒子分析装置中测定上述精度管理用物质,获得分布数据后,计算出该分布数据与上述分布数据的类似度δ,通过判断该类似度δ的值是否在规定值以上,就可以判断其分析精度是否足够,从而能对该粒子分析装置的精度进行管理。而且,当类似度δ比规定值小时,显示装置20即发出警报。
此外,每天都分析该精度管理物质,算出类似度δ,通过测定该δ随天数的变化,可以了解该粒子分析装置的精度随天数变化的情况。进而,数据库存储单元24中预先存储相对于精度管理物质的多个基本数据,例如G1,G2,G3,对各基本数据G1,G2,G3,算出与测定的分布数据的类似度δ1,δ2,δ3,然后取其平均值算出类似度δ,进而可得到精度高的δ值。于是,通过使用类似度δ,就可以对该粒子分析装置进行精度管理。
此外,将该粒子分析装置用于生产中的检查时,数据库存储单元24中,预先存储粒子分析装置正常动作时,测定某种基准管理用物质时的基准分布数据。而且,实际上对上述基准管理用物质进行测定,获得分布数据后,计算出相对于该上述基准分布数据的类似度δ。进行多次测定,观察每次获得的类似度δ的再现性(类似度δ的相同程度),就可以知道分布数据的再现性,继而知道该粒子分析装置的再现性。
此外,使温度变化,在各种温度下测定上述基准管理用物质,获得分布数据后,计算出类似度δ,通过观察该δ值随温度变化的变化,就可简单地了解分布数据的温度稳定性,进而可了解粒子分析装置的温度稳定性。于是,通过观察类似度δ,就可以评价该粒子分析装置的性能。
在上述实施例中,虽然是从1个粒子获得2种数据,但是也可以获得更多的数据,而且,还可以最低限度,从1个粒子获得1种数据。此外,数据变换单元22,也可以对1种数据划分成适当的序列,求出柱状统计图。
如上所述,按照本发明方法,可以根据测得的分布数据和预定的基准分布数据,以数值形式获得测定的分布数据相对于基准分布数据的类似度δ。因此,可对至今虽然包含有许多资料信息却不能对它加以利用、仅凭人类的主观意识加以判断的分布数据,进行定量地评价,本发明特别适合于对疾病的诊断和精度管理以及粒子分析装置进行定量地评价。
以下简单说明附图。
图1是本发明粒子分析装置的一个实施例的数据分析装置18的方块图。
图2是该实施例的方块图。
图3示出设置在该实施例数据变换单元22中的存储内容。
图4示出该实施例的数据库存储单元24的内容。
图5用于说明一元的2个数据的相关系数。
图6用于说明二元的2个数据的相关系数。
图7示出由该粒子分析装置分析而得的二元分布数据。
图8示出对由该实施例的粒子分析装置分析的白血球的二元分布数据,调查粒度分布相互间的类似度的结果。
图中所用符号说明如下。
18数据分析装置
22数据变化单元
26数据库存储单元
28类似度计算单元
〔数学式1〕〔数学式2〕〔数学式3〕〔数学式4〕的内容相同
〔数学式5〕
〔数学式6〕
2 4
∑ ∑ F(X1,X2)·G(X1,X2)=F(1,1)·G(1,1)+F(2,1)·G(2,1)
X1=1 X2=1 +F(1,2)·G(1,2)+F(2,2)·G(2,2)
+F(1,3)·G(1,3)+F(2,3)·G(2,3)
+F(1,4)·G(1,4)+F(2,4)·G(2,4)
=2·10+3·6
+1·9+4·8
+4·6+11·7
+7·5+10·5
=265
[数学式7]
2 4
∑ ∑|F(X1,X2)|2=|F(1,1)|2+|F(2,1)|2
X1=1 X2=1 +|F(1,2)|2+|F(2,2)|2
+|F(1,3)|2+|F(2,3)|2
+|F(1,4)|2+|F(2,4)|2
=22+32
+12+42
+42+112
+72+102
=316
[数学式8]
2 4
∑ ∑|G(X1,X2)|2=|G(1,1)|2+|G(2,1)|2
X1=1 X2=1 +|G(1,2)|2+|G(2,2)|2
+|G(1,3)|2+|G(2,3)|2
+|G(1,4)|2+|G(2,4)|2
=102+62
+92+82
+62+72
+52+52
=416
〔数学式9〕
〔数学式10〕
〔数学式11〕
数学式12〕
〔数学式13〕
〔数学式14〕
〔数学式15〕
〔数学式16〕
C=exp-((am+b-m)2/2(a2+1)σ2)
〔数学式17〕
〔数学式18〕
〔数学式19〕
〔数学式20〕
〔数学式21〕
〔数学式22〕
〔数学式23〕
〔数学式24〕
〔数学式25〕
f(x,y)=f1(x)f2(y)
〔数学式26〕
∫∫f(x,y)dxdy=∫f1(x)dx·∫f2(y)dy
〔数学式27〕
〔数学式28〕
〔数学式29〕
〔数学式30〕
〔数学式31〕
〔数学式32〕
〔数学式33〕
〔数学式34〕
〔数学式35〕
Claims (11)
1、一种粒子分析装置,该装置包括:
通过对多个粒子中每一个粒子进行测定来对上述每一粒子获得至少一种数据的粒子检测单元;
将这些数据变换成分布数据的变换单元;
存储着相对于上述粒子的基准分布数据的数据库存储单元;
根据上述分布数据及上述基准分布数据计算出上述分布数据相对于上述基准分布数据的类似度δ的类似度计算单元。
3、根据权利要求1所述的粒子分析装置,在该装置中,所述粒子是血球,所述基准分布数据是于某种疾病中的典型分布数据。
4、根据权利要求1所述的粒子分析装置,在该装置中,所述基准分布数据是粒子分析装置的精度管理用的基准分布数据。
5、根据权利要求1所述的粒子分析装置,在该装置中,所述基准分布数据是装置质量评价用的基准分布数据。
6、根据权利要求1所述的粒子分析装置,在该装置中,所述粒子检测单元包括一个粒子检测器、一个放大电路、一个波形处理电路和一个A/D转换器。
7、根据权利要求1所述的粒子分析装置,该装置还包括显示单元,用于显示由计算单元计算出的类似度δ。
8、根据权利要求7所述的粒子分析装置,在该装置中,由显示单元显示的类似度δ包括用数字表示的值。
9、根据权利要求7所述的粒子分析装置,在该装置中,由显示单元显示的类似度δ包括图象。
10、一种粒子分析方法,它包括以下步骤:
对多个粒子中每一个粒子进行测定从而对上述每一个粒子获得至少一种数据;
将这些数据变换成分布数据;以及
根据上述分布数据和基准数据计算出上述分布数据相对于上述基准分布数据的类似度δ。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP35989392A JP3213097B2 (ja) | 1992-12-28 | 1992-12-28 | 粒子分析装置及び方法 |
JP4-359893 | 1992-12-28 | ||
JP359893/92 | 1992-12-28 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1091830A true CN1091830A (zh) | 1994-09-07 |
CN1054203C CN1054203C (zh) | 2000-07-05 |
Family
ID=18466835
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN93121213A Expired - Fee Related CN1054203C (zh) | 1992-12-28 | 1993-12-25 | 粒子分析装置及方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US5532943A (zh) |
JP (1) | JP3213097B2 (zh) |
CN (1) | CN1054203C (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7493219B1 (en) | 2007-09-13 | 2009-02-17 | Shenzhen Mindray Bio-Medical Electronics Co. | Method of discriminating particle groups and particle analyzer |
CN102033031A (zh) * | 2009-09-29 | 2011-04-27 | 希森美康株式会社 | 粒子分析仪 |
Families Citing this family (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3504029B2 (ja) * | 1995-07-04 | 2004-03-08 | シスメックス株式会社 | 粒子分析装置 |
FR2743143B1 (fr) * | 1995-12-28 | 1998-02-27 | Elf Antar France | Procede de determination de la valeur d'une grandeur physique |
JP4136017B2 (ja) * | 1996-09-19 | 2008-08-20 | シスメックス株式会社 | 粒子分析装置 |
JP3620946B2 (ja) * | 1997-03-17 | 2005-02-16 | シスメックス株式会社 | スキャッタグラムの表示方法および粒子計測装置 |
US6122599A (en) | 1998-02-13 | 2000-09-19 | Mehta; Shailesh | Apparatus and method for analyzing particles |
JP4708605B2 (ja) * | 2000-07-24 | 2011-06-22 | シスメックス株式会社 | 粒子分析装置とその粒子分画方法 |
US6535836B1 (en) * | 2000-09-29 | 2003-03-18 | Coulter International Corp. | Method for the analysis of abnormal particle populations |
JP3572319B2 (ja) * | 2001-11-15 | 2004-09-29 | 独立行政法人理化学研究所 | 液体中微粒子分析装置 |
US20070250548A1 (en) * | 2006-04-21 | 2007-10-25 | Beckman Coulter, Inc. | Systems and methods for displaying a cellular abnormality |
JP5025371B2 (ja) | 2007-07-31 | 2012-09-12 | シスメックス株式会社 | 血液分析装置 |
JP4580976B2 (ja) * | 2007-12-26 | 2010-11-17 | シスメックス株式会社 | 精度管理機能を備えた生体試料分析装置および精度管理測定結果の表示方法 |
JP2009222610A (ja) * | 2008-03-18 | 2009-10-01 | Hitachi High-Technologies Corp | 自動分析装置 |
CN101672759B (zh) * | 2008-09-12 | 2013-01-23 | 深圳迈瑞生物医疗电子股份有限公司 | 一种粒子分类统计方法及装置 |
FI20105645A0 (fi) * | 2010-06-07 | 2010-06-07 | Environics Oy | Laite ja menetelmä biologisen materiaalin havaitsemiseksi |
JP5892967B2 (ja) | 2013-03-29 | 2016-03-23 | シスメックス株式会社 | 細胞分析装置、細胞測定装置の管理方法およびコンピュータプログラム |
KR20230156814A (ko) | 2017-10-26 | 2023-11-14 | 파티클 머슈어링 시스템즈, 인크. | 입자 측정을 위한 시스템 및 방법 |
JP6906477B2 (ja) * | 2018-05-22 | 2021-07-21 | 株式会社日立製作所 | データ処理装置およびデータ処理方法 |
WO2020219841A1 (en) | 2019-04-25 | 2020-10-29 | Particle Measuring Systems, Inc. | Particle detection systems and methods for on-axis particle detection and/or differential detection |
EP4062146A4 (en) | 2019-11-22 | 2024-02-28 | Particle Measuring Syst | ADVANCED SYSTEMS AND METHODS FOR INTERFEROMETRIC PARTICLE DETECTION AND DETECTION OF PARTICLES WITH SMALL DIMENSIONS |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4661913A (en) * | 1984-09-11 | 1987-04-28 | Becton, Dickinson And Company | Apparatus and method for the detection and classification of articles using flow cytometry techniques |
US4706207A (en) * | 1985-06-24 | 1987-11-10 | Nova Celltrak, Inc. | Count accuracy control means for a blood analyses system |
US4741043B1 (en) * | 1985-11-04 | 1994-08-09 | Cell Analysis Systems Inc | Method of and apparatus for image analyses of biological specimens |
JPH03131756A (ja) * | 1989-10-18 | 1991-06-05 | Hitachi Ltd | 血球自動分類装置 |
JP2643596B2 (ja) * | 1990-11-29 | 1997-08-20 | 株式会社日立製作所 | スカラ量分布表示方法 |
US5187673A (en) * | 1991-02-05 | 1993-02-16 | Edward L. Carver, Jr. | Method and apparatus for determining the distribution of constituent subpopulations within a population of particles having overlapping subpopulations |
-
1992
- 1992-12-28 JP JP35989392A patent/JP3213097B2/ja not_active Expired - Lifetime
-
1993
- 1993-12-25 CN CN93121213A patent/CN1054203C/zh not_active Expired - Fee Related
- 1993-12-28 US US08/174,249 patent/US5532943A/en not_active Expired - Lifetime
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7493219B1 (en) | 2007-09-13 | 2009-02-17 | Shenzhen Mindray Bio-Medical Electronics Co. | Method of discriminating particle groups and particle analyzer |
US8055471B2 (en) | 2007-09-13 | 2011-11-08 | Shenzhen Mindray Bio-Medical Electronics Co., Ltd. | Method of discriminating particle groups and particle analyzer |
CN102033031A (zh) * | 2009-09-29 | 2011-04-27 | 希森美康株式会社 | 粒子分析仪 |
Also Published As
Publication number | Publication date |
---|---|
JPH06201569A (ja) | 1994-07-19 |
US5532943A (en) | 1996-07-02 |
JP3213097B2 (ja) | 2001-09-25 |
CN1054203C (zh) | 2000-07-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1091830A (zh) | 粒子分析装置及方法 | |
Kim et al. | Demystifying “drop-outs” in single-cell UMI data | |
Raychaudhuri et al. | Principal components analysis to summarize microarray experiments: application to sporulation time series | |
CN1221911C (zh) | 用于组织数据的系统与方法 | |
Seita et al. | Gene Expression Commons: an open platform for absolute gene expression profiling | |
CN1615434A (zh) | 适用于处理电化学信号的方法和装置 | |
CN101069184A (zh) | 信息处理装置和方法、及程序 | |
CN101059824A (zh) | 支持基因相互作用网络分析的方法和设备以及计算机产品 | |
Jones et al. | Adjusting for cell type composition in DNA methylation data using a regression-based approach | |
CN1574269A (zh) | 不合格检测方法和不合格检测装置 | |
CN1961212A (zh) | 诊断系统性自身免疫疾病的模式识别法 | |
CN1940949A (zh) | 通过利用可变聚类端点的聚类分析的ct确定 | |
CN108664766B (zh) | 拷贝数变异的分析方法、分析装置、设备及存储介质 | |
CN1319212A (zh) | 用于统计过程控制的多维方法和系统 | |
CN1298486A (zh) | 评价化学和生物学分析的方法 | |
CN100351610C (zh) | 状态判断方法 | |
CN1577328A (zh) | 基于视觉的文档分割 | |
CN1871563A (zh) | 加工信息产生装置、程序和加工信息产生方法 | |
CN1838107A (zh) | 处理系列波谱尤其是nmr波谱的方法 | |
CN1923140A (zh) | 图像测量装置、方法以及肾小球滤过率的图像测量系统 | |
US20010049689A1 (en) | Molecular database for antibody characterization | |
Evin et al. | Optimizing digitalization effort in morphometrics | |
Sanderson et al. | Reconstructing past admixture processes from local genomic ancestry using wavelet transformation | |
CN101067833A (zh) | 将临床医学上各种检测或化验结果的正常范围参考值及其实际测量值统一标化的方法 | |
TW201248425A (en) | Comprehensive glaucoma determination method utilizing glaucoma diagnosis chip and deformed proteomics cluster analysis |
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 | ||
C19 | Lapse of patent right due to non-payment of the annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |