血液凝固时间测定方法

技术领域

本发明涉及血液凝固时间测定方法。

背景技术

血液凝固检查是用于通过向患者的血液检体添加规定的试剂来测定血液凝固时间等，从而诊断患者的血液凝固能力的检查。作为血液凝固时间的典型的例子，存在凝血酶原时间(PT)、活化部分凝血活酶时间(APTT)、凝血酶时间等。通过血液凝固检查，能够调查患者的止血能力、纤溶能力。血液凝固能力的异常主要引起凝固时间的延长。例如，凝固时间的延长的原因在于抗凝药剂的影响、参与凝固的成分的减少、先天性血液凝固因子的缺乏、抑制后天性凝固反应的自身抗体等。

近年，进行血液凝固检查的自动测量的自动分析装置是通用的，能够便于实施血液凝固检查。例如，在某种自动分析装置中，对向血液检体添加试剂而得的混合液照射光，并基于所得到的散射光量的变化来测量该血液检体的凝固反应。在通常的血液凝固反应中，在从添加试剂起经过了某种程度的时间的时刻，散射光量因凝固的进行而急剧上升，其后凝固反应接近结束，并且散射光量饱和，而达到平稳。能够基于这样的散射光量的时间的变化来测定血液凝固时间。作为基于自动分析装置的凝固时间的计算法，使用百分比检测法、微分法等几种方法。其中，在百分比检测法中，每单位时间的散射光变化量大时，例如通过检测最大散射光量的50％到达点，即使是低纤维蛋白原检体、乳糜检体、溶血检体等异常检体，也能够相当准确地进行凝固时间的测定。

另一方面，在分析装置中的测光数据中包含由装置、试剂、检体的状态等引起的各种噪声，这些噪声会导致凝固时间的误检测。在血液检体的自动分析中，要求去除噪声的负面影响来计算具有可靠性的凝固时间。考虑了避免由测光数据的噪声导致的凝固时间的误检测的方法。在专利文献1中记载有如下这样的血液凝固时间测定方法，即，对从分析装置实时得到的散射光量数据进行平滑化以及原点调整而作为基准数据X，根据该基准数据，进一步运算对该基准数据进行了积分后的基准积分数据Y、和作为上述基准数据的各相邻的微小时间的累计值之比的基准比数据Z，选出该基准比数据Z成为预先确定的一定的基准比数据值Zs的时刻中的、该基准比数据Z的峰值以后的时刻且该基准积分数据Y成为预先确定的阈值Ys以上的时刻的基准数据值Xd，将直到该Xd的1/N(N为1以上的一定整数)的值所对应的时刻为止的混合时刻起的时间作为凝固时间。

专利文献1：日本特开平6-249855号公报

发明内容

本发明涉及能够准确地测定表示各种血液凝固反应曲线的血液检体的凝固时间的血液凝固时间测定方法。

即，本发明提供以下内容。

＜1＞一种方法，该方法是血液凝固时间测定方法，

该方法包含：

〔1〕对关于血液检体的凝固反应的测量值P(i)进行平滑化以及零点调整，取得反应X(i)；

〔2〕取得该反应X(i)的累计比Z(i)，这里，该Z(i)是第一测量区间中的X的累计值、与和该第一区间相邻的第二测量区间中的X的累计值的比；

〔3〕计算参数te[k]、X(te[k])、tc[k]、te[k-1]、X(te[k-1])、以及tc[k-1]，

这里，

K为2以上的整数，

te[k]以及te[k-1]分别是满足Z(i)＜Zs[k]的i下的测量时间、以及满足Z(i)＜Zs[k-1]的i下的测量时间，

X(te[k])以及X(te[k-1])分别是te[k]以及te[k-1]下的反应X，

tc[k]以及tc[k-1]分别是满足X(i)＝{X(te[k])×Q％}的i下的测量时间、以及满足{X(te[k-1])×Q％}的i下的测量时间，

1＜Q＜100；

〔4〕计算指标R[k]、以及/或指标TR[k]，

这里，

R[k]＝X(te[k])/X(te[k-1]) (2)

TR[k]＝Δtc[k]/Δte[k] (5)

Δtc[k]＝tc[k]-tc[k-1] (3)

Δte[k]＝te[k]-te[k-1] (4)；以及

〔5〕在该R[k]以及/或TR[k]满足规定的条件的情况下，将该tc[k]或该tc[k-1]决定为凝固时间Tc。

＜2＞根据＜1＞记载的方法，在工序〔5〕中该R[k]或TR[k]不满足该规定的条件的情况下，设为k＝k+1，反复进行工序〔3〕～〔5〕。

＜3＞根据＜1＞或＜2＞记载的方法，Zs[k]＜Zs[k-1]，Zs[k]大于1，且Zs[1]为1.100以下。

＜4＞根据＜1＞～＜3＞中任一项记载的方法，k为10以下。

＜5＞根据＜1＞～＜4＞中任一项记载的方法，其中，Zs[k-1]与Zs[k]的差为0.050～0.001。

＜6＞根据＜1＞～＜5＞中任一项记载的方法，Z(i)由下述式表示：

Z(i)＝{X(i+1)+X(i+2)+...+X(i+m)}/{X(i-m)+X(i-m+1)+...+X(i-1)}(1)

这里，m＝10～30。

＜7＞根据＜1＞～＜6＞中任一项记载的方法，在i达到规定的计算开始点s以后进行工序〔3〕。

＜8＞根据＜7＞记载的方法，上述计算开始点s是比凝固反应的速度成为最大的时间靠后的测量点。

＜9＞根据＜7＞或＜8＞记载的方法，该方法还包括：在V(i)＝Vs的两点间的测量点数超过规定的值时，将成为V(i)＝Vs的较慢的测量点检测作为计算开始点s，这里，V(i)是上述反应X(i)的微分值。

＜10＞根据＜7＞记载的方法，上述计算开始点s是在工序〔2〕中取得的Z(i)达到峰值以后的测量点。

＜11＞根据＜1＞～＜10＞中任一项记载的方法，并行地进行工序〔1〕～〔2〕和工序〔3〕～〔5〕。

＜12＞根据＜2＞～＜11＞中任一项记载的方法，该方法还包含：通过反复进行上述工序〔3〕～〔5〕而没有决定Tc时，取得Tc未被正常地决定的结果。

根据本发明的方法，能够对表示包含正常检体、异常检体在内的各种血液凝固反应曲线的血液检体的凝固时间准确地进行测定。另外，在通过自动分析装置实时地分析大量的检体的情况，本发明的方法防止凝固时间的误检测并且使每一个检体的分析时间最佳化，使分析效率提高。

