JP2005201843A - 分析装置 - Google Patents

Abstract

【課題】
目的成分の濃度や活性値が分析装置の測定限界を超える高濃度・高活性試料に対して、不要な再検を削減し、試料中の目的成分の濃度や活性値が測定可能な分析技術を提供する。
【解決手段】
試料と１種類以上の試薬を混合した混合液を光学的に測定した測光値の時系列データをもとに前記試料中の目的成分の濃度又は活性値を測定する分析装置であって、前記試料の分析が正常に実施されたときの前記時系列データのモデルを決定付ける前記目的成分毎の化学反応に応じた1種類以上の係数と、前記試料の成分の濃度に応じた1種類以上の係数と、を含むパラメータを記憶する記憶手段１０２と、実際に測定された１点以上の測光値と前記パラメータを用いて、測定開始から所定の経過時間までの予測測光値を算出する算出手段１０６とを有することを特徴とする。
【選択図】 図１

Description

本発明は、試料中の目的成分の濃度や活性値を測定する分析技術に関する。

自動分析装置とは、試料と試薬を反応容器内で混合し、化学反応が進行するに連れて変化する吸光度の測光時系列データ（以降、反応過程データと呼ぶ）をもとに、一定時間（通常１０分程度）の吸光度変化量や吸光度変化率を計算し、試料中の目的成分の濃度や活性値を算出する装置である（例えば、特許文献１参照）。

このような自動分析装置は、主に医療機関において生化学検査や免疫検査等に用いられ、その検査結果は、医師が患者の病状把握や治療効果の判定、予後の経過観察など各種の診断を行う上で極めて重要な役割を担っている。

このため、患者の疾患の種類や重症度によっては、目的成分の濃度や活性値が極めて高値となり、装置の測定限界（正確さを保証できる範囲）を超えてしまうことがある。このような場合、混合液に対する試料の比率が小さくなるように設定し、再度同様の検査を実施する（以降、再検と呼ぶ）ことで、再検結果と試料混合比率をもとに、試料中の目的成分の濃度や活性値を算出することが可能となる。

特開平０９−３２５１５０号公報

この再検に関して、従来の自動分析装置では、試料を予め設定された量や比率で減量・増量・希釈して自動的に再検する自動分析装置がある。また、検査技師が手作業で再検する場合は、経験に基づいて試料を希釈し、測定結果をもとの濃度に換算する再検方法を実施する場合もある。

このため、再検時の試料混合比率が不適切だった場合、再検結果が測定限界を超えてしまい、再度減量・増量・希釈して再検しなければならないため、検査コストが上がるという問題があった。

本発明の目的は、目的成分の濃度や活性値が分析装置の測定限界を超える高濃度・高活性試料に対して、不要な再検を削減し、試料中の目的成分の濃度や活性値が測定可能な分析技術を提供することにある。

上記目的を達成するために、本発明による分析装置は、以下に示すような特徴を有する。

（１）試料と１種類以上の試薬とを混合した混合液を光学的に測定した測光値の時系列データをもとに前記試料中の目的成分の濃度または活性値を測定する分析装置であって、前記試料の分析が正常に実施されたときの前記時系列データのモデルを決定付ける前記目的成分毎の化学反応に応じた1種類以上の係数と、前記試料の成分の濃度に応じた1種類以上の係数とを含むパラメータを記憶する記憶手段と、実際に測定された１点以上の測光値と前記パラメータを用いて、測定開始から所定の経過時間までの予測測光値を算出する算出手段とを有することを特徴とする。

（２）前記（１）の分析装置であって、前記予測測光値が分析結果の正確さを保証できる範囲を示す測定限界を超えるか否かを判定する判定手段を有することを特徴とする。これにより、実際の分析で測定限界超過が判明するよりも早いタイミングで再検等の次の処理に移行することができるので、ＴＡＴ（Turn Around Time：検査依頼から結果の報告までの時間）を短縮しつつ、解決できる。

（３）前記（１）又は（２）の分析装置であって、前記パラメータと前記予測測光値を用いて、前記混合液に対する前記試料の最適な混合比率を算出する混合比率算出手段を有することを特徴とする。これにより、試料混合比率を適切に設定できるので、不要な再検を削減しつつ、解決できる。

（４）前記（１）、（２）又は（３）の分析装置であって、前記予測測光値をもとに、前記試料中の目的成分の濃度または活性値を算出して出力する出力手段を有することを特徴とする。これにより、分析結果が出力されるよりも早く予測分析結果を出力できるので、技師の結果報告や再検等の判断の迅速性を向上しつつ、解決できる。

（５）前記（４）の分析装置であって、前記出力手段が、実際に測定された測光値と前記予測測光値とを同時に異なる表現方法で出力することを特徴とする。これにより、実際に測定された測光値と前記予測測光値を比較することができるので、技師の再検の必要性の判断を効率化しつつ、解決できる。

（６）本発明の分析装置は、試料と試薬とを混合した混合溶液の化学反応の進行につれて変化する吸光度を所定の時間間隔で複数回測光して、前記吸光度の測光時系列データを測定する手段と、前記測光時系列データをもとに、前記吸光度の一定時間での変化量もしくは変化率を算出し、前記試料中の目的成分の濃度もしくは活性値を求める手段とを有し、かつ、前記試料の分析が正常に実施されたときの前記測光時系列データのモデルを決定付ける前記目的成分毎の化学反応に応じた1種類以上の係数と、前記試料の成分の濃度に応じた1種類以上の係数とを含むパラメータを記憶する記憶手段と、前記測光時系列データと前記パラメータを用いて、前記所定の時間間隔での予測吸光度を算出する算出手段とを設けることにより、前記予測吸光度をもとに前記試料中の目的成分の濃度もしくは活性値を求めるよう構成したことを特徴とする。

（７）前記（６）の分析装置において、前記予測吸光度が分析結果の正確さを保証できる範囲を示す測定限界を超えるか否かを判定する判定手段を有することを特徴とする。

（８）前記（７）の分析装置において、前記パラメータと前記予測吸光度とを用いて、前記混合液に対する前記試料の混合比率を算出する混合比率算出手段を有することを特徴とする。

（９）前記（６）、（７）又は（８）の分析装置において、前記予測吸光度をもとに、前記試料中の目的成分の濃度または活性値を算出して出力する出力手段を有することを特徴とする。

（１０）前記（９）の分析装置において、前記出力手段が、前記予測吸光度と共に、前記測光時系列データを出力し表示する構成されていることを特徴とする。

（１１）前記（９）の分析装置において、前記出力手段が、前記予測吸光度が分析結果の正確さを保証できる範囲を示す前記測定限界を超える場合における前記混合比率を、前記測光時系列データおよび前記予測吸光度のうち少なくとも一つと共に出力し表示するよう構成されていることを特徴とする。

本発明は、目的成分の濃度や活性値が自動分析装置の測定限界を超える高濃度・高活性試料に対して、化学反応モデルを表現するパラメータを用いて測定結果を予測し、不要な再検を削減して、試料中の目的成分の濃度や活性値の測定を可能にする分析技術を実現する。

以下、本発明の実施例について、図面を参照して詳述する。

図１に、本発明の一実施例における吸光度予測システムの構成図を示す。

本システム１００は、制御部１０１と、一時記憶装置１０２と、反応モデルパラメータＤＢ読み出し部１０３と、反応過程データ抽出部１０４と、測定パラメータＤＢ読み出し部１０５と、予測測光値算出部１０６と、予測分析結果算出部１０７と、試料混合比率算出部１０８と、判定部１０９と、出力部１１０と、反応モデルパラメータＤＢ１１１と、測定パラメータＤＢ１１２と、で構成される。

