(57)【要約】試料と参照データのそれぞれのガスクロマトグラフデータとの一致の程度を決定するためにガスクロマトグラフィを利用する方法であって、該方法は、クロマトグラフデータを獲得する手順と、モデル関数がクロマトグラフデータに実質的に適合するようにモデル関数における複数のパラメータを決定する手順であって、モデル関数は、確率分布関数の線形結合の和として画定される、決定する手順と、モデル関数とクロマトグラフデータとの適合度に従って、一致の程度を評価する手順と、を含む。【選択図】図１
2013502575.tif 000002
【特許請求の範囲】
【請求項１】
試料と参照データのそれぞれのガスクロマトグラフデータとの一致の程度を決定するためにガスクロマトグラフィを利用する方法であって、
前記ガスクロマトグラフデータを獲得する手順と、
モデル関数を前記ガスクロマトグラフデータに実質的に適合させるように、前記モデル関数における複数のパラメータを決定する手順であって、前記モデル関数は確率分布関数の線形結合の和として画定される、決定する手順と、
前記モデル関数と前記ガスクロマトグラフデータとの適合度に従って、一致の前記程度を評価する手順と、を含む方法。
【請求項２】
獲得する前記手順は、一次元ガスクロマトグラフ分離技法によって実行される、請求項１に記載の方法。
【請求項３】
前記試料を試料収集デバイスを介して収集する最初の手順をさらに含む、請求項１に記載の方法。
【請求項４】
獲得する前記手順は、検出器によって獲得される信号からのものである、請求項１に記載の方法。
【請求項５】
前記試料を実質的に構成する成分の不均衡な濃度の存在を占めるように、前記信号を正規化する手順をさらに含む、請求項４に記載の方法。
【請求項６】
メモリデバイス内の前記ガスクロマトグラフデータを記憶する手順をさらに含み、前記ガスクロマトグラフデータは、少なくとも１つのクロマトグラフピークを含む、請求項４に記載の方法。
【請求項７】
データベースを形成するように、前記参照データを記憶する手順をさらに含み、前記データベース内の前記参照データは、少なくとも１つの公知の化合物に実質的に対応する少なくとも１つのクロマトグラフ参照ピークを含む、請求項６に記載の方法。
【請求項８】
前記少なくとも１つの公知の化合物に対応する、質量分析データを有する前記参照データを較正する手順をさらに含む、請求項７に記載の方法。
【請求項９】
前記少なくとも１つの公知の化合物の特定のものは、そこから前記試料が獲得される生物の少なくとも１つの生物学的状態で実質的に表わされる少なくとも１つのバイオマーカを画定する、請求項７に記載の方法。
【請求項１０】
前記少なくとも１つのクロマトグラフピークの保持時間値を、前記少なくとも１つのクロマトグラフ参照ピークのそれぞれの参照保持時間値と比較することによって、前記少なくとも１つのクロマトグラフピークを、前記少なくともクロマトグラフ参照に登録する手順をさらに含む、請求項９に記載の方法。
【請求項１１】
登録する前記手順は、前記少なくとも１つのクロマトグラフピークの時間領域の状態位置を、前記参照データ内のそれぞれの前記少なくとも１つの参照クロマトグラフピークのそれぞれの位置と比較するものである、請求項１０に記載の方法。
【請求項１２】
登録する前記手順は、ｓ（ｆ（ｔ））が対応する参照データｒ（ｔ）と最大限に一致するように選択された、変形関数ｆ（ｔ）を利用することによって実行される、請求項１０に記載の方法。
【請求項１３】
前記変形関数は、現在現れている前記少なくとも１つのクロマトグラフピークを考慮に入れるように選択される、請求項１２に記載の方法。
【請求項１４】
登録の前記手順はさらに、それぞれの少なくとも１つの公知のクロマトグラフピークを発生するように、少なくとも１つの公知の保持時間を有する少なくとも１つの参照化学物質を前記試料に追加するものであり、それによって前記変形関数が前記少なくとも１つの公知のクロマトグラフピークに従って選択される、請求項１２に記載の方法。
【請求項１５】
登録する前記手順により、その一時的属性に従って前記少なくとも１つのクロマトグラフピークを分類する手順をさらに含む、請求項１２に記載の方法。
【請求項１６】
前記少なくとも１つのクロマトグラフピークが、それぞれの前記確率分布関数によってモデル化されるように、分類する前記手順に従って、前記モデル関数を構築する手順をさらに含む、請求項１５に記載の方法。
【請求項１７】
前記モデル関数は以下の一般の形であり、
【数１】
2013502575.tif 000013
ただし、Ｖｉはｉ個の前記確率分布関数を表し、αｉは前記線形結合における前記確率分布関数のｉ個の係数を表し、ｉは正の整数である、請求項１６に記載の方法。
【請求項１８】
前記モデル関数は、以下の特定の形であり、
【数２】
2013502575.tif 000014
ただし、ｘ（ｔ）は時間依存である前記モデル関数を表し、それぞれの前記確率分布関数Ｖｉは前記確率分布関数Ｄｊ（ｔ）、Ｈｋ（ｔ）、Ｏｌ（ｔ）およびＩｍ（ｔ）のうちの１つによって表わされ、βｊ、ηｋ、δｌおよびιｍはそれぞれの前記確率分布関数の線形結合におけるそれぞれの係数であり、ｊ、ｋ、ｌ、およびｍは正の整数であり、それぞれの前記確率分布関数は前記複数のパラメータの少なくとも１つによって特徴付けられ、
前記Ｄｊ（ｔ）は、前記少なくとも１つのクロマトグラフピークを、それぞれの前記少なくとも１つのバイオマーカの前記少なくとも１つのクロマトグラフ参照ピークにモデル化し、前記少なくとも１つのクロマトグラフピークは、前記ガスクロマトグラフデータ内で少なくとも部分的に合成され、前記少なくとも１つのバイオマーカは、悪い前記医学的状態を実質的に表し、前記Ｈｋ（ｔ）は、前記ガスクロマトグラフデータ内で少なくとも部分的に合成される、前記少なくとも１つのクロマトグラフピークを、それぞれの前記少なくとも１つのバイオマーカと前記悪い医学的状態を表すことが未知であるのいずれか１つである、前記少なくとも１つのクロマトグラフ参照ピークにモデル化し、前記Ｉｍ（ｔ）は、実質的に分解された前記少なくとも１つのクロマトグラフピークをモデル化し、前記Ｏｌ（ｔ）は、前記複数の参照データに実質的に適合しない前記少なくとも１つのクロマトグラフピークと、前記モデル関数のモデリングから得た剰余項のいずれか１つを前記ガスクロマトグラフデータにモデル化する、請求項１７に記載の方法。
【請求項１９】
前記確率分布関数は、以下の一覧、
指数関数変形ガウス（ＥＭＧ）関数、
ガンマ確率分布関数、
多項式変形ガウス関数、
非対称正規分布関数、
カイ分布関数、
正規化した分子速度のマクスウェル・ボルツマン分布関数、
保持時間に対して変形されたマクスウェル・ボルツマン分布関数、
ポアソン分布、
レイリー分布関数
から選択された、請求項１８に記載の方法。
【請求項２０】
前記少なくとも１つのクロマトグラフピークのどれが、前記確率分布関数のどれによってモデル化されるかを決定する手順をさらに含む、請求項１６に記載の方法。
【請求項２１】
前記モデル関数における前記複数のパラメータは、以下の一覧、
前記確率分布関数のモーメント、
前記線形結合における前記確率分布関数の前記係数、
から選択される、請求項１８に記載の方法。
【請求項２２】
前記複数のパラメータを決定する前記手順は、前記複数のパラメータに対する最初の値を評価する手順をさらに含む、請求項１８に記載の方法。
【請求項２３】
評価する前記手順は、勾配降下最適化方法を利用する、請求項２２に記載の方法。
【請求項２４】
前記複数のパラメータを決定する前記手順は、前記モデル関数と前記信号との差の二乗和が最小化されるように、前記複数のパラメータの候補解を反復して評価する手順を含む、請求項２１に記載の方法。
【請求項２５】
前記信号と前記モデル関数との差を取ることにより、時間依存モデルエラーを計算する手順をさらに含む、請求項２４に記載の方法。
【請求項２６】
どの前記少なくとも１つのクロマトグラフピークが合成物であるかを、それに対して前記時間依存モデルエラーが時間依存モデル閾値パラメータを超える時間値を評価することによって決定する手順をさらに含む、請求項２５に記載の方法。
【請求項２７】
前記合成物である前記少なくとも１つのクロマトグラフピークが、それによって少なくとも２つの前記確率分布関数に再モデル化されることを再決定する手順をさらに含む、請求項２６に記載の方法。
【請求項２８】
前記合成物を、前記合成物を少なくとも１つの実質的に分解されたクロマトグラフピークに実質的に分解するように、再決定する前記手順に従って決定された前記少なくとも１つのクロマトグラフピークに、再モデル化する手順さらに含む、請求項２７に記載の方法。
【請求項２９】
前記少なくとも１つの実質的に分解されたクロマトグラフピークを占めるために、前記モデル関数を精製する手順をさらに含む、請求項２８に記載の方法。
【請求項３０】
精製する前記手順は、前記時間依存モデルエラーが最小化されるように、反復して実行される、請求項２９に記載の方法。
【請求項３１】
それぞれの前記少なくとも１つのクロマトグラフピークをモデル化する前記確率分布関数の前記係数は、それぞれの前記少なくとも１つのクロマトグラフピークの最大値を四分位範囲（ＩＱＲ）で割ることによって正規化される、請求項３０に記載の方法。
【請求項３２】
各前記係数に対する閾値パラメータを画定することによって、またそれぞれの前記少なくとも１つのクロマトグラフピークをモデル化する、それぞれの前記確率分布関数である、前記正規化された前記係数のどれが、それぞれの前記閾値パラメータを超えるかを評価することによって、有意なピークを決定する手順をさらに含む、請求項３１に記載の方法。
【請求項３３】
一致の前記程度を前記評価することは、前記少なくとも１つの実質的に分解されたクロマトグラフピークと、前記少なくとも１つのクロマトグラフ参照ピークとの統計的距離測定に従う、請求項３２に記載の方法。
【請求項３４】
一致の前記程度を前記評価することは、それぞれの１つの前記有意ピークと、それぞれの前記少なくとも１つのクロマトグラフ参照ピークとの統計的距離測定に従う、請求項３３に記載の方法。
【請求項３５】
前記統計的距離測定は、カルバック・ライブラー発散に従って決定される、請求項３４に記載の方法。
【請求項３６】
一致の前記程度を前記評価することは、前記時間依存モデルエラーに従う、請求項３２に記載の方法。
【請求項３７】
ガスクロマトグラフデータを分析するためのシステムであって、
試料を複数の成分に分離するクロマトグラフ分離カラムであって、前記クロマトグラフ分離カラムはインレットおよびアウトレットを含む、クロマトグラフ分離カラムと、
前記試料を前記クロマトグラフ分離カラムに提供するための、前記クロマトグラフ分離カラムに前記インレットで結合された試料送達デバイスと、
前記複数の成分の少なくとも一部を検出するための、前記クロマトグラフ分離カラムの前記アウトレットと連通する検出器であって、前記検出器は、前記試料の検出された前記少なくとも一部のそれぞれの特徴の前記ガスクロマトグラフデータを含む、検出器と、
前記ガスクロマトグラフデータおよび複数の参照データを記憶するメモリデバイスと、
前記検出器と結合し、かつ前記メモリデバイスと結合する処理装置であって、前記処理装置は、モデル関数を前記ガスクロマトグラフデータに実質的に適合させるように、前記モデル関数における複数のパラメータを決定し、前記モデル関数は、確率分布関数の線形結合の和として画定され、前記処理装置は、前記ガスクロマトグラフデータと、前記モデル関数と前記ガスクロマトグラフデータとの適合度に従った前記複数の参照データとの一致の程度を評価する、処理装置と、を備えるシステム。
【請求項３８】
前記クロマトグラフ分離カラムのインレットと前記試料送達デバイスとの間を結合するインレット・チャンバをさらに備える、請求項３７に記載のシステム。
【請求項３９】
前記クロマトグラフ分離カラムのアウトレットと前記試料送達デバイスとの間を結合するアウトレット・チャンバをさらに備える、請求項３７に記載のシステム。
【請求項４０】
前記試料送達デバイスは、以下の一覧、
スプリット／スプリットレス（Ｓ／ＳＬ）インジェクタ、
パージトラップ（Ｐ／Ｔ）システム、
ガス源交換システム、
固定相マイクロ抽出（ＳＰＭＥ）システム、
プログラマブル温度気化注入、
熱脱着装置、
微量注射器直接注入、
から選択される、請求項３７に記載のシステム。
【請求項４１】
前記アウトレット・チャンバは、溶出ジェット・インターフェースと噴霧化液体導入システムのいずれか１つを利用する、請求項３９に記載のシステム。
【請求項４２】
前記クロマトグラフ分離カラムは、前記クロマトグラフ分離カラムを特定の温度に実質的に保持するための炉内に実質的に含まれる、請求項３７に記載のシステム。
【請求項４３】
前記検出器は、以下の一覧、
炎イオン化検出器（ＦＩＤ）、
熱伝導度検出器（ＴＣＤ）、
電子捕獲型検出器（ＥＤＣ）、
窒素リン検出器、
炎光光度検出器（ＦＰＤ）、
光イオン化検出器（ＰＩＤ）、
ホール電解伝導度検出器、
放電イオン化検出器（ＤＩＤ）、
パルス放電イオン化検出器（ＰＤＤ）、
質量選択検出器（ＭＳＤ）、
ヘリウムイオン化検出器（ＨＩＤ）、
熱エネルギー（伝導度）分析器／検出器（ＴＥＡ／ＴＣＤ）
から選択される、請求項３７に記載のシステム。
【請求項４４】
前記検出器によって発生された前記信号を増幅させるために、前記検出器と前記処理装置との間で作動可能な増幅器をさらに備える、請求項３７に記載のシステム。
【請求項４５】
前記処理装置は、前記試料を実質的に構成する前記複数の成分の不均衡な濃度の存在を占めるように、前記信号を正規化する、請求項３７に記載のシステム。
【請求項４６】
前記ガスクロマトグラフデータは、少なくとも１つのクロマトグラフピークを含む、請求項３７に記載のシステム。
【請求項４７】
前記複数の参照データはデータベースを形成し、前記データベース内の前記複数の参照データは、少なくとも１つの公知の化合物に実質的に対応する、少なくとも１つのクロマトグラフピークを含む、請求項４６に記載のシステム。
【請求項４８】
前記複数の参照データは、前記少なくとも１つの公知の化合物に対応する質量分析データと較正される、請求項４７に記載のシステム。
【請求項４９】
前記少なくとも１つの公知の化合物の特定のものは、それから前記試料が獲得される生物の少なくとも１つの生物学的状態に実質的に表わされる、少なくとも１つのバイオマーカを画定する、請求項４７に記載のシステム。
【請求項５０】
前記メモリデバイスは、前記少なくとも１つのバイオマーカを記憶する、請求項４９に記載のシステム。
【請求項５１】
前記処理装置は、前記少なくとも１つのクロマトグラフピークを、前記少なくとも１つのクロマトグラフピークの保持時間値を、前記少なくとも１つのクロマトグラフ参照ピークのそれぞれの参照保持時間値と比較することによって、前記少なくとも１つのクロマトグラフ参照ピークに登録する、請求項４９に記載のシステム。
【請求項５２】
前記処理装置は、前記少なくとも１つのクロマトグラフピークの時間領域における状態位置を、前記複数の参照データ内のそれぞれの前記少なくとも１つの参照クロマトグラフピークのそれぞれの位置と比較することによって前記登録する、請求項５１に記載のシステム。
【請求項５３】
前記処理装置は、ｓ（ｆ（ｔ））が対応する前記複数の参照データｒ（ｔ）と最大限に一致するように選択された、変形関数ｆ（ｔ）を利用することによって前記登録する、請求項５１に記載のシステム。
【請求項５４】
前記変形関数は、現在現れている前記少なくとも１つのクロマトグラフピークを考慮に入れるように選択される、請求項５３に記載のシステム。
【請求項５５】
前記処理装置は、それぞれの少なくとも１つの公知のクロマトグラフピークを発生するように、前記試料に追加され、少なくとも１つの公知の保持時間を有する、少なくとも１つの参照化学物質を前記登録し、それによって前記変形関数が前記少なくとも１つの公知のクロマトグラフピークに従って選択される、請求項５３に記載のシステム。
【請求項５６】
前記処理装置は、その一時的帰属に従った前記少なくとも１つのクロマトグラフピークの分類を、前記参照データと比較することによって実行する、請求項５３に記載のシステム。
【請求項５７】
前記処理装置は、前記少なくとも１つのクロマトグラフピークがそれぞれの前記確率分布関数によってモデル化されるように、前記分類に従って、前記モデル関数をさらに構築する、請求項５６に記載のシステム。
【請求項５８】
前記モデル関数は以下の一般の形であり、
【数３】
2013502575.tif 000015
ただし、Ｖｉはｉ個の前記確率分布関数を表し、αｉは前記線形結合における前記確率分布関数のｉ個の係数を表し、ｉは正の整数である、請求項５７に記載のシステム。
【請求項５９】
前記モデル関数は、以下の特定の形であり、
【数４】
2013502575.tif 000016
ただし、ｘ（ｔ）は時間依存である前記モデル関数を表し、各前記確率分布関数Ｖｉは前記確率分布関数Ｄｊ（ｔ）、Ｈｋ（ｔ）、Ｏｌ（ｔ）およびＩｍ（ｔ）のうちの１つによって表わされ、βｊ、ηｋ、δｌおよびιｍは各前記確率分布関数と線形結合におけるそれぞれの係数であり、ｊ、ｋ、ｌ、およびｍは正の整数であり、各前記確率分布関数は前記複数のパラメータの少なくとも１つによって特徴付けられ、
前記Ｄｊ（ｔ）は前記少なくとも１つのクロマトグラフピークをそれぞれの前記少なくとも１つのバイオマーカの前記少なくとも１つのクロマトグラフ参照ピークにモデル化し、前記少なくとも１つのクロマトグラフピークは、前記ガスクロマトグラフデータ内で少なくとも部分的に合成され、前記少なくとも１つのバイオマーカは、悪い前記医学的状態を実質的に表し、前記Ｈｋ（ｔ）は、前記ガスクロマトグラフデータ内で少なくとも部分的に合成される、前記少なくとも１つのクロマトグラフピークを、それぞれが前記少なくとも１つのバイオマーカであるか、または前記悪い医学的状態を表すことが未知であるかのうちのいずれか１つである、前記少なくとも１つのクロマトグラフ参照ピークにモデル化し、前記Ｉｍ（ｔ）は、実質的に分解された前記少なくとも１つのクロマトグラフピークをモデル化し、前記Ｏｌ（ｔ）は、前記複数の参照データに実質的に適合しない前記少なくとも１つのクロマトグラフピークか、または前記モデル関数のモデリングから獲得された剰余項のいずれか１つを、前記ガスクロマトグラフデータにモデル化する、請求項５８に記載のシステム。
【請求項６０】
前記確率分布関数は、以下の一覧、
指数関数変形ガウス（ＥＭＧ）関数、
ガンマ確率分布関数、
多項式変形ガウス関数、
非対称正規分布関数、
カイ分布関数、
正規化した分子速度のマクスウェル・ボルツマン分布関数、
保持時間に対して変形されたマクスウェル・ボルツマン分布関数、
ポアソン分布、
レイリー分布関数
から選択される、請求項５９に記載のシステム。
【請求項６１】
前記処理装置は、前記少なくとも１つのクロマトグラフピークの前記状態位置と、前記それぞれの前記少なくとも１つのクロマトグラフ参照ピークのそれぞれの位置との評価を含むように一致の前記程度を前記評価する、請求項５９に記載のシステム。
【請求項６２】
前記処理装置は、前記確率分布関数のどれかによってモデル化される、前記少なくとも１つのクロマトグラフピークのどれかを決定する、請求項５９に記載のシステム。
【請求項６３】
前記処理装置は、前記確率分布関数Ｄｊ（ｔ）、Ｈｋ（ｔ）、Ｏｌ（ｔ）およびＩｍ（ｔ）のどれが、それによって前記確率分布関数の前記一覧内の前記確率分布関数にモデル化されるかを決定する、請求項６０に記載のシステム。
【請求項６４】
前記モデル関数における前記複数のパラメータは、以下の一覧、
前記確率分布関数のモーメント、
前記線形結合における前記確率分布関数の前記係数、
から選択される、請求項５９に記載のシステム。
【請求項６５】
前記処理装置は、前記複数のパラメータを、前記複数のパラメータに対する最初の値を評価することによって前記決定する、請求項５９に記載のシステム。
【請求項６６】
前記処理装置は、前記複数のパラメータを、勾配降下最適化手法を利用することによって評価する、請求項６５に記載のシステム。
【請求項６７】
前記処理装置は、前記モデル関数と前記信号との差の二乗和が最小化されるように、前記複数のパラメータを、前記複数のパラメータの候補解を反復して評価することによって前記決定する、請求項６５に記載のシステム。
【請求項６８】
前記処理装置は、時間依存モデルエラーを、前記信号と前記モデル関数との差を取ることによって計算する、請求項６７に記載のシステム。
【請求項６９】
時間依存モデル閾値パラメータは画定され、前記処理装置は、それに対して前記時間依存モデルエラーが、前記時間依存モデルエラー閾値パラメータを超える時間値を評価することによって、前記少なくとも１つのクロマトグラフピークのどれが合成物であるかを決定する、請求項６８に記載のシステム。
【請求項７０】
前記処理装置は、前記合成物である前記少なくとも１つのクロマトグラフピークが再モデル化され、それによって少なくとも２つの前記確率分布関数を再決定する、請求項６９に記載のシステム。
【請求項７１】
前記処理装置は、前記合成物であり、それぞれの前記確率分布関数によって以前モデル化された前記少なくとも１つのクロマトグラフピークが、前記少なくとも２つの前記確率分布関数によって前記モデル化されるように、前記モデル関数を精製する、請求項７０に記載のシステム。
【請求項７２】
前記処理装置は、前記モデル関数を、前記時間依存モデルエラーが最小化されるように、各前記少なくとも１つのクロマトグラフピークに対して繰り返し精製する、請求項７１に記載のシステム。
【請求項７３】
前記処理装置は、前記モデル関数を、前記合成物である前記少なくとも１つのクロマトグラフピークを、少なくとも１つの実質的に分解されたクロマトグラフピークに実質的に分解するように繰り返し精製する、請求項７２に記載のシステム。
【請求項７４】
それぞれの前記少なくとも１つのクロマトグラフピークをモデル化する前記確率分布関数の前記係数は、前記少なくとも１つのクロマトグラフピークのそれぞれの最大値を四分位範囲（ＩＱＲ）で割ることにより、前記処理装置によって正規化される、請求項７３に記載のシステム。
【請求項７５】
少なくとも１つの閾値パラメータは、各前記係数に対して画定され、前記処理装置は、それぞれの前記少なくとも１つのクロマトグラフピークをモデル化する、それぞれの前記確率分布関数である、どの前記正規化された前記係数が、それぞれの前記少なくとも１つの閾値パラメータを超えるかを評価することによって、有意なクロマトグラフピークを決定する、請求項７４に記載のシステム。
【請求項７６】
前記処理装置は、前記少なくとも１つの実質的に分解されたクロマトグラフピークと、それぞれの少なくとも１つのクロマトグラフ参照ピークとの統計的距離測定に従って、一致の前記程度を評価する、請求項７５に記載のシステム。
【請求項７７】
前記処理装置は、それぞれの１つの前記有意なクロマトグラフピークと、それぞれの少なくとも１つのクロマトグラフ参照ピークとの統計的距離測定に従って、一致の前記程度を前記評価する、請求項７５に記載のシステム。
【請求項７８】
前記統計的距離測定は、前記カルバック・ライブラー発散に従って決定される、請求項７６に記載のシステム。
【請求項７９】
前記統計的距離測定は、前記カルバック・ライブラー発散に従って決定される、請求項７７に記載のシステム。
【請求項８０】
前記処理装置は、前記時間依存モデルエラーに従って一致の前記程度を評価する、請求項７５に記載のシステム。
【発明の詳細な説明】
【技術分野】
【０００１】
開示された技法は、一般にガスクロマトグラフィに関し、詳細には、ガスクロマトグラフデータを分析する方法およびシステムに関する。
【背景技術】
【０００２】
ガス液体分布クロマトグラフィ（ＧＬＰＣ）、気相クロマトグラフィ（ＶＰＣ）、ガス液体クロマトグラフィ、またより簡略にガスクロマトグラフィ（ＧＣ）としても公知であるものは、化学分解なしに気化することが可能な、化学合成液または化合物を分解および分析するために利用される分析化学技法の名前である。ＧＣは、ガス状合成物などの試料をその化合物に分解するために利用され、成分の相対量が決定されることがある。ＧＣはまた、物質の純度の試験のため、化合物同定の支援のため、および合成物から純粋化合物を抽出するために利用されることもある。ＧＣは、概してガスクロマトグラフィまたはガス分離器と呼ばれる機器によって実行される。概して、ＧＣ技法は、気化した形（たとえば、直接噴射、パージトラップ（Ｐ／Ｔ）法を介して）の試料を、異なる固体固定相または液体固定相（すなわち、吸収剤）で覆われた不活性固体支持体を有するように内部で構築された、ＧＣカラム（以下「カラム」と言う）の一端に導入するものである。移動相（すなわち、ヘリウムなどの搬送ガス）は、試料をカラムを介して掃引するために利用される。試料などの固定相と様々な相互作用をする試料の異種成分は、カラムを介して掃引され、それにより各成分が異なる時に溶出される（すなわち、成分の保持時間として公知である）。カラムを通る試料の異なる化合物の割合は、それらの化学的性質および物理的性質ならびに固定相との相互作用に依存する。成分のそれぞれの保持時間に依存して成分がカラムの他端から異なった時に出現するときに、様々な検出技術を利用する検出器によって検出される可能性がある。通常、検出器は、試料内の成分濃度に応答して電気信号を発生する。クロマトグラフデータは通常、時間関数（保持時間）として検出器応答（濃度）のグラフ（たとえば、スペクトル）の形で表わされる（クロマトグラムと呼ばれる）。したがって、各試料に対して、ＧＣは、ピークのスペクトルを有する対応するクロマトグラムを生成し、ピークは、異なった時にカラムから溶出する試料内の被分析物の存在を表す。試料のクロマトグラムにおけるスペクトルパターンの存在を定量的に分析することによって、それらを公知の被分析物濃度を含む一定基準と比較することによって、試料内の被分析物濃度を決定することが可能である。
【０００３】
したがって、ＧＣは、生物医学的応用、環境問題への応用、法医学的分析、石油化学の分析などの広い分野で利用される。たとえば、ＧＣは、揮発性有機化合物（ＶＯＣ）に対して人および動物の呼気の分析に利用される。概してＶＯＣは、様々な材料（たとえば、清掃用品、ペンキ、農薬、建築材料）によって排出される気体または蒸気であり、これらは生物に健康への悪影響をもたらすことがある。人は、吸入、経口、経皮吸収などを通してＶＯＣに必然的に曝されている。人の呼気を検査することによって、それには数百のＶＯＣが必然的に含まれており、体内に有害な化学物質が増大する可能性の示唆を提供することが可能である。したがって人の呼気内に検出されたＶＯＣは、肺癌、乳癌、糖尿病、および統合失調症などの疾患の存在の可能性の試験に生物学的マーカ（すなわち、バイオマーカ）として働く可能性がある。
【０００４】
しかし、クロマトグラフデータの分析、具体的には試料を成分に完全分離および完全分解することは、クロマトグラムに存在する重複ピークの現象が発生するため困難であり得ることは公知である。基本的には、この問題が起きるのは、試料の２つ以上の異なる成分が、実質的に同じ割合で溶出し（すなわち、それらは実質的に同様の保持時間を有する）、それらが単一成分であるかのように検出されるときである。
【０００５】
様々なタイプの装置およびクロマトグラフ分離方法は、当技術分野に公知である。重複クロマトフラフピークの検出を高めるそのような方法の１つは、多次元ガスクロマトグラフィ（本明細書ではＭＤＧＣと略す）を使用するものであり、試料の成分は、異なる特性を有する２つ以上のカラムを使用する２つ以上の分離ステップの影響を受ける。たとえば、二次元（２−Ｄ）ガスクロマトグラフィ（本明細書では２Ｄ−ＧＣと略す）では、さらなる分析を必要とするクロマトグラム内の領域が強化され（「中間留分」）、第２のカラム上で検査される。別の方法は、総合的に２Ｄ−ＧＣ（本明細書ではＧＣｘＧＣと略す）を使用するものであり、これは第１のカラムおよび流出物の一部の周期的な再注入から異なる特性を有する第２のカラムへの流出の収集に基づく。この方法では、第１のカラムからの流出は、試料全体が分離ステップのすべてに影響を及ぼす（すなわち、寸法）一方で、分離をそれぞれの前のステップから保護するように、複数回抽出される。この方法は、第１のカラムと第２のカラムを連結するインターフェースに依拠し、これによって周期的な注入を引き起こすことが可能になる。それにもかかわらず、これらの技法の使用は、追加機器ならびにスペクトルデータの複数のチャネルの分析を伴い、これらは個々の試料を含むすべての成分の完全な同定を最終的に保証しない。
【０００６】
ガスクロマトグラフデータを分析する方法およびシステムも、当技術分野では公知である。たとえば、クロマトグラフピークの形状を特徴付ける際に指数関数変形ガウス（ＥＭＧ）関数を利用することが当技術分野では公知であり、この理論的正当性はクロマトフラフピークが通常非対称の特性を示すという事実にある。他の方法は、デコンボリューション技法、反復ターゲット変換因子分析（ＩＴＴＦＡ）、パターン認識およびニューラルネットワーク技法などを含む。「Ｍｅｔｈｏｄ  ａｎｄ  Ａｐｐａｒａｔｕｓ  ｆｏｒ  Ｃｈｒｏｍａｔｏｇｒａｐｈｉｃ  Ｄａｔａ  Ｐｒｏｃｅｓｓｉｎｇ」という名称のイトウ（伊藤）らに対する米国特許第７，４０３，８５９Ｂ２号は、複数の重複ピークを含むクロマトグラムに対して、線形最小二乗法を利用することによって、曲線当てはめを促進するための液体クロマトグラフ分析器を対象とする。液体クロマトグラフ分析器は、カラム、試料供給部、液体ポンプ、制御装置、試料、および検出器を含む。試料供給部は、液体ポンプとカラムとの間に配置される。溶離液は、制御装置からの命令により、液体ポンプを使用してカラムに汲み出される。試料は、制御装置の命令によりサンプラから溶離液に供給される。試料は、カラムによって分離され、検出器によって検出される。検出されたデータのクロマトグラムは、分析されるために制御装置に送信される。
【０００７】
制御装置によって処理されるクロマトグラムのデータは、当てはめを実行するための時間間隔の指定、波形関数の選択、重み付けパターンの選択、当てはめ方向の選択、当てはめ実行ボタンのクリック、結果の表示および出力を含む手順によって実行される。最初は、特定の選択されたクロマトグラムに対して、クロマトグラムにおける時間間隔は、開始時刻および終了時刻を入力することによって当てはめのために選択される。したがって、ガウス関数すなわちＥＭＧ関数は、当てはめのための波形関数として使用される。重み付け関数の選択は、重み付け関数のグラフ表示をクロマトグラム上にポインティングデバイスを介して重ね合わせるものである。当てはめ方向の選択は、処理がクロマトグラムにおける選択された時間間隔の正面から実行されるか、または背面から実行されるかの方向の設定をするものである。当てはめ処理（実行）は、当てはめのために波形関数を利用し、この波形関数はガウス関数と基線（すなわち、線形一次方程式）の和である。当てはめ処理は、ガウス関数における当てはめパラメータが、計測されたクロマトグラムの信号強度における波形関数とそれぞれの点との差の二乗和を最小化するために決定されるように、最小二乗法を利用する。
【発明の概要】
【課題を解決するための手段】
【０００８】
開示された技法の目的は、ガスクロマトグラフィの新規のシステムおよび方法を提供することであり、これは先行技術の不都合を克服する。したがって、開示された技法に従って、それぞれの試料のガスクロマトグラフデータと参照データとの一致の程度を決定する方法が提供される。方法は、ガスクロマトグラフデータを獲得し、モデル関数がガスクロマトグラフデータに実質的に適合するように、モデル関数における複数のパラメータを決定し、モデル関数とガスクロマトグラフデータとの適合度に従って、一致の程度を評価する手順を含む。モデル関数は、確率分布関数の線形結合の和として画定される。
【０００９】
したがって、開示された技法の別の態様に従って、ガスクロマトグラフデータの分析システムが提供される。システムは、試料を複数の成分に分離するためのクロマトグラフ分離カラム、試料送達デバイス、検出器、メモリデバイス、および処理装置を含む。クロマトグラフ分離カラムは、インレットおよびアウトレットを含む。試料送達デバイスは、試料をクロマトグラフ分離カラムに提供するために、そのインレットでクロマトグラフ分離カラムと結合される。検出器は、クロマトグラフ分離カラムのアウトレットと連通し、複数の成分の少なくとも一部を検出し、試料の検出された部分の特性のそれぞれのガスクロマトグラフデータを含む信号を発生する。メモリデバイスは、処理装置と結合され、ガスクロマトグラフデータおよび複数の参照データを記憶する。処理装置は、検出器と結合され、モデル関数がガスクロマトグラフデータと実質的に適合するように、モデル関数における複数のパラメータを決定する。モデル関数は、確率分布関数の線形結合の和として画定される。処理装置は、ガスクロマトグラフデータと、モデル関数とガスクロマトグラフデータとの適合度に従った複数の参照データとの一致の程度を評価する。
【００１０】
本開示の技法は、図面に関連して以下の詳細な説明から、より完全に理解され認識されよう。
【図面の簡単な説明】
【００１１】
【図１】開示された技法の実施形態により構築され作動可能なガスクロマトグラフデータの分析システムの概略図である。
【図２Ａ】図１に示されたシステムによって獲得された、代表的なクロマトグラムの概略図である。
【図２Ｂ】図２Ａのクロマトグラムによりモデル化された時間依存モデル関数の最初の評価のグラフの概略図である。
【図２Ｃ】時間依存モデルエラー閾値関数のグラフと共に描かれた、図２Ｂの最初に評価されたモデル関数がもたらした、計算された時間依存モデルエラーのグラフの概略図である。
【図２Ｄ】図２Ａのクロマトグラムによりモデル化された、図２Ｂの時間依存モデル関数の精密評価の概略図である。
【図３Ａ】その異なる成分が、開示された技法の別の実施形態により構築され作動可能な、所与の試料を構成する重複クロマトグラフピーク内で成分を分解し同定する方法を示す概略構成図である。
【図３Ｂ】図３Ａの方法の続きを示す概略構成図である。
【発明を実施するための形態】
【００１２】
開示された技法は、その異なる成分が、確率密度関数の線形結合の和として画定されるモデル関数を利用することによって、所与の試料を構成する重複クロマトグラフピーク内で成分を分解し同定する方法およびシステムを提供することによって、先行技術の不都合を克服する。所与の試料を構成する化合物に関連したクロマトグラフデータは、一次元ＧＣ（本明細書では１Ｄ−ＧＣと略す）ガスクロマトグラフ分離技法（すなわち、ＭＤＧＣおよび２Ｄ−ＧＣなどの多次元ガスクロマトグラフ技法と対照的に）によって獲得される。試料のクロマトグラム内の顕著な特徴（たとえば、クロマトグラフピーク）は、分類され、その後、実施された分類に従った確率密度関数の特定のタイプによって表される（すなわち、モデル化される）可能性があるような方法で数学的に分解される。各確率密度関数を特徴付ける複数のパラメータは、最適化技法によって評価され、その後、確率密度関数の線形結合の和内の複数の線形係数パラメータが、最小二乗法によって決定される。時間依存モデルエラー関数およびモデルエラー閾値パラメータは画定される。合成されたと疑われるクロマトグラフピークは、それに対して時間依存モデルエラー閾値パラメータが時間依存モデルエラーを超える時間値を評価することによって実質的に決定される（すなわち、算定され、評価される）。精製されたモデル関数は、複数の確率密度関数によって合成されたと疑われるピークを再モデル化し、それぞれのピークの対応するモデルエラーを考慮に入れることによって構築され、それによって合成クロマトグラフピークを分解する。最適化技法は、最小二乗法を最小化するように、モデル関数がクロマトグラフデータに実質的に適合するために反復される。各反復で、精製されたモデル関数は、モデルエラーが最小化されるまで、前のモデル関数を置換する。開示された技法は、分析された試料にインタレストの特定のバイオマーカの有無を推測するために、参照ピーク間の一致の程度、データベース内に記憶されたその一致の程度の情報、ならびに試料の新発見され、分解されたピークを含む複数のピークを評価する。概して、開示された技法は通常、特定の悪い医学的状態（たとえば、癌）を有すると疑われる個体から収集された、呼気試料内のマルチ・バイオマーカの存在の確率判定表示を提供するために実施されることがある。
【００１３】
発明を実施するための形態、図、および特許請求の範囲の全体を通して使用される用語「確率密度関数」および「確率分布関数」は、交換可能である。ここで、開示された技法の一実施形態により、構築され作動可能なガスクロマトグラフデータ（概ね１００を参照）の分析システムの概略図である図１を参照する。システム１００は、クロマトグラフ分離カラム１０２、試料送達デバイス１０４、検出器１０６、処理装置１０８、およびメモリデバイス１１０を含む。システム１００は、恣意的にインレット・チャンバ１１２およびアウトレット・チャンバ１１４をさらに含むことがある。クロマトグラフ分離カラム１０２は、インレット１１６およびアウトレット１１８を含む。試料送達デバイス１０４は、インレット１１２を介してクロマトグラフ分離カラム１０２と連結される。あるいは、試料送達デバイス１０４は、インレット・チャンバ１１２（図１に示されている）を介してクロマトグラフ分離カラム１０２と連結されることがある。検出器１０６は、アウトレット１１４を介してクロマトグラフ分離カラム１０２と連結される。あるいは、検出器１０６は、アウトレット・チャンバ１１４（図１に示されている）を介してクロマトグラフ分離カラム１０２と連結される。検出器１０６は、処理装置１０８と連結され、次いで処理装置１０８は、メモリデバイス１１０と連結される。
【００１４】
まず、分析される試料（図示せず）（たとえば、呼気試料）は、試料送達デバイス１０４の中に提供される。あるいは、試料は、プローブ採取装置（ＰＳＤ）などの密閉吸着管（図示せず）内にまず収集され（すなわち、試料収集デバイスを介して）、その後試料送達デバイス１０４に分注される可能性がある。インレット・チャンバ１１２が利用されない場合は、試料送達デバイス１０４は、試料をヘリウム、窒素、アルゴン、および風乾などの搬送ガス（図示せず）の連続流の中に導入し、これは試料をクロマトグラフ分離カラム１０２のインレット１１６（「オンカラム・インレット」と呼ばれる）に掃引する。インレット１１６への試料の導入は、当技術分野で公知の自動採取装置および自動インジェクタの使用などを通じて自動的に達成される可能性がある。インレット・チャンバ１１２が利用される場合は、インレット・チャンバ１１２は概して、試料の揮発を促進するために、通常Ｓ／ＳＬ（スプリット／スプリットレス）インジェクタ（すなわち、試料送達デバイスのタイプ）を使用して、蒸発室（すなわち、蒸発室は温度制御されている）として機能する。たとえば、試料送達デバイスの他のタイプおよび技法は、Ｐ／Ｔ（パージトラップ）システム、ガス源交換システム、ＳＰＭＥ（固定相マイクロ抽出）、ＰＴＶ（プログラマブル温度気化）注入、微量注射器直接注入、熱脱着装置などに利用される。このような実装の一部に対して、システム１００は、搬送ガスを供給するための搬送ガスタンク（図示せず）をさらに含むことがあり、流れ制御装置、弁、圧力センサなどの、この目的のために他の様々な関連機器（図示せず）も利用されることがある。
【００１５】
試料がクロマトグラフ分離カラム１０２を通過すると、試料の様々な成分（図示せず）が吸収によって分離され、それらの成分がアウトレット１１８からアウトレット・チャンバ１１４の中に出現するのとは異なる割合で溶出する。たとえば、アウトレット・チャンバ１１４は、溶出ジェット・インターフェース、噴霧化液体導入システムなどを含むことがある。噴霧化液体導入システムでは、溶出ガス混合物が噴霧され（すなわち、エアロゾルとして）、検出器１０６の中にあるいはアウトレット・チャンバ１１４の一部の中に直接噴霧され、したがって均一に改善されたエアロゾルを生成する。たとえば、溶出ジェットまたは噴霧化液体導入システムを利用することによって、充填毛管カラムは炎イオン化、無炎熱電子イオン化、測光型検出器などに基づいて、直接検出器に結合されることがある。クロマトグラフ分離カラム１０２は、毛管型カラムであることが好ましく、概して充填カラム型の感度より比較的高い感度を提供する（すなわち、全体として検出されたクロマトグラフピークはより高くはるかに鋭いので、それによってより良い信号対雑音比が得られる）。しかし、開示された技法は、他のタイプのカラムが利用されることがある（たとえば、充填カラム、内部加熱マイクロＦＡＳＴカラム、マイクロ充填カラム）ように、特定のタイプのクロマトグラフカラムに限定されない。分子吸着および試料がクロマトグラフ分離カラム１０２を通じて発達する割合は、温度に依存するので、通常クロマトグラフ分離カラム１０２の温度を制御する必要がある。このような目的で、炉（図示せず）が通常クロマトグラフ分離カラム１０２を所望の温度で収納し保持するために利用される。炉の温度は、クロマトグラフ分離カラム１０２を、実行される各分析に対して特定の等温状態に通常維持するために電子的に制御される。
【００１６】
溶出（すなわち、流出物）がクロマトグラフ分離カラム１０２から現れると、試料を構成する成分の少なくとも画分が（アウトレット１１８と連通するように配置された）検出器１０６によって検出される。検出器の多くのタイプはＧＣに使用され得る。ＧＣ検出器は、それらの感度（すなわち、相対的に、他の要素または化合物より特定の要素または化合物に応答するための検出器の可能な測定）、および検出器が濃度依存検出器であるか、または量流検出器であるかなどの他の要因に従って分類され得る。たとえば、選択的検出器は、化学または物理に共通の特性を有する多様な成分に応答するが、非選択的（万能）検出器は、搬送ガスとは別に実質的にすべての成分に応答する。開示された技法によって利用され得る様々なタイプの検出器は、炎イオン化検出器（ＦＩＤ）、熱伝導度検出器（ＴＣＤ）、電子捕獲型検出器（ＥＤＣ）、窒素リン検出器、炎光光度検出器（ＦＰＤ）、光イオン化検出器（ＰＩＤ）、ホール電解伝導度検出器、放電イオン化検出器（ＤＩＤ）、パルス放電イオン化検出器（ＰＤＤ）、質量選択検出器（ＭＳＤ）、ヘリウムイオン化検出器（ＨＩＤ）、熱エネルギー（伝導度）分析器／検出器（ＴＥＡ／ＴＣＤ）などを含む。ＴＣＤは、万能選択性を有する濃度依存検出器の一例である。ＦＰＤは、流量タイプの選択検出器の一例であり、その選択性はリン、錫、ゲルマニウム、硫黄、セレンなどに関する。検出器１０６は、通常、試料内の成分の検出された濃度に応答して電気信号ｓ（ｔ）を時間関数として発生する。この電気信号は、処理および分析のために処理装置１０８に送信される。あるいは、システム１００は、検出器１０６と処理装置１０８との間の操作で、検出器１０６によって発生される電気信号を増幅するために、増幅段階（図示せず）をさらに含むことがある。増幅段階は、前置増幅器、増幅器、電子増幅器（ＥＭＡ）などによって実施され得る。
【００１７】
電気信号は、クロマトグラフデータ（図示せず）の表示であり、処理装置１０８はクロマトグラフデータを記憶および検索のためにメモリデバイス１１０に転送する。処理装置１０８によって分析される電気信号のそれぞれのクロマトグラフデータは、クロマトグラムの形に配置され提供され得る。次にさらに図２Ａおよび図２Ｂを参照する。図２Ａは、図１に示されたシステムによって獲得された、代表的クロマトグラム（概ね２００を参照）の概略図である。図２Ｂは、図２Ａのクロマトグラムに従って、モデル化された時間依存モデル関数の最初の評価のグラフの概略図である。クロマトグラム２００は、特定の試料のクロマトグラフ分離のグラフ記録を表し、デカルト座標系に表され、その垂直軸は検出された溶出物質の濃度の程度（すなわち、検出器応答）を時間関数として表す。クロマトグラム２００は、それぞれが出現した異なる成分（すなわち、ＧＣによって分離されなかった）の特定の成分または組合せを表す、複数のクロマトグラフピーク２０２、２０４、２０６、２０８、２１０、２１２、および２１４を表す。検出された電気信号ｓ（ｔ）は、所与の試料を構成する成分の不均衡な濃度の存在を占める（たとえば、補償する）ために正規化されることが可能であり、たとえばそれは、他の化学物質から、または利用される検出器の特有の事前選択からなどの外的影響に起因することがある。
【００１８】
メモリデバイス１１０は、公知の化合物に対応する複数の参照ＧＣデータのデータベース（図示せず）を記憶する。具体的には、データベースは、ピークのセットＤ’に対応するデータを記憶し、この場合、このセット内の各要素は、特定の悪い医学的状態（たとえば、疾病、感染症）に関連した公知の化合物のクロマトグラフピークを表す。データベース内の化合物の単一または結合に対応するデータは、バイオマーカ（図示せず）を画定するためにグループ化されることがある。たとえば、サブセット｛ｄ８’，ｄ３４’，ｄ３７１’｝⊂Ｄ’は、特定の疾病のバイオマーカを画定することがある。バイオマーカは概して、試料のクロマトグラフデータにその質的存在および量的存在の有無が、生物（たとえば、人、犬、猫）の特定の生物学的状態の指針である、１つの成分（または複数の成分）を指す。データベースはさらに、バイオマーカのセットＭ’を記憶し、この場合、各バイオマーカの成分はサブセットＤ’として画定される。本明細書に準備された指標は、参照データを示す。前述の例を鑑み、バイオマーカｍ１’⊂Ｍ’は、ｍ１’＝｛ｄ８’，ｄ３４’，ｄ３７１’｝と画定され得る。同様に、データベースは、ピークのセットＨ’に対応するデータを記憶し、この場合、このセット内の各要素は、特定の悪い医学的状態（たとえば、通常健康な個体に現れる）に関連すると未知であるか、または特定の悪い医学的状態に関連すると公知であるかのいずれか１つであるが、それにも関わらず検出のためのインタレストではない、化合物のクロマトグラフピークを表す。
【００１９】
データベースは、学習段階および較正段階で最初に構築される。この段階では、複数の公知の化合物および恐らく未知の化合物からのクロマトグラフデータ（すなわち、クロマトグラム）が獲得され、最終的にそこからクロマトグラフデータを参照として構成することになる。具体的には、複数のＶＯＣからのクロマトグラフデータ（たとえば、ピーク）は、（すなわち、検出において）インタレストの特定の医学的状態を診断された個体から（たとえば、呼気試料を通して）獲得され、インタレストの医学的状態を特徴付ける（すなわち、バイオマーカ）クロマトグラフデータを同定するためにインタレストの特定の医学的状態を有さないと診断された個体から獲得された複数のＶＯＣと比較する。質量分析法（ＭＳ）ならびに分光法は、較正方法としてこの段階に利用されることがあり、この場合、収集された各試料の元素組成は、試料の各成分のそれぞれの保持時間と比較され、それぞれの保持時間に関連する。概して、「健康」な個体および「不健康」な個体の両方からのＶＯＣのクロマトグラフデータが、収集され、分析され、データベースに記憶される。クロマトグラフ参照データの分析は、たとえば主成分分析（ＰＣＡ）などにより、クロマトグラフピークの検出によって実行され得る。それぞれの検出されたクロマトグラフピークは、本明細書の以下により詳細に説明される方法に従って、特定の確率密度関数によってモデル化される可能性がある。
【００２０】
開示された技法は、その異なる成分が、モデル関数を確率密度関数（確率分布関数とも呼ばれる）の線形結合として利用することにより、所与の試料を構成する重複クロマトグラフピーク内で成分を分解し同定する。以下でＶｉは一般形を有する。
【数１】
2013502575.tif 000003
上式で、αｉは確率密度関数の係数であり、ｉは正の整数である。具体的には、開示された技法に従って、式（１）における確率密度関数の線形結合は、以下の形を有する確率密度関数の線形結合に分解され得ると仮定される。
【数２】
2013502575.tif 000004
上式で、ｘ（ｔ）は検出器１０６によって獲得された、電気信号ｓ（ｔ）をモデル化するために利用された時間依存モデル関数を表す。電気信号ｓ（ｔ）は、修正（たとえば、増幅、事前処理）を施された可能性があることに留意されたい。Ｄｊ（ｔ）は、セットＤ’（すなわち、特定の悪い医学的状態に関連した）における特定のクロマトグラフピークに対応する可能性を有する、それぞれのクロマトグラフピーク（すなわち、実質的に分解されなかった）をモデル化する第ｊ番目の時間依存確率密度関数を表す。それぞれの第ｋ番目の時間依存確率密度関数Ｈｋ（ｔ）は、セットＨ’（すなわち、特定の悪い医学的状態に関連すると未知である、または特定の悪い医学的状態に関連すると公知であるかのいずれか１つであるが、それにも関わらず検出のためのインタレストではない）における特定のクロマトグラフピークに対応する可能性を有する、クロマトグラフピーク（すなわち、概ね部分的に分解された）をモデル化する。分離されたクロマトグラフピーク（すなわち、クロマトグラフピークが概ね分解された）が特定の医学的状態に関連すると公知であろうと未知であろうと、それらは第ｍ番目の時間依存確率密度関数Ｉｍ（ｔ）（すなわち、セットＨ’またはＤ’のどちらかにおける特定のクロマトグラフピークに対応する可能性を有する）によってモデル化される。Ｏｌ（ｔ）は、未知のクロマトグラフピーク（すなわち、データベースの一部ではない分類不能なクロマトグラフデータ）、またはモデリング手順から生じた剰余項をそれぞれにモデル化する第ｌ番目の時間依存確率密度関数を表す。スカラー重みβｊ、ηｋ、δｌおよびιｍは、それぞれの確率密度関数Ｄｊ（ｔ）、Ｈｋ（ｔ）、Ｏｌ（ｔ）およびＩｍ（ｔ）との線形結合における係数である。指数ｊ、ｋ、ｌおよびｍは正の整数である。
【００２１】
様々な確率密度関数は、ＥＭＧ、ガンマ分布（すなわち、その確率密度関数）、多項式変形ガウス分布、非対称正規分布、カイ分布、ポアソン分布、正規化した分子速度のマクスウェル・ボルツマン分布（すなわち、３個の自由度（ＤＯＦ）を有するカイ分布）、保持時間に対して修正されたマクスウェル・ボルツマン分布、レイリー分布（すなわち、２個のＤＯＦおよび標準偏差σ＝１を有するカイ分布）などのＤｊ（ｔ）、Ｈｋ（ｔ）、Ｏｌ（ｔ）およびＩｍ（ｔ）のために使用される可能性がある。
【００２２】
モデル処理は、最初に分離されたクロマトグラフピーク（すなわち、ピーク２０２および２１２）をモデル化することがあり、それらはクロマトグラム２００に現れる。これらのピークに対して、また概して分離されたピークであると疑われる各ピークｍに対して、処理装置１０８は、それぞれの時間依存確率密度関数Ｉｍ（ｔ）を見出し、Ｉｍ（ｔ）はそのピークに対する数学的モデルとして働く。使用され得る時間依存確率密度関数の具体的なパラメータ・ファミリは、形状パラメータζ≧０、（κ∈Ｒ）およびスケール・パラメータθ≧０（θ∈Ｒ）を単位としてパラメータ化された、ガンマ確率密度関数であり、以下の一般形を有する。
【数３】
2013502575.tif 000005
上式で、ｔ≧０であり、Γ（κ）はガンマ関数であり、以下のように表される。
【数４】
2013502575.tif 000006
したがって、モデル処理は、クロマトグラム２００に現れる他のピーク（すなわち、ピーク２０４、２０６、２１０、２１２、および２１４）をモデル化するために、ガンマ確率密度関数を利用する。時間軸に沿って各ピーク（たとえば、その最高ピーク高さ）の状態位置を、メモリデバイス１１０内に記憶されたセットＤ’およびＨ’における参照クロマトグラフピークの位置に対応するデータと比較することによって、処理装置１０６は、クロマトグラム２００における各ピークとそれぞれの参照クロマトグラフピークとの一致の可能性を評価する。クロマトグラム２００におけるピークは、参照クロマトグラフピークと実質的に一致し、このようにして、それらのタイプに従って分類される。したがって、各クロマトグラフピークは、分離されたピーク、未知のピーク、またはデータベース内に記憶されたセットＤ’およびＨ’のいずれかにおいて対応する参照クロマトグラフピークと実質的に一致するピークのいずれかとして分類される。たとえば、処理装置１０６は、ピーク２０４および２０８が、セットＤ’におけるそれぞれの参照クロマトグラフピークｄ１およびｄ２と実質的に一致すること、ピーク２０６が、セットＨ’における参照クロマトグラフピークｈ１と実質的に一致すること、ならびにピーク２１０および２１４が、未知であると分類されることを評価する。モデリング処理における少なくとも予備段階では、それらのクロマトグラフピークは、未知であると分類され、セットＤ’およびＨ’における参照クロマトグラフピークに実質的に応答しない。以前の同定不能なクロマトグラフピークが一旦同定されると、それに応じて分類される可能性がある。開示された技法を解明するために、ピーク２１０は合成されている（すなわち、ある定度重複する少なくとも２つの成分からなる）と想定される。処理装置１０６は、演繹的知識なしに、最初にピーク２１０を未知のピークとして分類し、それに応じて、確率密度関数Ｏｌ（ｔ）によってモデル化される。分離されたピークとして分類されたクロマトグラフピークはまた、セットＤ’またはＨ’における参照クロマトグラフピークにも応答する可能性があることに留意されたい。この場合、これらの分離されたピークは、上述された分離されたピークに対して時間依存確率密度関数Ｉｍ（ｔ）に従ってモデル化される。たとえば、ピーク２１２は、分離されたピークとして分類されて、モデル化されるが、このピークは、セットＨ’における参照クロマトグラフピークに帰属可能である。したがって、それぞれの分類されたクロマトグラフピークは、そのそれぞれの確率密度関数（すなわち、Ｄｊ（ｔ）、Ｈｋ（ｔ）、Ｏｌ（ｔ）およびＩｍ（ｔ））に従ってモデル化される。
【００２３】
処理装置１０６は、クロマトグラフピークのタイプに従って（たとえば、各クロマトグラフピークの一時的帰属に従って）、クロマトグラフピークの分類を促進するために、登録手順を利用し得る。具体的には、処理装置１０６は、データベースに記憶された参照クロマトグラフピークを有する、検出された電気信号ｓ（ｔ）のクロマトグラフデータ内のクロマトグラフピークを、クロマトグラフピークの保持時間値を参照クロマトグラフピークの対応する参照保持時間値と比較することによって登録する。処理装置１０６は、各クロマトグラフピークの時間領域（時間軸に沿って）における状態（手段）位置を、メモリデバイス１１０内に記憶された参照クロマトグラフピークの位置に対応するデータと比較し得る。登録は、ｓ（ｆ（ｔ））がデータベース・エントリｒ（ｔ）に一致するように、単調変換関数ｆ（ｔ）を利用するものである。好ましくは、変換関数は線形（すなわち、ｆ（ｔ）＝ａ・ｔ＋ｂであり、ａおよびｂはパラメータである）であるが、変換関数はまた、非線形であってもよい。変換関数は、一致得点（すなわち、対応するｒ（ｔ）’ｓと一致するｓ（ｆ（ｔ））から得られる）がａおよびｂの既定の範囲内で最大であるように選択される。これは、網羅的な検索技法を利用することによって、または好ましくはガウス・ニュートン法などの最適化手法を使用することによって達成され得る。あるいは、変換関数は、反復的に現れる（たとえば、２メチルウンデカンのクロマトグラフピーク）クロマトグラフピークを考慮に入れた方法で選択される。さらに別法として、登録は、その保持時間がそれぞれの公知の保持時間を有する、公知のクロマトグラフピークを発生するように公知の、特有の化学物質を（インレット１１２を介して）挿入するものである。変換関数は、登録を促進するために、これらの公知のクロマトグラフピークを占めるように構築される。
【００２４】
時間領域内に対応する参照クロマトグラフピークに登録されたクロマトグラフピークは、それらのタイプ（たとえば、分離されたクロマトグラフピーク、参照クロマトグラフピークに実質的に一致するクロマトグラフピーク、未知のクロマトグラフピーク）に従って分類される。それぞれの分類されたクロマトグラフピークをモデル化するガンマ確率密度関数は、時間軸（たとえば、平均値μ＝ζθ）、ζ、およびθに関してピークの位置によって特徴付けられる。処理装置１０６は、最初クロマトグラフピークをモデル化するために使用される各確率密度関数に対してこれらのパラメータを推定評価する。たとえば、セットＤ’における参照クロマトグラフピークに実質的に対応するクロマトグラフピークとして分類されたクロマトグラフピークは、確率密度関数Ｄｊ（ｔ；ζｊ，θｊ）によってモデル化される。最初の推定評価を最適化するために、処理装置１０６は、クロマトグラム２００におけるクロマトグラフピークをモデル化する、各確率密度関数（すなわち、評価関数）におけるパラメータの改善された解を検索するために、最急降下法（すなわち、勾配降下）などの最適化技法を利用する。ピーク位置周辺の重み付け平均を利用することは、確率密度関数が、最初の推定評価において少なくともその近傍で、十分に平滑であること、ならびに確率密度関数に対して方向導関数が存在することを実質的に確実にする。各確率密度関数に対してパラメータ・ベクトルｐを実数値パラメータｐ＝（μ、ζ、θ）の設定数のカラムベクトルと画定することによって、新しい解が次の反復規則に従って生成される。
【数５】
2013502575.tif 000007
上式で「ｐｄｆ」は確率密度関数ｒ≧１を意味し、∇ｐｄｆ（ｐｒ）はｐｒにおける特定の密度関数の勾配であり、ｓｒは選択されたステップサイズ・パラメータである。この方法に従って、パラメータ・ベクトルｐは、候補解の評価をおそらく低減させるはずであるわずかな方向によって、各確率密度関数におけるモーメント・パラメータに調節される（すなわち、摂動される）。概して、各反復は、モデルエラーを低減するので、勾配降下法によって生成された反復解は、実質的な最適値ｐ０＝（μ０，ζ０，θ０）を網羅する。勾配降下法によって生成された反復解がローカルミニマムに捕捉される場合は、開示された技法は、シミュレーテッド・アニーリング法などを使用し得ることに留意されたい。あるいは、パラメータ・ベクトルｐは、ｐ＝（μ，ｖａｒ，γ，κ）であるように、ガンマ分布関数（すなわち、またはその問題のための他の分布関数）の最初の４個のモーメントのカラムベクトルと画定されることがあり、その平均、分散、歪度、および尖度（特に、過剰尖度）は、それぞれμ＝ζθ、ｖａｒ＝ζθ２、γ＝２／√ζ、およびκ＝６／ζによって与えられる。通常、モーメントの１つ（たとえば、尖度）は、最初の推定評価値に固定される一方で、勾配降下最適化手法は、評価関数におけるその他のモーメントに対して候補解を見出すことに進む。勾配降下最適化手法から獲得された結果ｐ０＝（μ０，ｖａｒ０，γ０）の良好な定性的測定は、尖度に対する計算値を最適化手法から獲得された値から推定された尖度の値と比較することによって、実質的に検証される可能性がある。あるいは、開示された技法は、ニュートン法、準ニュートン法、ガウス・ニュートン法、レーベンバーグ・マルカート・アルゴリズム（ＬＭＡ）などの他の最適化法を利用してもよい。たとえば、ニュートン法では、ローカルミニマムへの収束は、勾配降下のローカルミニマムへの収束より大幅に速いが、確率分布関数のヘッシアン行列の逆数を計算する必要があり、これは問題になることが多い（たとえば、曖昧である）。
【００２５】
確率密度関数の候補パラメータは、勾配降下最適化手法から獲得され、モデル関数を特徴付けるために利用される。最小二乗法は、モデル関数を実験データ（電気信号の実験データｓ（ｔ））に適合させるために利用される。具体的には、時間依存モデル関数と、検出された電気信号の実験データｓ（ｔ）におけるそれぞれの点の任意の整数（たとえば、ｎ＞０）との差の二乗和Ｓが以下のように最小化される。
【数６】
2013502575.tif 000008
【００２６】
処理装置１０６は、未知数より多い方程式があり得るので、最小二乗法によって、線形係数パラメータ（すなわち、スカラー重み）βｊ、ηｋ、δｌおよびιｍをｎ個の方程式から決定する。モデル関数の最初の評価は、線形係数パラメータが実質的に公知であると画定される。時間依存モデル関数ｘ（ｔ）の最初の評価のグラフが、図２Ｂに示されている。モデル関数の恐らく改善された評価を獲得するために、勾配降下法が、確率密度関数のパラメータ（たとえば、μ、ζ、θ）の値を最適化するために、方程式（５）に従って再度適用され、小摂動はこれらのパラメータに導入される。各確率密度関数に対して以前計算されたパラメータ値ｐ０＝（μ０，ζ０，θ０）は、推奨されたローカルミニマムに対するそれぞれの候補推測として使用される。
【００２７】
モデルエラーに関する定量的評価は、観察されたデータ（すなわち、電気信号）とモデル関数との差を取ることによって（処理装置１０８を介して）計算される。具体的には以下の通りである。
Δ＝ｘ（ｔ）−ｓ（ｔ）  （７）
あるいは、モデルエラーは、時間依存モデルエラー関数Δ（ｔ）＝ｘ（ｔ）−ｓ（ｔ）として画定される可能性がある。（世界的規模の）モデルエラー閾値パラメータがεに画定されるのは、Δ＞εである場合は、モデル関数が観察されたデータに十分に適合しないと言われているためである。概して、モデルエラー閾値パラメータは、時間依存関数ε（ｔ）であり得、その結果、不等式Δ（ｔ）＞ε（ｔ）を満たすあらゆる時間価値に対して、モデル関数が観察されたデータにその時間価値において十分に適合しないと言われている。この場合、モデルエラーΔは、分解されなかった重複ピーク（たとえば、ピーク２１０）の状況などの、分解されなかった成分（たとえば、クロマトグラフピーク、ノイズ）に起因することが仮定される。示されたモデルエラーと分解されなかったクロマトグラフピークとの関係をさらに詳説するために、次にさらに図２Ｃを参照する。具体的には、図２Ｃは、時間依存モデルエラー閾値関数のグラフに関連して描かれた、図２Ｂの最初に評価されたモデル関数に起因する、計算された時間依存モデルエラーのグラフの概略図である。図２Ｃは、最大のモデルエラーはｔ２とｔ４との間、特にｔ３で発生し、これはピーク２１０の一時的近傍に対応することを示す。その近傍におけるモデルエラーが時間依存モデルエラー閾値パラメータに対する値を超えることを考えれば、したがってピーク２１０は合成物であると疑われる。したがって、このモデルエラーは、最初に評価されたモデル関数において未確認であり不明であった、分解されなかった、または隠蔽されていたクロマトグラフピークによって引き起こされる可能性がある。ピーク２１０の一時的近傍の分析は、モデルエラーがｔ１およびｔ６において実質的にごくわずかであり、ピーク２１０に対するモデル関数の最大値がｔ４で発生することを示す。疑わしい合成ピーク内に隠蔽されたピークの数を評価するために、処理装置１０６は、たとえば、その二次導関数（たとえば、変曲点）に含まれていた情報などの、時間依存モデルエラー（関数）の曲率を分析してもよい。最初に評価されたモデル関数における単一ピーク（たとえば、確率密度関数
【数７】
2013502575.tif 000009
によって）として実際にモデル化されたピーク２１０は、次に合成物（すなわち、複数のピークを含む）であると疑われ、複数ｑの確率密度関数（たとえば、
【数８】
2013502575.tif 000010
を使用して剰余モデルエラーを考慮することによって再モデル化される。精製された時間依存モデル関数ｘ２（ｔ）は、合成物であると疑われるピーク２１０（すなわち、または概ねその問題に対する他のピーク）に対して以下のような再モデル化された式を含むことによって画定される。
【数９】
2013502575.tif 000011
【００２８】
次に、精製された時間依存モデル関数が、現行のモデル関数として取られ、モデル処理は、精製されたモデル関数ｘ３、ｘ４、ｘ５…と方程式（７）におけるモデルエラーが最小化されるまで引き続いて取られることによって繰り返される。ピーク２１０が合成物である仮説の試験は、モデル処理における各反復に精製された時間依存モデル関数を引き続き使用することによって、モデルエラーが徐々に低減し、最小値に収束するかどうかを示すことによって、実質的に支持され得る。実際にモデルエラーが、具体的な数（たとえば、２）の確率密度関数をモデルピーク２１０に利用することによって、最小値に低減された場合は、ピーク２１０は合成物であり、また、ピーク２１０はその具体的な数の重複ピークから構成されていることを示すのにある程度役に立つ。それからピーク２１０が構成されたと画定される各ピークは、それぞれの確率密度関数によってモデル化される。例示を目的として、次に、図２Ａのクロマトグラムに従ってモデル化された、図２Ｂの時間依存モデル関数の精密評価の概略図である、図２Ｄをさらに参照する。示された例では、ピーク２１０（図２Ｂ）は、２つの明確なピーク２１６および２１８（図２Ｄ）に分解され、それらの最大値はそれぞれｔ２およびｔ５（図２Ｂおよび図２Ｃ）で発生する。ｔ２およびｔ５は、モデル処理の開始時では同定されなかった。この時点で、これらの分解されたピークが、後続のモデル関数においてデータベース（すなわち、セットＤ’およびＨ’のいずれかにおいて）と比較するとき参照ピークに実質的に一致する場合は、これらのピークは、それらのそれぞれの決定された分類に従って再分類され、再モデル化される。ガンマ確率分布関数に対してカルバック・ライブラー発散（すなわち、情報発散）などの統計距離測定（すなわち、統計発散）は、データベースに記憶された参照ピークと、それぞれの参照ピークに対応すると疑われる新しく同定された分解されたピークとの一致の程度あるいは差の程度を決定するための試験として利用されてもよい。以下の方程式（９）により与えられる。
【数１０】
2013502575.tif 000012
上式で、Γ（ρＲ，σＲ）は、（すなわち、データベースに記憶された特定の参照クロマトグラフピークの）参照（Ｒ）クロマトグラフデータに関連したガンマ確率密度関数であり、Γ（ρ，σ）は、試験される（たとえば、新しく分解されたクロマトグラフピークに対応する）ガンマ確率密度関数であり、Ψ（ρＲ）はディガンマ関数である。パラメータρは形状パラメータζであり、σは速度パラメータであり（すなわち、逆目盛りパラメータσ＝１／θとして画定される）、下付き文字「Ｒ」は参照データのパラメータを意味する。カルバック・ライブラー発散によって戻された最小値は、特定の対の確率分布関数すなわち、データベースに記憶された参照と、参照と実質的に一致する疑いで試験されるものに対する最良に達成した一致を示す。あるいは、カルバック・ライブラー発散は、他の対の参照と観察されたクロマトグラフピークとの差の程度を試験するために利用されてもよい。したがって、カルバック・ライブラー発散は、データベース内の多型マーカ（複数のマーカ）と所与の例（たとえば、多型比較試験など）の複数のそれぞれのピークとの差の程度を試験するために利用されてもよい。概して、多型マーカのライブラリを考えると、最大情報発散を有するマーカは最も検出される可能性がある。さらに別法として、（すなわち、ピークの）分布間の交点を評価するための他の統計距離測定は、カルバック・ライブラー発散基準の代わりに利用されてもよい。
【００２９】
一旦モデルエラーが最小化されると、モデル処理が終了し、精製されたモデル関数が再現性の実質的に妥当なレベルで実質的に決定される。精製されたモデル関数におけるそれぞれの決定された係数βｊ、ηｋ、δｌおよびιｍは、そのそれぞれの確率密度関数に対する重み付けされた項を表し、次いで確率密度関数はそれぞれのクロマトグラフピークをモデル化する。換言すると、各係数は、試料内の特定の化学物質に対する検出された濃度の相対値を表す。通常、所与の試料内の成分の不均衡な濃度の存在を占めるために、方程式（８）の係数は、四分位範囲（ＩＱＲ）などの統計分散の程度を評価することによって正規化される。ＩＱＲは、第３四分位値と第１四分位値の差と定義され（Ｑ３−Ｑ１）、それぞれの検出されたピークを正規化するために計算され使用される（すなわち、（そのそれぞれの検出された最大濃度に対応する）各ピークの最大値をＩＱＲで割る）。
【００３０】
それにもかかわらず、その検知された濃度のある種の化合物は、統計的に有意でない可能性があるように、既定値より下回ることがある。たとえば、ある種のバイオマーカを画定する特定の化学物質の検出された低い濃度は、このバイオマーカに起因する特定の疾患が無いことを示す可能性がある。したがって、方程式（８）における各係数に対して、最小値を設定するそれぞれの閾値パラメータ（図示せず）が画定される。これは、その閾値パラメータを越した場合、その係数に対応する確率密度関数が有意であると見なされるからである。したがって、分解されたピークの１つが、たとえば重複ピーク現象に起因して事前に検知されたなった特定のバイオマーカの同定を必要とする化合物に対応する場合は、ここで検出され得る。システム１００は、特定の試料が分析不可能（たとえば、モデル化の失敗）である場合に、表示（図示せず）を発生することができることに留意されたい。
【００３１】
次に図３Ａおよび図３Ｂを参照する。図３Ａは、その異なる成分が、開示された技法の別の実施形態に従って構築され作動可能な、所与の試料（概ね３００を参照）を構成する、重複クロマトグラフピーク内で成分を分解し同定する方法を示す概略構成図である。図３Ｂは、図３Ａからの方法の継続を示す概略構成図である。手順３０２では、それぞれの参照クロマトグラフデータのデータベースを構築するように、複数の化合物からクロマトグラフデータが獲得される。図１に関して、システム１００は、メモリ１１０内に記憶されるそれぞれの参照クロマトグラフデータのデータベースを構築するように、検出器１０６を介してクロマトグラフデータを複数の化合物（図示せず）から獲得する。
【００３２】
手順３０４では、分析される試料のクロマトグラフデータが獲得される。この場合、クロマトグラフデータは、複数のピークを有するクロマトグラムとして表される。図１および図２Ａに関して、システム１００（図１）は、検出器１０６を介して分析される試料のクロマトグラフデータを獲得する。試料の獲得されたクロマトグラフデータは、複数のクロマトグラフピーク２０２、２０４、２０６、２０８、２１０、２１２および２１４を有するクロマトグラム２００（図２Ａ）として表される。
【００３３】
手順３０６では、クロマトグラフデータ内の複数のピークは、データベース内に記憶された参照クロマトグラフデータ内の参照クロマトグラフピークに、各クロマトグラフピークの保持時間値を参照クロマトグラフピークの対応する参照保持時間値と比較することによって登録される。
【００３４】
手順３０８では、獲得されたクロマトグラフデータの各ピークは、対応する参照クロマトグラフデータと比較することによって、少なくともその一時的帰属に従って分類される。
【００３５】
手順３１０では、モデル関数は、各ピークが決定された分類に従ってそれぞれの確率密度関数によってモデル化されるように、確率密度関数の線形結合の和から構築される。ただし、各確率密度関数は少なくとも１つのパラメータによって特徴付けられる。方程式（２）に関して、モデル関数ｘ（ｔ）は、複数の確率密度関数Ｄｊ（ｔ）、Ｈｋ（ｔ）、Ｏｌ（ｔ）およびＩｍ（ｔ）を用いてモデル化される。
【００３６】
手順３１２では、各確率密度関数のパラメータは、勾配降下最適化手法によって評価される。方程式（５）に関して、各確率密度関数の実数値パラメータｐ＝（μ、ζ、θ）の既定数のカラムベクトルｐが評価される。
【００３７】
手順３１４では、確率密度関数の線形結合における線形係数パラメータは、モデル関数と対応するクロマトグラフデータとの差の二乗和Ｓを最小化するように決定される。方程式（６）に関して、線形係数パラメータβｊ、ηｋ、δｌおよびιｍは、方程式（６）に画定された和Ｓを最小化するように決定される。各確率密度関数のパラメータは、手順３１２において勾配降下最適化手法によって再度評価される。手順３１２および３１４は、和Ｓが最小化されるまで旋回される（すなわち、数回にわたって反復されることがある）。
【００３８】
手順３１６では、時間依存モデルエラーは、モデル関数からクロマトグラフデータを引くことによって計算される。図２Ｃおよび方程式（７）に関して、モデルエラーは、観察されたデータ（すなわち、電気信号）とモデル関数との差を取ることによって計算される。
【００３９】
手順３１８では、時間依存モデルエラー閾値パラメータが画定される。このパラメータは、時間依存関数として画定されてもよい。図２Ｃに関して、時間依存モデルエラー閾値パラメータεが描かれる。
【００４０】
手順３２０では、合成物であると疑われるピークは、それに対して時間依存モデルエラーが時間依存モデルエラー閾値パラメータを超える時間値を評価することによって決定される。図２Ａおよび図２Ｃに関して、ピーク２１０に一時的に対応する時間依存モデルエラーは、時間値ｔ２とｔ５との間のモデルエラー閾値パラメータを実質的に超える。
【００４１】
手順３２２では、精製されたモデル関数は、複数の確率密度関数によって合成物であると疑われるピークを再モデル化することによって構築され、それぞれのピークの対応するモデルエラーを考慮に入れ、それによって合成ピークを分解する。引き続き精製されたモデル関数は、手順３１６におけるモデルエラーが最小化されるまで、反復して手順３１０におけるモデル関数と置換される。図２Ａおよび方程式（８）に関して、ピーク２１０は、精製された時間依存モデル関数が画定されるように、複数の確率密度関数によって合成されたと疑われ、再モデル化される。精製された時間依存モデル関数は、方程式（２）における現行のモデル関数として取られ、モデル処理は、方程式（７）におけるモデルエラーが最小化されるまで、精製されたモデル関数を引き続いて取ることによって反復的に繰り返される（すなわち、ステップ３１０から）。
【００４２】
手順３２４では、ピークに関連する線形係数パラメータは、各ピークのそれぞれの最大ピーク値をＩＱＲで割ることによって正規化される。方程式（８）に関して、線形係数パラメータβｊ、ηｋ、δｌおよびιｍは、計算されたＩＱＲによって正規化される。
【００４３】
手順３２６では、有意なピークは、それぞれの確率密度関数の正規化された線形係数パラメータが、それぞれの閾値パラメータを超えるかどうかを評価することによって判定される。方程式（８）に関して、有意なピーク（図示せず）は、線形係数パラメータβｊ、ηｋ、δｌおよびιｍが、それぞれの閾値パラメータ（図示せず）を超えるかどうかを評価することによって判定される。
【００４４】
手順３２８では、参照ピークと分解されたピークを含む複数のピークとの一致の程度が試験される。図１および図２Ｄならびに方程式（９）に関して、分解されたピーク２１６および２１８は、それらとメモリ１１０（図１）のデータベース内に記憶されたクロマトグラフ参照ピークとの一致の程度（または差の程度）を試験するために、カルバック・ライブラー発散で試験される。
【００４５】
開示された技法は、上記に具体的に示され説明されたことに限定されないことは当業者には理解されよう。そうではなく、開示された技法の範囲は、以下の特許請求の範囲によってのみ画定される。
【図１】
2013502575.tif 000017
【図２Ａ】
2013502575.tif 000018
【図２Ｂ】
2013502575.tif 000019
【図２Ｃ】
2013502575.tif 000020
【図２Ｄ】
2013502575.tif 000021
【図３Ａ】
2013502575.tif 000022
【図３Ｂ】
2013502575.tif 000023
【国際調査報告】
2013502575.tif 000024
2013502575.tif 000025
2013502575.tif 000026