附图说明

图1是基于本发明的血液凝固时间测定方法的一实施方式的基本流程。

图2A是反应X以及微分值V；图2B是反应X以及累计比Z。图2A、图2B都是左侧为正常检体，右侧为FVIII缺乏检体。

图3是基于V的峰值宽度的计算开始点s的检测方法的概念图。

图4A是Zs[1]、和所对应的te[1]、X(te[1])、tc[1]以及X(tc[1])的例子。图4B是Zs[1]～Zs[3]、和所对应的te[1]～te[3]、X(te[1])～X(te[3])、tc[1]～tc[3]的例子。

图5是本发明的血液凝固时间测定方法的一实施方式的详细流程。

图6是关于试样1(左)以及试样10(右)的X的反应曲线。

图7是不同的Zs[k]下的tc[k](左)以及te[k](右)。

图8A～图8B是R[k](％)相对于APTT的标绘图以及其放大图。图8C～图8D是表示基于R[k](％)进行的Tc计算的准确性的标绘图以及其放大图。图8E是不同的Zs[k]下的X(te[k])/X(te[8])。

图9A是不同的Zs[k]下的X(te[k])/X(te[8])。图9B～图9D是TR[k](ΔC/ΔE)(％)、ΔE、以及ΔC相对于APTT的标绘图。

图10A～图10B是根据试样1～10计算出的tc[k]以及te[k]，其中本法：决定Tc时的值，1.001：Zs＝1.001时的值、1.031：Zs＝1.031时的值。图10C～图10D是本法或Zs＝1.031与Zs＝1.001的tc[k]以及te[k]的差。

图11是表示在测量时间内反应没有结束的APTT延长检体的反应曲线、以及所取得的参数和指标的表。

图12是表示低纤维蛋白原试样的反应曲线、以及所取得的参数和指标的表。

图13是不同的条件1～3下检测到的计算开始点s对Tc计算用参数的影响。

具体实施方式

1.血液凝固反应测量

在血液凝固检查中，向血液检体添加规定的试剂，对其后的血液凝固反应进行测量，并根据凝固反应来测定血液凝固时间。凝固反应的测量使用一般的手段，例如测量散射光量、透过度、吸光度等的光学手段、或测量血浆的粘度的力学手段等。正常的血液检体的凝固反应曲线取决于测量手段，但基本上显示S形。例如，基于正常检体的散射光量的凝固反应曲线通常在从试剂添加起经过了某种程度的时间的时刻因凝固的进行而急剧上升，其后，凝固反应接近结束并且达到平稳。另一方面，具有凝固异常的异常检体的凝固反应曲线取决于曲线的升高时间的延慢、缓慢的上升等异常原因而显示各种形状。异常检体的凝固反应曲线的多样性使得自动分析装置中的凝固时间的准确的测定变得困难。

在以往的一般的血液凝固时间测定中，至少取得到凝固反应结束为止的数据，基于所取得的数据来计算凝固时间。例如，存在如下等方法，即，在将散射光量饱和的时刻判定为凝固反应结束后，将从试剂添加时刻到凝固反应结束时刻为止的期间、凝固反应曲线达到最大速度的时刻决定为凝固时间的方法(微分法)、和将达到凝固反应结束时刻的散射光量的1/Q的时刻决定为凝固时间的方法(百分比检测法，参照专利文献1)等。然而，因上述那样的异常检体的凝固反应曲线的异常的形状、噪声，会引起凝固反应结束的误检测，例如有时在过早时刻检测到凝固反应结束、或未检测到凝固反应结束。这样的凝固反应结束的误检测的结果是，计算出不准确的凝固时间。

在自动分析装置中，为了高效地分析大量的检体，而希望当判断为一个检体的凝固反应结束时，迅速地结束测量，并开始下一个检体的测量。但是，在这样的方法中，存在上述的在过早时刻的凝固反应结束的误检测尽管在测量中却导致结束测量，并遗漏必要的数据的风险。另一方面，若将每一个检体的测定时间固定为足够长的时间，则能够防止由凝固反应结束的误检测导致的数据遗漏，但这样的方法由于对大量的检体而言测定时间长至必要以上，使整体的分析效率降低。

例如，在用专利文献1所记载的方法测定第VIII因子(FVIII)缺乏检体那样的凝固反应曲线的变化小的检体的凝固时间的情况下，基准比数据值Zs的适当的值比较小(例如1.01)，这意味着在该检体的凝固时间测定中需要比较长时间的凝固反应测量。另一方面，适合于正常检体的Zs是更大的值(1.05左右)，这意味着在正常检体中能够以更短的凝固反应测量进行凝固时间的计算。因此，在使用自动分析装置并按照专利文献1所记载的方法，进行包含正常检体、异常检体在内的所有的检体的凝固时间测定的情况下，必须将用于一部分的异常检体的比较长的测量时间也应用于大量的正常检体。这不能说是高效率的分析方法。

本发明的血液凝固时间测定方法能够防止由上述那样的凝固反应曲线的异常的形状引起的凝固时间的误检测，而进行准确的凝固时间测定。另外，根据本发明，能够使测量时间最佳化，使得对包含正常检体、异常检体在内的各种血液检体，应用各自的凝固时间测定所需的最低限度的凝固反应测量时间。

在本发明的血液凝固时间测定方法(以下，也称为本发明的方法)中，对混合了试剂的被检血液检体的凝固反应进行测量。基于由该测量而得到的凝固反应的时间序列数据，来测定血液凝固时间。作为本发明的方法中所测定的血液凝固时间的例子，可列举凝血酶原时间(PT)、活化部分凝血活酶时间(APTT)、纤维蛋白原浓度(Fbg)测定中的凝固时间等。

在本发明的方法中，作为被检血液检体，优选使用被检者的血浆。在以下的本说明书中，也将血液检体简称为检体。在该检体中可添加凝固检查通常使用的抗凝固剂。例如，在使用加入柠檬酸钠的采血管进行采血后，通过离心分离而得到血浆。