本システム１００はハードウェア構成として記載しているが、本システム１００の機能はソフトウェアで構成されていてもよい。

本システム１００は、ネットワーク１２０を介して、自動分析装置１３０と通信できる。また、本システム１００は、ネットワーク１２０を介して、入出力端末１４０と通信できる。

ネットワーク１２０は、検査施設内のネットワークを前提としているが、検査部門を有する医療施設内のネットワークでもよい。

自動分析装置１３０は、反応過程データを反応過程データ記憶装置１３１に格納する。また、反応過程データ記憶装置１３１に保存された反応過程データは、ネットワーク１２０を介して、反応過程データ抽出部１０４によって抽出できる。また、反応過程データ記憶装置１３１に保存された反応過程データは、ネットワーク１２０を介して、入出力端末１４０で閲覧できる。

自動分析装置１３０は、本システム１００とは別のハードウェアとして記載しているが、本システム１００の構成が自動分析装置１３０内に構成されていてもよい。

反応過程データ記憶装置１３１は、自動分析装置１３０内に構成されることを前提としているが、自動分析装置１３０と反応過程データ記憶装置１３１が別のハードウェアで構成されていてもよい。また、反応過程データ記憶装置１３１が、本システム１００内に構成されていてもよい。また、反応過程データ記憶装置１３１が、入出力端末１４０内に構成されていてもよい。

入出力端末１４０は、キーボードやマウス等を入力機能、ＣＲＴディスプレイを出力機能とするパソコン等の情報機器を想定しているが、他の入出力機能を有していてもよい。また、入出力端末１４０は、Ｗｅｂブラウザ機能を搭載した端末でもよい。また、入出力端末１４０は、本システム１００とは別のハードウェアとして記載しているが、入出力端末１４０の入出力機能が本システム１００に搭載されていてもよい。また、入出力端末１４０は、自動分析装置１３０とは別のハードウェアとして記載しているが、入出力端末１４０の入出力機能が自動分析装置１３０に搭載されていてもよい。また、入出力端末１４０は、パーソナルコンピュータを前提としているが、臨床検査システムでもよい。

また、入出力端末１４０のユーザは、検査技師等、自動分析装置１３０の操作者を前提としているが、自動分析装置１３０のメンテナンス担当者等、他のユーザでもよい。

また、本システム１００と、自動分析装置１３０と、入出力端末１４０は、１つのハードウェアとして構成されていてもよい。

図２に、自動分析装置１３０の構成例を示す。自動分析装置１３０は、光源ランプ２０１と、恒温槽２０２と、セル２０３と、試料分注ノズル２０４と、第１試薬分注ノズル２０５ａと、第２試薬分注ノズル２０５ｂと、撹拌棒２０６と、分光器２０７と、検知器２０８と、増幅器２０９と、Ａ／Ｄ変換器２１０と、で構成される。分析時には、光源ランプ２０１から発せられた白色光(全波長)が、恒温槽２０２につけられたセル２０３を透過して分光器２０７に入り、分析項目によって異なる特定の単波長成分が検知器２０８で受光され、増幅器２０９による増幅後、Ａ／Ｄ変換器２１０でＡ／Ｄ変換され、吸光度として出力される。

セル２０３は反応容器になっており、試料分注ノズル２０４から試料が、第１試薬分注ノズル２０５ａから第１試薬（Ｒｅａｇｅｎｔ１、以降Ｒ１と呼ぶ）が、第２試薬分注ノズル２０５ｂから第２試薬（Ｒｅａｇｅｎｔ２、以降Ｒ２と呼ぶ）が、それぞれ分注され、撹拌棒２０６によって撹拌されることで、セル内部で化学反応が起こる。このときの化学反応について、経時的に吸光度を測定（測光）することで、試料中の分析物の濃度や活性値に換算することが可能となる。

以下に、分析時の自動分析装置１３０の動作について説明する。
（１）まず初めに、セル２０３を水(純水)で満たし、水ブランク測光を行う。この値は、以降測定される吸光度の基準となる。
（２）次に、セル２０３内の水を排出し、セル２０３内に試料を分注する。
（３）次に、セル２０３内にＲ１を分注し、撹拌する。
（４）次に、試料とＲ１の混合溶液を一定時間間隔で複数回測光する。
（５）次に、セル２０３内にＲ２を分注し、撹拌する。
（６）次に、試料とＲ１とＲ２の混合溶液を一定時間間隔で複数回測光する。
（７）１０分間の反応時間において計３４回の測光を行った後、セル２０３を洗浄し、分析を終了する。

このとき、Ｒ１添加からＲ２添加までの反応を第１反応と呼び、Ｒ２添加から測光終了までの反応を第２反応と呼ぶ。また、分析項目によっては、Ｒ１のみを使用し、上記（５）、（６）の動作を実行しないものもある。このように、使用する試薬の数により、１試薬系分析、２試薬系分析に分類される。

測定した吸光度は、吸光度が物質の濃度に比例する法則（Ｌａｍｂｅｒｔ−Ｂｅｅｒの法則）に基づき、（数１）で示すような検量線と呼ばれる直線を示す換算式により濃度あるいは活性値Ｃに換算され、反応過程データ記憶装置１３１に保存される。

Ｃ＝ａ×Ａ＋ｂ …（数１）
（数１）において、ａ、ｂは予め設定された値で、分析項目によって異なる。また、Ａは、ある測光ポイントでの吸光度変化量、または、複数の測光ポイント間の吸光度変化率である。このように、吸光度を試料中の目的成分の濃度あるいは活性値に換算する方法としては、吸光度変化量を利用するエンド法、及び、吸光度変化率を利用するレート法がある。

図３に、エンド法が用いられる分析項目の例として、ＴＰ（総蛋白）の反応過程データの例３００を示す。横軸は、反応時間（Ｔｉｍｅ）を示し、縦軸は、吸光度（ＡＢＳ）を示す。

ＴＰは１試薬系分析であり、図３下部に示す化学反応式３０１のように、試料中のペプチド（−ＣＯ−ＮＨ−）と、Ｒ１中の銅イオン（Ｃｕ^２＋）が、キレート化合物（金属イオンに有機分子が結合したもの）に変化する化学反応を利用し、キレート化合物の濃度変化を測定することで、ＴＰの濃度を算出する。そのため、測光波長は生成物であるキレート化合物の吸収ピーク波長を用いる。

反応過程データ例３００に示すように、反応開始直後はキレート化合物が増加し吸光度が急激に上昇する。しかし、分析終了間際には、試料中のペプチドが全てキレート化合物に変化するため、反応がほとんど進行しなくなり、吸光度の上昇は見られなくなる。この分析終了時点での吸光度変化量（Ａ＝ΔＡＢＳ）を用いることで、ＴＰの濃度を算出する。

しかし、試料中のＴＰの濃度が高値の場合、ＴＰやキレート化合物による光の散乱の影響が無視できなくなるため、Ｌａｍｂｅｒｔ−Ｂｅｅｒの法則に従わなくなる。このため、（数１）で示される吸光度変化量Ａと濃度Ｃの比例関係が成り立たなくなる。

図４に、高濃度域で非直線となるＴＰの検量線の例４００を示す。検量線４００には、（数１）の直線性が保証されるＴＰの測定限界濃度Ｃ_{ＬＩＭＩＴ}と、濃度Ｃ_{ＬＩＭＩＴ}のときの吸光度変化量Ａ_{ＬＩＭＩＴ}が予め設定されている。このとき、吸光度変化量Ａ_ＬＯＷ（≦Ａ_{ＬＩＭＩＴ}）、Ａ_ＨＩＧＨ（≧Ａ_{ＬＩＭＩＴ}）のときの濃度をそれぞれＣ_ＬＯＷ（≦Ｃ_{ＬＩＭＩＴ}）、Ｃ_ＨＩＧＨ（≧Ｃ_{ＬＩＭＩＴ}）とすると、Ａ_ＬＯＷとＣ_ＬＯＷには（数１）の関係が成り立つが、Ａ_ＨＩＧＨとＣ_ＨＩＧＨには（数１）の関係は成り立たない。以降、このＡ_{ＬＩＭＩＴ}に対応した吸光度を、エンド法における測定限界吸光度と呼ぶ。

