(57)【要約】【課題】  混合物の成分濃度推定方法及び装置を提供する。【解決手段】  濃度を推定しようとする所定の特定成分を含む複数の独立変数により決定される前記濃度値を従属変数とする補正データセットに対して全域補正モデルを構成する段階と、前記補正データセットを前記従属変数値によって少なくとも２つ以上の小グループに分離し、前記分離された各小グループに属する補正データセットを用いて各小グループ別のローカル補正モデルを構成する段階と、与えられた混合物スペクトルが属する小グループを決定し、決定された小グループのローカル補正モデルを用いて前記特定成分の濃度を推定する段階とよりなる混合物の成分濃度推定方法。【選択図】図３
2005121663.tif 000002
【特許請求の範囲】
【請求項１】
（ａ）濃度を推定しようとする所定の特定成分を含む複数の独立変数により決定される前記濃度値を従属変数とする補正データセットを前記従属変数値によって少なくとも２つ以上の小グループに分離する段階と、
（ｂ）前記分離された各小グループに属する補正データセットを用いて各小グループ別のローカル補正モデルを構成する段階と、を含むことを特徴とするローカル補正モデル生成方法。
【請求項２】
前記（ａ）段階において、前記小グループが２つである場合、前記独立変数の平均に該当する従属変数値を基準として分離することを特徴とする請求項１に記載のローカル補正モデル生成方法。
【請求項３】
前記（ａ）段階において、前記小グループが２つ以上である場合、前記従属変数値の平均を基準として分離することを特徴とする請求項１に記載のローカル補正モデル生成方法。
【請求項４】
前記（ａ）段階において、濃度を推定しようとする特定成分の妨害成分を用いて前記補正データセットを分離することを特徴とする請求項１に記載のローカル補正モデル生成方法。
【請求項５】
（ａ）濃度を推定しようとする所定の特定成分を含む複数の独立変数により決定される前記濃度値を従属変数とする補正データセットに対して全域補正モデルを構成する段階と、
（ｂ）前記補正データセットを前記従属変数値によって少なくとも２つ以上の小グループに分離し、前記分離された各小グループに属する補正データセットを用いて各小グループ別のローカル補正モデルを構成する段階と、
（ｃ）与えられた混合物スペクトルが属する小グループを決定し、決定された小グループのローカル補正モデルを用いて前記特定成分の濃度を推定する段階と、を含むことを特徴とする混合物の成分濃度推定方法。
【請求項６】
前記（ｂ）段階において、前記小グループが２つである場合、前記独立変数の平均に該当する従属変数値を基準として分離することを特徴とする請求項５に記載の混合物の成分濃度推定方法。
【請求項７】
前記（ｂ）段階において、前記小グループが２つ以上である場合、前記従属変数値の平均を基準として分離することを特徴とする請求項５に記載の混合物の成分濃度推定方法。
【請求項８】
前記（ｂ）段階において、濃度を推定しようとする特定成分の妨害成分を用いて前記補正データセットを分離することを特徴とする請求項５に記載の混合物の成分濃度推定方法。
【請求項９】
前記（ｃ）段階において、前記与えられた混合物スペクトルに対して前記全域補正モデルを適用して得られる前記特定成分の濃度値を用いて前記混合物スペクトルが属する小グループを決定することを特徴とする請求項５に記載の混合物の成分濃度推定方法。
【請求項１０】
前記方法は、
（ｄ）前記（ｃ）段階で推定された濃度値が推定時に適用されるローカル補正モデルによる濃度範囲に属しない場合、前記全域補正モデルを適用して前記特定成分の濃度値を推定する段階をさらに含むことを特徴とする請求項５に記載の混合物の成分濃度推定方法。
【請求項１１】
（ａ）濃度を推定しようとする所定の特定成分を含む複数の独立変数により決定される前記濃度値を従属変数とする補正データセットを前記従属変数値によって少なくとも２つ以上の小グループに分離する段階と、
（ｂ）前記分離された各小グループに属する補正データセットを用いて各小グループ別のローカル補正モデルを構成する段階と、を含むローカル補正モデル生成方法を実行できるプログラムを記録したコンピュータ可読記録媒体。
【請求項１２】
（ａ）濃度を推定しようとする所定の特定成分を含む複数の独立変数により決定される前記濃度値を従属変数とする補正データセットに対して全域補正モデルを構成する段階と、
（ｂ）前記補正データセットを前記従属変数値によって少なくとも２つ以上の小グループに分離し、前記分離された各小グループに属する補正データセットを用いて各小グループ別のローカル補正モデルを構成する段階と、
（ｃ）与えられた混合物スペクトルが属する小グループを決定し、決定された小グループのローカル補正モデルを用いて前記特定成分の濃度を推定する段階と、を含む混合物の成分濃度推定方法を実行できるプログラムを記録したコンピュータ可読記録媒体。
【請求項１３】
濃度を推定しようとする所定の特定成分を含む複数の独立変数により決定される前記濃度値を従属変数とする補正データセットに対して全域補正モデルを構成し、前記全域補正モデルを保存する第１補正モデル形成手段と、
前記補正データセットを前記従属変数値によって少なくとも２つ以上の小グループに分離し、前記分離された各小グループに属する補正データセットを用いて各小グループ別のローカル補正モデルを構成し、小グループ別のローカル補正モデルを保存する第２補正モデル形成手段と、
与えられた混合物スペクトルが属する小グループを決定する小グループ決定部と、
前記決定された小グループのローカル補正モデルを用いて前記与えられた混合物スペクトルに含まれた前記特定成分の濃度を推定する濃度推定部と、を含むことを特徴とする混合物の成分濃度推定装置。
【請求項１４】
前記第２補正モデル形成手段は、前記小グループが２つである場合に前記独立変数の平均に該当する従属変数値を基準とし、前記小グループが２つ以上である場合に前記従属変数値の平均を基準として分離することを特徴とする請求項１３に記載の混合物の成分濃度推定装置。
【請求項１５】
前記小グループ決定部は、前記与えられた混合物スペクトルに対して前記全域補正モデルを適用して得られる前記特定成分の濃度値を用いて前記混合物スペクトルが属する小グループを決定することを特徴とする請求項１３に記載の混合物の成分濃度推定装置。
【請求項１６】
前記濃度推定部は、前記決定された小グループのローカル補正モデルを用いて推定された前記特定成分の濃度値が推定時に適用されるローカル補正モデルによる濃度範囲に属しない場合、前記全域補正モデルを適用して前記特定成分の濃度値を推定することを特徴とする請求項１３に記載の混合物の成分濃度推定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、混合物の成分濃度推定に係り、特に複数のローカル補正モデルを用いて混合物スペクトルに含まれた特定成分の濃度を推定する方法及び装置に関する。
【背景技術】
【０００２】
混合物のスペクトルから特定成分の濃度を推定するために一般的に使われる方法、例えばＰＣＲ(Ｐｒｉｎｃｉｐａｌ  Ｃｏｍｐｏｎｅｎｔ  Ｒｅｇｒｅｓｓｉｏｎ)のような多変量分析法では補正データセットから回帰モデルである全域補正モデルを構成し、全域補正モデルを用いて与えられたスペクトルから特定成分の濃度を予測する。ところが、全域補正モデルを構成するための補正データセットを得る過程では必然的にノイズが発生する。一般的に、ノイズの大きさや特徴は補正データセット全体にわたって同一でないために、ノイズの大きさや特徴によって補正データセットを概念的に多数の小グループに分離する。この場合、補正データセットから全域補正モデルを構成するようになれば、何れか一つの小グループのノイズ特性が他の小グループに伝播される可能性がある。例えば、補正データセットの何れか一つの小グループにはノイズが全くないとしても、全域補正モデルを構成する場合には他の小グループのノイズが伝播されて予測誤差を誘発させる恐れがある。すなわち、図１Ａ及び図１Ｂに示されたように一つの全域補正モデルを適用する場合、独立変数ｘ値を基準として任意の小グループには適切に適用できるが、他のグループには正確に適用され難い問題点がある。
【発明の開示】
【発明が解決しようとする課題】
【０００３】
本発明がな解決しようとする技術的課題は、多変量分析法を利用した特定成分の濃度予測時に発生可能なノイズの伝播問題を解決するために、所定の補正データセットを複数の小グループに分離し、各小グループ別にローカル補正モデルを構成するための方法を提供することである。
【０００４】
本発明が解決しようとする他の技術的課題は、多変量分析法を利用した特定成分の濃度予測時に発生できるノイズの伝播問題を解決するために、複数のローカル補正モデルを用いて混合物の成分濃度を予測するための方法及び装置を提供することである。
【課題を解決するための手段】
【０００５】
前記技術的課題を達成するために、本発明に係るローカル補正モデル生成方法は、（ａ）濃度を推定しようとする所定の特定成分を含む複数の独立変数により決定される前記濃度値を従属変数とする補正データセットを前記従属変数値によって少なくとも２つ以上の小グループに分離する段階と、（ｂ）前記分離された各小グループに属する補正データセットを用いて各小グループ別のローカル補正モデルを構成する段階と、を含むことが望ましい。
【０００６】
前記他の技術的課題を達成するために、本発明に係る混合物の成分濃度推定方法は、（ａ）濃度を推定しようとする所定の特定成分を含む複数の独立変数により決定される前記濃度値を従属変数とする補正データセットに対して全域補正モデルを構成する段階と、（ｂ）前記補正データセットを前記従属変数値によって少なくとも２つ以上の小グループに分離し、前記分離された各小グループに属する補正データセットを用いて各小グループ別のローカル補正モデルを構成する段階と、（ｃ）与えられた混合物スペクトルが属する小グループを決定し、決定された小グループのローカル補正モデルを用いて前記特定成分の濃度を推定する段階と、を含むことが望ましい。
【０００７】
前記他の技術的課題を達成するために、本発明に係る混合物の成分濃度推定装置は、濃度を推定しようとする所定の特定成分を含む複数の独立変数により決定される前記濃度値を従属変数とする補正データセットに対して全域補正モデルを構成し、前記全域補正モデルを保存する第１補正モデル形成手段と、前記補正データセットを前記従属変数値によって少なくとも２つ以上の小グループに分離し、前記分離された各小グループに属する補正データセットを用いて各小グループ別のローカル補正モデルを構成し、小グループ別のローカル補正モデルを保存する第２補正モデル形成手段と、与えられた混合物スペクトルが属する小グループを決定する小グループ決定部と、前記決定された小グループのローカル補正モデルを用いて前記与えられた混合物スペクトルに含まれた前記特定成分の濃度を推定する濃度推定部と、を含むことが望ましい。
【０００８】
前記方法は、望ましくはコンピュータで実行させるためのプログラムを記録したコンピュータ可読記録媒体に具現できる。
【発明の効果】
【０００９】
本発明によれば、補正データセットに存在するノイズの大きさや特徴が独立変数値の全区間にわたって同一でない場合、補正データセットを少なくとも２つ以上の小グループに分離して各小グループ別にローカル補正モデルを構成し、与えられた混合物スペクトルが属する小グループのローカル補正モデルを適用して特定成分の濃度を推定することによって、何れか一つの小グループによる予測誤差が他の小グループには影響を及ぼさなく全体的に予測誤差を減少させうる。
【発明を実施するための最良の形態】
【００１０】
以下、添付した図面に基づき、本発明の望ましい実施例について詳細に説明する。
一般的に補正モデル、すなわち、回帰モデルの構成にあたって、変数ｘとｙは次の数式１のような直線回帰模型で表現できると仮定する。
【００１１】
【数１】
2005121663.tif 000003
ここで、εは誤差項であって、次の数式２のような統計的な特性を有すると仮定する。
【００１２】
【数２】
2005121663.tif 000004
すなわち、あらゆるｘに対して誤差項εは等分散性を有すると仮定する。しかし、実際においては誤差項の非同質性によってあらゆるｘに対して等分散性の条件を充足させない場合が多い。特に、ｘ値が大きくなるほど誤差項の分散が大きくなる場合には大きい値を有するｘにおける誤差が小さな値を有するｘに伝播されてくるために、小さな値を有するｘでの予測力に悪影響をおよぼす。レベレッジ(ｌｅｖｅｒａｇｅ)が大きい値を有するｘでの誤差が小さな値を有するｘに伝播されてくる理由は、通常の場合、レベレッジが大きい値を有するｘが回帰モデルの構成時に小さな値を有するｘより大きい影響をおよぼすためである。したがって、この場合には与えられた補正データセットを多数の小グループに分離し、各小グループごとにローカル補正モデルを構成すれば、誤差の伝播による予測力の損傷を減少させることができる。以下、補正データセットを多数の小グループに分離する方法について説明する。
【００１３】
まず、回帰線ｙ＾が次の数式３のような形式を有すると推定する。
【数３】
2005121663.tif 000005
ここで、ｂ０、ｂ１、ｙ＾はそれぞれβ０、β１、ｙの推定値を表す。
【００１４】
ｂ０及びｂ１は、最小ニ乗法を用いて次の数式４及び数式５のように求められる。
【数４】
2005121663.tif 000006
【数５】
2005121663.tif 000007
前記数式４及び数式５において、ｘ−及びｙ−はそれぞれｘの平均とｙの平均を表す。
【００１５】
一方、前記数式３は次の数式６ないし８のように再配列する。
【数６】
2005121663.tif 000008
【数７】
2005121663.tif 000009
【数８】
2005121663.tif 000010
【００１６】
ここで、Ｓ(ｘｙ)及びＳ(ｘｘ)は、次の数式９のように表される。
【数９】
2005121663.tif 000011
【００１７】
ｙの推定値ｙ＾に対する期待値E（ｙ＾）は次の数式１０のように表される。
【数１０】
2005121663.tif 000012
【００１８】
次に、ｙの推定値ｙ＾に対する分散Var(ｙ＾)は次の数式１１のように表される。
【数１１】
2005121663.tif 000013
【００１９】
すなわち、ｙの推定値ｙ＾に対する分散Var（ｙ＾）は次の数式１２のように表される。
【数１２】
2005121663.tif 000014
【００２０】
一方、一つのｙ値の信頼区間Var(ｙｓ＾)は、次の数式１３のように求められる。
【数１３】
2005121663.tif 000015
前記数式１３において、ｙ値の信頼区間Var(ｙｓ＾)は２次グラフに近似化させることができ、図４において回帰線の上部に位置するグラフに該当する。回帰線の下部に位置するグラフは、上部に位置するグラフを回帰線に対して対称させたものである。任意のｘに対して２つの２次グラフでｙ値の差をｙ値の信頼区間Var(ｙｓ＾)と決定する。
前記数式１３からｙ値の信頼区間Var(ｙｓ＾)はｘの関数であって、ｘがｘの平均と同じである時、すなわち、ｘ＝ｘ−  である場合に最小値となり、この最小値から対称にｘ値がｘの平均から遠いほど大きくなる。したがって、補正データセットを多数の小グループに分離する方法は、補正データセットをｘの平均を基準とした濃度によって分離することが望ましい。その理由は、複数のローカル補正モデルを利用する場合、与えられたスペクトルからどの小グループに属しているかを決定しなければならないが、２つの小グループに分離する場合にはｘの平均を基準として補正データセットを分離する場合にエラー率が最小になれる。
【００２１】
図２は、本発明に係る混合物の成分濃度推定装置の一実施例の構成を示したブロック図であって、全域補正モデル生成部２１１、第１貯蔵部２１２、ローカル補正モデル生成部２１３、第２貯蔵部２１４、小グループ決定部２１５、及び濃度推定部２１６からなる。
【００２２】
図２を参照すれば、全域補正モデル生成部２１１では濃度を推定しようとする所定の特定成分を含む複数の独立変数により決定される濃度値を従属変数とする補正データセットに対して全域補正モデルを構成し、全域補正モデルを第１貯蔵部２１２に保存する。
【００２３】
ローカル補正モデル生成部２１３では補正データセットを前記従属変数値によって少なくとも２つ以上の小グループに分離し、前記分離された各小グループに属する補正データセットを用いて各小グループ別のローカル補正モデルを構成し、小グループ別のローカル補正モデルを第２貯蔵部２１４に保存する。分離しようとする小グループが２つである場合、独立変数の平均に該当する従属変数値を基準として分離することが望ましい。分離しようとする小グループが２つ以上である場合、従属変数値の平均を基準として分離することが望ましい。他の成分としては、濃度を推定しようとする特定成分、例えばグルコースの妨害成分、すなわちヘモグロビンを用いて補正データセットを分離できる。
【００２４】
小グループ決定部２１５では全域補正モデル生成部２１１から生成された全域補正モデルおよびローカル補正モデル生成部２１３から分離された小グループ別補正データセットを参照して、与えられた混合物スペクトルが属する小グループを決定する。この時、与えられた混合物スペクトルに対して全域補正モデルを適用して得られる特定成分の濃度値を用いて混合物スペクトルが属する小グループを決定する。
【００２５】
濃度推定部２１６では小グループ決定部２１５から決定された小グループのローカル補正モデルを用いて与えられた混合物スペクトルに含まれた特定成分の濃度を推定する。この時、推定された濃度値が推定時に適用されるローカル補正モデルによる濃度範囲に属しない場合、与えられた混合物スペクトルに全域補正モデルを適用して特定成分の濃度値を推定する。
【００２６】
図３は、本発明に係る混合物の成分濃度推定方法の一実施例の動作を説明したフローチャートである。
図３を参照すれば、３１１段階においては、一般的な多変量分析法を適用して補正データセット全体に対する全域回帰モデルを構成する。
【００２７】
３１２段階においては、補正データセットをn個(ここでｎは２以上の整数）の小グループに分離する。この時、ｎが２である場合、独立変数の平均に該当する従属変数値を基準とし、ｎが２以上である場合、従属変数値の平均を基準として分離する。一方、濃度を推定しようとする特定成分の妨害成分を用いて補正データセットを分離できる。例えば、特定成分がグルコースである場合ヘモグロビンが妨害成分に該当する。
【００２８】
３１３段階においては、全域補正モデルの構成に使われた多変量分析法を用いて各小グループに対するローカル補正モデルを構成する。
【００２９】
３１４段階においては、与えられた混合物スペクトルが属する小グループを決定する。このために、与えられた混合物スペクトルに全域補正モデルを適用して予備濃度値を算出し、与えられた混合物スペクトルの小グループを予備濃度値が属する小グループに決定する。
【００３０】
３１５段階においては、与えられた混合物スペクトルが属する小グループのローカル補正モデルを適用して特定成分の濃度を予測する。
【００３１】
３１６段階においては、３１５段階で予測された濃度が該当小グループに設定されている濃度範囲内に存在しているかを判断する。
【００３２】
３１７段階においては、３１６段階での判断結果、予測された濃度が該当小グループに設定されている濃度範囲内に存在していない場合には、３１１段階で生成された全域補正モデルを与えた混合物スペクトルに適用して特定成分の濃度を予測する。
【００３３】
３１８段階においては、３１５段階または３１７段階で予測された濃度を特定成分に対する推定濃度として決定する。
【００３４】
本発明は、またコンピュータ可読記録媒体にコンピュータ可読コードとして具現することができる。コンピュータ可読記録媒体は、コンピュータシステムによって読取り可能なデータが保存されるあらゆる種類の記録装置を含む。コンピュータ可読記録媒体の例としては、ＲＯＭ、ＲＡＭ、ＣＤ−ＲＯＭ、磁気テープ、フロッピー（登録商標）ディスク、光データ貯蔵装置などがあり、また、キャリアウェーブ(例えばインターネットを通じた伝送）の形態に具現されるものも含む。また、コンピュータ可読記録媒体は、ネットワークにより連結されたコンピュータシステムに分散され、分散方式によってコンピュータが読取り可能なコードが保存されて実行できる。そして、本発明を具現するための機能的なプログラム、コード及びコードセグメントは本発明が属する技術分野のプログラマーにより容易に推論できる。
【００３５】
一方、本発明による混合物の成分濃度推定方法の性能を評価するために１００個のインビトロ(ｉｎ−ｖｉｔｒｏ)補正データセットに対して複数のローカル補正モデルを適用した。この時、補正データセットのサンプル数は１００個であり、濃度の最小値は３１ｍｇ／ｄｌ、最大値は４８５ｍｇ／ｄｌ、平均は２２１ｍｇ／ｄｌであった。このような補正データセットに対して従来のように全域補正モデルを使用する場合、図５Ａに示されたように、予測誤差は１７.０１％であり、濃度を３区間に分離して予測誤差を分析する場合、図５Ｂに示されたように濃度が小さな区間であるほど相対誤差が大きくなることがわかる。特に、濃度が３１〜１４５ｍｇ／ｄｌ区間では相対誤差が３１％に至った。
【００３６】
一方、本発明に係る複数のローカル補正モデルを使用する場合、図６Ａに示すように補正データセットを２つの小グループ６１０、６２０に分離し、各小グループに対してローカル補正モデルを構成し、予測誤差を分析した結果は、図６Ｂに示すように全域補正モデルを使用する場合に比べて特に低濃度での相対誤差が２５.７１％に減少した。また、複数のローカル補正モデルを使用する場合、予測誤差は１５.２９％に全域補正モデルに比べて１.７１％減少した。
【００３７】
以上、図面と明細書において最適の実施例が開示された。ここで特定の用語が使われたが、これは単に本発明を説明するために使われたものであって、意味限定や特許請求の範囲上に記載された本発明の範囲を制限するために使われたものではない。したがって、当業者であれば、これより多様な変形及び均等な他の実施例が可能である。したがって、本発明の真の技術的な保護範囲は特許請求の範囲上の技術的な思想によって決まるべきである。
【産業上の利用可能性】
【００３８】
本発明は、混合物で特定成分の濃度予測のために使われる多変量分析法に適用されて、補正データセットを複数の小グループに分離して各小グループ別にローカル補正モデルを構成し、与えられた混合物スペクトルが属する小グループのローカル補正モデルを適用して特定成分の濃度を推定することによって、特定成分の濃度予測時に発生可能なノイズの伝播問題を解決する。
【図面の簡単な説明】
【００３９】
【図１Ａ】一般的な多変量分析法による濃度推定時に予測誤差を説明する図面である。
【図１Ｂ】一般的な多変量分析法による濃度推定時に予測誤差を説明する図面である。
【図２】本発明に係る混合物の成分濃度推定装置の一実施例の構成を示したブロック図である。
【図３】本発明に係る混合物の成分濃度推定方法の一実施例の動作を説明したフローチャートである。
【図４】本発明に係るローカル補正モデル構成時に使われる信頼区間を説明する図面である。
【図５Ａ】一般的な全域補正モデル適用時の実測濃度と予測濃度とを示すグラフである。
【図５Ｂ】一般的な全域補正モデル適用時の実測濃度と予測濃度とを示すグラフである。
【図６Ａ】本発明に係るローカル補正モデル適用時の実測濃度と予測濃度とを示すグラフである。
【図６Ｂ】本発明に係るローカル補正モデル適用時の実測濃度と予測濃度とを示すグラフである。
【符号の説明】
【００４０】
２１１・・・全域補正モデル生成部
２１２・・・第１貯蔵部
２１３・・・ローカル補正モデル生成部
２１４・・・第２貯蔵部
２１５・・・小グループ決定部
２１６・・・濃度推定部
【図１Ａ】
2005121663.tif 000016
【図１Ｂ】
2005121663.tif 000017
【図２】
2005121663.tif 000018
【図３】
2005121663.tif 000019
【図４】
2005121663.tif 000020
【図５Ａ】
2005121663.tif 000021
【図５Ｂ】
2005121663.tif 000022
【図６Ａ】
2005121663.tif 000023
【図６Ｂ】
2005121663.tif 000024