向该被检检体添加凝固时间测定试剂，使血液凝固反应开始。可测量试剂添加后的混合液的凝固反应。所使用的凝固时间测定试剂能够根据测定目标而任意地选择。用于各种凝固时间测定的试剂在市场中出售(例如，APTT测定试剂Coagpia APTT-N；积水医疗株式会社制)。凝固反应的测量使用一般的手段，例如测量散射光量、透过度、吸光度等的光学手段、或测量血浆的粘度的力学手段等即可。凝固反应的反应开始时刻典型地可以定义为向检体混合试剂而使凝固反应开始的时刻，但也可以将其他时机定义为反应开始时刻。持续凝固反应的测量的时间例如可以是从检体与试剂的混合的时刻起数十秒～7分钟左右。该测量时间也可以是任意确定的固定的值，但也可以设为直到检测到各检体的凝固反应的结束的时刻为止。该测量时间期间，可以以规定的间隔反复进行凝固反应的进行状况的测量(在光学检测的情况下为测光)。例如，可以以0.1秒间隔进行测量。该测量中的混合液的温度在通常的条件下例如为30℃以上40℃以下，优选为35℃以上39℃以下。另外，测量的各种条件可根据被检检体、试剂、测量手段等而适当设定。

上述凝固反应测量中的一系列的操作能够使用自动分析装置来进行。作为自动分析装置的一个例子，可列举血液凝固自动分析装置CP3000(积水医疗株式会社制)。或者，也可以通过人工作业来进行一部分的操作。例如，由人进行被检检体的调制，此后的操作通过自动分析装置进行。

1.血液凝固时间的测定

以下，参照图1所示的基于本发明的血液凝固时间测定方法的一实施方式的基本流程，对本发明进行说明。

1.1反应X、微分值V、累计比Z的取得

在本发明的方法中，通过分析装置，逐次取得针对从反应开始时刻起的凝固反应的测量值P(i)(例如散射光量的测光值)。这里“i”表示测量点编号。例如，若测量(测光)间隔为0.1秒，则P(i)表示测量开始0.1×i秒后的测量值。

由于在该测量值P(i)中包含测光时的噪声、与紧接在测光开始之后出现的反应无关的变动，所以用公知的方法对该测量值进行平滑化处理。另外，在用散射光量对凝固反应进行测光的情况下，进行减去出自反应前的反应液的散射光量的零点调整处理。因此，所取得的时间序列的测量值P(i)被逐次平滑化以及零点调整，而取得反应X(i)(工序S1)。在测量值的平滑化处理中，可使用与噪声去除相关的各种公知的方法中的任一种。例如，作为平滑化处理，可列举滤波处理、或在通过差值、后述的区间内平均斜率的运算等而求出微分值后对其进行积分的处理等。在零点调整中，例如只要以测量开始时刻的值成为0的方式调整平滑化后的测量值即可。并且，也可以对测量值P(i)进行初始变动去除处理。初始变动去除处理只要使得从测光开始时刻到预先确定的初始变动去除时间为止的所有的值成为0即可。基本上，如图2所示，反应X(i)构成进行了平滑化以及零点调整的凝固反应曲线。

接着，根据求出的反应X(i)，取得该反应X(i)的累计比Z(i)(工序S2)。并且，也可以取得该反应X(i)的微分值V(i)。微分值V(i)能够用于后述的计算开始点s的检测。

微分值V(i)是对反应X(i)进行一次微分而得到的。微分处理能够通过任意的方法实施，例如可通过区间内平均斜率值的计算来进行。在区间内平均斜率值的计算中，可利用各测量点i的前后的一定数量的测量点，例如从i-K到i+K为止的2K+1个测量点。例如，可利用第i-2、i-1、i、i+1、i+2这5个点的测量点。平均斜率值是指对这些多个测量点进行了直线近似时的斜率值。直线近似的运算方法可利用最小二乘法等常规方法。这些测量点的平均斜率值可视为在测量点i的一次微分值。

累计比Z(i)是在第一测量区间的X的累计值、与和该第一区间相邻的第二测量区间中的反应X的累计值的比。更详细而言，Z(i)是比测量点i靠前的第一测量区间中的反应X(例如在测量点[i-m]～[i-1]的X(i-m)～X(i-1))的累计值、与比该测量点i靠后的第二测量区间中的反应X(例如在测量点[i+1]～[i+m]的X(i+1)～X(i+m))的累计值之比。这里，第一测量区间与该第二测量区间优选为相同的长度。即，在将测量点i的反应X设为X(i)，将m设为常数时，测量点i处的累计比Z(i)通过下式(1)计算。

Z(i)＝{X(i+1)+X(i+2)+...+X(i+m)}/{X(i-m)+X(i-m+1)+...+X(i-1)} (1)

在式(1)中，m优选设定为该微小测量区间的时间宽度成为1秒至3秒。即，若测量(测光)时间间隔为0.1秒，则m优选为10至30，若测量(测光)时间间隔为0.2秒，则m优选为5至15。只要考虑测量条件来选择m的最佳值，以便能够计算适当的凝固时间即可。

因此，在作为关于X的第2m+1个数据，取得测量点i+m处的X(i+m)以后，能够执行累计比Z(i)的计算。但是，为了避免在反应初始的X的变动的影响，优选将从X(i-m)到X(i-1)为止全部为规定的阈值Xs以上作为条件，开始Z(i)的计算。由于基于散射光测量的反应曲线成为S形，所以通常当X(i)一旦超过Xs时，则其后也成为X(i)＞Xs。如果在成为X(i)＞Xs后成为X(i)≤Xs的情况下，数据可被视为异常。Xs优选为预测的反应X的最大值的10％至50％。例如，在估算反应X的最大值为约700的情况下，Xs优选为70～350。

通过以上过程，能够求出各测量点i的反应X(i)、微分值V(i)、以及累计比Z(i)。参照图2，对关于APTT为正常的正常检体、和作为APTT延长的例子的FVIII缺乏检体(以下，称为延长检体)的反应X、微分值V、以及累计比Z的每一个的形状的相对的差异进行说明。图2A表示反应X和微分值V，图2B表示反应X以及累计比Z。图2A、B均是左为正常检体，右为延长检体。横轴换算为从反应开始时刻起的时间。

反应X)根据图2A所示，X构成进行了平滑化以及零点调整的凝固反应曲线。在正常检体(图2A左)中，X的升高点快，上升的斜率大。但是，在延长检体(图2A右)中，X的升高点慢，上升的斜率小。

微分值V)在正常检体(图2A左)中，V的峰值的顶部位置高，形状大致左右对称。但是，在延长检体(图2A右)中，V的峰值的顶部位置低(在该例子中，小于正常检体的六分之一)，形状为左右非对称，峰值并不是一个，而是呈双峰性形状。

累计比Z)在正常检体(图2B左)中，在反应曲线的升高时Z的最大值约为8，其后朝向1急剧减少。另一方面，在延长检体(图2B右)中Z的最大值约为2.5，减少也缓慢。无论在哪一种情况下，累计比Z都在凝固反应曲线接近平稳(血液凝固反应接近结束)的同时接近1。

1.2凝固时间Tc的计算