このように、エンド法が用いられる項目で、分析が終了する前に測定限界吸光度に達してしまう高濃度試料を分析する場合、通常の分析方法では濃度を算出することは困難である。したがって、測定限界吸光度に達した場合は、分析終了時点での吸光度が測定限界吸光度よりも小さくなるように試料混合比率を変更して、再検を実施する必要がある。

また、図５に、レート法が用いられる分析項目の例として、ＬＤＨ（乳酸脱水素酵素）の反応過程データの例５００を示す。

ＬＤＨは２試薬系分析であり、図５下部に示す化学反応式５０１のように、試料中の酵素ＬＤＨを触媒として、Ｒ２中のピルビン酸が乳酸に、Ｒ１中のＮＡＤＨ（β−ニコチンアミドアデニンジヌクレオチド還元型）がＮＡＤ（β−ニコチンアミドアデニンジヌクレオチド酸化型）に変化する化学反応を利用し、ＮＡＤＨの減少率を測定することで、ＬＤＨの活性値を算出する。そのため、測光波長はＮＡＤＨの吸収ピーク波長を用いる。

反応過程データ例５００に示すように、第１反応では吸光度はほぼ一定であるが、第２反応では分析終了まで吸光度が一定の傾きで低下している。これは、試料中のＬＤＨの活性に応じて、ＮＡＤＨが徐々に減少している様子を示している。この第２反応時の吸光度変化率Ａ＝ΔＡＢＳ／ΔＴｉｍｅ（予め設定された２点以上の測光ポイントにおける吸光度の傾き）を用いることで、試料中のＬＤＨの活性値を算出する。

しかし、試料中のＬＤＨの活性値が高く、吸光度変化率を計算するために予め設定された測光ポイントよりも前に全てのＮＡＤＨがＮＡＤに変化してしまうと、適切な吸光度変化率を計算できず、試料中のＬＤＨの活性値を適切に算出できない。

図６に、高活性ＬＤＨの反応過程データの例６００を示す。反応過程データ例６００では、時間Ｔ_Ｌ〜Ｔ_Ｒ間の吸光度変化率６０１を用いている。しかし、ＬＤＨの活性が高く、時間Ｔ（＜Ｔ_Ｒ）で反応が終了してしまい、吸光度６０２に達すると、実際の吸光度変化率６０３を計算できない。以降、この吸光度６０２を、レート法における測定限界吸光度と呼ぶ。

このように、レート法が用いられる項目で、吸光度変化率を計算するのに用いる測光ポイント（図６中のＴ_ＬやＴ_Ｒ）よりも前に測定限界吸光度に達して、反応が終了してしまう高活性試料を分析する場合、通常の分析方法では活性値を算出することは困難である。したがって、測定限界吸光度に達した場合は、吸光度変化率を計算するのに用いる測光ポイントよりも前に測定限界吸光度に達しないように、試料混合比率を変更して、再検を実施する必要がある。

ここで、エンド法において、分析終了まで光の散乱が無視できると仮定した場合、高濃度域の検量線も直線として考えることができる。これを利用して、本システムでは、化学反応をモデル化し、検量線が直線である吸光度域（図４中のＡ_{ＬＩＭＩＴ}以下の範囲）の反応過程データをもとに、分析終了まで光の散乱が無視できると仮定したときの各測光ポイントにおける予測吸光度を算出することで、高濃度試料の測定を実現する。

また、レート法において、吸光度変化率を計算するのに用いる測光ポイント（図６中のＴ_ＬやＴ_Ｒ）よりも前に反応が終了しないように、試薬中の反応物質が混合液中に十分な量存在すると仮定した場合、吸光度変化率を計算することは可能となる。これを利用して、本システムでは、化学反応をモデル化し、反応が終了する測光ポイント（図６中のＴ）よりも前の反応過程データをもとに、試薬中の反応物質が十分な量存在すると仮定したときの各測光ポイントにおける予測吸光度を算出することで、高活性試料の測定を実現する。

このように、本システムでは、化学反応をモデル化し、反応過程データをもとに、最適条件下で分析できると仮定したときの各測光ポイントにおける予測吸光度を算出することで、高濃度・高活性試料の測定を実現する。

以下、その手順を示す。

本実施例では、まず、セル中で起こる化学反応の過程を、試料、試薬、装置全てに異常が無い最適条件下で測定されると仮定し、（数２）のような時間ｔを変数とする指数関数でモデル化する。この（数２）で表現される反応過程データを、反応モデルデータと呼ぶ。

ＡＢＳ_{ＭＯＤＥＬ}（ｔ）＝Ａ_０＋Ａ_１×（１−ｅｘｐ（−ｋ×ｔ））
…（数２）
このとき、ｋは化学反応の種類（分析項目と試薬の組み合わせ）によって一意に決まる反応の速度、Ａ_０は反応開始時の初期吸光度、Ａ_１は目的成分の濃度、を示すパラメータ（以降反応モデルパラメータと呼ぶ）である。

図７（Ａ）に、反応モデルパラメータＤＢ１１１の例７００を示す。また、図７（Ｂ）に、測定パラメータＤＢ１１２の例７１０を示す。

反応モデルパラメータＤＢ１１１は、精度管理試料（Quality Control、以降、ＱＣと呼ぶ）の種類を識別するＱＣＩＤを格納するフィールド７０１と、分析項目名を格納するフィールド７０２と、測定パラメータを識別する測定パラメータＩＤを格納するフィールド７０３と、分析項目の濃度あるいは活性値の測定値を格納するフィールド７０４と、レート法の分析項目に関して、正規化（後述）後の第１反応の吸光度の基準値を格納するフィールド７０５と、反応モデルパラメータｋ、Ａ_０、Ａ_１をそれぞれ格納するフィールド７０６、７０７、７０８と、で構成されている。

測定パラメータＤＢ１１２は、測定パラメータを識別する測定パラメータＩＤを格納するフィールド７１１と、分析項目名を格納するフィールド７１２と、Ｒ１とＲ２の種類を識別する試薬ＩＤを格納するフィールド７１３、７１４と、分析に使用した波長を格納するフィールド７１５と、測定限界吸光度を格納するフィールド７１６と、レート法において吸光度変化率を計算するのに利用する測光ポイント（図６中のＴ_Ｌ、Ｔ_Ｒ）を示す番号を格納するフィールド７１７、７１８と、濃度または活性値への換算式を決定する（数１）におけるａ、ｂをそれぞれ格納するフィールド７１９、７２０と、で構成されている。

反応モデルパラメータ例７００及び測定パラメータ例７１０では、測定パラメータＩＤ「００１」で識別される試薬「Ｒ_Ａ１」と「Ｒ_Ａ２」を使用して波長「λ_Ａ」で分析したＬＤＨの、活性値が「１８５」のＱＣＩＤ「Ａ」で識別されるＱＣについて、正規化基準値「１００００」のとき、「ｋ＝０.０００６１」、「Ａ_０＝１４３５０」、「Ａ_１＝−３４５０」であることを示している。

