自動分析装置及び自動分析方法
【課題】試料と試薬を混合し、混合液の経時的な変化を測定する自動分析装置において、反応過程に現れる検体固有の直線範囲を自動的に決定又は予測する技術を提供する。
【解決手段】本発明は、反応過程データを関数で近似し、当該関数に基づいて反応初期及び又は後期の曲線部分を自動的に判定する。曲線部分を除いた直線範囲を検体ごとに決定し、決定した直線範囲の吸光度データを用いて検査値を算出する。また、本発明は、反応過程の途中までに得られた吸光度データに基づいて反応初期の直線開始時刻を自動的に判定し、当該直線終了時刻と事前に予定した直線終了時刻とに基づいて直線範囲を予測し、その結果に基づいて予測値を算出する。
【解決手段】本発明は、反応過程データを関数で近似し、当該関数に基づいて反応初期及び又は後期の曲線部分を自動的に判定する。曲線部分を除いた直線範囲を検体ごとに決定し、決定した直線範囲の吸光度データを用いて検査値を算出する。また、本発明は、反応過程の途中までに得られた吸光度データに基づいて反応初期の直線開始時刻を自動的に判定し、当該直線終了時刻と事前に予定した直線終了時刻とに基づいて直線範囲を予測し、その結果に基づいて予測値を算出する。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、血液、尿その他の生体サンプルを定性又は定量分析する自動分析装置及び方法に関し、例えば臨床検査時の反応をモニタリングする機能を有する自動分析装置及び方法に関する。
【背景技術】
【０００２】
臨床検査用の自動分析装置は、試料と試薬を一定量分注して反応させた後、一定時間にわたり反応液の吸光度を測定し、測定結果に基づき測定対象物質の検査値（濃度や活性値）等を求める。
【０００３】
臨床検査の測定法の一つであるレート法は、主に試料に含まれる酵素成分の活性値を測定するために使用される。測定方法は、試薬として一定量の基質を添加し、酵素が基質を消費して変化する要素を測定する。通常、試薬には十分量の基質が含まれている。このため、試料と試薬の反応が正常に行なわれていれば、一般に反応は時間変化に対して吸光度が一定量ずつ直線的に変化する。レート法ではこの直線の傾き（反応速度）から測定対象物質の活性値を求めている。
【０００４】
しかしながら、測定項目、検体の濃度、撹拌の状況、反応温度などの理由から、反応速度が一定（直線）になるまでにラグタイム部と称される曲線部をもつことがある。また、高活性検体等の場合には反応が早く進むことがある。この場合、測定中に試薬の基質が十分量でなくなってしまい、後半の反応速度が一定でなくなることがある。すなわち、後半の反応速度が非直線状（曲線）に変化することがある。
【０００５】
このように反応過程データには、反応の開始部と後部に曲線部分が出現することがある。このため、従来装置の場合には、検査値の換算に用いる時間範囲と吸光度限界値(Absリミット)を項目ごとに定める手法を採用する。また、従来装置の中には、リニアリティチェックという直線性のチェック機能を備えるものもある。リニアリティチェックでは、一定の測光範囲について、その前半と後半における吸光度変化量の差が算出され、その差が指定したリニアリティチェック値を超える場合に、直線でないと判断される。直線でないと判断された場合、従来装置は「異常」を通知する。
【０００６】
検査値の換算に用いる時間範囲（直線範囲）の決定方法には、特許文献１及び２に開示された方法がある。特許文献１には、隣接する２つの測光ポイント間毎の吸光度差を算出し、それらのうちで一番大きな吸光度差(dmax)を示す測光ポイント間を基準点に定め、その前後の測定範囲の中で吸光度差が基準点の吸光度差(dmax)の80％以上となる区間を直線範囲に決定する方法が開示されている。特許文献２には、換算に用いる範囲幅を一定に設定し、先頭点の測光ポイントから１ポイントずつずらしながら、各範囲内に対応する吸光度の相関係数を最小二乗計算法により算出し、一番相関の良かった範囲を直線範囲に決定する方法が示されている。
【０００７】
一方、近年、病院における臨床検査の分野では、患者検体の測定結果をできる限り早急に報告することが求められている。自動分析装置の測定時の反応時間は１０分ほどであるが、単に測定時間を短くするだけでは、正確な結果が得られなくなってしまう。例えば反応が不完全な時点で測定を終了し、その時点の吸光度から検査値を測定することになった場合には、正確な結果が得られなくなってしまう。そこで、反応終了時の吸光度を予測して検査値に換算する方法が考えられる。反応終了時の吸光度を予測する方法としては、特許文献３が公知である。特許文献３には、（数１）に示す近似式を用いて反応過程データを近似し、レート法に関しては、予め定めた２点間又は近似式を時間で微分した後に予め定めた１点から直線の傾きを求める方法が示されている。なお、（数１）において、ｔは時間、ｙは吸光度を表し、Ａ、Ｂ、Ｋはパラメータである。
（数１） y＝A ＋ (B−A)／eKt
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】特開昭６４−６８６４２号公報
【特許文献２】特開平１０−１８５６８７号公報
【特許文献３】特開平６−１９４３１３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００９】
例えばレート法の主試薬添加後の反応過程データの例を図２の（ａ）及び(ｂ）に示す。横軸１１０は時間の経過を示し、縦軸１２０は吸光度を示す。前述したように、従来装置は、検査項目毎に時間範囲（直線範囲）を定めている。すなわち、直線開始時刻と直線終了時刻を事前に定めている。図中の範囲１６０は、ある検査項目について予め定められた時間範囲を示す。
【００１０】
図２の（ａ）の符号１３０、１４０、１５０は、同じ項目の異なる検体の吸光度データを示し、範囲１３５、１４５、１５５は各吸光度データのラグタイム部を示す。なお、符号１７０は予め定められた吸光度の限界値を示す。
【００１１】
符号１４０で示す吸光度データのように、本来の直線部（本来の検査対象）だけが予め定められた時間範囲１６０に含まれる場合がある一方で、符号１５０で示す吸光度データのように、ラグタイム部１５５の一部が予め定められた時間範囲１６０内に含まれる場合もある。
【００１２】
この符号１５０で示す吸光度データの場合、従来装置のリニアリティチェック機能では、反応過程データに異常があると判断され、その検体は再検査の対象になる。しかし、符号１５０で示す吸光度データの場合でも、その直線部を正確に判断できさえすれば、再検査は必要ない。また、符号１５０で示す吸光度データのように、その直線部とラグタイム部１５５の直線性にそれほどの差がない場合には、リニアリティチェックにおいて異常（ラグタイム部１５５が直線範囲の測定期間に含まれること）を検出することができず、曲線部の情報を含んだまま直線の傾きが算出されてしまう場合もある。この場合、直線の傾きが精度良く算出されないので検査値の誤差が大きくなる可能性がある。
【００１３】
また、符号１３０で示す吸光度データの場合、予め定められた範囲１６０で測定結果は直線を示している。しかし、符号１３０で示す吸光度データは、範囲１６０の後半部分で吸光度の限界値を超える。このため、符号１３０で示す吸光度データの検体も再検査の対象となる。しかし、吸光度の限界値を超えた後も、当該検体の吸光度が直線的に変化することが保証できれば、吸光度の限界値を超えた以降のデータも検査値の換算に用いることが可能となり、再検査を減らすことができる。
【００１４】
図２の（ｂ）に示す符号１８０も、ある検体のある検査項目に関する吸光度データである。この検体の場合、範囲１６０の後半に曲線部分が現れる。このため、この検体も、従来装置のリニアリティチェック機能により、再検査の対象に分類されてしまう。しかし、この検体の場合でも、符号１８０で示す吸光度データの直線範囲だけを抽出することができれば検査値の換算が可能となり、再検査を減らすことができる。
【００１５】
因みに、特許文献１のリニアリティチェック機能を符号１５０で示す吸光度データに適用する場合、一番大きな吸光度差を示す測光ポイントは、図中に円で囲んで示す時点となる。ところが、この時点は、ラグタイム部１５５と直線部の境である。特許文献１のリニアリティチェック機能では、円で囲んだ測光ポイントを基準として、前後の期間を直線範囲に定める。このことは、本来の直線部でないラグタイム部１５５が、検査値の換算に用いられることを意味する。すなわち、検査値の誤差が大きくなる可能性がある。
【００１６】
また、特許文献２のリニアリティチェック機能の場合にも、換算に使用する範囲幅が一定であるので、範囲幅が長く設定されると直線部でない範囲が検査値の換算に用いられる可能性がある。このため、検査値の誤差がやはり大きくなる可能性がある。その一方で、範囲幅が短く設定された場合、検査値の計算に用いる吸光度データが少なくなり、やはり誤差が大きくなる可能性がある。
【００１７】
同様に、反応終了時の吸光度を予測する特許文献３の場合にも、予め定めた範囲（開始時点と終了時点で与えられる）について、反応過程データの直線の傾きを算出することになるため、予め定めた範囲が実測定データの直線範囲の外に位置する場合には、正確な予測値が得られない可能性がある。
【課題を解決するための手段】
【００１８】
そこで、本発明は、試料と試薬を混合し、混合液の経時的な変化を測定する自動分析装置において、(a) 試料と試薬の反応過程から複数の測定点データを取得する測定点データ取得部と、(b) 測定点データを処理するデータ処理部と、(c) データ処理部で使用する第１の近似式を格納した記憶部と、(d) データ処理部の処理結果を出力する出力部とを有するものを提案する。ここで、データ処理部では、記憶部に格納された第１の近似式を複数の測定点データに近似する処理と、近似処理の結果得られる第２の近似式に基づいて前記反応過程の直線範囲を決定する処理とを実行する。
【００１９】
なお、データ処理部は、所定時刻までに取得された複数の測定点データを用いて第２の近似式を算出し、当該第２の近似式に基づいて反応過程の直線範囲を推定又は予測する機能を有することが望ましい。
【００２０】
また、本発明は、試料と試薬を混合した混合液の経時的な変化を、自動分析装置を用いて測定する自動分析方法において、(a) 自動分析装置が、試料と試薬の反応過程から複数の測定点データを取得する処理と、(b) 自動分析装置が、記憶部から読み出した第１の近似式を複数の測定点データに近似する処理と、(c) 自動分析装置が、近似処理の結果得られる第２の近似式に基づいて反応過程の直線範囲を決定する処理とを実行するものを提案する。やはり、所定時刻までに取得された複数の測定点データを用いて第２の近似式を算出し、当該第２の近似式に基づいて反応過程の直線範囲を推定又は予測する処理を有することが望ましい。
【００２１】
なお、本発明は、直線範囲の精度指標を算出し、処理結果として出力部に出力する機能を有することが望ましい。
【発明の効果】
【００２２】
本発明により、検査値の換算に用いる直線範囲を検体毎に定めることが可能となる。結果的に検体毎に正確な検査値が得られる。また、従来装置に比して再検査の対象となる検体数を減らすことができる。さらに、所定時刻までに取得された複数の測定点データを用いて第２の近似式を算出し、当該第２の近似式に基づいて反応過程の直線範囲を推定する場合には、直線範囲の決定までに要する時間を短縮することができる。また、直線範囲の精度指標を出力する機能を有する場合には、検査値や予測値の精度保証が可能となる。
【図面の簡単な説明】
【００２３】
【図１】第１の実施例に係る処理フローチャート。
【図２】レート法の反応過程データ例を説明する図。
【図３】自動分析装置の概略構成例を示す図。
【図４】制御部の内部構成例を示す図。
【図５】検査項目と、試薬の種類と、最適な近似式の関係を記述したテーブル例を示す図。
【図６】ラグタイム部を数値化する方法を説明する図。
【図７】反応後期の曲線部分を数値化する方法を説明する図。
【図８】反応後期の曲線部分を判定する処理フローチャート。
【図９】検体と検査値との関係を記述したテーブル例を示す図。
【図１０】第１の実施例を用いて決定される直線範囲例を説明する図。
【図１１】第２の実施例に係る処理フローチャート。
【図１２】近似式パラメータの分布例を示す図。
【図１３】近似式パラメータ比の分布例を示す図。
【図１４】第２の実施例を用いて決定される直線範囲例を説明する図。
【図１５】第３の実施例に係る処理フローチャート。
【図１６】予測値と検査値と精度指標を記述したテーブル例を示す図。
【発明を実施するための形態】
【００２４】
以下、図面に基づいて、本発明の実施例を説明する。なお、後述する装置構成や処理動作の内容は発明を説明するための一例である。本発明には、後述する装置構成や処理動作に既知の技術を組み合わせた発明や後述する装置構成や処理動作の一部を既知の技術と置換した発明も含まれる。
【実施例】
【００２５】
［実施例１］
以下、第１の実施例に係る自動分析装置の装置構成及び処理動作を、図面を参照して詳細に説明する。図３に、本発明に係る分析機能を実装する生化学自動分析装置の概略構成を示す。
【００２６】
生化学自動分析装置は、サンプルディスク１、試薬ディスク２、反応ディスク３、反応槽４、サンプリング機構５、ピペッティング機構６、攪拌機構７、測光機構８、洗浄機構９、コンピュータ（PC）１０、記憶装置１２、制御部１３、圧電素子ドライバ１４、攪拌機構コントローラ１５、試料容器１６、円形ディスク１７，１９、試薬ボトル１８、保冷庫２０、反応容器２１、反応容器ホルダ２２、駆動機構２３、プローブ２４，２７、支承軸２５，２８、アーム２６，２９、固定部３１、ノズル３３、上下駆動機構３４を有している。記憶装置１２には、分析パラメータ、各試薬ボトルの分析可能回数、最大分析可能回数、キャリブレーション結果、分析結果等が記憶される。
【００２７】
生化学自動分析装置における試料の分析は、サンプリング、試薬分注、撹拌、測光、反応容器の洗浄、濃度換算等のデータ処理の順番に実行される。
【００２８】
サンプルディスク１は、制御部１３によりコンピュータ１０を介して制御される。サンプルディスク１上には、複数の試料容器１６が円周上に並んで設置される。試料容器１６は、分析順に従ってサンプリングプローブ２４の下まで移動される。試料容器１６中の検体は、検体サンプリング機構５に連結された試料用ポンプにより反応容器２１の中に所定量ずつ分注される。
【００２９】
試料が分注された反応容器２１は、反応槽４の中を第１試薬の添加位置まで移動される。移動された反応容器１６には、試薬分注プローブ２７に連結された試薬用ポンプ（図示せず）により試薬容器１８から吸引された試薬が所定量だけ加えられる。第１試薬が添加された後、反応容器２１は撹拌機構７の位置まで移動され、最初の撹拌が行われる。このような試薬の添加と撹拌が、例えば第１〜第４試薬について行われる。
【００３０】
内容物が撹拌された反応容器２１は、光源から発せられる光束中に配置される。光束の一部は反応容器２１を通過し、一部は内容物で吸収される。吸収の度合いは、例えば多波長光度計で構成される測光機構８により検知される。測光機構８は、時間の経過と共に検出された検体に関する吸収の度合いを測定点データ（吸光度信号）として制御部１３に出力する。制御部１３は、後述するデータ処理を通じ、測定点データの直線範囲を決定し、その後、検体の検査値（直線範囲の傾き）を算出する。算出された検査値（データ）は記憶装置１２に記憶され、コンピュータ１０に付属する表示装置に表示される。測光の終了した反応容器２１は、洗浄機構９の位置まで搬送され、洗浄の後、次の分析に使用される。
【００３１】
次に、制御部１３において実行される直線範囲の決定処理手順を説明する。図１に、当該処理手順を示す。なお、図１は、制御部１３内で実行される直線範囲の決定プログラムに関わる部分のみを示している。また、図４に、直線範囲の決定プログラムを通じて実現される制御部１３の機能ブロックを示す。なお、図４は、当該決定プログラムをハードウェアの観点から表している。図４に示す制御部１３は、吸光度データ取得部４１０、近似式パラメータ算出部４２０、ラグタイム部判定部４２５、後半曲線部判定部４３５、直線範囲内吸光度データ数判定部４４０、検査値算出部４４５、これらを相互に接続するデータバス４５０で構成される。データバス４５０を通じ、各部は、相互にデータを受け渡すことができる。なお、制御部１３を構成する機能ブロックは、制御部１３以外のハードウェアやＣＰＵで構成しても良い。勿論、制御部１３を構成する機能ブロックは、同一のＣＰＵ内にソフトウェアモジュールとして実装しても良い。
【００３２】
まず、ある検体に対し、ある検査項目の測定が開始される。この測定開始と同時に、ステップＳ５において、近似式パラメータ算出部４２０が、直線範囲の決定に適した近似式（第１の近似式）を選択する。選択対象としての吸光度の時間変化を表す複数の近似式は、記憶装置１２に予め記憶されているものとする。近似式パラメータ算出部４２０は、これら複数の近似式の中から検査項目と試薬の組み合わせに対応した最適な近似式を選択的に読み出す。
【００３３】
この実施例の場合、（数２）〜（数５）に示す関数を選択可能な近似式とする。なお、各関数において、ｔは時刻、ｘは吸光度を表すものとする。また、ａ、ｂ、ｃ、ｄ、ｅ、ｋ、ｐ、ｑ、ｒ、ｕ、ｖ、ｗはパラメータである。
（数２） x＝a * t ＋ b ＋ c * exp(-k * t)
（数３） x＝a * t ＋ b ＋ e ／ (t ＋ d)
（数４） x＝a * t ＋ b ＋ w ／ {exp(u * t) ＋ v}
（数５） x＝a * t ＋ b ＋ p * log{1 ＋ q * exp(r * t)}
【００３４】
例示した４つの関数は、いずれも、時間に比例した成分と、定数成分と、時間変化率が異なる非直線成分で構成される。勿論、これ以外の関数を近似式として用意することも可能である。
【００３５】
なお、近似式パラメータ算出部４２０による近似式の選択は、検査項目と試薬の組み合わせに基づいて自動的に実行される手法を採用しても良いし、ユーザーが自由に選択する手法を採用しても良い。前者の機能の実現には、例えば検査項目と試薬の組み合わせ毎に最も適した近似式をテーブルとして記憶装置１２に記憶しておけば良い。近似式パラメータ算出部４２０は、このテーブルを検査項目と試薬の組み合わせによって検索し、組み合わせに対応した最適な近似式を選択する。図５に、この種のテーブル例を示す。テーブル５００は、列５１０、５２０及び５３０で構成される。列５１０には検査項目が記述され、列５２０には試薬の種類が記述され、列５３０には検査項目と試薬の種類に対応付けられた最適な近似式の種類が記述されている。
【００３６】
この実施例の場合、近似式パラメータ算出部４２０は、検査項目と試薬の組み合わせに基づいてテーブル５００を検索し、直線範囲の決定に最適な近似式を選択する。なお、テーブル５００に記憶された対応関係の内容は、ユーザーが変更可能な構成としても良い。
【００３７】
吸光度は、時間の経過と共に複数回測定される。次のステップＳ１０では、吸光度データ取得部４１０が、１回の測定又は複数回の測定平均の吸光度データを、測光機構８から入力する。すなわち、吸光度データは制御部１３に入力される。試薬と検体との反応に伴う色調変化に吸光度が大きく変化する波長（主波長）の光と、吸光度が殆ど変化しない波長（副波長）の光の２波長光を用いる測定方式においては、主波長光の吸光度と、副波長光の吸光度との差を、吸光度データとして入力する。
【００３８】
ステップＳ１５では、吸光度データ取得部４１０が、主反応試薬が添加されたかどうかを判定する。主反応試薬が添加されていないと判定された場合、処理プロセスは、ステップＳ１０に戻り、次の吸光度データの入力が行われる。主反応試薬の添加が判定されるまで、以上の判定動作が繰り返し実行される。ここで、主反応試薬とは、複数の試薬を用いる反応において、主たる吸光度変化をもたらす試薬（通常は最終の試薬）を指す。主反応試薬が添加された場合には、処理プロセスはステップＳ２０に移行する。
【００３９】
ステップＳ２０では、入力された吸光度データを、吸光度データ取得部４１０が記憶装置１２に記憶する。
【００４０】
ステップＳ２５では、吸光度データ取得部４１０が、最後の吸光度データが記憶されたか否かを判定する。すなわち、この実施例では、直線範囲の決定前に、予め定めた反応期間に対応する全ての吸光度データが記憶装置１２に記憶されることが必要とされる。ここで、最後の吸光度データが記憶装置１２に記憶されていないと判定された場合、処理プロセスはステップＳ１０に戻る。このループ動作（吸光度データの入力と記憶）は、必要なデータ数が記憶装置１２に記憶されるまで繰り返し実行される。必要なデータ数の蓄積が吸光度データ取得部４１０により判定された場合、処理プロセスはステップＳ３０に進む。
【００４１】
ステップＳ３０では、ステップＳ５において選択された近似式（第１の近似式）によって表される吸光度の時間変化と、実際の吸光度の時間変化との差がなるべく小さくなるように、近似式パラメータ算出部４２０が、数式中のパラメータの値を算出する。具体的には、測定結果として記憶された吸光度データと、近似式により算出される対応時点の吸光度データとの二乗誤差がなるべく小さくなるように、数式中のパラメータ値が算出される。その際、各近似式の最終項（（数２）〜（数５）の場合であれば、（a＊t ＋ b）以降の非直線成分）においては、時刻ｔが増大するとパラメータ値が“０（ゼロ）”に漸近するように制限する。さらに、各近似式の最終項においては、反応初期の吸光度データに近似式をフィッティングさせるようにパラメータの初期値を設定する。
【００４２】
パラメータ値の算出には、既存の最小二乗計算方法が利用可能である。この他、様々な形式の数式に対応可能な方法として、例えば最急降下法により二乗誤差が最小となるパラメータ値を算出する方法もある。この明細書では、各項のパラメータ値を測定結果（吸光度データ）に応じて最適化した後の近似式を第２の近似式と呼び、検査項目と試薬の組み合わせに応じて選択される第１の近似式と区別する。
【００４３】
更に、ステップＳ３０において、近似式パラメータ算出部４２０は、各時刻について、第２の近似式により算出される吸光度（近似値）と実際に測定された吸光度（実測値）との差（誤差）を算出する。
【００４４】
ステップＳ３５では、近似式パラメータ算出部４２０が、第２の近似式を規定するパラメータ値と、各時刻について算出された近似値と実測値との差（誤差）を記憶装置１２に記憶する。
【００４５】
次のステップＳ４０では、ラグタイム部判定部４２５が、反応初期の吸光度が曲線的に変化するラグタイム部分を第２の近似式に基づいて計算する。図６に、ラグタイム部の判定方法を説明する。
【００４６】
図６において横軸１１０は反応開始からの経過時間を示し、縦軸１２０は吸光度を示す。曲線２１０は、第２の近似式により求められる吸光度変化の近似曲線を示す。直線２２０は曲線２１０が漸近する直線である。また、横軸１１０上の点２３０は、曲線２１０が直線２２０に十分漸近する時刻(Ｔｌとする)を示し、横軸１１０の０〜Ｔlの範囲がラグタイム部に相当する。ラグタイム部の終了時刻が直線開始時刻に相当するので、点２３０に対応する時刻が直線開始時刻となる。
【００４７】
ここで、十分漸近した時刻とは、例えば微小な値εを予め定めておき、曲線２１０と直線２２０との差がε以下になった時刻として定義する。εは一定値としても良いし、初期吸光度や吸光度の変化幅に応じて設定しても良い。例えば初期吸光度に定数を乗じた値、又は、初期吸光度と最終吸光度の差に定数を乗じた値をεとしても良い。また、十分漸近した時刻は微小な値δを定めておき、曲線２１０と直線２２０の傾きの差がδ以下になった時刻として定義しても良い。この場合、δは一定値としても良いし、直線２２０の傾きに応じて設定しても良い。例えば直線２２０の傾きに定数を乗じた値をδとしても良い。
【００４８】
ステップＳ４５では、ラグタイム部判定部４２５が、算出した時刻を直線開始時刻として記憶装置１２に記憶する。
【００４９】
次のステップＳ５５では、後半曲線部判定部４３５が、反応過程データの後半に現れる曲線部（以下、「後半曲線部」という。）を判定する。図７と図８を用い、後半曲線部の判定方法を説明する。
【００５０】
図７の上図の横軸１１０は反応開始からの経過時間を表し、縦軸１２０は吸光度を示す。曲線２１０は第２の近似式により求められる吸光度変化の近似曲線である。記号１８０は検体の吸光度データを示す。図７の下図の横軸１１０は反応開始からの経過時間を示し、縦軸７００は誤差を示す。ここでの誤差とは、実測値と近似値の差分値として与えられる。例えば曲線７１０は実測値から近似値を引いた値の時間変化を示す。反応過程データが減少曲線の場合は、例えば求める誤差は近似値から実測値を引いた値にする。この場合も図７の下図に示すような図となる。（数２）〜（数５）の近似式には後半曲線部を表す式が含まれていないので、記号１８０のように後半が曲線になった場合、近似曲線２１０は図７の上図に示すように記号１８０の反応過程曲線と交差する。よって、図７の下図に示すように、近似値と実測値の差（誤差）の時間変化を図示すると、誤差は曲線の開始時刻で極大を示すことになる。この時刻を、この実施例では、直線終了時刻とする。
【００５１】
図８を用い、後半曲線部判定部４３５で実行される処理動作を説明する。まず、ステップＳ８００において、後半曲線部判定部４３５は、記憶装置１２に記憶された誤差を読み出し、誤差が予め定めた許容範囲内か否かを判定する。判定に使用する許容範囲は、記憶装置１２に予め記憶されているものとする。
【００５２】
誤差が許容範囲内の場合、後半曲線部判定部４３５は、ステップＳ８１０に進む。ステップＳ８１０において、後半曲線部判定部４３５は、反応過程データの後半部分に曲線部がないと判定する。この場合、直線終了時刻は反応終了時刻となる。ステップＳ８００で誤差が許容範囲外と判定された場合、後半曲線部判定部４３５は、ステップＳ８２０に進む。
【００５３】
ステップＳ８２０において、後半曲線部判定部４３５は、図７の下図に示すような誤差の時間変化を監視し、誤差が極大値を示す時刻７２０（Ｔｅとする）を検出する。
【００５４】
次のステップＳ８３０において、後半曲線部判定部４３５は、時刻７２０以降の誤差の分布が単調減少か否かを判定する。単調減少するか否かは、例えば隣り合う時間の誤差の差分の符号が一定か否かで判定することができる。単調減少を示していた場合、図７の上図に示すように後半部分に曲線部があることを示す。従って、後半曲線部判定部４３５はステップＳ８４０に進み、反応過程データの後半に曲線部があると判定する。この場合、直線終了時刻は曲線開始時刻となる。すなわち、図７の時刻７２０となる。ステップＳ８３０で単調減少とならないと判定された場合、後半曲線部判定部４３５は、ステップＳ８５０に進んで測定の「異常」を示す信号を出力する。「異常」を示す信号の出力方法には、アラームを鳴らす方法、図３のコンピュータ１０に付属する表示装置に表示する方法等が考えられる。
【００５５】
図１の説明に戻る。ステップＳ５５の処理が終了すると、後半曲線部判定部４３５は、ステップＳ６０に進む。ステップＳ６０において、後半曲線部判定部４３５は、判定結果として得られた直線終了時刻を記憶装置１２に記憶する。
【００５６】
ステップＳ７０では、直線範囲内吸光度データ数判定部４４０が、記憶装置１２に記憶された直線開始時刻と直線終了時刻とから直線範囲を決定する。
【００５７】
続くステップＳ７５では、直線範囲内吸光度データ数判定部４４０が、直線範囲内に出現する吸光度データの数が予め定めた閾値以上か否かを判定する。なお、直線式は、例えば（数６）に示すような１次式で表される。（数６）を用いて直線の傾きを求める場合に必要な吸光度データ数の範囲は２個以上である。ただし、（数６）において、ｔは時刻、ｘは吸光度、ａ、ｂはパラメータを示す。
（数６）x＝a * t ＋ b
【００５８】
また、信頼できる検査値を算出するために必要な最低限の吸光度データ数を予め実験から定めておいても良い。
【００５９】
ところで、ステップＳ７５の判定処理において、直線範囲内の吸光度データ数が閾値より小さいと判定された場合、直線範囲内吸光度データ数判定部４４０は、ステップＳ９０に進み、測定が「異常」であると判定する。「異常」を示す信号の出力方法には、アラームを鳴らす方法、図３のコンピュータ１０に付属する表示装置に表示する方法等が考えられる。
【００６０】
一方、ステップＳ７５の判定処理において、直線範囲内の吸光度データ数が閾値以上と判定された場合、直線範囲内吸光度データ数判定部４４０は、処理プロセスをステップＳ８０に進める。ステップＳ８０においては、検査値算出部４４５が、直線範囲内の吸光度データと近似する直線式を算出し、その直線式から検査値換算に用いる傾きを算出する。
【００６１】
ここでの直線式の算出（傾きを算出する処理に対応する）には、例えば既存の最小二乗計算方法等を利用する。検査値算出部４４５は、直線範囲内の吸光度データに近似する（数６）で与えられる直線を算出する。すなわち、直線の傾きを算出する。また、検査値算出部４４５は、傾きから検量線に基づいて検査値を算出することもできる。検量線データは予め記憶装置１２に記憶しておく。直線式を算出した直線範囲と直線式パラメータ、算出された検査値は記憶装置１２に記憶される。
【００６２】
次のステップＳ８５において、検査値算出部４４５は、算出された検査値を出力する。検査値の出力方法には、図３のコンピュータ１０に付属する表示装置に表示する等が考えられる。表示例を図９に示す。図９に示す測定結果の表示画面は、列９００から９６０で構成される。列９００には検体IDが記述され、列９１０には検査値が記述される。列９０５には直線範囲の決定に使用した近似式（第１の近似式）が記述される。列９２０には検査値の換算に用いた直線範囲の開始時刻（直線開示時刻）が記述される。列９３０には検査値の換算に用いた直線範囲の終了時刻（直線終了時刻）が記述される。列９４０には検査値の換算に用いた直線範囲内に出現する吸光度データ数が記述される。列９５０，列９６０には検査値の換算に用いた直線範囲とその吸光度データ数とで検査値を算出した場合の精度指標が記述されている。精度指標には、例えばCV値を用いる。CV値はキャリブレーション用の検体等、検査値が明確になっている検体を用いて実験を行い、予め定めておけば良い。実験では、様々な直線開始時刻と直線終了時刻の組み合わせ(吸光度データ数)についてCV値(分散や標準偏差等ばらつきを表す値)を算出しておけば良い。精度指標の他の例としては，例えば吸光度データ数の平均値とばらつき(分散や標準偏差等)を用いる。吸光度データ数の平均値とばらつきは，直近に測定した一般(患者)検体やキャリブレーション用検体などの結果から，項目ごとに算出すれば良い。なお，算出に用いる検体，検体数，期間はユーザーが自由に選択できるようにすれば良い。通常、直線範囲内の吸光度データ数が多いほど、検出値の信頼性は一般的に高くなる。
【００６３】
以上説明した第１の実施例においては、図１に示す処理の全てを制御部１３で実行する例を説明した。しかし、制御部１３以外の処理装置を使用して、同様の処理を実行することも可能である。例えばコンピュータ(PC)１０内で実行されるソフトウェア処理として図１の処理を実行することも可能である。また、記憶装置１２としてコンピュータ(PC)１０内部の記憶装置を使用することも可能である。
【００６４】
以上説明した第１の実施例により直線範囲を決定した例を、図１０の（ａ）及び(ｂ）に示す。（ａ）及び（ｂ）の横軸１１０は反応開始からの経過時間を示し、縦軸１２０は吸光度を示す。時刻１０００は、本実施例を通じて決定された直線開始時刻を示し、時刻１０２０は本実施例を通じて決定された直線終了時刻を示す。範囲１０３０は本実施例により決定された直線範囲を表す。なお、図１０の（ａ）はラグタイムのある反応過程データの例であり、図１０の(ｂ）は反応過程の後半に曲線のある反応過程データの例を示している。どちらの場合も、曲線部を除いた直線部だけを直線範囲と決定できることが確認できた。
【００６５】
以上の通り、本実施例で説明した直線範囲の決定機能の採用により、検査値の換算に用いる直線範囲を検体ごとに決定することができる。そのため、従来装置に比して、精度の高い検査値を得ることが可能となる。同時に、再検査を減らすことが可能となる。更に、検査値換算に用いた直線範囲、吸光度データ数、参考CV値や参考吸光度データ平均値等の精度指標も出力することができるため、検体ごとの検査値の精度保証が可能となる。
【００６６】
［実施例２］
次に、第２の実施例に係る自動分析装置を、図面を参照しながら説明する。本実施例の場合も、自動分析装置は生化学自動分析装置であるものとする。従って、その装置構成は、第１の実施例と同じである。すなわち、図３に示す装置構成を有している。また、制御部１３以外の動作は、第１の実施例と同じである。従って、制御部１３以外の詳細な説明は省略する。
【００６７】
以下、制御部１３の処理動作を、本実施例に特有の処理動作を中心に説明する。本実施例の場合、制御部１３は、図１１に示す処理手順に基づいて直線範囲を決定する。なお、図１１のうち図１に示した処理と同じ処理を行う処理ステップには同じ符号を付して示す。
【００６８】
第２の実施例の場合、制御部１３は、反応過程の全ての吸光度データが測定されるのを待つことなく、吸光度データの測定動作と並行的に直線範囲の開始点の決定処理を開始する。
【００６９】
図１１に示す処理プロセスのうちステップＳ５からステップＳ２０までの処理は、図１におけるステップＳ５からステップＳ２０までの処理と同一である。
【００７０】
次のステップＳ１１００では、近似式パラメータ算出部４２０が、現在までに取得された吸光度データに基づいて近似式パラメータとパラメータ変化値を算出する。近似式パラメータの算出方法は、図１に示すステップＳ３０と同一である。ただし、第１の実施例では、反応過程の最後の吸光度データが取得されてから、取得した全ての吸光度データを用いて直線範囲を算出したのに対し、この第２の実施例の場合には、吸光度データが取得される度に近似式パラメータとパラメータ変化値を順次算出する。
【００７１】
ここで、パラメータ変化値とは、近似式パラメータの時間変化を数値化した値のことであり、算出には直線部の傾きを表す近似式パラメータを使用する。例えば第２の近似式が（数２）〜（数５）のいずれかの形式で与えられる場合、近似式パラメータａを使用する。パラメータの変化を数値化する方法には、様々な方法が利用可能である。例えば一回前に算出されたパラメータ値との差を算出する方法、一回前に算出されたパラメータとの比を算出する方法等を利用できる。なお、この実施例では直線部の傾きを与えるパラメータにのみ着目しているが、その他のパラメータ、例えば非直線部のパラメータの変化に着目しても良い。また、直線部と非直線部の両方のパラメータの変化に着目しても良い。
【００７２】
ステップＳ１１１０では、近似式パラメータ算出部４２０が、算出した近似式パラメータとパラメータ変化値を記憶装置１２に記憶する。
【００７３】
次のステップＳ４０とステップＳ４５は、図１に示すステップＳ４０とステップＳ４５と同じである。この実施例の場合、パラメータ変化値が許容範囲以下となるタイミングをラグタイム部の終了時点と判定する。
【００７４】
次のステップＳ１１２０では、後半曲線部判定部４３５が、算出されたパラメータ変化値が予め定めた範囲内か否かを判定し、反応過程データの後半曲線部の判定を行う。図１２と図１３を用いて後半曲線部の判定方法を説明する。
【００７５】
図１２は、直線の傾きを表す近似式パラメータとして、例えば（数２）〜（数５）のパラメータａを用いた場合の経時変化を示す。横軸１１０は反応開始からの経過時間を示し、縦軸１２００は近似式パラメータａの値を示す。記号１２１０は、各時刻に算出された近似式パラメータａの値を示している。近似式パラメータａがほぼ一定の間は反応過程データが直線上に位置することが分かり、近似式パラメータａが一定値から外れると反応過程データが曲線に沿って変化していることが分かる。
【００７６】
例えばパラメータ変化値に近似式パラメータａを用いる場合、後半曲線部判定部４３５は、近似式パラメータａが一定と判断された場合、後半曲線部判定部４３５は、直線と判断してステップＳ２５に進む。一方、近似式パラメータａが変化したと判定された場合、後半曲線部判定部４３５は、直線範囲が終了したと判断してステップＳ１１３０に進む。
【００７７】
図１３は、図１２に示すパラメータの分布に基づいて、一回前に算出されたパラメータとの比を、パラメータ変化値として数値化する場合を示す。横軸１１０は反応開始からの経過時間を示し、縦軸１３００はパラメータ比を示す。記号１３１０は、各時刻の一回前に算出されたパラメータとの比の値を示す。図１３に示すように、反応過程データが直線を示す場合は、パラメータ比は１に近い値となり、直線を示さなくなった場合にはパラメータ比は１から離れた値となる。
【００７８】
例えばパラメータ変化値にパラメータ比を用いる場合、後半曲線部判定部４３５は、予め定められた許容範囲１３２０を記憶装置１２から読み出して、ステップＳ１１２０の判定処理に使用する。ここで、パラメータ比が許容範囲内であると判定された場合、後半曲線部判定部４３５は、直線と判断してステップＳ２５に進む。一方、パラメータ比が許容範囲外と判定された場合、後半曲線部判定部４３５は、直線範囲が終了したと判断してステップＳ１１３０に進む。
【００７９】
処理プロセスがステップＳ２５に進んだ場合、吸光度データ取得部４１０は、入力された吸光度データが最後の吸光度データか否かを判定する。最後の吸光度データでないと判定された場合、吸光度データ取得部４１０は、処理プロセスをステップＳ１０に戻す。一方、ステップＳ２５において、必要なデータ数が蓄積されたと判定された場合、吸光度データ取得部４１０は、処理プロセスをステップＳ７０に進める。
【００８０】
なお、ステップＳ１１３０に進んでいた場合、後半曲線部判定部４３５は、判定された直線終了時刻を記憶装置１２に記憶し、その後、ステップＳ７０に処理プロセスを進める。
【００８１】
ステップＳ７０〜ステップＳ９０は、図１に示すステップＳ７０〜ステップＳ９０の処理と同じである。
【００８２】
以上説明した第２の実施例では、ラグタイム部の判定（すなわち、直線開始時刻の判定）に、第１の実施例と同じ手法を適用する場合について説明した。しかし、本実施例のステップＳ１１２０で説明した後半曲線判定方法を、ラグタイム部の判定に用いることも可能である。その場合、パラメータ変化値が、予め定めた許容範囲外から許容範囲内に転じた時刻を直線開始時刻と判定すれば良い。
【００８３】
以上述べた第２の実施例では、図１１に示す処理を制御部１３で行う例を説明したが、装置の他の部分で処理を行うことも可能である。例えばコンピュータ(PC)１０内においてソフトウェアとして図１１の処理を実行することも可能である。また、記憶装置１２としてコンピュータ(PC)１０内部の記憶装置を使用することも可能である。
【００８４】
図１４に、本実施例で説明した処理手法を用いて直線範囲を決定した例を示す。横軸１１０は反応開始からの経過時間を示し、縦軸１２０は吸光度を示す。時刻１４００は本実施例により決定された直線開始時刻を示し、時刻１４１０は本実施例により決定された直線終了時刻を示し、範囲１４２０は本実施例により決定された直線範囲を示す。図１４に示すように、曲線部を除いた直線部だけを直線範囲と決定できることが確認できた。
【００８５】
以上説明したように、第２の実施例を用いることにより、検査値の換算に使用する直線範囲を検体ごとに決定することができる。そのため、従来装置に比して、精度の高い検査値を得ることが可能となる。同時に、再検査を減らすことが可能となる。また、最終吸光度データを得る以前に後半の曲線部を判定できるので処理時間の短縮が可能となる。
【００８６】
また、第１の実施例と同様に、検査値換算に用いた直線範囲、吸光度データ数、参考CV値や参考吸光度データ平均値等の精度指標も出力することができるため、検体ごとの検査値の精度保証が可能となる。
【００８７】
［実施例３］
次に、第３の実施例に係る自動分析装置を、図面を参照しながら説明する。本実施例の場合も、自動分析装置は生化学自動分析装置であるものとする。従って、その装置構成は、第１の実施例と同じである。すなわち、図３に示す装置構成を有している。また、制御部１３以外の動作は、第１の実施例と同じである。従って、制御部１３以外の詳細な説明は省略する。
【００８８】
以下、制御部１３の処理動作を、本実施例に特有の処理動作を中心に説明する。本実施例の場合、制御部１３は、図１５に示す処理手順に基づいて直線範囲を決定する。なお、図１５のうち図１に示した処理と同じ処理を行う処理ステップには同じ符号を付して示す。
【００８９】
第３の実施例の場合、制御部１３は、反応過程が終了するまでに（全ての吸光度データが測定される前に）、検査値を予測的に算出することを特徴とする。
【００９０】
図１５に示す処理プロセスのうちステップＳ５からステップＳ２０までの処理は、図１におけるステップＳ５からステップＳ２０までの処理と同一である。
【００９１】
次のステップＳ１５００では、吸光度データ取得部４１０が、予測時点の吸光度データが記憶されたか否かを判定する。予測時点は、検査点の予測値（以下、「予測検出点」という。）を算出するために反応過程の中間に定める点である。ここでの予測時点は、検査項目に関係なく１つの時点に定めても良いし、検査項目毎に最適な予測時点を定めても良い。予測時点は、記憶装置１２に予め記憶されている。
【００９２】
予測時点の吸光度データが記憶されていない場合、吸光度データ取得部４１０は、処理プロセスをステップＳ１０に戻し、予測時点に達するまで、吸光度データの入力と記憶を繰り返す。一方、予測時点の吸光度データが記憶された場合、吸光度データ取得部４１０は、処理プロセスをステップＳ１５１０に進める。
【００９３】
ステップＳ１５１０では、近似式パラメータ算出部４２０が、図１のステップＳ３０と同じ方法で近似式パラメータ値を算出する。すなわち、本実施例では、反応過程の開始から予測時点までに取得された吸光度データだけに基づいて第２の近似式を算出する。
【００９４】
次のステップＳ１５２０では、近似式パラメータ算出部４２０が、算出した近似式パラメータ値を記憶装置１２に記憶する。
【００９５】
次のステップＳ４０は、図１のステップＳ４０と同一である。すなわち、算出された第２の近似式に基づいて、ラグタイム部判定部４２５が、反応過程データのラグタイム部を算出する。
【００９６】
次のステップＳ１５４０では、ラグタイム部判定部４２５が、算出されたラグタイム部の長さが予め定めた範囲内か否かを判定する。この処理ステップは、反応過程データが、検出値の予測に適するか否かを判定するために実行される。ラグタイム部の範囲は項目毎に予め実験的に定めておき、記憶装置１２に記憶させておく。判定されたラグタイム部が予め定めた範囲内に存在する場合、ラグタイム部判定部４２５は、処理プロセスをステップＳ１５６０に進める。判定されたラグタイム部が予め定めた範囲外の場合、ラグタイム部判定部４２５は、処理プロセスをステップＳ１５５０に進める。
【００９７】
ステップＳ１５５０に処理プロセスが進んだ場合、ラグタイム部判定部４２５は、予測処理による検査値の算出を停止し、反応過程の全ての吸光度データが取得されるまで、吸光度データの蓄積を継続する。すなわち、予測時点以降もバックグラウンドで進められていた吸光度データの取得を、全ての反応過程が終了するまで継続する。
【００９８】
ステップＳ１５５０に続くステップＳ１５５５では、前述した第１の実施例や第２の実施例と同じ処理手法にて、検体の検査値が実測データに基づいて算出され、出力される。
【００９９】
一方、ステップＳ１５４０からステップＳ１５６０に進んだ場合、ラグタイム部判定部４２５は、算出されたラグタイム部の終了時刻、すなわち直線開始時刻を記憶装置１２に記憶する。
【０１００】
次のステップＳ７０とステップＳ７５は、図１のステップＳ７０及びステップＳ７５と同じである。ただし、ステップＳ７０では、直線終了時刻として記憶装置１２に記憶されている予測時点を用いる。従って、ステップＳ１５６０において直線開始時刻が確定すると、検体の直線範囲も自動的に確定する。
【０１０１】
ステップＳ７５では、直線範囲内吸光度データ数判定部４４０が、吸光度データ数が予め定めた数値の範囲内か否かを判定する。ここで、吸光度データ数が予め定めた数値の範囲外と判定された場合、直線範囲内吸光度データ数判定部４４０は、処理プロセスをステップＳ１５５０に移す。すなわち、予測処理を停止する。一方、吸光度データ数が範囲内の場合、直線範囲内吸光度データ数判定部４４０は、処理をステップＳ１５７０に移す。
【０１０２】
ステップＳ１５７０では、検査値算出部４４５が、直線範囲内の吸光度データに近似する直線式を算出し、当該直線式から検査値換算に用いる傾きを算出する。例えば既存の最小二乗計算方法等を利用し、直線範囲内の吸光度データに近似するように、（数６）で与えられる直線のパラメータを算出する。すなわち、直線の傾きを算出する。
【０１０３】
ステップＳ１５８０では、検査値算出部４４５が、近似直線の傾きから検量線に基づいて検査値の予測値を算出する。検量線データは記憶装置１２から読み出す。その後、検査値算出部４４５は、直線式の算出に用いた直線範囲と直線式パラメータ、算出された予測値を記憶装置１２に記憶する。
【０１０４】
次のステップＳ１５９０では、検査値算出部４４５は、算出された予測値を出力する。検査値の出力方法には、図３のコンピュータ１０に付属する表示装置に表示する方法等が考えられる。表示例を図１６に示す。図１６に示す測定結果の表示画面は、列９００、９０５、９１０、９２０、９３０、９４０、９５０、９６０、１６０５、１６２０、１６３０、１６４０、１６５０、１６６０、１６７０で構成される。列９００には検体IDが記述され、列９０５には予測や検査値の出力に用いた近似式（第１の近似式）が記述される。列１６０５には予測値が記述され、列１６２０には予測に用いた直線の開始時刻（直線開始時刻）が記述される。列１６３０には予測に用いた直線の終了時刻（直線終了時刻）が記述され、列１６４０には予測に用いた直線範囲内の吸光度データ数が記述される。列１６５０，列１６６０，列１６７０には予測に用いた直線範囲の吸光度データで予測値を算出した場合の精度指標が記述される。精度指標には、例えばCV値を用いる。CV値はキャリブレーション用の検体等、検査値が明確になっている検体を用いて実験を行い、予め定めておけば良い。実験では、様々な直線開始時刻と直線終了時刻の組み合わせ(吸光度データ数)で予測した時のCV値(分散や標準偏差等のばらつきを表す値)を算出しておけば良い。精度指標の他の例としては、例えば吸光度データ数の平均値とばらつき(分散や標準偏差等)を用いる。吸光度データ数の平均値とばらつきは、直近に測定した一般(患者)検体やキャリブレーション用検体などの予測結果から、項目ごとに算出すれば良い。なお、算出に用いる検体，検体数，期間はユーザーが自由に選択できるようにすれば良い。更に、精度指標の他の例としては、例えば正答率(一致率)を用いる。正答率は、直近に測定した一般検体やキャリブレーション用検体などから得られた予測値と検査値の結果から、項目ごと、吸光度データ数ごとに算出すれば良い。なお、算出に用いる検体，検体数，期間はユーザーが自由に選択できるようにすれば良い。
【０１０５】
なお、図１５のステップＳ１５５０に処理が進んだ場合、予測値ではなく検査値の出力となる。その場合は、結果を図１６の列９１０〜９６０に記述するようにすれば良い。列９１０〜９６０は、第１の実施例と同じであるので詳細な説明は省略する。
【０１０６】
以上説明した第３の実施例においては、予測値を出力した場合には通常の検査値を出力しない場合について説明した。しかしながら、予測と通常の測定を同時に行い、通常の測定から得られた検査値と予測値の両方を出力する方式を採用しても良い。ここで、通常の測定には、第１の実施例や第２の実施例で説明した処理を用いれば良い。この場合の出力例は図１６の検体ID３に示すように、予測値に関する列１６０５〜列１６７０と検査値に関する列９１０〜列９６０の両方を記述するようにすれば良い。
【０１０７】
以上説明した第３の実施例では、図１５に示す処理の全てを制御部１３で実行する例を説明した。しかし、制御部１３以外の処理装置を使用して、同様の処理を実行することも可能である。例えばコンピュータ(PC)１０内においてソフトウェアとして図１５の処理を実行することも可能である。また、記憶装置１２としてコンピュータ(PC)１０内部の記憶装置を使用することも可能である。
【０１０８】
以上説明したように、第３の実施例では、予測に用いる吸光度データは直線部のみとなるので近似精度が向上し、精度の良い予測が可能となる。また、直線範囲を検体ごとに決定するので、検体ごとにより高い精度で予測値を出力することができる。また、予測処理では反応過程の終了を待つ必要がない。従って、予測値が出力されるまでの時間を短縮することができる。さらに、予測に用いた直線範囲、吸光度データ数、参考CV値や参考正答率等の精度指標を出力するので、検体ごとの予測値の精度保証が可能となる。
【産業上の利用可能性】
【０１０９】
第１及び第２の実施例で説明した通り、本発明を適用した自動分析装置では、レート法により測定される検査項目に関し、精度の良い検査値を提供することが可能になり、検査の信頼性向上に貢献することができる。また、第３の実施例で説明した通り、本発明を適用した自動分析装置では、レート法により測定される検査項目に関し、精度の良い予測が可能になり、検査時間短縮に貢献することができる。
【符号の説明】
【０１１０】
１：サンプルディスク、２：試薬ディスク、３：反応ディスク、４：反応槽、５：サンプリング機構、６：ピペッティング機構、７：攪拌機構、８：測光機構、９：洗浄機構、１０：コンピュータ（PC）、１２：記憶装置、１３：制御部、１４：圧電素子ドライバ、１５：攪拌機構コントローラ、１６：試料容器、１７：円形ディスク、１８：試薬ボトル、１９：円形ディスク、２０：保冷庫、２１：反応容器、２２：反応容器ホルダ、２３：駆動機構、２４：プローブ、２５：支承軸、２６：アーム、２７：プローブ、２８：支承軸、２９：アーム、３１：固定部、３３：ノズル、３４：上下駆動機構、１１０：横軸（時間の経過）、１２０：縦軸（吸光度）、１３０：符号（各時点の吸光度）、１３５：範囲（ラグタイム部）、１４０：符号（各時点の吸光度）、１４５：範囲（ラグタイム部）、１５０：符号（各時点の吸光度）、１５５：範囲（ラグタイム部）、１６０：時間範囲（検査値換算に用いる範囲）、１８０：符号（各時点の吸光度）、２１０：曲線（近似式より算出された吸光度を表す曲線）、２２０：直線（反応過程データを近似した曲線が漸近する直線）、２３０：点（反応過程データを近似した曲線２１０が直線２２０に十分漸近した時間）、４１０：吸光度データ取得部、４２０：近似式パラメータ算出部、４２５：ラグタイム部判定部、４３５：後半曲線部判定部、４４０：直線範囲内吸光度データ数判定部、４４５：検査値算出部、４５０：データバス、５００：テーブル（検査項目／試薬／最適な近似式）、５１０：列（検査項目）、５２０：列（試薬の種類）、５３０：列（近似式の種類）、７００：縦軸（誤差）、７１０：曲線（各時点の誤差を表す曲線）、７２０：時刻（誤差の極大値を示す時刻）、９００：列（検体ID）、９０５：列（近似式の種類）、９１０：列（検査値）、９２０：列（直線開始時刻）、９３０：列（直線終了時刻）、９４０：列（直線範囲内の吸光度データ数）、９５０：列（精度指標１）、９６０：列（精度指標２）、１０００：時刻（直線開始時刻）、１０２０：時刻（直線終了時刻）、１０３０：範囲（直線範囲）、１２００：縦軸（近似式パラメータの値）、１２１０：記号（各時点における近似式パラメータの算出値）、１３００：縦軸（近似式パラメータの比）、１３１０：記号（各時点における近似式パラメータの比の算出値）、１３２０：許容範囲、１４００：時刻（直線開始時刻）、１４１０：時刻（直線終了時刻）、１４２０：範囲（直線範囲）、１６１０：列（予測値）、１６２０：列（予測に用いた直線開始時刻）、１６３０：列（予測に用いた直線終了時刻）、１６４０：列（予測に用いた直線範囲内の吸光度データ数）、１６５０：列（精度指標１）、１６６０：列（精度指標２）、１６７０：列（精度指標３）。
【技術分野】
【０００１】
本発明は、血液、尿その他の生体サンプルを定性又は定量分析する自動分析装置及び方法に関し、例えば臨床検査時の反応をモニタリングする機能を有する自動分析装置及び方法に関する。
【背景技術】
【０００２】
臨床検査用の自動分析装置は、試料と試薬を一定量分注して反応させた後、一定時間にわたり反応液の吸光度を測定し、測定結果に基づき測定対象物質の検査値（濃度や活性値）等を求める。
【０００３】
臨床検査の測定法の一つであるレート法は、主に試料に含まれる酵素成分の活性値を測定するために使用される。測定方法は、試薬として一定量の基質を添加し、酵素が基質を消費して変化する要素を測定する。通常、試薬には十分量の基質が含まれている。このため、試料と試薬の反応が正常に行なわれていれば、一般に反応は時間変化に対して吸光度が一定量ずつ直線的に変化する。レート法ではこの直線の傾き（反応速度）から測定対象物質の活性値を求めている。
【０００４】
しかしながら、測定項目、検体の濃度、撹拌の状況、反応温度などの理由から、反応速度が一定（直線）になるまでにラグタイム部と称される曲線部をもつことがある。また、高活性検体等の場合には反応が早く進むことがある。この場合、測定中に試薬の基質が十分量でなくなってしまい、後半の反応速度が一定でなくなることがある。すなわち、後半の反応速度が非直線状（曲線）に変化することがある。
【０００５】
このように反応過程データには、反応の開始部と後部に曲線部分が出現することがある。このため、従来装置の場合には、検査値の換算に用いる時間範囲と吸光度限界値(Absリミット)を項目ごとに定める手法を採用する。また、従来装置の中には、リニアリティチェックという直線性のチェック機能を備えるものもある。リニアリティチェックでは、一定の測光範囲について、その前半と後半における吸光度変化量の差が算出され、その差が指定したリニアリティチェック値を超える場合に、直線でないと判断される。直線でないと判断された場合、従来装置は「異常」を通知する。
【０００６】
検査値の換算に用いる時間範囲（直線範囲）の決定方法には、特許文献１及び２に開示された方法がある。特許文献１には、隣接する２つの測光ポイント間毎の吸光度差を算出し、それらのうちで一番大きな吸光度差(dmax)を示す測光ポイント間を基準点に定め、その前後の測定範囲の中で吸光度差が基準点の吸光度差(dmax)の80％以上となる区間を直線範囲に決定する方法が開示されている。特許文献２には、換算に用いる範囲幅を一定に設定し、先頭点の測光ポイントから１ポイントずつずらしながら、各範囲内に対応する吸光度の相関係数を最小二乗計算法により算出し、一番相関の良かった範囲を直線範囲に決定する方法が示されている。
【０００７】
一方、近年、病院における臨床検査の分野では、患者検体の測定結果をできる限り早急に報告することが求められている。自動分析装置の測定時の反応時間は１０分ほどであるが、単に測定時間を短くするだけでは、正確な結果が得られなくなってしまう。例えば反応が不完全な時点で測定を終了し、その時点の吸光度から検査値を測定することになった場合には、正確な結果が得られなくなってしまう。そこで、反応終了時の吸光度を予測して検査値に換算する方法が考えられる。反応終了時の吸光度を予測する方法としては、特許文献３が公知である。特許文献３には、（数１）に示す近似式を用いて反応過程データを近似し、レート法に関しては、予め定めた２点間又は近似式を時間で微分した後に予め定めた１点から直線の傾きを求める方法が示されている。なお、（数１）において、ｔは時間、ｙは吸光度を表し、Ａ、Ｂ、Ｋはパラメータである。
（数１） y＝A ＋ (B−A)／eKt
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】特開昭６４−６８６４２号公報
【特許文献２】特開平１０−１８５６８７号公報
【特許文献３】特開平６−１９４３１３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００９】
例えばレート法の主試薬添加後の反応過程データの例を図２の（ａ）及び(ｂ）に示す。横軸１１０は時間の経過を示し、縦軸１２０は吸光度を示す。前述したように、従来装置は、検査項目毎に時間範囲（直線範囲）を定めている。すなわち、直線開始時刻と直線終了時刻を事前に定めている。図中の範囲１６０は、ある検査項目について予め定められた時間範囲を示す。
【００１０】
図２の（ａ）の符号１３０、１４０、１５０は、同じ項目の異なる検体の吸光度データを示し、範囲１３５、１４５、１５５は各吸光度データのラグタイム部を示す。なお、符号１７０は予め定められた吸光度の限界値を示す。
【００１１】
符号１４０で示す吸光度データのように、本来の直線部（本来の検査対象）だけが予め定められた時間範囲１６０に含まれる場合がある一方で、符号１５０で示す吸光度データのように、ラグタイム部１５５の一部が予め定められた時間範囲１６０内に含まれる場合もある。
【００１２】
この符号１５０で示す吸光度データの場合、従来装置のリニアリティチェック機能では、反応過程データに異常があると判断され、その検体は再検査の対象になる。しかし、符号１５０で示す吸光度データの場合でも、その直線部を正確に判断できさえすれば、再検査は必要ない。また、符号１５０で示す吸光度データのように、その直線部とラグタイム部１５５の直線性にそれほどの差がない場合には、リニアリティチェックにおいて異常（ラグタイム部１５５が直線範囲の測定期間に含まれること）を検出することができず、曲線部の情報を含んだまま直線の傾きが算出されてしまう場合もある。この場合、直線の傾きが精度良く算出されないので検査値の誤差が大きくなる可能性がある。
【００１３】
また、符号１３０で示す吸光度データの場合、予め定められた範囲１６０で測定結果は直線を示している。しかし、符号１３０で示す吸光度データは、範囲１６０の後半部分で吸光度の限界値を超える。このため、符号１３０で示す吸光度データの検体も再検査の対象となる。しかし、吸光度の限界値を超えた後も、当該検体の吸光度が直線的に変化することが保証できれば、吸光度の限界値を超えた以降のデータも検査値の換算に用いることが可能となり、再検査を減らすことができる。
【００１４】
図２の（ｂ）に示す符号１８０も、ある検体のある検査項目に関する吸光度データである。この検体の場合、範囲１６０の後半に曲線部分が現れる。このため、この検体も、従来装置のリニアリティチェック機能により、再検査の対象に分類されてしまう。しかし、この検体の場合でも、符号１８０で示す吸光度データの直線範囲だけを抽出することができれば検査値の換算が可能となり、再検査を減らすことができる。
【００１５】
因みに、特許文献１のリニアリティチェック機能を符号１５０で示す吸光度データに適用する場合、一番大きな吸光度差を示す測光ポイントは、図中に円で囲んで示す時点となる。ところが、この時点は、ラグタイム部１５５と直線部の境である。特許文献１のリニアリティチェック機能では、円で囲んだ測光ポイントを基準として、前後の期間を直線範囲に定める。このことは、本来の直線部でないラグタイム部１５５が、検査値の換算に用いられることを意味する。すなわち、検査値の誤差が大きくなる可能性がある。
【００１６】
また、特許文献２のリニアリティチェック機能の場合にも、換算に使用する範囲幅が一定であるので、範囲幅が長く設定されると直線部でない範囲が検査値の換算に用いられる可能性がある。このため、検査値の誤差がやはり大きくなる可能性がある。その一方で、範囲幅が短く設定された場合、検査値の計算に用いる吸光度データが少なくなり、やはり誤差が大きくなる可能性がある。
【００１７】
同様に、反応終了時の吸光度を予測する特許文献３の場合にも、予め定めた範囲（開始時点と終了時点で与えられる）について、反応過程データの直線の傾きを算出することになるため、予め定めた範囲が実測定データの直線範囲の外に位置する場合には、正確な予測値が得られない可能性がある。
【課題を解決するための手段】
【００１８】
そこで、本発明は、試料と試薬を混合し、混合液の経時的な変化を測定する自動分析装置において、(a) 試料と試薬の反応過程から複数の測定点データを取得する測定点データ取得部と、(b) 測定点データを処理するデータ処理部と、(c) データ処理部で使用する第１の近似式を格納した記憶部と、(d) データ処理部の処理結果を出力する出力部とを有するものを提案する。ここで、データ処理部では、記憶部に格納された第１の近似式を複数の測定点データに近似する処理と、近似処理の結果得られる第２の近似式に基づいて前記反応過程の直線範囲を決定する処理とを実行する。
【００１９】
なお、データ処理部は、所定時刻までに取得された複数の測定点データを用いて第２の近似式を算出し、当該第２の近似式に基づいて反応過程の直線範囲を推定又は予測する機能を有することが望ましい。
【００２０】
また、本発明は、試料と試薬を混合した混合液の経時的な変化を、自動分析装置を用いて測定する自動分析方法において、(a) 自動分析装置が、試料と試薬の反応過程から複数の測定点データを取得する処理と、(b) 自動分析装置が、記憶部から読み出した第１の近似式を複数の測定点データに近似する処理と、(c) 自動分析装置が、近似処理の結果得られる第２の近似式に基づいて反応過程の直線範囲を決定する処理とを実行するものを提案する。やはり、所定時刻までに取得された複数の測定点データを用いて第２の近似式を算出し、当該第２の近似式に基づいて反応過程の直線範囲を推定又は予測する処理を有することが望ましい。
【００２１】
なお、本発明は、直線範囲の精度指標を算出し、処理結果として出力部に出力する機能を有することが望ましい。
【発明の効果】
【００２２】
本発明により、検査値の換算に用いる直線範囲を検体毎に定めることが可能となる。結果的に検体毎に正確な検査値が得られる。また、従来装置に比して再検査の対象となる検体数を減らすことができる。さらに、所定時刻までに取得された複数の測定点データを用いて第２の近似式を算出し、当該第２の近似式に基づいて反応過程の直線範囲を推定する場合には、直線範囲の決定までに要する時間を短縮することができる。また、直線範囲の精度指標を出力する機能を有する場合には、検査値や予測値の精度保証が可能となる。
【図面の簡単な説明】
【００２３】
【図１】第１の実施例に係る処理フローチャート。
【図２】レート法の反応過程データ例を説明する図。
【図３】自動分析装置の概略構成例を示す図。
【図４】制御部の内部構成例を示す図。
【図５】検査項目と、試薬の種類と、最適な近似式の関係を記述したテーブル例を示す図。
【図６】ラグタイム部を数値化する方法を説明する図。
【図７】反応後期の曲線部分を数値化する方法を説明する図。
【図８】反応後期の曲線部分を判定する処理フローチャート。
【図９】検体と検査値との関係を記述したテーブル例を示す図。
【図１０】第１の実施例を用いて決定される直線範囲例を説明する図。
【図１１】第２の実施例に係る処理フローチャート。
【図１２】近似式パラメータの分布例を示す図。
【図１３】近似式パラメータ比の分布例を示す図。
【図１４】第２の実施例を用いて決定される直線範囲例を説明する図。
【図１５】第３の実施例に係る処理フローチャート。
【図１６】予測値と検査値と精度指標を記述したテーブル例を示す図。
【発明を実施するための形態】
【００２４】
以下、図面に基づいて、本発明の実施例を説明する。なお、後述する装置構成や処理動作の内容は発明を説明するための一例である。本発明には、後述する装置構成や処理動作に既知の技術を組み合わせた発明や後述する装置構成や処理動作の一部を既知の技術と置換した発明も含まれる。
【実施例】
【００２５】
［実施例１］
以下、第１の実施例に係る自動分析装置の装置構成及び処理動作を、図面を参照して詳細に説明する。図３に、本発明に係る分析機能を実装する生化学自動分析装置の概略構成を示す。
【００２６】
生化学自動分析装置は、サンプルディスク１、試薬ディスク２、反応ディスク３、反応槽４、サンプリング機構５、ピペッティング機構６、攪拌機構７、測光機構８、洗浄機構９、コンピュータ（PC）１０、記憶装置１２、制御部１３、圧電素子ドライバ１４、攪拌機構コントローラ１５、試料容器１６、円形ディスク１７，１９、試薬ボトル１８、保冷庫２０、反応容器２１、反応容器ホルダ２２、駆動機構２３、プローブ２４，２７、支承軸２５，２８、アーム２６，２９、固定部３１、ノズル３３、上下駆動機構３４を有している。記憶装置１２には、分析パラメータ、各試薬ボトルの分析可能回数、最大分析可能回数、キャリブレーション結果、分析結果等が記憶される。
【００２７】
生化学自動分析装置における試料の分析は、サンプリング、試薬分注、撹拌、測光、反応容器の洗浄、濃度換算等のデータ処理の順番に実行される。
【００２８】
サンプルディスク１は、制御部１３によりコンピュータ１０を介して制御される。サンプルディスク１上には、複数の試料容器１６が円周上に並んで設置される。試料容器１６は、分析順に従ってサンプリングプローブ２４の下まで移動される。試料容器１６中の検体は、検体サンプリング機構５に連結された試料用ポンプにより反応容器２１の中に所定量ずつ分注される。
【００２９】
試料が分注された反応容器２１は、反応槽４の中を第１試薬の添加位置まで移動される。移動された反応容器１６には、試薬分注プローブ２７に連結された試薬用ポンプ（図示せず）により試薬容器１８から吸引された試薬が所定量だけ加えられる。第１試薬が添加された後、反応容器２１は撹拌機構７の位置まで移動され、最初の撹拌が行われる。このような試薬の添加と撹拌が、例えば第１〜第４試薬について行われる。
【００３０】
内容物が撹拌された反応容器２１は、光源から発せられる光束中に配置される。光束の一部は反応容器２１を通過し、一部は内容物で吸収される。吸収の度合いは、例えば多波長光度計で構成される測光機構８により検知される。測光機構８は、時間の経過と共に検出された検体に関する吸収の度合いを測定点データ（吸光度信号）として制御部１３に出力する。制御部１３は、後述するデータ処理を通じ、測定点データの直線範囲を決定し、その後、検体の検査値（直線範囲の傾き）を算出する。算出された検査値（データ）は記憶装置１２に記憶され、コンピュータ１０に付属する表示装置に表示される。測光の終了した反応容器２１は、洗浄機構９の位置まで搬送され、洗浄の後、次の分析に使用される。
【００３１】
次に、制御部１３において実行される直線範囲の決定処理手順を説明する。図１に、当該処理手順を示す。なお、図１は、制御部１３内で実行される直線範囲の決定プログラムに関わる部分のみを示している。また、図４に、直線範囲の決定プログラムを通じて実現される制御部１３の機能ブロックを示す。なお、図４は、当該決定プログラムをハードウェアの観点から表している。図４に示す制御部１３は、吸光度データ取得部４１０、近似式パラメータ算出部４２０、ラグタイム部判定部４２５、後半曲線部判定部４３５、直線範囲内吸光度データ数判定部４４０、検査値算出部４４５、これらを相互に接続するデータバス４５０で構成される。データバス４５０を通じ、各部は、相互にデータを受け渡すことができる。なお、制御部１３を構成する機能ブロックは、制御部１３以外のハードウェアやＣＰＵで構成しても良い。勿論、制御部１３を構成する機能ブロックは、同一のＣＰＵ内にソフトウェアモジュールとして実装しても良い。
【００３２】
まず、ある検体に対し、ある検査項目の測定が開始される。この測定開始と同時に、ステップＳ５において、近似式パラメータ算出部４２０が、直線範囲の決定に適した近似式（第１の近似式）を選択する。選択対象としての吸光度の時間変化を表す複数の近似式は、記憶装置１２に予め記憶されているものとする。近似式パラメータ算出部４２０は、これら複数の近似式の中から検査項目と試薬の組み合わせに対応した最適な近似式を選択的に読み出す。
【００３３】
この実施例の場合、（数２）〜（数５）に示す関数を選択可能な近似式とする。なお、各関数において、ｔは時刻、ｘは吸光度を表すものとする。また、ａ、ｂ、ｃ、ｄ、ｅ、ｋ、ｐ、ｑ、ｒ、ｕ、ｖ、ｗはパラメータである。
（数２） x＝a * t ＋ b ＋ c * exp(-k * t)
（数３） x＝a * t ＋ b ＋ e ／ (t ＋ d)
（数４） x＝a * t ＋ b ＋ w ／ {exp(u * t) ＋ v}
（数５） x＝a * t ＋ b ＋ p * log{1 ＋ q * exp(r * t)}
【００３４】
例示した４つの関数は、いずれも、時間に比例した成分と、定数成分と、時間変化率が異なる非直線成分で構成される。勿論、これ以外の関数を近似式として用意することも可能である。
【００３５】
なお、近似式パラメータ算出部４２０による近似式の選択は、検査項目と試薬の組み合わせに基づいて自動的に実行される手法を採用しても良いし、ユーザーが自由に選択する手法を採用しても良い。前者の機能の実現には、例えば検査項目と試薬の組み合わせ毎に最も適した近似式をテーブルとして記憶装置１２に記憶しておけば良い。近似式パラメータ算出部４２０は、このテーブルを検査項目と試薬の組み合わせによって検索し、組み合わせに対応した最適な近似式を選択する。図５に、この種のテーブル例を示す。テーブル５００は、列５１０、５２０及び５３０で構成される。列５１０には検査項目が記述され、列５２０には試薬の種類が記述され、列５３０には検査項目と試薬の種類に対応付けられた最適な近似式の種類が記述されている。
【００３６】
この実施例の場合、近似式パラメータ算出部４２０は、検査項目と試薬の組み合わせに基づいてテーブル５００を検索し、直線範囲の決定に最適な近似式を選択する。なお、テーブル５００に記憶された対応関係の内容は、ユーザーが変更可能な構成としても良い。
【００３７】
吸光度は、時間の経過と共に複数回測定される。次のステップＳ１０では、吸光度データ取得部４１０が、１回の測定又は複数回の測定平均の吸光度データを、測光機構８から入力する。すなわち、吸光度データは制御部１３に入力される。試薬と検体との反応に伴う色調変化に吸光度が大きく変化する波長（主波長）の光と、吸光度が殆ど変化しない波長（副波長）の光の２波長光を用いる測定方式においては、主波長光の吸光度と、副波長光の吸光度との差を、吸光度データとして入力する。
【００３８】
ステップＳ１５では、吸光度データ取得部４１０が、主反応試薬が添加されたかどうかを判定する。主反応試薬が添加されていないと判定された場合、処理プロセスは、ステップＳ１０に戻り、次の吸光度データの入力が行われる。主反応試薬の添加が判定されるまで、以上の判定動作が繰り返し実行される。ここで、主反応試薬とは、複数の試薬を用いる反応において、主たる吸光度変化をもたらす試薬（通常は最終の試薬）を指す。主反応試薬が添加された場合には、処理プロセスはステップＳ２０に移行する。
【００３９】
ステップＳ２０では、入力された吸光度データを、吸光度データ取得部４１０が記憶装置１２に記憶する。
【００４０】
ステップＳ２５では、吸光度データ取得部４１０が、最後の吸光度データが記憶されたか否かを判定する。すなわち、この実施例では、直線範囲の決定前に、予め定めた反応期間に対応する全ての吸光度データが記憶装置１２に記憶されることが必要とされる。ここで、最後の吸光度データが記憶装置１２に記憶されていないと判定された場合、処理プロセスはステップＳ１０に戻る。このループ動作（吸光度データの入力と記憶）は、必要なデータ数が記憶装置１２に記憶されるまで繰り返し実行される。必要なデータ数の蓄積が吸光度データ取得部４１０により判定された場合、処理プロセスはステップＳ３０に進む。
【００４１】
ステップＳ３０では、ステップＳ５において選択された近似式（第１の近似式）によって表される吸光度の時間変化と、実際の吸光度の時間変化との差がなるべく小さくなるように、近似式パラメータ算出部４２０が、数式中のパラメータの値を算出する。具体的には、測定結果として記憶された吸光度データと、近似式により算出される対応時点の吸光度データとの二乗誤差がなるべく小さくなるように、数式中のパラメータ値が算出される。その際、各近似式の最終項（（数２）〜（数５）の場合であれば、（a＊t ＋ b）以降の非直線成分）においては、時刻ｔが増大するとパラメータ値が“０（ゼロ）”に漸近するように制限する。さらに、各近似式の最終項においては、反応初期の吸光度データに近似式をフィッティングさせるようにパラメータの初期値を設定する。
【００４２】
パラメータ値の算出には、既存の最小二乗計算方法が利用可能である。この他、様々な形式の数式に対応可能な方法として、例えば最急降下法により二乗誤差が最小となるパラメータ値を算出する方法もある。この明細書では、各項のパラメータ値を測定結果（吸光度データ）に応じて最適化した後の近似式を第２の近似式と呼び、検査項目と試薬の組み合わせに応じて選択される第１の近似式と区別する。
【００４３】
更に、ステップＳ３０において、近似式パラメータ算出部４２０は、各時刻について、第２の近似式により算出される吸光度（近似値）と実際に測定された吸光度（実測値）との差（誤差）を算出する。
【００４４】
ステップＳ３５では、近似式パラメータ算出部４２０が、第２の近似式を規定するパラメータ値と、各時刻について算出された近似値と実測値との差（誤差）を記憶装置１２に記憶する。
【００４５】
次のステップＳ４０では、ラグタイム部判定部４２５が、反応初期の吸光度が曲線的に変化するラグタイム部分を第２の近似式に基づいて計算する。図６に、ラグタイム部の判定方法を説明する。
【００４６】
図６において横軸１１０は反応開始からの経過時間を示し、縦軸１２０は吸光度を示す。曲線２１０は、第２の近似式により求められる吸光度変化の近似曲線を示す。直線２２０は曲線２１０が漸近する直線である。また、横軸１１０上の点２３０は、曲線２１０が直線２２０に十分漸近する時刻(Ｔｌとする)を示し、横軸１１０の０〜Ｔlの範囲がラグタイム部に相当する。ラグタイム部の終了時刻が直線開始時刻に相当するので、点２３０に対応する時刻が直線開始時刻となる。
【００４７】
ここで、十分漸近した時刻とは、例えば微小な値εを予め定めておき、曲線２１０と直線２２０との差がε以下になった時刻として定義する。εは一定値としても良いし、初期吸光度や吸光度の変化幅に応じて設定しても良い。例えば初期吸光度に定数を乗じた値、又は、初期吸光度と最終吸光度の差に定数を乗じた値をεとしても良い。また、十分漸近した時刻は微小な値δを定めておき、曲線２１０と直線２２０の傾きの差がδ以下になった時刻として定義しても良い。この場合、δは一定値としても良いし、直線２２０の傾きに応じて設定しても良い。例えば直線２２０の傾きに定数を乗じた値をδとしても良い。
【００４８】
ステップＳ４５では、ラグタイム部判定部４２５が、算出した時刻を直線開始時刻として記憶装置１２に記憶する。
【００４９】
次のステップＳ５５では、後半曲線部判定部４３５が、反応過程データの後半に現れる曲線部（以下、「後半曲線部」という。）を判定する。図７と図８を用い、後半曲線部の判定方法を説明する。
【００５０】
図７の上図の横軸１１０は反応開始からの経過時間を表し、縦軸１２０は吸光度を示す。曲線２１０は第２の近似式により求められる吸光度変化の近似曲線である。記号１８０は検体の吸光度データを示す。図７の下図の横軸１１０は反応開始からの経過時間を示し、縦軸７００は誤差を示す。ここでの誤差とは、実測値と近似値の差分値として与えられる。例えば曲線７１０は実測値から近似値を引いた値の時間変化を示す。反応過程データが減少曲線の場合は、例えば求める誤差は近似値から実測値を引いた値にする。この場合も図７の下図に示すような図となる。（数２）〜（数５）の近似式には後半曲線部を表す式が含まれていないので、記号１８０のように後半が曲線になった場合、近似曲線２１０は図７の上図に示すように記号１８０の反応過程曲線と交差する。よって、図７の下図に示すように、近似値と実測値の差（誤差）の時間変化を図示すると、誤差は曲線の開始時刻で極大を示すことになる。この時刻を、この実施例では、直線終了時刻とする。
【００５１】
図８を用い、後半曲線部判定部４３５で実行される処理動作を説明する。まず、ステップＳ８００において、後半曲線部判定部４３５は、記憶装置１２に記憶された誤差を読み出し、誤差が予め定めた許容範囲内か否かを判定する。判定に使用する許容範囲は、記憶装置１２に予め記憶されているものとする。
【００５２】
誤差が許容範囲内の場合、後半曲線部判定部４３５は、ステップＳ８１０に進む。ステップＳ８１０において、後半曲線部判定部４３５は、反応過程データの後半部分に曲線部がないと判定する。この場合、直線終了時刻は反応終了時刻となる。ステップＳ８００で誤差が許容範囲外と判定された場合、後半曲線部判定部４３５は、ステップＳ８２０に進む。
【００５３】
ステップＳ８２０において、後半曲線部判定部４３５は、図７の下図に示すような誤差の時間変化を監視し、誤差が極大値を示す時刻７２０（Ｔｅとする）を検出する。
【００５４】
次のステップＳ８３０において、後半曲線部判定部４３５は、時刻７２０以降の誤差の分布が単調減少か否かを判定する。単調減少するか否かは、例えば隣り合う時間の誤差の差分の符号が一定か否かで判定することができる。単調減少を示していた場合、図７の上図に示すように後半部分に曲線部があることを示す。従って、後半曲線部判定部４３５はステップＳ８４０に進み、反応過程データの後半に曲線部があると判定する。この場合、直線終了時刻は曲線開始時刻となる。すなわち、図７の時刻７２０となる。ステップＳ８３０で単調減少とならないと判定された場合、後半曲線部判定部４３５は、ステップＳ８５０に進んで測定の「異常」を示す信号を出力する。「異常」を示す信号の出力方法には、アラームを鳴らす方法、図３のコンピュータ１０に付属する表示装置に表示する方法等が考えられる。
【００５５】
図１の説明に戻る。ステップＳ５５の処理が終了すると、後半曲線部判定部４３５は、ステップＳ６０に進む。ステップＳ６０において、後半曲線部判定部４３５は、判定結果として得られた直線終了時刻を記憶装置１２に記憶する。
【００５６】
ステップＳ７０では、直線範囲内吸光度データ数判定部４４０が、記憶装置１２に記憶された直線開始時刻と直線終了時刻とから直線範囲を決定する。
【００５７】
続くステップＳ７５では、直線範囲内吸光度データ数判定部４４０が、直線範囲内に出現する吸光度データの数が予め定めた閾値以上か否かを判定する。なお、直線式は、例えば（数６）に示すような１次式で表される。（数６）を用いて直線の傾きを求める場合に必要な吸光度データ数の範囲は２個以上である。ただし、（数６）において、ｔは時刻、ｘは吸光度、ａ、ｂはパラメータを示す。
（数６）x＝a * t ＋ b
【００５８】
また、信頼できる検査値を算出するために必要な最低限の吸光度データ数を予め実験から定めておいても良い。
【００５９】
ところで、ステップＳ７５の判定処理において、直線範囲内の吸光度データ数が閾値より小さいと判定された場合、直線範囲内吸光度データ数判定部４４０は、ステップＳ９０に進み、測定が「異常」であると判定する。「異常」を示す信号の出力方法には、アラームを鳴らす方法、図３のコンピュータ１０に付属する表示装置に表示する方法等が考えられる。
【００６０】
一方、ステップＳ７５の判定処理において、直線範囲内の吸光度データ数が閾値以上と判定された場合、直線範囲内吸光度データ数判定部４４０は、処理プロセスをステップＳ８０に進める。ステップＳ８０においては、検査値算出部４４５が、直線範囲内の吸光度データと近似する直線式を算出し、その直線式から検査値換算に用いる傾きを算出する。
【００６１】
ここでの直線式の算出（傾きを算出する処理に対応する）には、例えば既存の最小二乗計算方法等を利用する。検査値算出部４４５は、直線範囲内の吸光度データに近似する（数６）で与えられる直線を算出する。すなわち、直線の傾きを算出する。また、検査値算出部４４５は、傾きから検量線に基づいて検査値を算出することもできる。検量線データは予め記憶装置１２に記憶しておく。直線式を算出した直線範囲と直線式パラメータ、算出された検査値は記憶装置１２に記憶される。
【００６２】
次のステップＳ８５において、検査値算出部４４５は、算出された検査値を出力する。検査値の出力方法には、図３のコンピュータ１０に付属する表示装置に表示する等が考えられる。表示例を図９に示す。図９に示す測定結果の表示画面は、列９００から９６０で構成される。列９００には検体IDが記述され、列９１０には検査値が記述される。列９０５には直線範囲の決定に使用した近似式（第１の近似式）が記述される。列９２０には検査値の換算に用いた直線範囲の開始時刻（直線開示時刻）が記述される。列９３０には検査値の換算に用いた直線範囲の終了時刻（直線終了時刻）が記述される。列９４０には検査値の換算に用いた直線範囲内に出現する吸光度データ数が記述される。列９５０，列９６０には検査値の換算に用いた直線範囲とその吸光度データ数とで検査値を算出した場合の精度指標が記述されている。精度指標には、例えばCV値を用いる。CV値はキャリブレーション用の検体等、検査値が明確になっている検体を用いて実験を行い、予め定めておけば良い。実験では、様々な直線開始時刻と直線終了時刻の組み合わせ(吸光度データ数)についてCV値(分散や標準偏差等ばらつきを表す値)を算出しておけば良い。精度指標の他の例としては，例えば吸光度データ数の平均値とばらつき(分散や標準偏差等)を用いる。吸光度データ数の平均値とばらつきは，直近に測定した一般(患者)検体やキャリブレーション用検体などの結果から，項目ごとに算出すれば良い。なお，算出に用いる検体，検体数，期間はユーザーが自由に選択できるようにすれば良い。通常、直線範囲内の吸光度データ数が多いほど、検出値の信頼性は一般的に高くなる。
【００６３】
以上説明した第１の実施例においては、図１に示す処理の全てを制御部１３で実行する例を説明した。しかし、制御部１３以外の処理装置を使用して、同様の処理を実行することも可能である。例えばコンピュータ(PC)１０内で実行されるソフトウェア処理として図１の処理を実行することも可能である。また、記憶装置１２としてコンピュータ(PC)１０内部の記憶装置を使用することも可能である。
【００６４】
以上説明した第１の実施例により直線範囲を決定した例を、図１０の（ａ）及び(ｂ）に示す。（ａ）及び（ｂ）の横軸１１０は反応開始からの経過時間を示し、縦軸１２０は吸光度を示す。時刻１０００は、本実施例を通じて決定された直線開始時刻を示し、時刻１０２０は本実施例を通じて決定された直線終了時刻を示す。範囲１０３０は本実施例により決定された直線範囲を表す。なお、図１０の（ａ）はラグタイムのある反応過程データの例であり、図１０の(ｂ）は反応過程の後半に曲線のある反応過程データの例を示している。どちらの場合も、曲線部を除いた直線部だけを直線範囲と決定できることが確認できた。
【００６５】
以上の通り、本実施例で説明した直線範囲の決定機能の採用により、検査値の換算に用いる直線範囲を検体ごとに決定することができる。そのため、従来装置に比して、精度の高い検査値を得ることが可能となる。同時に、再検査を減らすことが可能となる。更に、検査値換算に用いた直線範囲、吸光度データ数、参考CV値や参考吸光度データ平均値等の精度指標も出力することができるため、検体ごとの検査値の精度保証が可能となる。
【００６６】
［実施例２］
次に、第２の実施例に係る自動分析装置を、図面を参照しながら説明する。本実施例の場合も、自動分析装置は生化学自動分析装置であるものとする。従って、その装置構成は、第１の実施例と同じである。すなわち、図３に示す装置構成を有している。また、制御部１３以外の動作は、第１の実施例と同じである。従って、制御部１３以外の詳細な説明は省略する。
【００６７】
以下、制御部１３の処理動作を、本実施例に特有の処理動作を中心に説明する。本実施例の場合、制御部１３は、図１１に示す処理手順に基づいて直線範囲を決定する。なお、図１１のうち図１に示した処理と同じ処理を行う処理ステップには同じ符号を付して示す。
【００６８】
第２の実施例の場合、制御部１３は、反応過程の全ての吸光度データが測定されるのを待つことなく、吸光度データの測定動作と並行的に直線範囲の開始点の決定処理を開始する。
【００６９】
図１１に示す処理プロセスのうちステップＳ５からステップＳ２０までの処理は、図１におけるステップＳ５からステップＳ２０までの処理と同一である。
【００７０】
次のステップＳ１１００では、近似式パラメータ算出部４２０が、現在までに取得された吸光度データに基づいて近似式パラメータとパラメータ変化値を算出する。近似式パラメータの算出方法は、図１に示すステップＳ３０と同一である。ただし、第１の実施例では、反応過程の最後の吸光度データが取得されてから、取得した全ての吸光度データを用いて直線範囲を算出したのに対し、この第２の実施例の場合には、吸光度データが取得される度に近似式パラメータとパラメータ変化値を順次算出する。
【００７１】
ここで、パラメータ変化値とは、近似式パラメータの時間変化を数値化した値のことであり、算出には直線部の傾きを表す近似式パラメータを使用する。例えば第２の近似式が（数２）〜（数５）のいずれかの形式で与えられる場合、近似式パラメータａを使用する。パラメータの変化を数値化する方法には、様々な方法が利用可能である。例えば一回前に算出されたパラメータ値との差を算出する方法、一回前に算出されたパラメータとの比を算出する方法等を利用できる。なお、この実施例では直線部の傾きを与えるパラメータにのみ着目しているが、その他のパラメータ、例えば非直線部のパラメータの変化に着目しても良い。また、直線部と非直線部の両方のパラメータの変化に着目しても良い。
【００７２】
ステップＳ１１１０では、近似式パラメータ算出部４２０が、算出した近似式パラメータとパラメータ変化値を記憶装置１２に記憶する。
【００７３】
次のステップＳ４０とステップＳ４５は、図１に示すステップＳ４０とステップＳ４５と同じである。この実施例の場合、パラメータ変化値が許容範囲以下となるタイミングをラグタイム部の終了時点と判定する。
【００７４】
次のステップＳ１１２０では、後半曲線部判定部４３５が、算出されたパラメータ変化値が予め定めた範囲内か否かを判定し、反応過程データの後半曲線部の判定を行う。図１２と図１３を用いて後半曲線部の判定方法を説明する。
【００７５】
図１２は、直線の傾きを表す近似式パラメータとして、例えば（数２）〜（数５）のパラメータａを用いた場合の経時変化を示す。横軸１１０は反応開始からの経過時間を示し、縦軸１２００は近似式パラメータａの値を示す。記号１２１０は、各時刻に算出された近似式パラメータａの値を示している。近似式パラメータａがほぼ一定の間は反応過程データが直線上に位置することが分かり、近似式パラメータａが一定値から外れると反応過程データが曲線に沿って変化していることが分かる。
【００７６】
例えばパラメータ変化値に近似式パラメータａを用いる場合、後半曲線部判定部４３５は、近似式パラメータａが一定と判断された場合、後半曲線部判定部４３５は、直線と判断してステップＳ２５に進む。一方、近似式パラメータａが変化したと判定された場合、後半曲線部判定部４３５は、直線範囲が終了したと判断してステップＳ１１３０に進む。
【００７７】
図１３は、図１２に示すパラメータの分布に基づいて、一回前に算出されたパラメータとの比を、パラメータ変化値として数値化する場合を示す。横軸１１０は反応開始からの経過時間を示し、縦軸１３００はパラメータ比を示す。記号１３１０は、各時刻の一回前に算出されたパラメータとの比の値を示す。図１３に示すように、反応過程データが直線を示す場合は、パラメータ比は１に近い値となり、直線を示さなくなった場合にはパラメータ比は１から離れた値となる。
【００７８】
例えばパラメータ変化値にパラメータ比を用いる場合、後半曲線部判定部４３５は、予め定められた許容範囲１３２０を記憶装置１２から読み出して、ステップＳ１１２０の判定処理に使用する。ここで、パラメータ比が許容範囲内であると判定された場合、後半曲線部判定部４３５は、直線と判断してステップＳ２５に進む。一方、パラメータ比が許容範囲外と判定された場合、後半曲線部判定部４３５は、直線範囲が終了したと判断してステップＳ１１３０に進む。
【００７９】
処理プロセスがステップＳ２５に進んだ場合、吸光度データ取得部４１０は、入力された吸光度データが最後の吸光度データか否かを判定する。最後の吸光度データでないと判定された場合、吸光度データ取得部４１０は、処理プロセスをステップＳ１０に戻す。一方、ステップＳ２５において、必要なデータ数が蓄積されたと判定された場合、吸光度データ取得部４１０は、処理プロセスをステップＳ７０に進める。
【００８０】
なお、ステップＳ１１３０に進んでいた場合、後半曲線部判定部４３５は、判定された直線終了時刻を記憶装置１２に記憶し、その後、ステップＳ７０に処理プロセスを進める。
【００８１】
ステップＳ７０〜ステップＳ９０は、図１に示すステップＳ７０〜ステップＳ９０の処理と同じである。
【００８２】
以上説明した第２の実施例では、ラグタイム部の判定（すなわち、直線開始時刻の判定）に、第１の実施例と同じ手法を適用する場合について説明した。しかし、本実施例のステップＳ１１２０で説明した後半曲線判定方法を、ラグタイム部の判定に用いることも可能である。その場合、パラメータ変化値が、予め定めた許容範囲外から許容範囲内に転じた時刻を直線開始時刻と判定すれば良い。
【００８３】
以上述べた第２の実施例では、図１１に示す処理を制御部１３で行う例を説明したが、装置の他の部分で処理を行うことも可能である。例えばコンピュータ(PC)１０内においてソフトウェアとして図１１の処理を実行することも可能である。また、記憶装置１２としてコンピュータ(PC)１０内部の記憶装置を使用することも可能である。
【００８４】
図１４に、本実施例で説明した処理手法を用いて直線範囲を決定した例を示す。横軸１１０は反応開始からの経過時間を示し、縦軸１２０は吸光度を示す。時刻１４００は本実施例により決定された直線開始時刻を示し、時刻１４１０は本実施例により決定された直線終了時刻を示し、範囲１４２０は本実施例により決定された直線範囲を示す。図１４に示すように、曲線部を除いた直線部だけを直線範囲と決定できることが確認できた。
【００８５】
以上説明したように、第２の実施例を用いることにより、検査値の換算に使用する直線範囲を検体ごとに決定することができる。そのため、従来装置に比して、精度の高い検査値を得ることが可能となる。同時に、再検査を減らすことが可能となる。また、最終吸光度データを得る以前に後半の曲線部を判定できるので処理時間の短縮が可能となる。
【００８６】
また、第１の実施例と同様に、検査値換算に用いた直線範囲、吸光度データ数、参考CV値や参考吸光度データ平均値等の精度指標も出力することができるため、検体ごとの検査値の精度保証が可能となる。
【００８７】
［実施例３］
次に、第３の実施例に係る自動分析装置を、図面を参照しながら説明する。本実施例の場合も、自動分析装置は生化学自動分析装置であるものとする。従って、その装置構成は、第１の実施例と同じである。すなわち、図３に示す装置構成を有している。また、制御部１３以外の動作は、第１の実施例と同じである。従って、制御部１３以外の詳細な説明は省略する。
【００８８】
以下、制御部１３の処理動作を、本実施例に特有の処理動作を中心に説明する。本実施例の場合、制御部１３は、図１５に示す処理手順に基づいて直線範囲を決定する。なお、図１５のうち図１に示した処理と同じ処理を行う処理ステップには同じ符号を付して示す。
【００８９】
第３の実施例の場合、制御部１３は、反応過程が終了するまでに（全ての吸光度データが測定される前に）、検査値を予測的に算出することを特徴とする。
【００９０】
図１５に示す処理プロセスのうちステップＳ５からステップＳ２０までの処理は、図１におけるステップＳ５からステップＳ２０までの処理と同一である。
【００９１】
次のステップＳ１５００では、吸光度データ取得部４１０が、予測時点の吸光度データが記憶されたか否かを判定する。予測時点は、検査点の予測値（以下、「予測検出点」という。）を算出するために反応過程の中間に定める点である。ここでの予測時点は、検査項目に関係なく１つの時点に定めても良いし、検査項目毎に最適な予測時点を定めても良い。予測時点は、記憶装置１２に予め記憶されている。
【００９２】
予測時点の吸光度データが記憶されていない場合、吸光度データ取得部４１０は、処理プロセスをステップＳ１０に戻し、予測時点に達するまで、吸光度データの入力と記憶を繰り返す。一方、予測時点の吸光度データが記憶された場合、吸光度データ取得部４１０は、処理プロセスをステップＳ１５１０に進める。
【００９３】
ステップＳ１５１０では、近似式パラメータ算出部４２０が、図１のステップＳ３０と同じ方法で近似式パラメータ値を算出する。すなわち、本実施例では、反応過程の開始から予測時点までに取得された吸光度データだけに基づいて第２の近似式を算出する。
【００９４】
次のステップＳ１５２０では、近似式パラメータ算出部４２０が、算出した近似式パラメータ値を記憶装置１２に記憶する。
【００９５】
次のステップＳ４０は、図１のステップＳ４０と同一である。すなわち、算出された第２の近似式に基づいて、ラグタイム部判定部４２５が、反応過程データのラグタイム部を算出する。
【００９６】
次のステップＳ１５４０では、ラグタイム部判定部４２５が、算出されたラグタイム部の長さが予め定めた範囲内か否かを判定する。この処理ステップは、反応過程データが、検出値の予測に適するか否かを判定するために実行される。ラグタイム部の範囲は項目毎に予め実験的に定めておき、記憶装置１２に記憶させておく。判定されたラグタイム部が予め定めた範囲内に存在する場合、ラグタイム部判定部４２５は、処理プロセスをステップＳ１５６０に進める。判定されたラグタイム部が予め定めた範囲外の場合、ラグタイム部判定部４２５は、処理プロセスをステップＳ１５５０に進める。
【００９７】
ステップＳ１５５０に処理プロセスが進んだ場合、ラグタイム部判定部４２５は、予測処理による検査値の算出を停止し、反応過程の全ての吸光度データが取得されるまで、吸光度データの蓄積を継続する。すなわち、予測時点以降もバックグラウンドで進められていた吸光度データの取得を、全ての反応過程が終了するまで継続する。
【００９８】
ステップＳ１５５０に続くステップＳ１５５５では、前述した第１の実施例や第２の実施例と同じ処理手法にて、検体の検査値が実測データに基づいて算出され、出力される。
【００９９】
一方、ステップＳ１５４０からステップＳ１５６０に進んだ場合、ラグタイム部判定部４２５は、算出されたラグタイム部の終了時刻、すなわち直線開始時刻を記憶装置１２に記憶する。
【０１００】
次のステップＳ７０とステップＳ７５は、図１のステップＳ７０及びステップＳ７５と同じである。ただし、ステップＳ７０では、直線終了時刻として記憶装置１２に記憶されている予測時点を用いる。従って、ステップＳ１５６０において直線開始時刻が確定すると、検体の直線範囲も自動的に確定する。
【０１０１】
ステップＳ７５では、直線範囲内吸光度データ数判定部４４０が、吸光度データ数が予め定めた数値の範囲内か否かを判定する。ここで、吸光度データ数が予め定めた数値の範囲外と判定された場合、直線範囲内吸光度データ数判定部４４０は、処理プロセスをステップＳ１５５０に移す。すなわち、予測処理を停止する。一方、吸光度データ数が範囲内の場合、直線範囲内吸光度データ数判定部４４０は、処理をステップＳ１５７０に移す。
【０１０２】
ステップＳ１５７０では、検査値算出部４４５が、直線範囲内の吸光度データに近似する直線式を算出し、当該直線式から検査値換算に用いる傾きを算出する。例えば既存の最小二乗計算方法等を利用し、直線範囲内の吸光度データに近似するように、（数６）で与えられる直線のパラメータを算出する。すなわち、直線の傾きを算出する。
【０１０３】
ステップＳ１５８０では、検査値算出部４４５が、近似直線の傾きから検量線に基づいて検査値の予測値を算出する。検量線データは記憶装置１２から読み出す。その後、検査値算出部４４５は、直線式の算出に用いた直線範囲と直線式パラメータ、算出された予測値を記憶装置１２に記憶する。
【０１０４】
次のステップＳ１５９０では、検査値算出部４４５は、算出された予測値を出力する。検査値の出力方法には、図３のコンピュータ１０に付属する表示装置に表示する方法等が考えられる。表示例を図１６に示す。図１６に示す測定結果の表示画面は、列９００、９０５、９１０、９２０、９３０、９４０、９５０、９６０、１６０５、１６２０、１６３０、１６４０、１６５０、１６６０、１６７０で構成される。列９００には検体IDが記述され、列９０５には予測や検査値の出力に用いた近似式（第１の近似式）が記述される。列１６０５には予測値が記述され、列１６２０には予測に用いた直線の開始時刻（直線開始時刻）が記述される。列１６３０には予測に用いた直線の終了時刻（直線終了時刻）が記述され、列１６４０には予測に用いた直線範囲内の吸光度データ数が記述される。列１６５０，列１６６０，列１６７０には予測に用いた直線範囲の吸光度データで予測値を算出した場合の精度指標が記述される。精度指標には、例えばCV値を用いる。CV値はキャリブレーション用の検体等、検査値が明確になっている検体を用いて実験を行い、予め定めておけば良い。実験では、様々な直線開始時刻と直線終了時刻の組み合わせ(吸光度データ数)で予測した時のCV値(分散や標準偏差等のばらつきを表す値)を算出しておけば良い。精度指標の他の例としては、例えば吸光度データ数の平均値とばらつき(分散や標準偏差等)を用いる。吸光度データ数の平均値とばらつきは、直近に測定した一般(患者)検体やキャリブレーション用検体などの予測結果から、項目ごとに算出すれば良い。なお、算出に用いる検体，検体数，期間はユーザーが自由に選択できるようにすれば良い。更に、精度指標の他の例としては、例えば正答率(一致率)を用いる。正答率は、直近に測定した一般検体やキャリブレーション用検体などから得られた予測値と検査値の結果から、項目ごと、吸光度データ数ごとに算出すれば良い。なお、算出に用いる検体，検体数，期間はユーザーが自由に選択できるようにすれば良い。
【０１０５】
なお、図１５のステップＳ１５５０に処理が進んだ場合、予測値ではなく検査値の出力となる。その場合は、結果を図１６の列９１０〜９６０に記述するようにすれば良い。列９１０〜９６０は、第１の実施例と同じであるので詳細な説明は省略する。
【０１０６】
以上説明した第３の実施例においては、予測値を出力した場合には通常の検査値を出力しない場合について説明した。しかしながら、予測と通常の測定を同時に行い、通常の測定から得られた検査値と予測値の両方を出力する方式を採用しても良い。ここで、通常の測定には、第１の実施例や第２の実施例で説明した処理を用いれば良い。この場合の出力例は図１６の検体ID３に示すように、予測値に関する列１６０５〜列１６７０と検査値に関する列９１０〜列９６０の両方を記述するようにすれば良い。
【０１０７】
以上説明した第３の実施例では、図１５に示す処理の全てを制御部１３で実行する例を説明した。しかし、制御部１３以外の処理装置を使用して、同様の処理を実行することも可能である。例えばコンピュータ(PC)１０内においてソフトウェアとして図１５の処理を実行することも可能である。また、記憶装置１２としてコンピュータ(PC)１０内部の記憶装置を使用することも可能である。
【０１０８】
以上説明したように、第３の実施例では、予測に用いる吸光度データは直線部のみとなるので近似精度が向上し、精度の良い予測が可能となる。また、直線範囲を検体ごとに決定するので、検体ごとにより高い精度で予測値を出力することができる。また、予測処理では反応過程の終了を待つ必要がない。従って、予測値が出力されるまでの時間を短縮することができる。さらに、予測に用いた直線範囲、吸光度データ数、参考CV値や参考正答率等の精度指標を出力するので、検体ごとの予測値の精度保証が可能となる。
【産業上の利用可能性】
【０１０９】
第１及び第２の実施例で説明した通り、本発明を適用した自動分析装置では、レート法により測定される検査項目に関し、精度の良い検査値を提供することが可能になり、検査の信頼性向上に貢献することができる。また、第３の実施例で説明した通り、本発明を適用した自動分析装置では、レート法により測定される検査項目に関し、精度の良い予測が可能になり、検査時間短縮に貢献することができる。
【符号の説明】
【０１１０】
１：サンプルディスク、２：試薬ディスク、３：反応ディスク、４：反応槽、５：サンプリング機構、６：ピペッティング機構、７：攪拌機構、８：測光機構、９：洗浄機構、１０：コンピュータ（PC）、１２：記憶装置、１３：制御部、１４：圧電素子ドライバ、１５：攪拌機構コントローラ、１６：試料容器、１７：円形ディスク、１８：試薬ボトル、１９：円形ディスク、２０：保冷庫、２１：反応容器、２２：反応容器ホルダ、２３：駆動機構、２４：プローブ、２５：支承軸、２６：アーム、２７：プローブ、２８：支承軸、２９：アーム、３１：固定部、３３：ノズル、３４：上下駆動機構、１１０：横軸（時間の経過）、１２０：縦軸（吸光度）、１３０：符号（各時点の吸光度）、１３５：範囲（ラグタイム部）、１４０：符号（各時点の吸光度）、１４５：範囲（ラグタイム部）、１５０：符号（各時点の吸光度）、１５５：範囲（ラグタイム部）、１６０：時間範囲（検査値換算に用いる範囲）、１８０：符号（各時点の吸光度）、２１０：曲線（近似式より算出された吸光度を表す曲線）、２２０：直線（反応過程データを近似した曲線が漸近する直線）、２３０：点（反応過程データを近似した曲線２１０が直線２２０に十分漸近した時間）、４１０：吸光度データ取得部、４２０：近似式パラメータ算出部、４２５：ラグタイム部判定部、４３５：後半曲線部判定部、４４０：直線範囲内吸光度データ数判定部、４４５：検査値算出部、４５０：データバス、５００：テーブル（検査項目／試薬／最適な近似式）、５１０：列（検査項目）、５２０：列（試薬の種類）、５３０：列（近似式の種類）、７００：縦軸（誤差）、７１０：曲線（各時点の誤差を表す曲線）、７２０：時刻（誤差の極大値を示す時刻）、９００：列（検体ID）、９０５：列（近似式の種類）、９１０：列（検査値）、９２０：列（直線開始時刻）、９３０：列（直線終了時刻）、９４０：列（直線範囲内の吸光度データ数）、９５０：列（精度指標１）、９６０：列（精度指標２）、１０００：時刻（直線開始時刻）、１０２０：時刻（直線終了時刻）、１０３０：範囲（直線範囲）、１２００：縦軸（近似式パラメータの値）、１２１０：記号（各時点における近似式パラメータの算出値）、１３００：縦軸（近似式パラメータの比）、１３１０：記号（各時点における近似式パラメータの比の算出値）、１３２０：許容範囲、１４００：時刻（直線開始時刻）、１４１０：時刻（直線終了時刻）、１４２０：範囲（直線範囲）、１６１０：列（予測値）、１６２０：列（予測に用いた直線開始時刻）、１６３０：列（予測に用いた直線終了時刻）、１６４０：列（予測に用いた直線範囲内の吸光度データ数）、１６５０：列（精度指標１）、１６６０：列（精度指標２）、１６７０：列（精度指標３）。
【特許請求の範囲】
【請求項１】
試料と試薬を混合し、混合液の経時的な変化を測定する自動分析装置において、
前記試料と前記試薬の反応過程から複数の測定点データを取得する測定点データ取得部と、
前記測定点データを処理するデータ処理部と、
前記データ処理部で使用する第１の近似式を格納した記憶部と、
前記データ処理部の処理結果を出力する出力部と、を有し、
前記データ処理部は、前記記憶部に格納された第１の近似式を前記複数の測定点データに近似する処理と、近似処理の結果得られる第２の近似式に基づいて前記反応過程の直線範囲を決定する処理とを実行する
ことを特徴とする自動分析装置。
【請求項２】
請求項１に記載の自動分析装置において、
前記データ処理部は、所定時刻までに取得された前記複数の測定点データを用いて前記第２の近似式を算出し、当該第２の近似式に基づいて前記反応過程の直線範囲を推定する
ことを特徴とする自動分析装置。
【請求項３】
請求項２に記載の自動分析装置において、
前記データ処理部は、前記第２の近似式が漸近する直線を算出し、当該直線に前記第２の近似式が基準値以下に達した時刻を直線開始時刻とし、当該直線開始時刻と事前に定めた直線終了時刻との間の期間を直線範囲として推定する
ことを特徴とする自動分析装置。
【請求項４】
請求項２に記載の自動分析装置において、
前記記憶部は、前記直線開始時刻の基準範囲を記憶し、
前記データ処理部は、前記直線範囲の推定値と前記基準範囲を比較し、比較結果に基づいて直線範囲の推定が可能か否かを判定する
ことを特徴とする自動分析装置。
【請求項５】
請求項２に記載の自動分析装置において、
前記記憶部は、前記直線範囲が含むべき前記測定点データの数の基準値を記憶し、
前記データ処理部は、信号処理を通じて推定した前記直線範囲に含まれる測定点データの数をカウントする処理と、カウントされた前記測定点データの数と前記基準値とを比較し、直線範囲の推定が可能か否かを判定する処理とを有する
ことを特徴とする自動分析装置。
【請求項６】
請求項１に記載の自動分析装置において、
前記第１の近似式は、以下に示す数１〜数４のいずれか１つである
ことを特徴とする自動分析装置。
数１：x ＝ a * t + b + c * exp( -k * t )
数２：x ＝ a * t + b + e / ( t + d )
数３：x ＝ a * t + b + w / { exp( u * t ) + v }
数４：x ＝ a * t + b + p * log{ 1 + q * exp( r * t ) }
ただし、tは時刻、xは吸光度、a、b、c、d、e、k、p、q、r、u、v、wはパラメータである。
【請求項７】
請求項１に記載の自動分析装置において、
前記データ処理部は、前記第２の近似式が漸近する直線を算出し、当該直線に前記第２の近似式が基準値以下に達した時刻を直線開始時刻に決定する
ことを特徴とする自動分析装置。
【請求項８】
請求項１に記載の自動分析装置において、
前記データ処理部は、直線開始時刻以降の測定点データと前記第２の近似式が漸近する直線との差分値を算出し、当該差分値の経時変化に基づいて直線終了時刻を決定する
ことを特徴とする自動分析装置。
【請求項９】
請求項１に記載の自動分析装置において、
前記データ処理部は、各時点までに取得された複数の測定点データを用いて前記第２の近似式を逐次算出する処理と、逐次算出される前記第２の近似式のうち少なくとも１つのパラメータの変化が許容範囲以下に小さくなった時点を直線開始時刻に決定する
ことを特徴とする自動分析装置。
【請求項１０】
請求項１に記載の自動分析装置において、
前記データ処理部は、直線開始時刻の決定後も、各時点までに取得された複数の測定点データを用いて前記第２の近似式を逐次算出する処理と、逐次算出される前記第２の近似式のうち少なくとも１つのパラメータに閾値を越える変化が検出された時点を直線終了時刻に決定する
ことを特徴とする自動分析装置。
【請求項１１】
請求項１０に記載の自動分析装置において、
前記１つのパラメータは、前記第２の近似式を規定する直線成分の傾きを与えるパラメータである
ことを特徴とする自動分析装置。
【請求項１２】
請求項１に記載の自動分析装置において、
前記記憶部は、前記直線範囲が含むべき前記測定点データの数の基準値を記憶し、
前記データ処理部は、信号処理を通じて決定した直線範囲に含まれる測定点データの数をカウントする処理と、カウントされた前記測定点データの数と前記基準値とを比較して決定された直線範囲に関する異常の有無を判定する処理とを有する
ことを特徴とする自動分析装置。
【請求項１３】
請求項１に記載の自動分析装置において、
前記データ処理部は、決定された前記直線範囲の精度指標を算出する処理と、算出された前記精度指標を前記出力部に出力する処理とを有する
ことを特徴とする自動分析装置。
【請求項１４】
試料と試薬を混合した混合液の経時的な変化を、自動分析装置を用いて測定する自動分析方法において、
前記自動分析装置が、前記試料と前記試薬の反応過程から複数の測定点データを取得する処理と、
前記自動分析装置が、記憶部から読み出した第１の近似式を前記複数の測定点データに近似する処理と、
前記自動分析装置が、近似処理の結果得られる第２の近似式に基づいて前記反応過程の直線範囲を決定する処理と
を有することを特徴とする自動分析方法。
【請求項１】
試料と試薬を混合し、混合液の経時的な変化を測定する自動分析装置において、
前記試料と前記試薬の反応過程から複数の測定点データを取得する測定点データ取得部と、
前記測定点データを処理するデータ処理部と、
前記データ処理部で使用する第１の近似式を格納した記憶部と、
前記データ処理部の処理結果を出力する出力部と、を有し、
前記データ処理部は、前記記憶部に格納された第１の近似式を前記複数の測定点データに近似する処理と、近似処理の結果得られる第２の近似式に基づいて前記反応過程の直線範囲を決定する処理とを実行する
ことを特徴とする自動分析装置。
【請求項２】
請求項１に記載の自動分析装置において、
前記データ処理部は、所定時刻までに取得された前記複数の測定点データを用いて前記第２の近似式を算出し、当該第２の近似式に基づいて前記反応過程の直線範囲を推定する
ことを特徴とする自動分析装置。
【請求項３】
請求項２に記載の自動分析装置において、
前記データ処理部は、前記第２の近似式が漸近する直線を算出し、当該直線に前記第２の近似式が基準値以下に達した時刻を直線開始時刻とし、当該直線開始時刻と事前に定めた直線終了時刻との間の期間を直線範囲として推定する
ことを特徴とする自動分析装置。
【請求項４】
請求項２に記載の自動分析装置において、
前記記憶部は、前記直線開始時刻の基準範囲を記憶し、
前記データ処理部は、前記直線範囲の推定値と前記基準範囲を比較し、比較結果に基づいて直線範囲の推定が可能か否かを判定する
ことを特徴とする自動分析装置。
【請求項５】
請求項２に記載の自動分析装置において、
前記記憶部は、前記直線範囲が含むべき前記測定点データの数の基準値を記憶し、
前記データ処理部は、信号処理を通じて推定した前記直線範囲に含まれる測定点データの数をカウントする処理と、カウントされた前記測定点データの数と前記基準値とを比較し、直線範囲の推定が可能か否かを判定する処理とを有する
ことを特徴とする自動分析装置。
【請求項６】
請求項１に記載の自動分析装置において、
前記第１の近似式は、以下に示す数１〜数４のいずれか１つである
ことを特徴とする自動分析装置。
数１：x ＝ a * t + b + c * exp( -k * t )
数２：x ＝ a * t + b + e / ( t + d )
数３：x ＝ a * t + b + w / { exp( u * t ) + v }
数４：x ＝ a * t + b + p * log{ 1 + q * exp( r * t ) }
ただし、tは時刻、xは吸光度、a、b、c、d、e、k、p、q、r、u、v、wはパラメータである。
【請求項７】
請求項１に記載の自動分析装置において、
前記データ処理部は、前記第２の近似式が漸近する直線を算出し、当該直線に前記第２の近似式が基準値以下に達した時刻を直線開始時刻に決定する
ことを特徴とする自動分析装置。
【請求項８】
請求項１に記載の自動分析装置において、
前記データ処理部は、直線開始時刻以降の測定点データと前記第２の近似式が漸近する直線との差分値を算出し、当該差分値の経時変化に基づいて直線終了時刻を決定する
ことを特徴とする自動分析装置。
【請求項９】
請求項１に記載の自動分析装置において、
前記データ処理部は、各時点までに取得された複数の測定点データを用いて前記第２の近似式を逐次算出する処理と、逐次算出される前記第２の近似式のうち少なくとも１つのパラメータの変化が許容範囲以下に小さくなった時点を直線開始時刻に決定する
ことを特徴とする自動分析装置。
【請求項１０】
請求項１に記載の自動分析装置において、
前記データ処理部は、直線開始時刻の決定後も、各時点までに取得された複数の測定点データを用いて前記第２の近似式を逐次算出する処理と、逐次算出される前記第２の近似式のうち少なくとも１つのパラメータに閾値を越える変化が検出された時点を直線終了時刻に決定する
ことを特徴とする自動分析装置。
【請求項１１】
請求項１０に記載の自動分析装置において、
前記１つのパラメータは、前記第２の近似式を規定する直線成分の傾きを与えるパラメータである
ことを特徴とする自動分析装置。
【請求項１２】
請求項１に記載の自動分析装置において、
前記記憶部は、前記直線範囲が含むべき前記測定点データの数の基準値を記憶し、
前記データ処理部は、信号処理を通じて決定した直線範囲に含まれる測定点データの数をカウントする処理と、カウントされた前記測定点データの数と前記基準値とを比較して決定された直線範囲に関する異常の有無を判定する処理とを有する
ことを特徴とする自動分析装置。
【請求項１３】
請求項１に記載の自動分析装置において、
前記データ処理部は、決定された前記直線範囲の精度指標を算出する処理と、算出された前記精度指標を前記出力部に出力する処理とを有する
ことを特徴とする自動分析装置。
【請求項１４】
試料と試薬を混合した混合液の経時的な変化を、自動分析装置を用いて測定する自動分析方法において、
前記自動分析装置が、前記試料と前記試薬の反応過程から複数の測定点データを取得する処理と、
前記自動分析装置が、記憶部から読み出した第１の近似式を前記複数の測定点データに近似する処理と、
前記自動分析装置が、近似処理の結果得られる第２の近似式に基づいて前記反応過程の直線範囲を決定する処理と
を有することを特徴とする自動分析方法。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【図１３】
【図１４】
【図１５】
【図１６】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【図１３】
【図１４】
【図１５】
【図１６】
【公開番号】特開２０１１−２２６９０９（Ｐ２０１１−２２６９０９Ａ）
【公開日】平成２３年１１月１０日（２０１１．１１．１０）
【国際特許分類】
【出願番号】特願２０１０−９６７７０（Ｐ２０１０−９６７７０）
【出願日】平成２２年４月２０日（２０１０．４．２０）
【出願人】（５０１３８７８３９）株式会社日立ハイテクノロジーズ (4,325)
【Ｆターム（参考）】
【公開日】平成２３年１１月１０日（２０１１．１１．１０）
【国際特許分類】
【出願日】平成２２年４月２０日（２０１０．４．２０）
【出願人】（５０１３８７８３９）株式会社日立ハイテクノロジーズ (4,325)
【Ｆターム（参考）】
[ Back to top ]