1.2.1计算开始点的决定

在本发明的方法中，使用上述的反应X(i)、微分值V(i)、以及累计比Z(i)的值来计算凝固时间Tc。在本发明的方法中，X、V、Z的取得可以与后述的凝固反应测量过程并行地进行。首先，基于目前为止取得的X、V、Z，来决定开始凝固时间Tc的计算处理的时机，即计算开始点s。

在一个例子中，能够将Z达到峰值以后的任意的测量点设为计算开始点s。例如，能够将Z达到阈值后，经过峰值再次达到该阈值的测量点设为计算开始点s。以图2B为例，能够将累计比Z达到2后，经过峰值再次达到2的测量点设为计算开始点s。

在另一个例子中，求出凝固反应速度成为最大值Vmax的测量点，并将该测量点设为计算开始点s。最简单来说，能够将与由逐次取得的V(i)构成的曲线中的最初的峰值顶部对应的测量点设为与Vmax对应的计算开始点s。

另一方面，在计算开始点s的决定时，希望排除在凝固因子缺乏检体中常见的凝固反应速度曲线成为双峰性的情况下(例如图2A右)可能引起的Vmax的误检测、由噪声引起的变动的影响。因此，在优选的例子中，能够将逐次取得的V(i)超过规定的值之后，经过峰值而达成规定的条件时的测量点设为计算开始点s。在图3示出有计算开始点s的决定方法的一个例子的概念图。在图3示出有V(i)、和用于计算开始点s的决定的参数。V(i)形成具有峰值Vmax的曲线。Vs是将Vmax设定为100％时的S％的值。V(i)一旦达到Vs后，经过峰值，在V(s)中再次达到Vs。此时，在峰值宽度Ws超过了规定的值时，检测该V(s)时的测量点s来作为计算开始点s，该峰值宽度Ws是V(i)＝Vs的两点间的宽度(测量点数)。

1.2.2参数取得

接下来，取得Tc计算所需的参数(工序S3)。

首先，检测在计算开始点s以后Z(i)初次达到不足预先确定的第一个阈值Zs[1](Z(i)＜Zs[1])时的测量点处的测量时间(从反应开始时刻起的时间)te[1]。根据需要，继续取得Z(i)，直到满足Z(i)＜Zs[1]为止。

若检测到te[1]，则求出在te[1]处的反应X(te[1])。接下来，检测与具有或接近X(te[1])的Q％的值的X(i)对应的测量点处的测量时间tc[1](这里tc[1]＜te[1])。例如，能够检测满足X(tc[1])＝{X(te[1])×Q％}的时间tc[1]。另外，例如能够将最接近{X(te[1])×Q％}的值设为X(tc[1])。另外，例如在X(tc’[1])＜{X(te[1])×Q％}＜X(tc”[1])(tc’[1]与tc”[1]是相邻的测量点的时间)时，能够检测tc’[1]和tc”[1]中的任一个来作为tc[1]。在测量(测光)间隔足够短的(例如0.1秒等)的情况下，最终的Tc的计算值不会产生大的差异，因此也可以选择tc’[1]和tc”[1]中的任一个。Q只要是0＜Q＜100的任意的值即可，优选为10～80，更优选为20～70(例如50)。

在图4A中示出有表示Zs[1]、和凝固反应曲线(由X(i)构成的曲线)上的te[1]及X(te[1])、和tc[1]及X(tc[1])的关系的概念图。

并且，检测在Z(i)初次达到不足预先确定第二个阈值Zs[2](Z(i)＜Zs[2])时的测量点处的测量时间te[2](这里Zs[2]＜Zs[1])。根据需要，继续取得Z(i)，直到满足Z(i)＜Zs[2]为止。在te[2]处的X为X(te[2])。接下来，检测与具有或接近X(te[2])的Q％的值的X(i)对应的测量点处的时间tc[2](此时，tc[2]＜te[2])。

上述所求出的te是Z(i)达到阈值Zs的时刻，且表示临时的凝固结束点。另一方面，tc是X的值成为临时的凝固结束点处的X的Q％的时刻，且表示临时的凝固点。图4B是表示Zs[1]～Zs[3]、以及与这些对应的反应曲线X上的临时的凝固结束点te[1]～te[3]、以及此时的X(te[1])～X(te[3])、和临时的凝固点tc[1]～tc[3]的概念图。如图4B所示，随着反应曲线接近平稳(凝固反应接近结束)，te的增加幅度变大，另一方面，X(te)的增加幅度逐渐变小。另一方面，随着凝固反应曲线接近平稳，tc也逐渐增加，但其变化与te相比较小。

1.2.3Tc的计算

在本发明的方法中，进行Tc计算步骤，该Tc计算步骤利用上述那样的随着凝固反应的进行而产生的临时的凝固结束点te以及临时的凝固点tc的变化，来决定临时的凝固点tc是否为真正的凝固时间Tc。即，在本发明的方法中，利用出自上述所求出的Zs的参数te、tc、X(te)等，取得用于Tc计算的指标(工序S4)。

因此，在本发明的方法中，检测Z(i)初次达到不足第k个阈值Zs[k](Z(i)＜Zs[k])时的测量点(这里为Zs[k]＜Zs[k-1]、k≥2的整数)。根据需要，继续取得Z(i)，直到满足Z(i)＜Zs[k]为止。接下来，取得出自Zs[k]的参数(例如te[k]、tc[k]、X(te[k])等)以及出自在其之前取得的第k-1个阈值Zs[k-1](这里Zs[k]＜Zs[k-1])的参数(例如te[k-1]、tc[k-1]、X(te[k-1])等)。

优选该指标可以基于te[k]、tc[k]以及X(te[k])中的任一个以上、和te[k-1]、tc[k-1]以及X(te[k-1])中的任一个以上而取得。

在Tc计算步骤中，在该指标满足预先确定的基准的情况下，将临时的凝固点tc[k]或tc[k-1]决定为凝固时间Tc。

因此，该指标的取得可以在取得出自第二个阈值Zs[2]的te[2]、tc[2]、X(te[2])以后进行。优选为在k＝2时，最初的指标可基于出自Zs[1]的te[1]、tc[1]或X(te[1])、和出自Zs[2]的te[2]、tc[2]或X(te[2])而取得。

在一实施方式中，该指标是出自第k个阈值Zs[k]的X(te[k])、与出自在其之前的阈值Zs[k-1]的X(te[k-1])的差或比。优选为，按照下述式(2)得出的R[k]被用作用于Tc计算的指标。

R[k]＝X(te[k])/X(te[k-1]) (2)