また、反応モデルパラメータ例７００及び測定パラメータ例７１０では、測定パラメータＩＤ「００２」で識別される試薬「Ｒ_Ｂ１」を使用して波長「λ_Ｂ」で分析したＴＰの、濃度が「６．８」のＱＣＩＤ「Ｂ」で識別されるＱＣについて、反応モデルパラメータは、「ｋ＝０.００４８５」、「Ａ_０＝１２５５」、「Ａ_１＝８１０」であることを示している。

また、測定パラメータ例７１０では、測定パラメータＩＤ「００１」で測定されるＬＤＨについて、測定限界吸光度が「≦２０００」で、吸光度２０００以下のとき測定限界に達したことを示し、吸光度変化率を計算するのに利用する測光ポイントが「１９」ポイント目と「３３」ポイント目で、吸光度変化率Ａから活性値Ｃへの換算式がＣ＝−１４０×Ａであることを示している。

また、例７１０では、測定パラメータＩＤ「００２」で測定されるＴＰについて、測定限界吸光度が「≧３５７０」で、吸光度３５７０以上のとき測定限界に達したことを示し、吸光度変化量Ａから濃度Ｃへの換算式がＣ＝０.００３５×Ａ＋０．５であることを示している。

本システムは、反応モデルパラメータをもとに、高濃度・高活性試料の反応過程データＡＢＳ_Ｈ（ｔ）に対して、最適条件下で分析できると仮定したとき、（数２）と同様の（数３）で示すような反応モデルデータＡＢＳ_{Ｈ＿ＭＯＤＥＬ}（ｔ）を算出する。

ＡＢＳ_{Ｈ＿ＭＯＤＥＬ}（ｔ）＝
Ａ_Ｈ０＋Ａ_Ｈ１×（１−ｅｘｐ（−ｋ×ｔ）） …（数３）
このとき、ｋは反応の速度、Ａ_Ｈ０は反応開始時の初期吸光度、Ａ_Ｈ１は目的成分の濃度、を示すパラメータである。また、この反応モデルデータＡＢＳ_{Ｈ＿ＭＯＤＥＬ}（ｔ）は、最適条件下で分析できると仮定したときの各測光ポイントにおける予測測光値と言い換えることができる。

本システムは、予測測光値をもとに算出した濃度・活性値を、予測分析結果として出力する。

図８に、反応モデルパラメータをもとに、反応過程データＡＢＳ_Ｈ（ｔ）に対して、最適条件下で分析できると仮定したとき、（数２）と同様の（数３）で示すような反応モデルデータＡＢＳ_{Ｈ＿ＭＯＤＥＬ}（ｔ）を算出し、入出力端末１４０に予測分析結果を出力するときのフローチャートを示す。

まず、制御部１０１が、測定パラメータＤＢ読み出し部１０５を起動し、測定パラメータＤＢ１１２の各フィールドの情報を、一時記憶装置１０２に保存するステップ８０１を実行する。

次に、制御部１０１が、予測測光値算出部１０６を起動し、反応開始時の測光ポイントＴ＿Ｓを設定して、一時記憶装置１０２に保存するステップ８０２を実行する。ステップ８０２では、１試薬系の分析項目の場合、Ｒ１投入直後の測光ポイント、２試薬系の分析項目の場合、Ｒ２投入直後の測光ポイントを、Ｔ_{ＳＴＡＲＴ}とする。

次に、制御部１０１が、反応過程データ抽出部１０４を起動し、反応過程データ記憶装置１３１から、ｔ＝Ｔ_１〜Ｔ_Ｎにおける反応過程データＡＢＳ_Ｈ（ｔ）を抽出し、一時記憶装置１０２に保存するステップ８０３を実行する。ただし、Ｔ_１は最初に測定された測光ポイント（１ポイント目）、Ｔ_Ｎはステップ８０３実行直前に測定された測光ポイント（Ｎポイント目）、とする。

次に、制御部１０１が、予測測光値算出部１０６を起動し、「Ｔ_{ＳＴＡＲＴ}＜Ｔ_Ｎ」の真偽を判定するステップ８０４を実行する。

ステップ８０４において偽と判定された場合、本システムは動作を終了する。

ステップ８０４において真と判定された場合、制御部１０１が、予測測光値算出部１０６を起動し、測定限界吸光度に達する直前の測光ポイントＴ_ＥＮＤを設定して、一時記憶装置１０２に保存するステップ８０５を実行する
図９（Ａ）、（Ｂ）に、測定限界吸光度「≧３５７０」であるエンド法を用いる分析項目で、ステップ８０５でＴ_ＥＮＤを設定するときの反応過程データの例９００及び例９１０を示す。図９（Ａ）では、Ｔ_Ｎの時点で吸光度ＡＢＳ_Ｈ（Ｔ_Ｎ）が３５７０よりも小さいため、Ｔ_ＥＮＤ＝Ｔ_Ｎと設定する。図９（Ｂ）では、Ｔ_Ｎの時点で吸光度ＡＢＳ_Ｈ（Ｔ_Ｎ）が３５７０よりも大きいため、反応過程データＡＢＳ_Ｈ（ｔ）が３５７０よりも大きくなる直前の測光ポイント（Ｎ−２ポイント目）を検出して、Ｔ_ＥＮＤ＝Ｔ_Ｎ−２と設定する。

また、図１０（Ａ）、（Ｂ）に、測定限界吸光度「≦２０００」であるレート法を用いる分析項目で、ステップ８０５でＴ_ＥＮＤを設定するときの反応過程データの例１０００及び例１０１０を示す。図１０（Ａ）では、Ｔ_Ｎの時点で吸光度ＡＢＳ_Ｈ（Ｔ_Ｎ）が２０００よりも大きいため、Ｔ_ＥＮＤ＝Ｔ_Ｎと設定する。図１０（Ｂ）では、Ｔ_Ｎの時点で吸光度ＡＢＳ_Ｈ（Ｔ_Ｎ）が２０００よりも小さいため、反応過程データＡＢＳ_Ｈ（ｔ）が２０００よりも小さくなる直前の測光ポイント（Ｎ−１ポイント目）を検出して、Ｔ_ＥＮＤ＝Ｔ_Ｎ−１と設定する。

ステップ８０５において、レート法を用いる分析項目で、隣り合う測光ポイント間の傾きの差がある一定の値よりも大きくなる測光ポイントを、反応が終了した時間として検出し、この測光ポイントをＴ_ＥＮＤとして設定することもできる。これにより、反応に無関係な物質の吸光が含まれる場合など、測定限界吸光度に達する前に反応が終了してしまう場合でも、Ｔ_ＥＮＤを適切に設定することができる。

次に、制御部１０１が、予測測光値算出部１０６を起動し、Ｔ_{ＳＴＡＲＴ}〜Ｔ_ＥＮＤ間に、一定数以上の測光ポイントが含まれているかどうかを判定するステップ８０６を実行する。

ステップ８０６により、利用する測光ポイント数をできるだけ多くすることで、以降の処理で算出するパラメータＡ_Ｈ１、Ａ_Ｈ０を、統計的な量として算出することができる。

ステップ８０６において偽と判定された場合、本システムは動作を終了する。

ステップ８０６において真と判定された場合、制御部１０１が、反応モデルパラメータＤＢ読み出し部１０３を起動し、反応モデルパラメータＤＢ１１１の各フィールドの情報を、一時記憶装置１０２に保存するステップ８０７を実行する。

次に、制御部１０１が、予測測光値算出部１０６を起動し、一時記憶装置１０２上の反応過程データＡＢＳ_Ｈ（ｔ）が、エンド法を用いる分析項目か、レート法を用いる分析項目かを判断するステップ８０８を実行する。