该R[k]可以用比率表示，也可以用百分率表示。若R[k]是预先确定的基准值Rs或比其小，则能够将该tc[k]或该tc[k-1]决定为凝固时间Tc。此时该te[k]或该te[k-1]能够被决定为真正的凝固结束点。

在另一实施方式中，首先，按照以下式(3)～(5)，计算Δtc[k]、Δte[k]以及TR[k]。该TR[k]被用作用于Tc计算的指标。

Δtc[k]＝tc[k]-tc[k-1] (3)

Δte[k]＝te[k]-te[k-1] (4)

TR[k]＝Δtc[k]/Δte[k] (5)

该TR[k]可以用比率表示，也可以用百分率表示。若TR[k]是预先确定的基准值TRs或比其小，则能够将该tc[k]或该tc[k-1]决定为凝固时间Tc。此时该te[k]或该te[k-1]能够被决定为真正的凝固结束点。

基准值Rs以及TRs能够适当设定，以使根据本发明而决定的Tc的准确性成为所希望的值。例如，优选调整Rs以及TRs，使得对于各种检体，通过本发明的方法决定的Tc成为实测的凝固时间的90％以上，优选为95％以上。实测的凝固时间能够从目测的凝固反应曲线的凝固反应结束点推断。或者，也可以将Zs接近1的情况下(例如Zs＝1.001)得到的Tc推断为实测的凝固时间。例如，Rs(％)优选为110％以下，更优选为105％以下，TRs(％)优选为15％以下，更优选为10％以下。

在优选的实施方式中，在上述两个基准都满足的情况下，例如在R[k]≤Rs、且TR[k]≤TRs时，能够将该tc[k]或该tc[k-1]决定为凝固时间Tc。

另一方面，在该指标不满足预先确定的基准的情况下，设为k＝k+1，并反复进行上述过程。边使k递增边反复进行该过程，直到取得满足基准的指标为止。

以下对指标的取得和Tc决定的过程的具体例进行说明。优选为，首先，使用与Zs[1]相关的参数te[1]、tc[1]、或X(te[1])、和与Zs[2]相关的参数te[2]、tc[2]、或X(te[2])来计算用于Tc计算的指标。

在一实施方式中，该指标是X(te[1])与X(te[2])的差或比。优选为，按照下述式得出的R[2]被用作用于Tc计算的指标。

R[2](％)＝{X(te[2])/X(te[1])}(％) (2)’

在优选的实施方式中，首先，按照以下式，计算Δtc[2]、Δte[2]、以及TR[2]。该TR[2]被用作用于Tc计算的指标。

Δtc[2]＝tc[2]-tc[1] (3)’

Δte[2]＝te[2]-te[1] (4)’

TR[2](％)＝{Δtc[2]/Δte[2]}(％) (5)’

在该指标满足预先确定的基准的情况下，将tc[1]或tc[2]决定为凝固时间Tc。优选为将tc[2]决定为Tc。例如，若上述R[2]是预先确定的基准值Rs或比其小，则将该tc[2]或该tc[1]决定为Tc。另外，例如，若上述TR[2]是预先确定的基准值TRs或比其小，则将该tc[2]或该tc[1]决定为Tc。优选为，在上述两个基准均满足的情况下，例如在R[2]≤Rs、且TR[2]≤TRs时，将该tc[2]或该tc[1]决定为Tc。

另一方面，在该指标不满足预先确定的基准的情况下，针对第三个阈值Zs[3](这里Zs[3])＜Zs[2])，与上述同样地求出te[3]、tc[3]或X(te[3])，并根据这些值和te[2]、tc[2]或X(te[2])，按同样的过程计算指标。

例如，在R[2]大于Rs、且TR[2]大于TRs的情况下，针对第三个阈值Zs[3]，与上述同样地求出X(te[3])，接下来根据上述式(2)计算R[3]、以及/或求出te[3]以及tc[3]，接下来根据上述式(3)～(5)计算TR[3]。

若计算出的指标满足基准，则将该tc[2]或该tc[3]决定为Tc。优选为在R[3]≤Rs、以及/或TR[3]≤TRs的情况下，将该tc[3]或该tc[2]决定为Tc。

在指标不满足基准的情况下，进一步针对第四个阈值Zs[4](这里Zs[4]＜Zs[3])反复进行相同的过程。

这样，一边并行地逐次进行凝固反应的测量和X、Y、Z的取得，直到决定Tc为止，一边针对第N个阈值Zs[N](这里Zs[N]＜Zs[N-1]；N≥2)反复进行相同的过程。最终，该tc[N]或该tc[N-1]被决定为Tc。此外，若将Zs[N]与Zs[N-1]的差设定在适当的范围，则满足R[N]≤Rs以及/或TR[N]≤TRs时的tc[N]和tc(N-1)是表示比较接近的值，因此对于Tc，无论选择tc[N]，还是选择该tc(N-1)，在使用Tc的临床的判断中都不会产生实际的问题。但是，从提高Tc的计算值的可靠性的观点出发，优选对Tc选择tc[N]。

本发明的方法基本上是应用百分比法，该百分比法是以往用作血液凝固时间的一般的计算法。在百分比法中，将凝固反应进行到反应结束时的Q％(作为一个例子为50％)为止的时刻决定为凝固时间Tc。因此，在本发明的方法中，理想的也是凝固反应来到累计比Z为1或接近1时(血液凝固反应大致结束时)的Q％的时刻作为凝固时间Tc来进行计算。若计算条件适当，则te[k]成为S状的凝固反应曲线接近平稳之前的时刻，因此通过本发明的方法而决定的凝固时间Tc不会与理想值相差较大。另外，在百分比法中，在将凝固反应从反应结束时的50％偏离的比率C％(作为一个例子，为45％以下或55％以上)设为凝固时间的情况下，也可以在求出作为Q％＝50％的50％凝固时间后，基于相同的反应结束时刻来计算C％时的凝固时间。

1.2.4阈值Zs

如图2所示，Z(i)在由X(i)构成的凝固反应曲线接近平稳的同时接近1。因此，Zs[k]是大于1的值。在Zs[k]达到预先设定的最低值Zsmin时，结束上述过程的反复进行，将基于该Zs[k]计算出的该tc[k]或该tc[k-1](优选为tc[k])决定为凝固时间Tc。Zsmin只要大于1即可，只要根据测定条件(例如试剂、分析装置以及分析项目)适当设定即可。例如在APTT测定的情况下，Zsmin优选为1.0001～1.05，更优选为1.0002～1.01，进一步优选为1.0005～1.002(例如1.001)。

或者，在递增的整数k达到了预先设定的最大数N时，结束上述过程的反复进行，将基于Zs[N]计算出的tc[N]或tc[N-1](优选为tc[N])决定为凝固时间Tc。从凝固时间计算的效率的观点出发，N优选为20以下，更优选为10以下，进一步优选为8以下。

第一个阈值Zs[1]只要根据测定条件(例如试剂、分析装置以及分析项目)适当设定即可。例如在APTT测定的情况下，优选为1.100以下，更优选为1.080以下，进一步优选为1.050以下。例如，在FVIII缺乏检体的APTT测定中，凝固反应速度成为最大时的累计比Z大多为1.03至1.04，因此对于Zs[1]，选择约1.04即可。

阈值Zs的刻度宽度(Zs[k-1]与Zs[k]的差；ΔZs))可取决于Zs[1]、Zsmin、k的最大数N、或凝固反应的测量频度(取样周期)等而决定，但从凝固时间计算的精度以及效率的观点出发，优选为0.050～0.001的范围，更优选为0.020～0.002的范围，进一步优选为0.010～0.005的范围。

优选以Zs[1]、Zsmin、k的最大数N、以及ΔZs分别成为上述的适当的范围的方式，设定各个值。例如，能够决定优选的Zsmin、ΔZs、以及N，并与其对应地决定Zs[1]。或者，能够决定优选的Zs[1]、Zsmin、以及N，并与其对应地决定ΔZs。

例如，在Zsmin＝1.001、k≤8、ΔZs＝0.005的情况下，Zs[1]＝1.036。但是，本发明的方法中的Zs[1]、Zsmin、ΔZs、以及重复次数的组合的例子并不限定于这些，本领域技术人员能够适当设定适当的值。

1.2.5错误处理

在通过反复进行上述那样的过程而没有发现满足Z(i)＜Zs[k]的te[k]、或指标不满足基准的情况下，上述Tc计算步骤在没有决定Tc的状态下不完全地结束。在这样的情况下，在本发明的方法中，也可以将最后取得的tc[k]决定为Tc，或能够视为Tc不可计算。在任一种情况下，在没有决定Tc时，都能够取得Tc未被正常地决定的结果。

2.Tc决定的例示的过程

以下，参照图5所示的详细流程，对基于本发明进行的凝固时间Tc决定的更具体的过程的例子更详细地进行说明。图5所示的过程是例示，能够由本领域技术人员进行通常的变更，例如运算的顺序的更换等。

S001：如下述表1那样，预先阶段性地以从高的值到低的值的方式设定Zs，以能够通过变量k(计数器)依次选择各Zs的方式制成数组Zs[k]。k为整数，设定1，来作为初始值。

[表1]

S002：在与凝固反应相关的测量值P(i)(i表示测量点编号)的数据数量超过平滑化处理、初始变动去除处理以及零点调整处理所需的数据数量(例如50个)以后，开始反应X(i)、微分值V(i)、以及累计比Z(i)的取得。

S003：边继续取得测光数据、和X(i)V(i)以及Z(i)，边搜索在计算开始点s以后满足Z(i)＜Zs[1]的最早的i下的时间te[1]。根据搜索到的te[1]，如上述1.2.2中说明的那样，求出X(te[1])、tc[1]等参数。

S004：在S003中，不可取得满足Z(i)＜Zs[1]的i下的时间te[1]的情况下，在i达到上限值为止时分支到S020，判定为“未检测到反应”，凝固时间不可计算。例如，在以0.1秒间隔将测量(测光)时间设定为360秒的情况下，最大测量点数为3600，这是i的上限值。另一方面，在能够取得te[1]的情况下分支到S005。

S005：对k进行向上计数。

S006：在前处理为S005的情况下k＝2，在前处理为S010的情况下，k为3以上的整数。接下来，进一步边继续取得P(i)、X(i)、V(i)以及Z(i)，边搜索满足Z(i)＜Zs[k]的i下的时间te[k]。

S007：在S006中，在不可取得te[k]的情况下，在i达到上限值为止时分支到S030，判断为“反应中途”，不能计算出凝固时间。i的上限值与S004相同。在能够取得的情况下分支到S008。

S008：在S006中能够取得te[k]的情况下，如上述1.2.2中说明的那样，求出X(te[k])、tc[k]等参数。根据满足Z(i)＜Zs[k-1]的te[k-1]求出X(te[k-1])、tc[k-1]等参数，计算用于Tc决定的指标。作为该指标，例如可列举上述式(2)所示的R[k]、以及上述式(5)所示的TR[k]等。R[k]和TR[k]只要任一个用作指标即可，也可以将双方用作指标。

接下来，调查该指标是否满足基准。作为基准，可列举下述(a)或(b)，优选要求满足下述(a)以及(b)双方。

基准(a)：R[k]≤Rs

基准(b)：TR[k]≤TRs

在满足上述基准的情况下，分支到S040，在其他情况下分支到S009。

S009：在S008中基准(a)和(b)均不成立的情况下，将k向上计数。

S010：在向上计数的k超过预先设定的最大值时分支到S011。K为最大值以下时分支到S006，将k向上计数并继续进行新的te[k]的搜索。

S011：在S008中基准(a)和(b)不成立，而k超过最大值，因此Tc决定过程不完全结束。在进入该步骤的情况下，也可以不输出Tc。在输出Tc的情况下，将最后取得的tc(即tc[k-1])决定为Tc。Tc的输出可以自动输出，也可以用按钮操作等通过手动输出。

S020：在未发现满足Z(i)＜Zs[1]的te[1]而结束凝固反应测量的情况下，Tc不可计算。在进入该步骤的情况下，能够判定为因某种理由而没有凝固反应，因此也可以输出“未检测到反应”标志。

S030：未发现满足Z(i)＜Zs[k]的te[k]，因此Tc不可计算。在进入该步骤的情况下，能够判定为在凝固反应的中途阶段结束测量，因此也可以输出“反应中途”标志。

S040：能够确认凝固结束点，并且将tc[k-1]或tc[k]作为Tc输出。优选将tc(k)作为Tc输出。也可以没有标志、或输出“正常结束”标志。

3.针对其他凝固反应测量法的应用

本发明的方法基本上是使用基于散射光量的血液凝固反应测量的血液凝固时间测定法，但本领域技术人员可以将本发明的方法应用于其他测量方法，例如使用基于透过度、吸光度等的血液凝固反应测量的血液凝固时间测定法。

实施例

以下，列举实施例，对本发明进一步详细地进行说明，但本发明并不限定于这些实施例。

实施例1

1.方法