ステップ８０８でレート法と判断された場合、制御部１０１が、予測測光値算出部１０６を起動し、一時記憶装置１０２上の反応過程データＡＢＳ_Ｈ（ｔ）を、第１反応の吸光度の平均値が、反応モデルパラメータＤＢ１１１のフィールド７０５に格納された値になるように、第１反応、及び第２反応の測光ポイントの吸光度をそれぞれ正規化し、一時記憶装置１０２に保存するステップ８０９を実行する。

以下に、ステップ８０９の具体的な例を示す。

正規化後の第１反応の吸光度の平均値を「ｃ」とし、一時記憶装置１０２上の反応過程データＡＢＳ_Ｈ（ｔ）について、第１反応の吸光度の平均値「ｃ’」、第１反応時の反応液量「Ｖ_１」、第２反応時の反応液量「Ｖ_２」であったとすると、ステップ８０９では、予測測光値算出部１０６が、第１反応の測光ポイントの吸光度に値「ｃ−ｃ’」を加算し、第２反応の測光ポイントの吸光度に値「(ｃ−ｃ’)×Ｖ_１／Ｖ_２」を加算する。

ステップ８０９により、レート分析法において、第１反応時の吸光度のずれの影響による第２反応時の吸光度のずれを補正することが可能となる。

次に、制御部１０１が、予測測光値算出部１０６を起動し、一時記憶装置１０２上の反応過程データＡＢＳ_Ｈ（ｔ）から、（数４）に示す測光ポイントｔ〜ｔ＋１（ｔ≧Ｔ_{ＳＴＡＲＴ}、ｔ＋１≦Ｔ_ＥＮＤ）の反応過程データＡＢＳ_Ｈ（ｔ）の変化量ΔＡＢＳ_Ｈ（ｔ）を算出し、一時記憶装置１０２に保存するステップ８１０を実行する。

ΔＡＢＳ_Ｈ（ｔ）＝ＡＢＳ_Ｈ（ｔ＋１）−ＡＢＳ_Ｈ（ｔ）
…（数４）
次に、制御部１０１が、予測測光値算出部１０６を起動し、一時記憶装置１０２上の反応過程データＡＢＳ_Ｈ（ｔ）に対して、（数５）に示すＡ_Ｈ１（ｔ）と、そのときのＡ_Ｈ１（ｔ）の平均を算出して、一時記憶装置１０２に保存するステップ８１１を実行する。このとき算出したＡ_Ｈ１（ｔ）の平均が、反応モデルデータＡＢＳ_{Ｈ＿ＭＯＤＥＬ}（ｔ）の濃度を示すパラメータＡ_Ｈ１となる。尚、Ａ_Ｈ１の符号は、反応過程データＡＢＳ_Ｈ（ｔ）が単調増加の場合は＋、単調減少の場合は−となる。

Ａ_Ｈ１（ｔ）＝±√（｛ΔＡＢＳ_Ｈ（ｔ）｝^２／
｛（−ｅｘｐ（−ｋ×（ｔ＋１）））−（−ｅｘｐ（−ｋ×ｔ））｝^２）
…（数５）
これにより、第１反応の測光ポイントと、測定限界吸光度に達した測光ポイントを除外し、適切な条件下で測定された測光ポイントのみを利用できるので、反応モデルデータＡＢＳ_{Ｈ＿ＭＯＤＥＬ}（ｔ）の濃度を示すパラメータＡ_Ｈ１を正確に算出することが可能となる。

次に、制御部１０１が、予測測光値算出部１０６を起動し、ステップ８１１で算出したＡ_Ｈ１を（数３）に代入し、時間ｔ＝Ｔ_{ＳＴＡＲＴ}〜Ｔ_ＥＮＤ毎に、（数３）の値と一時記憶装置１０２上の反応過程データＡ_Ｈ１（ｔ）との二乗誤差総和が最小となるＡ_Ｈ０を算出して、一時記憶装置１０２に保存するステップ８１２を実行する。このとき算出したＡ_Ｈ０が、反応モデルデータＡＢＳ_{Ｈ＿ＭＯＤＥＬ}（ｔ）の初期吸光度を示すパラメータＡ_Ｈ０となる。

次に、制御部１０１が、予測測光値算出部１０６を起動し、反応モデルパラメータＤＢ１１１のフィールド７０６のｋ、及びステップ８１１、ステップ８１２で算出したＡ_Ｈ１、Ａ_Ｈ０を（数３）に代入し、時間ｔ毎の反応モデルデータＡＢＳ_{Ｈ＿ＭＯＤＥＬ}（ｔ）を算出して、一時記憶装置１０２に保存するステップ８１３を実行する。

次に、制御部１０１が、予測分析結果算出部１０７を起動し、（数１）及び反応モデルデータＡＢＳ_{Ｈ＿ＭＯＤＥＬ}（ｔ）をもとに、予測分析結果を算出して、一時記憶装置１０２に保存するステップ８１４を実行する。

図１１（Ａ）に、エンド法を用いる分析項目の反応過程データＡＢＳ_Ｈ（ｔ）の例１１００とその反応モデルデータＡＢＳ_{Ｈ＿ＭＯＤＥＬ}（ｔ）の例１１０１を示す。また、図１１（Ｂ）に、そのときの検量線の例１１０２を示す。図１１（Ａ）では、ＡＢＳ_Ｈ（ｔ）をもとに、分析終了時間まで最適条件下で分析できると仮定したときの予測測光値であるＡＢＳ_{Ｈ＿ＭＯＤＥＬ}（ｔ）を算出している。通常の分析では、ＡＢＳ_Ｈ（ｔ）が測定限界吸光度１１０３に達するため、検量線１１０２により濃度１１０４が得られるが、ＡＢＳ_{Ｈ＿ＭＯＤＥＬ}（ｔ）を利用することで、適切な濃度１１０５を得ることができる。

また、図１２に、レート法を用いる分析項目の反応過程データＡＢＳ_Ｈ（ｔ）の例１２００とその反応モデルデータＡＢＳ_{Ｈ＿ＭＯＤＥＬ}（ｔ）の例１２０１を示す。図１２では、ＡＢＳ_Ｈ（ｔ）をもとに、吸光度変化率を計算するのに利用する測光ポイント１２０２〜１２０３間まで最適条件下で分析できると仮定したときの予測測光値であるＡＢＳ_{Ｈ＿ＭＯＤＥＬ}（ｔ）を算出している。通常の分析では、ＡＢＳ_Ｈ（ｔ）が測定限界吸光度１２０４に達するため吸光度変化率１２０５が得られるが、ＡＢＳ_{Ｈ＿ＭＯＤＥＬ}（ｔ）を利用することで、活性値算出に必要な適切な吸光度変化率１２０６を得ることができる。

次に、制御部１０１が、判定部１０９を起動し、反応モデルデータＡＢＳ_{Ｈ＿ＭＯＤＥＬ}（ｔ）が、反応終了時間までに測定限界に達するかどうかを判定して、その判定結果を一時記憶装置１０２に保存するステップ８１５を実行する。

次に、制御部１０１が、出力部１１０を起動し、一時記憶装置１０２上の反応過程データＡＢＳ_Ｈ（ｔ）と、反応モデルデータＡＢＳ_{Ｈ＿ＭＯＤＥＬ}（ｔ）と、予測分析結果と、ステップ８１５の判定結果を、入出力端末１４０に出力するステップ８１６を実行する。