对将出自存在因血液凝固因子而引起凝固反应异常的受检者的血液检体(异常检体)、和正常的血液检体(正常检体)以分别不同的比例混合而得到的多个混合试样，实施了本发明的方法根据凝固时间的测定。

1)试样

作为异常检体，使用第VIII因子(FVIII)浓度为0.1％以下的Factor VIIIDeficient Plasma(George King Bio-Medical，Inc.制)(以下，称为FVIII缺乏血浆)。作为正常检体，使用FVIII浓度被视为100％的正常混合血浆(以下，称为正常血浆)。改变FVIII缺乏血浆与正常血浆的混合比率而调制FVIII浓度分别为50％、25％、10％、5％、2.5％、1％、0.75％、0.5％、0.25％的混合血浆(依次为试样1～试样9)。另外，仅将FVIII缺乏血浆(FVIII浓度0.1％以下)作为试样10。

2)试剂

作为APTT测定用试剂，使用Coagpia APTT-N(积水医疗株式会社制)。作为氯化钙液，使用Coagpia APTT-N氯化钙液(积水医疗株式会社制)。

3)凝固反应测量

凝固反应测量使用血液凝固自动分析装置CP3000(积水医疗株式会社制)来进行。将试样50μL分注到反应杯(反应容器)后，以37℃加热45秒，接着向反应杯添加加热至约37℃的APTT测定用试剂50μL，进一步经过171后添加25mM氯化钙液50μL而使凝固反应开始。反应在维持为约37℃的状态下进行。凝固反应的测量(测光)通过向反应杯照射将波长660nm的LED等作为光源的光，以0.1秒间隔对90度侧方散射光的散射光量进行测光而进行。最大测量时间为360秒(数据数量3600个，0.1秒间隔)。

4)凝固反应曲线的制作

根据测光数据求出测光曲线P，在对该P进行了包含噪声去除在内的平滑化处理后，进行零点调整处理以使测光开始时刻的散射光量成为0，从而计算出反应X。

5)微分值V、累计比Z的取得

根据反应X，按照一次微分值V以及下述式(1)计算出累计比Z。

Z(i)＝{X(i+1)+X(i+2)+...+X(i+m)}/{X(i-m)+X(i-m+1)+...+X(i-1)} (1)

m＝20

以X(i-m)到X(i-1)为止全部Xs在以上为条件而开始Z的计算。Xs设定为350。

6)计算开始点

计算开始点s设为在Z(i)超过规定值后，低于该规定值的点。该规定值设为1.2。

7)阈值Zs的设定

按照下述表2设定了Zs[k]。

[表2]

8)参数

te[k]、tc[k]、X(te[k])、te[k-1]、tc[k-1]、X(te[k-1])(k＝2～8的整数)。

9)凝固时间Tc计算用指标

指标(a)：R[k](％)(＝X(te[k])/X(te[k-1])×100)

指标(b)：TR[k](％)(＝ΔC/ΔE×100)

ΔC＝tc[k]-tc[k-1]

ΔE＝te[k]-te[k-1]

10)计算工序

在根据测光数据而逐次取得反应X、累计比Z以及根据需要而取得的微分值V，并在计算开始点s以后，使用表1所示的Zs[k]来计算参数：te[k]、tc[k]、X(te[k])、te[k-1]、tc[k-1]、以及X(te[k-1])，进一步计算指标(a)：R[k](％)以及指标(b)：TR[k](％)。将满足下述基准(a)以及(b)时的tc[k]设为凝固时间Tc。K设为初始值k＝1，在每次增加1最大增加到8为止的同时反复进行上述工序，直到满足下述基准(a)以及(b)为止。在k＝8且不满足该基准的情况下，将tc[8]设为Tc。

基准(a)：R[k](％)≤Rs，Rs＝105％

基准(b)：TR[k](％)≤TRs，TRs＝10％

此外，在本实施例以及以下的实施例中，为了评价所决定的Tc，在满足Tc决定基准(同时满足了下述基准(a)和基准(b))后，仍继续凝固反应测量直到最大测量时间为止，从而取得了参数和指标。

2.Tc计算基准的评价

1)临时凝固结束点te、临时凝固点tc的时间的变化

在图6示出有相对于测量时间的试样1(左，FVIII浓度50％)以及试样10(右，FVIII浓度0.1％以下)的X的反应曲线。对于曲线上的标识符，白色方块表示临时凝固结束点te[k]，白色圆圈表示临时凝固点tc[k](k＝2～8，Zs[k]＝1.031～1.001)。对于试样1、10的任意一个，tc[2]和te[2]最小，tc[8]和te[8]最大。对于试样1、10的任意一个，临时凝固结束点te[8](Zs[8]＝1.001)是能够在视觉上判定为真正的凝固反应结束点的值。tc[k]的值在试样1中大致一致，另一方面，在试样10中不同。

将满足Z(i)＜Zs[8](＝1.001)的时刻设为凝固结束点，将X为凝固结束点的50％的时刻设为APTT。在图7示出有关于试样1～10的Zs[k](k＝1～8，Zs[k]＝1.036～1.001)的tc[k](左)以及te[k](右)相对于APTT的标绘图。APTT越延长，临时凝固点tc[k]的分布宽度越宽。临时凝固结束点te[k]也确认到同样的趋势。其中，最慢的点te[8]处于远离其他点的位置。

2)关于指标(a)的Tc决定基准

图8A、B是针对试样1～10求出的指标(a)：R[k](％)相对于APTT的标绘图(A)以及其放大图(B)。R[k]即使出自相同的Zs值，其分布也根据APTT而不同。由图8B可知，通过将用于Tc决定的基准值Rs设为105％，从而在约96％的情况下成为R[k]≤Rs，能够准确地计算凝固时间Tc(APTT)。

图8C、D是将基于来自试样1～10的R[k](％)进行的Tc计算的准确性按对应的Zs[k]区分地表示的标绘图(C)以及其放大图(D)。纵轴是将各试样中的Zs＝1.001(Zs[8])时的Tc计算值设为100％时的Zs[1]～[7]的相对Tc计算值，且表示计算值的准确性(％)。对于准确性成为96％以上的Zs值而言，试样10(APTT＝约130秒)不符合(即使Zs[7]＝1.006，也达不到96％)，在试样9(APTT＝约90秒)中Zs[7]＝1.006和Zs[6]＝1.011符合。图8E是将X(te[k])/X(te[8])(％)按对应的Zs[k]区分地表示的标绘图。纵轴是将各试样中的Zs＝1.001(Zs[8])的X(te[8])设为100％时的Zs[1]～[7]的X(te[1])～X(te[7])的相对值，且表示反应曲线X的相对高度(％)。相对高度成为与APTT的大小对应的分布。与图8A对比可知，随着反应曲线的相对高度接近100％，R[k]也朝向100％降低。