図１３に、ｔ＝Ｔ_４（４ポイント目まで測定が終了した時点）において、図８に示すフローチャートを実行した場合で、ステップ８１６実行後の入出力端末１４０の画面例１３００を示す。画面例１３００では、ｔ＝Ｔ_４において、実際の測定済み反応過程データ１３０１と、反応モデルデータ１３０２と、反応モデルデータ１３０２を利用して算出した予測分析結果１３０３とを、将来測定限界に達することを知らせる表示とともに出力している。ｔ＝Ｔ_４の時点では、測定済み反応過程データ１３０１が測定限界吸光度１３０４に達していないため、反応開始時間１３０５からＴ_４までの全ての測光ポイントのデータを用いて、反応モデルデータ１３０２を算出している。

このように、吸光度を予測することによって、実際に測定された吸光度が測定限界に達するよりも早いタイミングで、測定限界に達することを技師が判断でき、速やかに再検等の次の処理に移行できるので、ＴＡＴを短縮することが可能となる。

また、吸光度を予測することによって、現在測定限界に達しているか否か、または将来測定限界に達するか否かに関わらず、技師は予測分析結果を参照できるので、迅速な結果報告の必要性の有無を、早いタイミングで判断できる。

また、吸光度を予測して予測分析結果を出力することにより、短時間での迅速分析が可能となる。

また、図１４に、ｔ＝Ｔ_８（８ポイント目まで測定が終了した時点）において、図８に示すフローチャートを実行した場合で、ステップ８１６実行後の入出力端末１４０の画面例１４００を示す。画面例１４００では、ｔ＝Ｔ_８において、実際の測定済み反応過程データ１４０１と、反応モデルデータ１４０２と、反応モデルデータ１４０２を利用して算出した予測分析結果１４０３とを、既に測定限界に達したことを知らせる表示とともに出力している。ｔ＝Ｔ_８の時点では、測定済み反応過程データ１４０１が測定限界吸光度１４０４に達したため、反応開始時間１３０５から測定限界内最終測光ポイント１４０６までの全ての測光ポイントのデータを用いて、反応モデルデータ１４０２を算出している。

また、図１５に、測定終了後に、図８に示すフローチャートを実行した場合で、ステップ８１６実行後の入出力端末１４０の画面例１５００を示す。画面例１５００では、実際の測定済み反応過程データ１５０１と、反応モデルデータ１５０２と、実際の分析結果１５０３と、反応モデルデータ１５０２を利用して算出した予測分析結果１５０４とを、既に測定限界に達したことを知らせる表示とともに出力している。

このように、測定終了後に実際の分析結果と予測分析結果を同時に出力することで、技師は再検の必要性の有無を効率的且つ容易に判断することができる。

本実施例では、ステップ８０２において、Ｔ_{ＳＴＡＲＴ}をＲ１またはＲ２投入直後の測光ポイントに設定したが、Ｔ_{ＳＴＡＲＴ}はユーザが任意に設定できる。これにより、試薬投入直後に温度が一時的に変化する場合、Ｔ_{ＳＴＡＲＴ}を反応容器内の温度が一定となる時間に設定することで、より適切なＡ_Ｈ１、Ａ_Ｈ０を算出することが可能となる。

また、本システムは、画面例１３００のように、実際に測定された吸光度が測定限界吸光度に達していない場合でも、反応モデルデータを予測することで、将来的に測定限界吸光度に達することが想定される場合、その時点で、再検のための試料希釈倍率を入出力端末１４０に出力することもできる。

図１６に、入出力端末１４０に試料希釈倍率を出力するときのフローチャート１６００を示す。

まず、ステップ８０１〜８１４までと同様の処理をするステップ１６０１を実行する。

次に、制御部１０１が、判定部１０９を起動し、反応モデルデータＡＢＳ_{Ｈ＿ＭＯＤＥＬ}（ｔ）が、反応終了時間までに測定限界に達するかどうかを判定して、その判定結果を一時記憶装置１０２に保存するステップ１６０２を実行する。

ステップ１６０２において偽と判定された場合、本システムは動作を終了する。

ステップ１６０２において真と判定された場合、制御部１０１が、試料混合比率算出部１０８を起動し、以下の手順で試料の希釈倍率ｒを設定して、一時記憶装置１０２に保存するステップ１６０３を実行する。

以下に、ステップ１６０３の具体的な手順を示す。

まず、希釈倍率ｒで希釈した希釈試料に対して、（数２）をもとに、（数６）で示すような希釈試料の反応モデルデータＡＢＳ_{Ｄ＿ＭＯＤＥＬ}（ｔ）を考える。

ＡＢＳ_{Ｄ＿ＭＯＤＥＬ}（ｔ）＝Ａ_Ｄ０＋Ａ_Ｄ１×（１−ｅｘｐ（−ｋ×ｔ））
…（数６）
このとき、ｋは反応の速度、Ａ_Ｄ０は反応開始時の初期吸光度、Ａ_Ｄ１は目的成分の濃度、を示すパラメータである。また、この反応モデルデータＡＢＳ_{Ｄ＿ＭＯＤＥＬ}（ｔ）は、最適条件下で分析できると仮定したときの各測光ポイントにおける予測測光値と言い換えることができる。このとき、濃度・活性値算出に用いる測光ポイントをＴとすると、（数１）より、希釈倍率ｒに対して、（数７）のような関係が成り立つ。

ａ×ＡＢＳ_{Ｄ＿ＭＯＤＥＬ}（Ｔ）＋ｂ＝
（ａ×ＡＢＳ_{Ｈ＿ＭＯＤＥＬ}（Ｔ）＋ｂ）／ｒ
…（数７）
ここで、ａ、ｂ、ＡＢＳ_{Ｈ＿ＭＯＤＥＬ}（Ｔ）を代入し、ＡＢＳ_{Ｄ＿ＭＯＤＥＬ}（Ｔ）が予め設定された吸光度（ただし、測定限界吸光度を超えない値）に一致するようにｒを設定する。

次に、制御部１０１が、予測測光値算出部１０６を起動し、（数６）のパラメータＡ_Ｄ０、Ａ_Ｄ１を算出して、一時記憶装置１０２に保存するステップ１６０４を実行する。

以下に、ステップ１６０４の具体的な手順を示す。

（数７）より、全測光ポイントのｔに対して、（数８）が得られる。

ＡＢＳ_{Ｄ＿ＭＯＤＥＬ}（ｔ）＝
ＡＢＳ_{Ｈ＿ＭＯＤＥＬ}（ｔ）／ｒ−（１−１／ｒ）×ｂ／ａ
…（数８）
このとき、ＡＢＳ_{Ｈ＿ＭＯＤＥＬ}（ｔ）はｋ、Ａ_Ｈ０、Ａ_Ｈ１をパラメータとした関数であるから、（数８）は（数９）のように示すことができる。

ＡＢＳ_{Ｄ＿ＭＯＤＥＬ}（ｔ）＝Ａ_Ｈ０／ｒ−（１−１／ｒ）×ｂ／ａ
＋(Ａ_Ｈ１／ｒ)×（１−ｅｘｐ（−ｋ×ｔ））
…（数９）
以上から、（数６）と（数９）の右辺同士を比較することにより、Ａ_Ｄ０＝Ａ_Ｈ０／ｒ−（１−１／ｒ）×ｂ／ａ、Ａ_Ｄ１＝Ａ_Ｈ１／ｒを算出することができる。

このように、ステップ１６０４により、希釈試料の反応モデルデータを算出することができる。

次に、制御部１０１が、出力部１１０を起動し、一時記憶装置１０２上の反応過程データＡＢＳ_Ｈ（ｔ）と、反応モデルデータＡＢＳ_{Ｈ＿ＭＯＤＥＬ}（ｔ）と、予測分析結果と、希釈倍率ｒと、希釈倍率ｒで希釈した希釈試料の反応モデルデータＡＢＳ_{Ｄ＿ＭＯＤＥＬ}（ｔ）と、ステップ１６０２の判定結果を、入出力端末１４０に出力するステップ１６０５を実行する。

図１７に、ｔ＝Ｔ_４（４ポイント目まで測定が終了した時点）において、図１６に示すフローチャートを実行した場合で、ステップ１６０５実行後の入出力端末１４０の画面例１７００を示す。

画面例１７００では、ｔ＝Ｔ_４において、実際の測定済み反応過程データ１７０１と、反応モデルデータ１７０２と、反応モデルデータ１７０２を利用して算出した予測分析結果１７０３と、再検時に推奨する希釈倍率１７０４と、希釈倍率１７０４で希釈した試料の反応モデルデータ１７０５とを、将来測定限界に達することを知らせる表示とともに出力している。また、希釈倍率１７０４は、希釈試料の反応モデルデータ１７０５が、濃度・活性値算出に用いる測光ポイント１７０６において測定限界１７０８に達しない吸光度１７０７になるように、設定されている。

このように、吸光度を予測することによって、実際に測定された吸光度が測定限界に達するよりも早いタイミングで、測定限界に達することを技師が視覚的に判断でき、希釈倍率を提示することによって、速やかに適切な再検を実施することができるので、ＴＡＴを短縮し、且つ不要な再検を削減し、再検コストを低減することが可能となる。

本実施例では、測定限界を、測定できる最大濃度を示す吸光度として扱ったが、装置が測定できる最低濃度を示す吸光度として扱ってもよい。

また、本実施例では、ステップ１６０３において、希釈倍率を設定したが、試料混合量の増減量を設定してもよい。また、希釈倍率と試料混合量の増減量を同時に設定してもよい。

また、本実施例では、ステップ１６０３において、ＡＢＳ_{Ｄ＿ＭＯＤＥＬ}（Ｔ）が一致する吸光度（ただし、測定限界吸光度を超えない値）が予め設定されていたが、この値は、分析装置の測定精度が最も高い吸光度または吸光度範囲が設定されていてもよい。これにより、希釈試料の測定精度を向上させることが可能となる。

また、本実施例では、演算によって希釈倍率を算出したが、予測分析結果に対応する希釈倍率を予めテーブルに保持していても良い。

以上詳述したように、本発明における吸光度予測システムにより、化学反応モデルを表現するパラメータを用いて、測定結果を予測し、測定限界を超える高濃度・高活性試料の再検時の希釈倍率をユーザに提供することができるので、不要な再検を削減し、検査コストの低減、ＴＡＴの短縮が可能となる。

本発明の一実施例における吸光度予測システムの構成例を説明する図。 図１における自動分析装置の構成例を説明する図。 ＴＰ（総蛋白）の反応過程データの例を説明する図。 ＴＰの検量線の例を説明する図。 ＬＤＨ（乳酸脱水素酵素）の反応過程データの例を説明する図。 高活性ＬＤＨの反応過程データの例を説明する図。 反応モデルパラメータＤＢの例（Ａ）と測定パラメータＤＢの例（Ｂ）を説明する図。 本発明における入出力端末に予測分析結果を出力するときのフローチャートを説明する図。 エンド法の項目においてステップ８０５（図８中）でＴ_ＥＮＤを設定するときの例を説明する図。 レート法の項目においてステップ８０５（図８中）でＴ_ＥＮＤを設定するときの例を説明する図。 エンド法の項目におけるＡＢＳ_Ｈ（ｔ）とＡＢＳ_{Ｈ＿ＭＯＤＥＬ}（ｔ）の例（Ａ）とそのときの検量線の例（Ｂ）を説明する図。 レート法の項目におけるＡＢＳ_Ｈ（ｔ）とＡＢＳ_{Ｈ＿ＭＯＤＥＬ}（ｔ）の例を説明する図。 図８におけるステップ８１６実行後の入出力端末の画面例（ｔ＝Ｔ_４）を説明する図。 図８におけるステップ８１６実行後の入出力端末の画面例（ｔ＝Ｔ_８）を説明する図。 図８におけるステップ８１６実行後の入出力端末の画面例（測定終了時）を説明する図。 本発明における入出力端末に試料希釈倍率を出力するときのフローチャートを説明する図。 図１６におけるステップ１６０５実行後の入出力端末の画面例（ｔ＝Ｔ_４）を説明する図。

符号の説明