根据以上，示出有在决定凝固时间Tc的情况下，作为表示反应曲线的相对高度达到了90％(即凝固反应处于结束阶段)的基准，能够设定[R[k](％)≤105％]。

3)对于指标(b)的Tc决定基准

图9A是将X(te[k])/X(te[8])(％)按对应的Zs[k]区分表示的标绘图(与图8E相同)，图9B～D是指标(b)：TR[k](＝ΔC/ΔE)(％)、ΔE以及ΔC相对于APTT的标绘图。ΔE和ΔC具有根据APTT值而增加的趋势，ΔE的变化量比ΔC大(图9C～D)。其结果，如图9B所示，作为两者的变化量的比的TR[k](＝ΔC/ΔE)也根据APTT而增加。TR[k]相对于APTT的分布示出与将图9A(X(te[k]))的纵轴反转时相同的行为。因此，TR[k]间接表示X(te[k])中的反应曲线的相对高度。

3.Tc决定

在表3～表6中示出有针对试样1～10所取得的各参数以及指标。在任一个试样均没有达到最大测量时间的情况下结束Tc计算工序。Tc被决定以后的参数和指标的栏用斜线表示。灰色的栏表示指标满足基准(a)或(b)。

[表3]

试样1

试样2

试样3

[表4]

试样4

试样5

试样6

[表5]

试样7

试样8

试样9

[表6]

试样10

在图10A、图10B中，将凝固时间Tc被决定时的临时凝固点tc[k](A)和临时凝固结束点te[k](B)相对于APTT进行了标绘(在图中，本法)。在图中，对作为对照Zs＝1.001(k＝8)、以及作为比较例Zs＝1.031(k＝2)时的tc[k]、te[k]也一同进行了标绘。另外，在图10C以及D中分别示出从对照(Zs＝1.001、k＝8)的tc[k]以及te[k]减去本法或比较例1中的值而得的差值。

在本法中Tc被决定时的tc[k](＝被决定的Tc)与对照中的tc[k](与APTT大致相等)重叠。本法与对照的tc[k]的差为1秒以内。另一方面，本法中的te[k]在APTT短的情况下接近比较例，但在APTT90秒以后与对照重叠。te[k]最短为33秒(右图)。这表示在APTT小于90秒的情况下，本法所得到Tc所需的凝固反应测量时间比对照短约20秒至约30秒。另一方面，在比较例与对照之间，tc[k](＝被决定的Tc)的差在APTT到50秒为止时为1秒以内，但随着APTT的延长而变大，最大为29秒，另一方面te[k]的差最小为26秒，最大为142秒。即，在比较例中，凝固反应测量时间与对照相比变短，但在APTT小于50秒的情况下可得到与对照同等的凝固时间，但当APTT大于50秒时，与对照相比，凝固时间被计算得较短。根据以上表明，通过本发明的方法能够兼顾更准确的凝固时间的计算和测量时间的缩短化。

实施例2

作为检体，使用XII因子浓度为0.1％以下的Factor XII Deficient Plasma(George King Bio-Medical，Inc.制)(以下，称为FXII缺乏血浆、FXII浓度0％)、和FXII浓度被视为100％的正常混合血浆(以下，称为正常血浆)，将FXII缺乏血浆与正常血浆的混合比率设为1：39，而调制了FXII浓度为0.25％的混合血浆(试样11)。

针对所调制的试样11，按照与实施例1的1.10)同样的过程搜索满足基准(a)以及(b)的Tc。但是，在该试样中，虽然到k＝6(Zs＝1.011)为止得到了te[k]、te[k]，但由于在取得k＝7(Zs＝1.006)之前达到了最大测量时间所以Tc计算工序结束。结果，在本例中，不输出凝固时间Tc而输出表示“反应中途”的标志。在图11中示出有针对试样11的反应曲线X(图示有te[k]、te[k])、以及计算出的参数以及指标。

实施例3

针对纤维蛋白原浓度为低值(62mg/dL)的检体(试样12)，按照与实施例1的1.10)同样的过程搜索满足基准(a)以及(b)的Tc。在k＝2(Zs＝1.031)时满足基准(a)以及(b)，凝固时间Tc决定为30.9秒。在图12中示出有试样12的反应曲线X(对k＝1～8的te[k]、te[k]进行了图示)、以及计算出的参数以及指标。图12的表中灰色的栏表示指标满足基准(a)或(b)。与纤维蛋白原浓度为正常水平的图6左边相比较，在图12中的凝固反应的最大值约为八分之一，这样能够确认即使在纤维蛋白原浓度低的检体中也能够通过本发明的方法来计算凝固时间Tc。

实施例4

对Tc计算的计算开始点s的条件进行评价。试样使用了在反应初始反应曲线持续缓慢的上升的检体(试样13)。按照以下条件来确定计算开始点s，按照与实施例1的1.10)同样的过程，对满足基准(a)以及(b)的Tc进行了搜索。

条件1(与以往技术相同)：在累计比Z越过峰值以后，即Z超过规定值后，将低于该规定值的点设为计算开始点s。Z的规定值设为1.2。

条件2：在夹着峰值而成为V(i)＝Vs的两点间的测量点数超过Ws时，将成为V(i)＝Vs的较慢的测量点设为计算开始点s(参照图3)。Vs设定为90％Vmax，Ws设定为10(1秒)。

条件3：与条件2相同，但Vs设定为50％Vmax，Ws设定为10(1秒)。

在图13示出有在条件1～3下取得的各参数以及指标。左图表示反应曲线X和微分值V，右图表示X和累计比Z。在图中的标识符中，白色方块表示Te[k]，白色圆圈表示Tc[k]。右端是各参数以及指标的表。表中的灰色的栏表示指标满足基准(a)或(b)。

在条件1中，由于在反应初始的反应曲线X的缓慢的上升部中，Z低于规定值，因此在V(i)成为最大之前开始Tc计算，而得到k＝1～4的Te[k]、Tc[k]。但在，在不满足基准(a)以及(b)的状态下进入Tc计算工序，最终在k＝8(Zs＝1.001)时决定了Tc163.6秒。

在条件2、3中，在V(i)超过峰值后开始Tc计算。图中的黑色方块为计算开始点s。条件2、3都是在k＝8(Zs＝1.001)时满足基准(a)以及(b)，决定了Tc163.6秒。但是在条件3中，k＝1～4的Te[k]是接近的，因此对k＝1～4时的指标以不采用到Tc计算中的方式进行了处理。