１００…吸光度予測システム、１０１…制御部、１０２…一時記憶装置、１０３…反応モデルパラメータＤＢ読み出し部、１０４…反応過程データ抽出部、１０５…測定パラメータＤＢ読み出し部、１０６…予測測光値算出部、１０７…予測分析結果算出部、１０８…試料混合比率算出部、１０９…判定部、１１０…出力部、１１１…反応モデルパラメータＤＢ、１１２…測定パラメータＤＢ、１２０…検査施設内のネットワーク、１３０…自動分析装置、１３１…反応過程データ記憶装置、１４０…入出力端末、２０１…光源ランプ、２０２…恒温槽、２０３…セル、２０４…試料分注ノズル、２０５ａ…第１試薬分注ノズル、２０５ｂ…第２試薬分注ノズル、２０６…撹拌棒、２０７…分光器、２０８…検知器、２０９…増幅器、２１０…Ａ／Ｄ変換器、３００…ＴＰの反応過程データの例、３０１…ＴＰの濃度を算出するのに利用する化学反応を示す化学反応式、４００…高濃度域で非直線となるＴＰの検量線の例、５００…ＬＤＨの反応過程データの例、５０１…ＬＤＨの活性値を算出するのに利用する化学反応を示す化学反応式、６００…高活性ＬＤＨの反応過程データの例、６０１…時間Ｔ_Ｌ〜Ｔ_Ｒ間の吸光度変化率、６０２…ＬＤＨの測定限界吸光度、６０３…実際の吸光度変化率、７００…反応モデルパラメータＤＢの例、７０１…ＱＣの種類を識別するＱＣＩＤを格納するフィールド、７０２…分析項目名を格納するフィールド、７０３…測定パラメータを識別する測定パラメータＩＤを格納するフィールド、７０４…分析項目の濃度あるいは活性値の測定値を格納するフィールド、７０５…正規化基準値を格納するフィールド（レート法のみ）、７０６…反応モデルパラメータｋを格納するフィールド、７０７…反応モデルパラメータＡ_０を格納するフィールド、７０８…反応モデルパラメータＡ_１を格納するフィールド、７１０…測定パラメータＤＢの例、７１１…測定パラメータを識別する測定パラメータＩＤを格納するフィールド、７１２…分析項目名を格納するフィールド、７１３…Ｒ１の種類を識別する試薬ＩＤを格納するフィールド、７１４…Ｒ２の種類を識別する試薬ＩＤを格納するフィールド、７１５…分析に使用した波長を格納するフィールド、７１６…測定限界吸光度、７１７…レート法において吸光度変化率を計算するのに利用する左側の測光ポイントを示す番号を格納するフィールド、７１８…レート法において吸光度変化率を計算するのに利用する右側の測光ポイントを示す番号を格納するフィールド、７１９…濃度あるいは活性値への換算式を決定する（数１）におけるａを格納するフィールド、７２０…濃度あるいは活性値への換算式を決定する（数１）におけるｂを格納するフィールド、８０１…測定パラメータＤＢを読み出すステップ、８０２…反応開始時の測光ポイントＴ_{ＳＴＡＲＴ}を設定するステップ、８０３… 反応過程データ記憶装置からｔ＝Ｔ_１〜Ｔ_Ｎにおける反応過程データＡＢＳ_Ｈ（ｔ）を抽出するステップ、８０４…「Ｔ_{ＳＴＡＲＴ}＜Ｔ_Ｎ」の真偽を判定するステップ、８０５…定限界吸光度に達する直前の測光ポイントＴ_ＥＮＤを設定するステップ、８０６…Ｔ_{ＳＴＡＲＴ}〜Ｔ_ＥＮＤ間に一定数以上の測光ポイントが含まれているかどうかを判定するステップ、８０７…反応モデルパラメータＤＢを読み出すステップ、８０８…エンド法かレート法かを判断するステップ、８０９…レート法の反応過程データを正規化するステップ、８１０…ΔＡＢＳ_Ｈ（ｔ）を算出するステップ、８１１…ＡＢＳ_{Ｈ＿ＭＯＤＥＬ}（ｔ）の濃度を示すパラメータＡ_Ｈ１を算出するステップ、８１２…ＡＢＳ_{Ｈ＿ＭＯＤＥＬ}（ｔ）の初期吸光度を示すパラメータＡ_Ｈ０を算出するステップ、８１３…ＡＢＳ_{Ｈ＿ＭＯＤＥＬ}（ｔ）を算出するステップ、８１４…予測分析結果を算出するステップ、８１５…ＡＢＳ_{Ｈ＿ＭＯＤＥＬ}（ｔ）が反応終了時間までに測定限界に達するかどうかを判定するステップ、８１６…算出した予測分析結果などを入出力端末に出力するステップ、９００…エンド法の項目においてステップ８０５でＴ_ＥＮＤを設定するときの反応過程データの例１、９１０…エンド法の項目においてステップ８０５でＴ_ＥＮＤを設定するときの反応過程データの例２、１０００…レート法の項目においてステップ８０５でＴ_ＥＮＤを設定するときの反応過程データの例１、１０１０…レート法の項目においてステップ８０５でＴ_ＥＮＤを設定するときの反応過程データの例２、１１００…エンド法の項目におけるＡＢＳ_Ｈ（ｔ）の例、１１０１…１１００の反応モデルデータＡＢＳ_{Ｈ＿ＭＯＤＥＬ}（ｔ）の例、１１０２…検量線、１１０３…測定限界吸光度、１１０４…実際に算出される濃度、１１０５…適切な濃度、１２００…レート法の項目におけるＡＢＳ_Ｈ（ｔ）の例、１２０１…１２００の反応モデルデータＡＢＳ_{Ｈ＿ＭＯＤＥＬ}（ｔ）の例、１２０２…吸光度変化率を計算するのに利用する左側の測光ポイント、１２０３…吸光度変化率を計算するのに利用する右側の測光ポイント、１２０４…測定限界吸光度、１２０５…実際の分析に利用される吸光度変化率、１２０６…適切な吸光度変化率、１３００…ステップ８１６実行後の入出力端末の画面例（ｔ＝Ｔ_４）、１３０１…測定済み反応過程データ、１３０２…反応モデルデータ、１３０３…予測分析結果、１３０４…測定限界吸光度、１３０５…反応開始時間、１４００…ステップ８１６実行後の入出力端末の画面例（ｔ＝Ｔ_８）、１４０１…測定済み反応過程データ、１４０２…反応モデルデータ、１４０３…予測分析結果、１４０４…測定限界吸光度、１４０５…反応開始時間、１４０６…測定限界内最終測光ポイント、１５００…ステップ８１６実行後の入出力端末の画面例（測定終了時）、１５０１…測定済み反応過程データ、１５０２…反応モデルデータ、１５０３…実際の分析結果、１５０４…予測分析結果、１６０１…ステップ８０１〜８１４までと同様の処理をするステップ、１６０２…ＡＢＳ_{Ｈ＿ＭＯＤＥＬ}（ｔ）が反応終了時間までに測定限界に達するかどうかを判定するステップ、１６０３…試料希釈倍率を設定するステップ、１６０４…ＡＢＳ_{Ｄ＿ＭＯＤＥＬ}（ｔ）を算出するステップ、１６０５…設定した希釈倍率などを入出力端末に出力するステップ、１７００…ステップ１６０５実行後の入出力端末の画面例（ｔ＝Ｔ_４）、１７０１…測定済み反応過程データ、１７０２…反応モデルデータ、１７０３…予測分析結果、１７０４…希釈倍率、１７０５…希釈試料の反応モデルデータ、１７０６…濃度・活性値算出に用いる測光ポイント、１７０７…測定限界吸光度に達しない吸光度。

Claims (11)

1. 試料と１種類以上の試薬とを混合した混合液を光学的に測定した測光値の時系列データをもとに前記試料中の目的成分の濃度または活性値を測定する分析装置であって、前記試料の分析が正常に実施されたときの前記時系列データのモデルを決定付ける前記目的成分毎の化学反応に応じた1種類以上の係数と、前記試料の成分の濃度に応じた1種類以上の係数とを含むパラメータを記憶する記憶手段と、実際に測定された１点以上の測光値と前記パラメータを用いて、測定開始から所定の経過時間までの予測測光値を算出する算出手段とを有することを特徴とする分析装置。
2. 請求項１に記載の分析装置であって、前記予測測光値が分析結果の正確さを保証できる範囲を示す測定限界を超えるか否かを判定する判定手段を有することを特徴とする分析装置。
3. 請求項２に記載の分析装置であって、前記パラメータと前記予測測光値を用いて、前記混合液に対する前記試料の最適な混合比率を算出する混合比率算出手段を有することを特徴とする分析装置。
4. 請求項１、２又は３に記載の分析装置であって、前記予測測光値をもとに、前記試料中の目的成分の濃度または活性値を算出して出力する出力手段を有することを特徴とする分析装置。
5. 請求項４に記載の分析装置であって、前記出力手段が、実際に測定された測光値と前記予測測光値とを同時に異なる表現方法で出力することを特徴とする分析装置。
6. 試料と試薬とを混合した混合溶液の化学反応の進行につれて変化する吸光度を所定の時間間隔で複数回測光して、前記吸光度の測光時系列データを測定する手段と、前記測光時系列データをもとに、前記吸光度の一定時間での変化量もしくは変化率を算出し、前記試料中の目的成分の濃度もしくは活性値を求める手段とを有し、かつ、前記試料の分析が正常に実施されたときの前記測光時系列データのモデルを決定付ける前記目的成分毎の化学反応に応じた1種類以上の係数と、前記試料の成分の濃度に応じた1種類以上の係数とを含むパラメータを記憶する記憶手段と、前記測光時系列データと前記パラメータを用いて、前記所定の時間間隔での予測吸光度を算出する算出手段とを設けることにより、前記予測吸光度をもとに前記試料中の目的成分の濃度もしくは活性値を求めるよう構成したことを特徴とする分析装置。
7. 請求項６に記載の分析装置において、前記予測吸光度が分析結果の正確さを保証できる範囲を示す測定限界を超えるか否かを判定する判定手段を有することを特徴とする分析装置。
8. 請求項７に記載の分析装置において、前記パラメータと前記予測吸光度とを用いて、前記混合液に対する前記試料の混合比率を算出する混合比率算出手段を有することを特徴とする分析装置。
9. 請求項６、７又は８に記載の分析装置において、前記予測吸光度をもとに、前記試料中の目的成分の濃度または活性値を算出して出力する出力手段を有することを特徴とする分析装置。
10. 請求項９に記載の分析装置において、前記出力手段が、前記予測吸光度と共に、前記測光時系列データを出力し表示する構成されていることを特徴とする分析装置。
11. 請求項９に記載の分析装置において、前記出力手段が、前記予測吸光度が分析結果の正確さを保証できる範囲を示す前記測定限界を超える場合における前記混合比率を、前記測光時系列データおよび前記予測吸光度のうち少なくとも一つと共に出力し表示するよう構成されていることを特徴とする分析装置。

Images