【課題】多成分物質の評価の精度及び効率を向上させることを可能とする。
【解決手段】対象ＦＰ作成工程１７３と、対象ＦＰピーク帰属工程１４９と、対象ＦＰピ
ーク特徴量作成工程１５１と、対象ＦＰタイプ２作成工程１５３と、対象ＦＰ領域分割特
徴量作成工程１５５と、対象ＦＰ特徴量統合工程１５７と、評価工程１７１とを備え、対
象ＦＰピーク特徴量及び対象ＦＰ領域分割特徴量を統合して対象ＦＰ統合特徴量を作成し
、対象ＦＰ統合特徴量とこの対象ＦＰ統合特徴量に対応し評価基準となる多成分物質の複
数の基準ＦＰに基づく基準ＦＰ統合特徴量とを比較評価することを特徴とする。

【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、パターンの特徴量作成方法、多成分物質、例えば多成分薬剤である漢方薬の
品質評価を行うための多成分物質のＦＰの特徴量作成方法、作成プログラム、及び作成装
置に関する。
【背景技術】
【０００２】
多成分物質として、例えば多成分で構成される薬剤（以下、多成分薬剤）である漢方薬
などの天然物由来の薬剤がある。これら薬剤の定量的、定性的プロフィールは、使用する
原料生薬の、地質学的要因、生態学的要因、収集時期、収集場所、収集年代、生育期の天
候等が原因で変化する。
【０００３】
そのため、これら多成分薬剤などについては、その安全性および有効性を担保するため
の品質として一定の基準を規定し、その基準に基づいて、国家の監督機関、化学的組織、
製造業者等が品質評価を行っている。
【０００４】
しかしながら、多成分薬剤の品質等の判定基準は、多成分薬剤中のある特徴的な１ない
しは数成分を選択して、その含有量などに基づいて設定することが一般的であった。
【０００５】
例えば、非特許文献１では、多成分薬剤において有効成分の同定ができていない場合、
定量分析が可能、水に溶けやすい、熱水中で分解しない、他の成分と化学反応をしない等
の物性を持った複数の成分を選択し、化学分析により得られるそれら成分の含有量を評価
の基準としている。
【０００６】
また、多成分薬剤にクロマトグラフィーを適用し、保持時間ごとに紫外可視吸収スペク
トルを得て、その中の一部の成分情報から評価の基準を設定することも知られている。
【０００７】
例えば、特許文献１では、ＨＰＬＣクロマトグラム・データ中の一部のピークを選択し
、バーコード化することによって多成分薬剤を評価している。
【０００８】
しかしながら、これらの方法は、評価の対象が「特定成分の含有量」又は「特定成分の
クロマト・ピーク」に限定されており、多成分薬剤が含有する成分の一部が評価対象にな
っているに過ぎない。このため、多成分薬剤については、評価の対象外となっている成分
が多数存在することから、多成分薬剤の評価方法として精度的に不十分である。
【０００９】
多成分薬剤の品質を正確に評価するには、全ピーク情報或いは数％の細かい情報を除外した全ピーク情報に近いピーク情報を網羅した評価が必要であり、そのため、多成分薬剤間で全ピーク或いはこれに近い個数のピークを対応させる必要がある。
【００１０】
ところが、複数のピークを高精度で効率的に対応させることが困難であり、このことが
多成分薬剤の高精度で効率的な評価の妨げとなっていた。
【００１１】
さらに説明すると、同じ製品名の多成分薬剤であっても生薬は天然物であるがために構
成成分が僅かに異なる場合がある。このため、同じ品質の薬剤であっても構成成分の含量
比率が異なったり、ある薬剤に存在する成分が他の薬剤には存在しないこと（以下、薬剤
間誤差）がある。さらに、クロマトにおけるピーク強度やピークの溶出時間には厳密な再
現性がない（以下、分析誤差）などの要因もある。これらにより、多成分薬剤間の全ピー
ク或いはこれに近い個数のピークについて同じ成分に由来するピークに対応させること（以下、ピーク帰属）ができない
ことが原因となって高精度で効率的な評価の妨げとなっていた。
【先行技術文献】
【特許文献】
【００１２】
【特許文献１】特開２００２−２１４２１５号公報
【非特許文献】
【００１３】
【非特許文献１】月刊薬事 ｖｏｌ．２８， Ｎｏ．３， ６７-７１ (１９８６)
【発明の概要】
【発明が解決しようとする課題】
【００１４】
解決しようとする問題点は、既存の評価方法では多成分物質の品質等を高精度で効率的
に評価することに限界があった点である。
【課題を解決するための手段】
【００１５】
本発明は、評価の精度及び効率を向上させることに寄与可能とするために、時系列でピ
ークが変化するパターンを複数の領域に分割し各領域に存在するピークの存在率又は存在
量からパターン領域分割特徴量を作成するパターン領域分割特徴量作成工程を備えたこと
をパターンの特徴量作成方法の特徴とする。
【００１６】
本発明は、多成分物質のクロマトから検出されたピークとそのリテンション・タイムと
で構成されるＦＰを複数の領域に分割し各領域に存在するピークの存在率又は存在量から
ＦＰ領域分割特徴量を作成するＦＰ領域分割特徴量作成工程を備えたことをＦＰの特徴量
作成方法の特徴とする。
【００１７】
本発明は、時系列でピークが変化するパターンを複数の領域に分割し各領域に存在する
ピークの存在率又は存在量からパターン領域分割特徴量を作成するパターン領域分割特徴
量作成機能をコンピュータに実現させることをパターンの特徴量作成プログラムの特徴と
する。
【００１８】
本発明は、多成分物質のクロマトから検出されたピークとそのリテンション・タイムと
で構成されるＦＰを複数の領域に分割し各領域に存在するピークの存在率又は存在量から
ＦＰ領域分割特徴量を作成するＦＰ領域分割特徴量作成機能を備えたことをＦＰの特徴量
作成プログラムの特徴とする。
【００１９】
本発明は、時系列でピークが変化するパターンを複数の領域に分割し各領域に存在する
ピークの存在率又は存在量からパターン領域分割特徴量を作成するパターン領域分割特徴
量作成部を備えたことをパターンの特徴量作成装置の特徴とする。
本発明は、多成分物質のクロマトから検出されたピークとそのリテンション・タイムとで
構成されるＦＰを複数の領域に分割し各領域に存在するピークの存在率又は存在量からＦ
Ｐ領域分割特徴量を作成するＦＰ領域分割特徴量作成部を備えたことをＦＰの特徴量作成
装置の特徴とする。
【発明の効果】
【００２０】
本発明のパターン又はＦＰの特徴量作成方法は、上記構成であるから、パターン又はＦ
Ｐの特徴量を領域分割により簡単に得ることが可能となる。したがって、例えば細かいピ
ークをも捉えて特徴量を作成することができる。
【００２１】
本発明のパターン又はＦＰの特徴量作成プログラムは、上記構成であるから、各機能を
コンピュータに実現させ、パターン又はＦＰの特徴量を簡単に得ることが可能となる。本
発明のパターン又はＦＰの特徴量作成装置は、上記構成であるから、各部を作用させ、パ
ターン又はＦＰの特徴量を簡単に得ることが可能となる。
【図面の簡単な説明】
【００２２】
【図１】多成分薬剤の評価装置のブロック図である。（実施例１）
【図２】多成分薬剤の評価手順を示すブロック図である。（実施例１）
【図３】三次元クロマトグラム・データ（以下、３Ｄクロマト）から作成したＦＰの説明図である。（実施例１）
【図４】（Ａ）は、薬剤Ａ、（Ｂ）は、薬剤Ｂ、（Ｃ）は、薬剤ＣのＦＰである。（実施例１）
【図５】対象ＦＰ及び基準ＦＰのリテンション・タイムを示す図である。（実施例１）
【図６】対象ＦＰのリテンション・タイム・出現パターンを示す図である。（実施例１）
【図７】基準ＦＰのリテンション・タイム・出現パターンを示す図である。（実施例１）
【図８】対象ＦＰと基準ＦＰのリテンション・タイム・出現距離の一致数を示す図である。（実施例１）
【図９】対象ＦＰと基準ＦＰのリテンション・タイム・出現パターンの一致度を示す図である。（実施例１）
【図１０】対象ＦＰの帰属対象ピークを示す図である。（実施例１）
【図１１】帰属対象ピークを含めたピーク３本によるピーク・パターン図である。（実施例１）
【図１２】帰属対象ピークを含めたピーク５本によるピーク・パターン図である。（実施例１）
【図１３】帰属対象ピークの許容幅を示す図である。（実施例１）
【図１４】帰属対象ピークに対する基準ＦＰの帰属候補ピークを示す図である。（実施例１）
【図１５】帰属対象ピークと帰属候補ピークのピーク３本によるピーク・パターン図である。（実施例１）
【図１６】帰属対象ピークと別な帰属候補ピークのピーク３本によるピーク・パターン図である。（実施例１）
【図１７】帰属対象ピークと別な帰属候補ピークのピーク３本によるピーク・パターン図である。（実施例１））
【図１８】帰属対象ピークと別な帰属候補ピークのピーク３本によるピーク・パターン図である。（実施例１）
【図１９】帰属対象ピークと帰属候補ピークのピーク５本によるピーク・パターン図である。（実施例１）
【図２０】帰属対象ピークと別な帰属候補ピークのピーク５本によるピーク・パターン図である。（実施例１）
【図２１】帰属対象ピークと別な帰属候補ピークのピーク５本によるピーク・パターン図である。（実施例１））
【図２２】帰属対象ピークと別な帰属候補ピークのピーク５本によるピーク・パターン図である。（実施例１）
【図２３】帰属対象ピークと帰属候補ピークのピーク・パターン構成候補ピークを示す図である。（実施例１）
【図２４】ピーク・パターン構成候補ピークを４本とした場合の帰属対象ピークの全ピーク・パターン数を示す図である。（実施例１）
【図２５】ピーク・パターン構成候補ピークを４本とした場合の帰属候補ピークの全ピーク・パターン数を示す図である。（実施例１）
【図２６】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図２７】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図２８】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図２９】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図３０】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図３１】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図３２】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図３３】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図３４】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図３５】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図３６】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図３７】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図３８】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図３９】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図４０】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図４１】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図４２】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図４３】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図４４】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図４５】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図４６】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図４７】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図４８】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図４９】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図５０】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図５１】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図５２】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図５３】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図５４】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図５５】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図５６】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図５７】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図５８】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図５９】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図６０】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図６１】帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）
【図６２】帰属対象ピークと帰属候補ピークのピーク３本によるピーク・パターンの一致度の算出方法を示す図である。（実施例１）
【図６３】帰属対象ピークと帰属候補ピークのピーク３本によるピーク・パターンの一致度の算出方法を示す図である。（実施例１）
【図６４】帰属対象ピークと帰属候補ピークのピーク５本によるピーク・パターンの一致度の算出方法を示す図である。（実施例１）
【図６５】帰属対象ピーク及び帰属候補ピークのＵＶスペクトルを示す図である。（実施例１）
【図６６】帰属対象ピークと帰属候補ピークのＵＶスペクトルの一致度の説明図である。（実施例１）
【図６７】ピーク・パターンとＵＶスペクトルの両方の比較による帰属候補ピークの一致度計算の説明図である。（実施例１）
【図６８】対象ＦＰの基準群ＦＰへの帰属を示す説明図である。（実施例１）
【図６９】対象ＦＰが基準群ＦＰへ帰属された状況を示す図である。（実施例１）
【図７０】領域分割による数量化を示す説明図である。（実施例１）
【図７１】リテンション・タイム等の変動との関係を示す説明図である。（実施例１）
【図７２】領域の位置を変更して数量化する説明図である。（実施例１）
【図７３】ＦＰタイプ２のデータを示す図表である。（実施例１）
【図７４】ＦＰタイプ２のパターンを示す説明図である。（実施例１）
【図７５】縦・横分割線での領域分割による領域ごとの特徴量化を示す説明図である。（実施例１）
【図７６】縦分割線（1本目）の設定を示す説明図である。（実施例１）
【図７７】横分割線（1本目）の設定を示す説明図である。（実施例１）
【図７８】縦・横分割線による領域分割を示す説明図である。（実施例１）
【図７９】特徴量化する領域の数を示す説明図である。（実施例１）
【図８０】領域１の特定を示す説明図である。（実施例１）
【図８１】全ピークの高さ及び合計を示す図表である。（実施例１）
【図８２】領域１のピーク高さの合計を示す説明図である。（実施例１）
【図８３】全領域の特徴量を示す図表である。（実施例１）
【図８４】縦１本目の位置を順次変更してできた各領域での特徴量を示す図表である。（実施例１）
【図８５】横１本目の位置を順次変更してできた各領域での特徴量を示す図表である。（実施例１）
【図８６】各縦・横分割線の位置を変更しない１通りでの特徴量を示す図表である。（実施例１）
【図８７】各種対象ＦＰとその評価値（ＭＤ値）を示す図である。（実施例１）
【図８８】各種対象ＦＰとその評価値（ＭＤ値）を示す図である。（実施例１）
【図８９】各種対象ＦＰとその評価値（ＭＤ値）を示す図である。（実施例１）
【図９０】各種対象ＦＰとその評価値（ＭＤ値）を示す図である。（実施例１）
【図９１】各種対象ＦＰとその評価値（ＭＤ値）を示す図である。（実施例１）
【図９２】多成分薬剤の評価方法を示す工程図である。（実施例１）
【図９３】多成分薬剤の品質評価フロー図である。（実施例１）
【図９４】多成分薬剤の品質評価フロー図である。（実施例１）
【図９５】単波長によるＦＰ作成機能におけるデータ処理フローチャートである。（実施例１）
【図９６】複数波長によるＦＰ作成機能におけるデータ処理フローチャートである。（実施例１）
【図９７】複数波長によるＦＰ作成機能におけるデータ処理フローチャートである。（実施例１）
【図９８】ピーク帰属処理１（基準ＦＰの選定）におけるデータ処理フローチャートである。（実施例１）
【図９９】ピーク帰属処理２（帰属スコアの算出）におけるデータ処理フローチャートである。（実施例１）
【図１００】ピーク帰属処理３（対応ピークの特定）におけるデータ処理フローチャートである。（実施例１）
【図１０１】ピーク帰属処理４（基準群ＦＰへの帰属）におけるデータ処理フローチャートである。（実施例１）
【図１０２】ピーク帰属処理４（基準群ＦＰへの帰属）におけるデータ処理フローチャートである。（実施例１）
【図１０３】ピーク帰属処理１（基準ＦＰの選定）におけるリテンション・タイム・出現パターンの一致度計算処理のフローチャートである。（実施例１）
【図１０４】ピーク帰属処理２（帰属スコアの算出）におけるＵＶスペクトルの一致度計算処理のフローチャートである。（実施例１）
【図１０５】ピーク帰属処理２（帰属スコアの算出）におけるピーク・パターンの一致度計算処理のフローチャートである。（実施例１）
【図１０６】「ＦＰ_type２作成」の詳細を示すフローチャートである。（実施例１）
【図１０７】「領域分割による対象ＦＰ_type２の特徴量化処理」の詳細を示すフローチャートである。（実施例１）
【図１０８】「対象ＦＰのピーク特徴量と領域分割特徴量の統合」の詳細を示すフローチャートである。（実施例１）
【図１０９】基準ＦＰ特徴量統合ファイルを作成するためのフローチャートである。（実施例１）
【図１１０】基準ＦＰ特徴量統合ファイルを作成するためのフローチャートである。（実施例１）
【図１１１】「基準ＦＰ帰属結果統合処理（基準ＦＰ対応表の作成）」の詳細を示すフローチャートである。（実施例１）
【図１１２】「基準ＦＰ帰属結果統合処理（基準ＦＰ対応表の作成）」の詳細を示すフローチャートである。（実施例１）
【図１１３】「ピーク特徴量化処理（基準群ＦＰの作成）」の詳細を示すフローチャートである。（実施例１）
【図１１４】「基準ＦＰ_type２の作成処理」の詳細を示すフローチャートである。（実施例１）
【図１１５】「領域分割による基準ＦＰの特徴量化処理」の詳細を示すフローチャートである。（実施例１）
【図１１６】基準ＦＰの特徴量統合処理に係るフローチャートである。（実施例１）
【図１１７】３Ｄクロマトのデータ例を示す図表である。（実施例１）
【図１１８】ピーク情報のデータ例を示す図表である。（実施例１）
【図１１９】ＦＰのデータ例を示す図表である。（実施例１）
【図１２０】対象ＦＰの基準ＦＰへの帰属スコア計算結果（判定結果ファイル例）を示す図表である。（実施例１）
【図１２１】対象ＦＰと基準ＦＰで対応するピークの照合過程を示す図表である。（実施例１）
【図１２２】対象ＦＰと基準ＦＰで対応するピークを特定した結果（照合結果ファイル）例を示す図表である。（実施例１）
【図１２３】基準群ＦＰのデータ例を示す図表である。（実施例１）
【図１２４】対象ＦＰピーク特徴量ファイル例を示す図表である。（実施例１）
【図１２５】対象及び基準ＦＰタイプ２のデータ例を示す図表である。（実施例１）
【図１２６】対象ＦＰ領域分割特徴量ファイル例を示す図表である。（実施例１）
【図１２７】対象ＦＰ統合特徴量ファイル例を示す図表である。（実施例１）
【図１２８】基準ｔｙｐｅ２群ＦＰ例を示す図表である。（実施例１）
【図１２９】基準群統合データ例を示す図表である。（実施例１）
【図１３０】図１０４に代えて適用するサブルーチン２の変形例の詳細を示すフローチャートである。（実施例１）
【図１３１】移動平均及び移動傾きの計算例を示す図表である。（実施例１）
【発明を実施するための形態】
【００２３】
評価の精度及び効率を向上させることに寄与可能にするという目的を、多成分物質のク
ロマトから検出されたピークとそのリテンション・タイムとで構成されるＦＰを複数の領
域に分割し各領域に存在するピークの存在率又は存在量からＦＰ領域分割特徴量を作成す
ることにより実現した。
【実施例１】
【００２４】
本発明の実施例１は、多成分物質、例えば多成分薬剤を評価する多成分薬剤の評価方法
、評価プログラム、パターンであるＦＰの特徴量作成方法、作成プログラム、作成装置で
ある。
【００２５】
多成分薬剤は、複数の有効化学成分を含有する薬剤と定義され、限定はされないが、生
薬、生薬の組合せ、それらの抽出物、漢方薬等が含まれる。また剤形も特に限定されず、
例えば、第１５改正日本薬局方の製剤総則で規定されている液剤、エキス剤、カプセル剤
、顆粒剤、丸剤、懸濁剤・乳剤、散剤、酒精剤、錠剤、浸剤・煎剤、チンキ剤、トローチ
剤、芳香水剤、流エキス剤等が含まれる。多成分物質としては、薬剤以外のものも含まれ
る。
【００２６】
漢方薬の具体例は、医療用漢方製剤１４８処方「使用上の注意」の業界統一と自主改訂
、一般用漢方処方の手引き(１９７８年)に記載されている。
【００２７】
多成分薬剤の評価では、評価対象薬剤が正常品と定めた複数の薬剤と同等であるかどう
かを評価するため、まず、評価対象薬剤の三次元クロマトグラム・データ（以下、３Ｄク
ロマト）からその薬剤特有の情報を抽出した対象ＦＰを作成する。
【００２８】
次に、対象ＦＰの各ピークを、全基準ＦＰをピーク帰属処理し作成した全基準ＦＰのピーク対応データ（以下、基準群ＦＰ）に帰属し、ピーク特徴量を得る。
【００２９】
さらに、対象ＦＰから帰属されたピークを除き、残ったピークでＦＰタイプ２を作成し
、このＦＰタイプ２を領域分割し領域分割特徴量を得る。
【００３０】
これら２つの特徴量を統合し、対象ＦＰ統合特徴量を得る。
【００３１】
この対象ＦＰ統合特徴量と全基準ＦＰから得た基準ＦＰ統合特徴量により、基準群ＦＰ
と対象ＦＰの同等性をＭＴ法で評価する。最後に、得られた評価値（以下、ＭＤ値）と予
め設定しておいた判定値（ＭＤ値の上限値）を比較し、評価対象薬剤が正常品と同等であ
るかどうかを判定する。
［多成分薬剤の評価装置］
図１は、多成分薬剤の評価装置のブロック図、図２は、多成分薬剤の評価手順を示すブ
ロック図、図３は、３Ｄクロマトから作成したＦＰの説明図、図４（Ａ）は、薬剤Ａ、（
Ｂ）は、薬剤Ｂ、（Ｃ）は、薬剤ＣのＦＰである。
【００３２】
図1のように、多成分薬剤の評価装置１は、ＦＰ作成部３と、対象ＦＰピーク帰属部５
と、対象ＦＰピーク特徴量作成部７と、対象ＦＰタイプ２作成部９と、対象ＦＰ領域分割
特徴量作成部１１と、対象ＦＰ特徴量統合部１３と、基準ＦＰピーク帰属部１５と、基準
ＦＰ帰属結果統合部１７と、基準ＦＰピーク特徴量作成部１９と、基準ＦＰタイプ２作成
部２１と、基準ＦＰ領域分割特徴量作成部２３と、基準ＦＰ特徴量統合部２５と、評価部
２７とを備えている。多成分薬剤の評価装置１は、パターンであるＦＰの特徴量作成装置
を含む。
【００３３】
ＦＰ作成部３は、対象ＦＰ作成部２９と、基準ＦＰ作成部３１とを備えている。
【００３４】
対象ＦＰピーク帰属部５は、基準ＦＰ選定部３３と、ピーク・パターン作成部３５と、
ピーク帰属部３７とを備えている。
【００３５】
この多成分薬剤の評価装置１は、例えば、一つのコンピュータで構成され、図示はしな
いが、ＣＰＵ、ＲＯＭ、ＲＡＭなどを備えている。多成分薬剤の評価装置１は、コンピュータにインストールされたパターンの特徴量作成プログラムとしてＦＰの特徴量作成プログラムを実現させ、ＦＰの特徴量を得ることができる。但し、ＦＰの特徴量作成プログラムは、これを記録したＦＰの特徴量作成プログラム記録媒体を用い、コンピュータで構成された多成分薬剤の評価装置１にこれを読み取らせることで、ＦＰの特徴量を得ることを実現させることもできる。
【００３６】
多成分薬剤の評価装置１は、各部をそれぞれ別々のコンピュータで構成することができ、例えば対象ＦＰピーク帰属部５と、対象ＦＰピーク特徴量作成部７と、対象ＦＰタイプ２作成部９と、対象ＦＰ領域分割特徴量作成部１１と、対象ＦＰ特徴量統合部１３と、評価部２７とを一つのコンピュータで構成し、基準ＦＰ作成部３１と、基準ＦＰピーク帰属部１５と、基準ＦＰ帰属結果統合部１７と、基準ＦＰピーク特徴量作成部１９と、基準ＦＰタイプ２作成部２１と、基準ＦＰ領域分割特徴量作成部２３と、基準ＦＰ特徴量統合部２５を他のコンピュータで構成することなどもできる。
【００３７】
この場合、基準ＦＰ統合特徴量は、他のコンピュータで作成し、評価装置１の評価部２
７に入力されることになる。
【００３８】
そして、対象ＦＰ作成部２９と、対象ＦＰピーク帰属部５と、対象ＦＰピーク特徴量作
成部７と、対象ＦＰタイプ２作成部９と、対象ＦＰ領域分割特徴量作成部１１と、対象Ｆ
Ｐ特徴量統合部１３とにより対象ＦＰ統合特徴量が作成され、基準ＦＰ作成部３１と、基
準ＦＰピーク帰属部１５と、基準ＦＰ帰属結果統合部１７と、基準ＦＰピーク特徴量作成
部１９と、基準ＦＰタイプ２作成部２１と、基準ＦＰ領域分割特徴量作成部２３と、基準
ＦＰ特徴量統合部２５とにより基準ＦＰ統合特徴量が作成され、これらが比較評価されて
対象ＦＰ４３と基準群ＦＰ４５との同等性を評価する。
【００３９】
ＦＰ作成部３の対象ＦＰ作成部２９は、例えば、図２、図３のように、漢方薬３９のク
ロマトとして三次元クロマトグラム・データである３Ｄクロマト４１から、特定の検出波
長における複数のピークとそのリテンション・タイムならびにＵＶスペクトルを抽出した
対象ＦＰ４３（以下単に、「ＦＰ４３」ということもある。）を作成する機能部である。
【００４０】
このＦＰ４３は、３Ｄクロマト４１と同様に三次元の情報（ピーク、リテンション・タ
イム及びＵＶスペクトル）で構成している。
【００４１】
このため、ＦＰ４３は、その薬剤特有の情報をそのまま継承したデータである。それに
も係らず、データ容量は約１／７０に圧縮されているため、３Ｄクロマト４１に比較して
、処理すべき情報量を大幅に減少させることができ処理速度を速めることができる。
【００４２】
３Ｄクロマト４１は、漢方薬３９に高速液体クロマトグラフィー（ＨＰＬＣ）を適用し
た結果である。この３Ｄクロマト４１は、各成分の移動速度として現れ、それを特定時間
における移動距離として表し、或いはカラム末端から時系列に現れる様をチャートに表し
たものである。ＨＰＬＣにおいては、時間軸に対する検出器応答をプロットしたもので、
ピークの出現時間を保持時間（リテンション・タイム）と呼んでいる。
【００４３】
検出器としては、特に限定されないが、光学的性質を利用した吸光度検出器(Absorbanc
e Detector) が使用され、ピークは、紫外線（ＵＶ）の検出波長に応じたシグナル強度として三次元的に得られたものである。光学的性質を利用したものとしては、透過率検出器（Transmittance Detector）を用いることもできる。
【００４４】
検出波長に限定はなく、好ましくは１５０ｎｍ〜９００ｎｍの範囲であり、特に好まし
くは２００ｎｍ〜４００ｎｍの紫外可視吸領域、更に好ましくは２００ｎｍ〜３００ｎｍ
から選ばれた複数の波長である。
【００４５】
そして、３Ｄクロマト４１は、少なくとも、漢方薬の番号（ロット番号）、リテンショ
ン・タイム、検出波長、及びピークをデータとして有するものである。
【００４６】
なお、３Ｄクロマト４１は、市販の装置によっても得ることができ、かかる市販の装置
としては、Ａｇｉｌｅｎｔ１１００システム等が挙げられる。また、クロマトグラフィー
は、ＨＰＬＣに限定されず、種々のものを採用することができる。
【００４７】
３Ｄクロマト４１は、図２、図３のようにｘ軸をリテンション・タイム、ｙ軸を検出波
長、ｚ軸をシグナル強度として表示する。
【００４８】
ＦＰ４３は、少なくとも、漢方薬の番号（ロット番号）、リテンション・タイム、特定
の波長におけるピーク、及びＵＶスペクトルをデータとして有するものである。
【００４９】
ＦＰ４３は、図２、図３のようにｘ軸をリテンション・タイム、ｙ軸を特定の検出波長
におけるピークとした二次元で表示するが、図３のように１ピークで示すＵＶスペクトル
４２と同様なＵＶスペクトル情報をピークごとに有したデータである。 ＦＰ４３を作成
する特定の検出波長は、特に限定されず種々選択できる。但し、ＦＰ４３には３Ｄクロマ
ト中の全てのピークを含めることが情報を継承するという点において重要である。このた
め、本実施例１では検出波長を３Ｄクロマト中の全てのピークを含んでいる２０３ｎｍと
した。
【００５０】
一方、単独の波長では全てのピークを含めることができないことがある。そのような場
合は検出波長を複数とし、後述のように複数の波長を組み合わせて全てのピークを含んだ
ＦＰを作成する。
【００５１】
本実施例１では、ピークをシグナル強度（ピーク高さ）の極大値としたが、ピークとし
て面積値を採用することもできる。また、ＦＰにＵＶスペクトルを含めず、ｘ軸をリテン
ション・タイム、ｙ軸を特定の検出波長におけるピークとした二次元の情報のみにするこ
とも可能である。この場合は、ＦＰを漢方薬の番号（ロット番号）、リテンション・タイ
ムをデータとして有するクロマトとしての２Ｄクロマトから作成することもできる。
【００５２】
図４の（Ａ）は、薬剤Ａ、（Ｂ）は、薬剤Ｂ、（Ｃ）は、薬剤ＣのＦＰ５５、５７、５
９である。
【００５３】
対象ＦＰピーク帰属部５は、前記対象ＦＰとこの対象ＦＰに対応し評価基準となる多成
分物質の基準ＦＰとのピークを比較し対応するピークを特定する機能部である。この対象
ＦＰピーク帰属部５は、基準ＦＰ選定部３３と、ピーク・パターン作成部３５と、ピーク
帰属部３７とからなっている。
【００５４】
基準ＦＰ選定部３３は、対象ＦＰのピーク帰属に適した多成分物質のＦＰを複数の基準
ＦＰから選定する機能部である。すなわち、対象ＦＰの各ピークのピーク帰属を高い精度
で行うため、図５〜図９（後述）のように対象ＦＰと基準ＦＰ間でピークのリテンション
・タイム・出現パターンの一致度を算出し、この一致度が最小となる基準ＦＰを全基準Ｆ
Ｐから選定する。
【００５５】
ピーク・パターン作成部３５は、図１０〜図１２（後述）のように対象ＦＰ６１におけ
る帰属の対象となるピーク（以下、帰属対象ピーク）に対し、時間軸方向前後の少なくと
も一方に存在するｎ本のピークを含めた計ｎ＋１本のピークで構成されるピーク・パター
ンを帰属対象ピークのピーク・パターンとして作成する機能部である。ｎは、自然数であ
る。
【００５６】
なお、図１１（後述）では時間軸方向前後の少なくとも一方に存在する２本のピークを
含めた計３本のピークで構成されるピーク・パターンを、図１２（後述）では時間軸方向
前後の少なくとも一方に存在する４本のピークを含めた計５本のピークで構成されるピー
ク・パターンを示す。
【００５７】
また、ピーク・パターン作成部３５は、図１３〜図２２（後述）のように基準ＦＰ８３
において帰属対象ピークのリテンション・タイムとの差が設定した範囲（許容幅）内の全
てのピーク（以下、帰属候補ピーク）に対し、時間軸方向前後の少なくとも一方に存在す
るｎ本のピークを含めた計ｎ＋１本のピークで構成されるピーク・パターンを帰属候補ピ
ークのピーク・パターンとして作成する機能部である。なお、図１５〜図１８（後述）で
は時間軸方向前後の少なくとも一方に存在する２本のピークを含めた計３本のピークで構
成されるピーク・パターンを示す。図１９〜図２２（後述）では時間軸方向前後の少なく
とも一方に存在する４本のピークを含めた計５本のピークで構成されるピーク・パターン
を示す。
【００５８】
許容幅に限定はなく、精度と効率の点から０．５分〜２分が好ましい。実施例１では、
１分とした。
【００５９】
さらに、ピーク・パターン作成部３５では、対象ＦＰ６１と基準ＦＰ８３のピーク数に
違いがある（つまり、どちらか一方に存在しないピークがある）場合に対しても柔軟に対
応できるようにする。このため、図２３〜図６１（後述）のように、帰属対象ピークなら
びに帰属候補ピークの両方でピーク・パターンを構成するピーク（以下、ピーク・パター
ン構成ピーク）を変化させて網羅的にピーク・パターンを作成する。なお、図２３〜図６
１では時間軸方向前後の少なくとも一方に存在する２本のピークを含めた計３本のピーク
で構成されるピーク・パターンの場合について示す。
【００６０】
ピーク帰属部３７は、対象ＦＰ、基準ＦＰそれぞれのピーク・パターンを比較し対応す
るピークを特定する機能部である。実施例では、帰属対象ピークのピーク・パターンと帰
属候補ピークのピーク・パターンとの一致度及びＵＶスペクトルの一致度を算出して対応
するピークを特定する。
【００６１】
また、この２つの一致度を統合した帰属候補ピークの一致度を算出し、この一致度に基
づき、対象ＦＰ６１の各ピークを基準ＦＰ８３の各ピークに帰属する機能部である。
【００６２】
ピーク帰属部３７において、ピーク・パターンの一致度は、図６２〜図６４（後述）の
ように帰属対象ピークと帰属候補ピークのピーク・パターン間の対応するピーク及びリテ
ンション・タイムの差をもとに算出する。また、ＵＶスペクトルの一致度は、図６５、図
６６（後述）のように帰属対象ピーク７３のＵＶスペクトル１３５と帰属候補ピーク９５
のＵＶスペクトル１３９の各波長における吸光度の差をもとに算出する。さらに、図６７
（後述）のようにこれら２つの一致度を乗じて帰属候補ピーク９５の一致度を算出する。
【００６３】
対象ＦＰピーク特徴量作成部７は、対象ＦＰピーク帰属部５で特定されて帰属したピー
クと複数の基準ＦＰである基準群ＦＰ４５のピークとを比較評価して特徴量化された対象
ＦＰピーク特徴量として作成する機能部である。複数の基準ＦＰは、評価基準となる多成
分物質である複数の漢方薬に対応して作成されており、この複数の漢方薬は、正常品とさ
れているものである。
【００６４】
すなわち、対象ＦＰピーク特徴量作成部７は、対象ＦＰ６１と基準ＦＰ８３の帰属結果
に基づき、最終的に図２、図６８、図６９（後述）のように対象ＦＰ４３の各ピークを基
準群ＦＰ４５の各ピークに帰属して特徴量化された対象ＦＰピーク特徴量４７を作成する
機能部である。
【００６５】
対象ＦＰタイプ２作成部９は、対象パターンから特徴量化されたピークを除いて残った
ピークで構成されるパターンを対象パターン・タイプ２として作成する。例えば、対象Ｆ
Ｐピーク特徴量作成部７において特定されたピーク４７を元の対象ＦＰ４３から除いて残
ったピークとそのリテンション・タイムとで構成されるＦＰを図２の対象ＦＰタイプ２（
４９）をパターンとして作成する機能部である。
【００６６】
この対象ＦＰタイプ２（４９）は、対象ＦＰピーク特徴量作成部７において特徴量化さ
れなかったピークを集めてＦＰとしたものである。この対象ＦＰタイプ２（４９）を特徴
量化し評価に加えることで、より正確な評価を行わせることができる。
【００６７】
対象ＦＰ領域分割特徴量作成部１１は、対象パターン・タイプ２を複数の領域に分割し
各領域に存在するピークの存在率から対象パターン領域分割特徴量を作成するＦＰ領域分
割特徴量作成部を構成し、対象ＦＰタイプ２（４９）を複数の領域に分割し各領域に存在
するピークの存在率から対象ＦＰ領域分割特徴量を対象パターン領域分割特徴量として作
成する機能部である。
【００６８】
なお、対象ＦＰ領域分割特徴量作成部１１は、存在率に代えて存在量を用いることもで
きる。存在率は、後述のように各領域のピーク高さの存在量を全体のピーク高さの合計（
つまり全体のピーク高さの存在量）で割った値である。したがって、各領域のピーク高さ
の存在量そのものを用いて領域分割特徴量を作成する構成にすることもできる。
【００６９】
この対象ＦＰ領域分割特徴量作成部１１は、例えば図７０（後述）のように、対象ＦＰ
タイプ２（４９）を、シグナル強度軸に平行な複数の縦分割線と時間軸に平行な複数の
横分割線とにより格子状の領域に分割し、図２の対象ＦＰ領域分割特徴量５１を作成する
。
【００７０】
対象ＦＰ特徴量統合部１３は、対象ＦＰピーク特徴量作成部７で作成された対象ＦＰピ
ーク特徴量４７と対象ＦＰ領域分割特徴量作成工程１１で作成された対象ＦＰ領域分割特
徴量５１とを統合して対象ＦＰ統合特徴量を作成する機能部である。
【００７１】
一方、ＦＰ作成部３の基準ＦＰ作成部３１は、対象ＦＰ作成部２９と同様にして複数の
基準ＦＰを作成する機能部である。例えば、正常品と判定されている複数の漢方薬（基準
漢方薬）の三次元クロマトグラム・データである各３Ｄクロマトから、特定の検出波長に
おける複数のピークとそのリテンション・タイムならびにＵＶスペクトルを抽出した基準
ＦＰを基準漢方薬ごとに作成する。
【００７２】
基準ＦＰピーク帰属部１５も、対象ピーク帰属部５と同様にパターン認識により帰属す
べきピークを特定する機能部である。但し、この基準ＦＰピーク帰属部１５では、全基準
ＦＰを対象とし、選択した組み合わせ且つ順番で帰属スコアを算出することでピークを特
定する。
【００７３】
基準ＦＰ帰属結果統合部１７は、基準ピーク帰属部１５で特定されて帰属したピークを
統合して基準ピーク対応表（後述）を作成する機能部である。
【００７４】
基準ＦＰピーク特徴量作成部１９は、基準ＦＰ帰属結果統合部１７で作成された基準ピ
ーク対応表に基づき前記複数の基準ＦＰを特徴量化した基準ＦＰピーク特徴量を作成する
機能部である。
【００７５】
基準ＦＰタイプ２作成部２１は、対象ＦＰタイプ２作成部９と同様に機能し、複数の各
基準ＦＰから前記特徴量化されたピークを除いて残ったピークとそのリテンション・タイ
ムとで構成されるＦＰを基準ＦＰタイプ２をパターンとして作成する機能部である。
【００７６】
基準ＦＰ領域分割特徴量作成部２３は、対象ＦＰ領域分割特徴量作成部１１と同様に機
能し、ＦＰ領域分割特徴量作成部として基準ＦＰタイプ２を複数の領域に分割し各領域に
存在するピークの存在率から基準ＦＰ領域分割特徴量を作成する機能部である。
【００７７】
但し、基準ＦＰ領域分割特徴量作成部２３では、分割された各領域の位置を変更し変更
前後で基準ＦＰ領域分割特徴量を作成する。すなわち、各縦・横分割線を設定範囲内で平
行移動させるように位置を変更設定することで前記各領域の位置を変更する。
【００７８】
基準ＦＰ特徴量統合部２５は、対象ＦＰ特徴量統合部１３と同様に機能し、基準ＦＰピ
ーク特徴量と基準ＦＰ領域分割特徴量とを統合して基準ＦＰ統合特徴量を作成する機能部
である。
【００７９】
評価部２７は、対象パターン統合特徴量とこの対象パターン統合特徴量に対応し評価基
準となる複数の基準パターンに基づく基準パターン統合特徴量とを比較評価する。すなわ
ち、評価部２７は、対象パターン統合特徴量としての対象ＦＰ統合特徴量と基準パターン
統合特徴量としての基準ＦＰ統合特徴量とを比較評価する機能部である。実施例では、Ｍ
Ｔ法で対象ＦＰ統合特徴量と基準ＦＰ統合特徴量との同等性を評価する。
【００８０】
ＭＴ法は、品質工学で一般に知られている計算手法を意味する。例えば、「品質工学の
数理」日本規格協会発行（２０００）第１３６−１３８頁、品質工学応用講座「化学・薬
学・生物学の技術開発」日本規格協会編（１９９９）第４５４−４５６頁及び品質工学
１１（５），７８−８４（２００３）、入門 ＭＴシステム（２００８）に記載がある。
【００８１】
また一般に市販されているＭＴ法プログラムソフトも使用できる。市販のＭＴ法プログ
ラムソフトとしては、アングルトライ（株）のＡＴＭＴＳ； （財）日本規格協会のＴＭ
−ＡＮＯＶＡ；（株）オーケンのＭＴ法 ｆｏｒ ｗｉｎｄｏｗｓ等が挙げられる。
【００８２】
評価部２７は、対象ＦＰ４３のうち、漢方薬のロット番号と、リテンション・タイム又
はＵＶ検出波長の一方とに対して、ＭＴ法における変数軸を割り振り、ピークをＭＴ法に
おける特徴量とする。
【００８３】
変数軸の割り振りには特に限定はないが、ＭＴ法におけるいわゆる項目軸にリテンショ
ン・タイムを割り振り、いわゆる番号列軸に多成分系薬剤の番号を割り振り、ＭＴ法にお
けるいわゆる特徴量にピークを割り振ることが好ましい。
【００８４】
ここで、上記項目軸と番号列軸は、以下のように定義される。すなわちＭＴ法において
は、データセットＸijについて、平均値mjと標準偏差σjを求め、Ｘijを規準化した値で
あるｘij＝（Ｘij−ｍ_j）／σ_jから、iとjの相関係数rを求めて、単位空間やマハラノビ
スの距離を得る。この時、「平均値mjと標準偏差σjは、項目軸の値ごとに、番号列軸の
値を変化させて求める」というように項目軸と番号列軸とは定義される。
【００８５】
軸が割り振られたデータと特徴量とから、ＭＴ法を用いて、基準点と単位量（以下、「
単位空間」と略記することがある）を得る。ここで、基準点、単位量及び単位空間は、上
記ＭＴ法の文献の記載に従い定義される。
【００８６】
ＭＴ法により、評価すべき薬剤の単位空間との相違の程度を表す値としてＭＤ値を得る
。ここでＭＤ値は、ＭＴ法の文献の説明と同様に定義され、またＭＤ値は、文献に記載の
方法で求められる。
【００８７】
このようにして得られたＭＤ値を用いて、評価すべき薬剤は、正常品として定めた複数
の薬剤からの相違の程度を判定し評価することができる。
【００８８】
例えば、図８７〜図９１の各対象ＦＰを前記のように帰属処理することで、上記ＭＴ法
によりＭＤ値（ＭＤ値：０．２６、２．２０等）を求めることができる。
【００８９】
このＭＤ値を正常品のＭＤ値に対して評価する場合、正常品と定めた複数の薬剤で同様
にＭＤ値を求める。この正常品のＭＤ値から閾値を設定し、図２の評価部２７の評価結果
５３に示すように評価対象薬剤のＭＤ値をプロットし、正常品か非正常品かの判定を行う
ことができる。図２の評価部２７の評価結果５３では、例えばＭＤ値１０以下を正常品と
している。
【００９０】
なお、評価部２７は、対象ＦＰ統合特徴量と基準ＦＰ統合特徴量との同等性を比較評価
できれば良いため、ＭＴ法以外のパターン認識手法等を適用することも可能である。
［ピーク・パターン処理の動作原理］
図５〜図６９は、前記基準ＦＰ選定部３３、ピーク・パターン作成部３５、ピーク帰属
部３７、及び対象ＦＰピーク特徴量作成部７の動作原理を説明するものである。
【００９１】
図５〜図９は、基準ＦＰ選定部３３に係る対象ＦＰと基準ＦＰとのリテンション・タイ
ム・出現パターンの一致度について説明した図である。図５は、対象ＦＰ及び基準ＦＰの
リテンション・タイムを示す図、図６は、対象ＦＰのリテンション・タイム・出現パター
ンを示す図、図７は、基準ＦＰのリテンション・タイム・出現パターンを示す図である。
図８は、対象ＦＰと基準ＦＰのリテンション・タイム・出現距離の一致数を示す図、図９
は、対象ＦＰと基準ＦＰのリテンション・タイム・出現パターンの一致度を示す図である
。
【００９２】
図５では、対象ＦＰ６１及び基準ＦＰ８３それぞれのリテンション・タイムを示す。図
６、図７では、対象ＦＰ６１及び基準ＦＰ８３それぞれのリテンション・タイムから全て
のリテンション・タイム間距離を算出し、それら距離を表形式にまとめたリテンション・
タイム・出現パターンを示す。図８では、これらの出現パターンからリテンション・タイ
ム・出現距離の一致数を算出し、それら一致数を表形式にまとめたリテンション・タイム
・出現距離の一致数を示す。図９では、この一致数をもとにリテンション・タイム・出現
パターンの一致度を算出し、それら一致度を表形式にまとめたリテンション・タイム・出
現パターンの一致度を示す。
【００９３】
図１０〜図１２は、ピーク・パターン作成部３５に係る帰属対象ピークとその周辺ピー
クで作成したピーク・パターンについて説明した図である。図１０は、対象ＦＰの帰属対
象ピークを示す図、図１１は、周辺ピーク２本を含めたピーク３本で作成したピーク・パ
ターンについて、図１２は、周辺ピーク４本を含めたピーク５本で作成したピーク・パタ
ーンについて説明した図である。
【００９４】
図１３、図１４は、ピーク・パターン作成部３５に係る帰属対象ピークと帰属候補ピー
クの関係について説明し、図１３は、帰属対象ピークの許容幅を示す図、図１４は、帰属
対象ピークに対する基準ＦＰの帰属候補ピークを示す図である。
【００９５】
図１５〜図１８は、ピーク・パターン作成部３５に係るピーク３本で作成した帰属対象
ピーク及び帰属候補ピークのピーク・パターン例である。図１５は、帰属対象ピークと帰
属候補ピークのピーク３本によるピーク・パターン図、図１６は、帰属対象ピークと別な
帰属候補ピークのピーク３本によるピーク・パターン図、図１７は、帰属対象ピークと別
な帰属候補ピークのピーク３本によるピーク・パターン図、図１８は、帰属対象ピークと
別な帰属候補ピークのピーク３本によるピーク・パターン図である。
【００９６】
図１９〜図２２は、ピーク・パターン作成部３５に係るピーク５本で作成した帰属対象
ピーク及び帰属候補ピークのピーク・パターン図である。
【００９７】
図２３〜図６１は、ピーク・パターン作成部３５に係る帰属対象ピーク及び帰属候補ピ
ークのピーク・パターンを網羅的に作成し、比較する網羅的比較の原理を説明した図であ
る。
【００９８】
図６２、図６３は、ピーク帰属部３７に係るピーク３本で作成したピーク・パターンの
一致度の算出方法について説明した図である。
【００９９】
図６４は、ピーク帰属部３７に係るピーク５本で作成したピーク・パターンの一致度の
算出方法について説明した図である。
【０１００】
図６５は、ピーク帰属部３７に係る帰属対象ピーク７３及び帰属候補ピーク９５のＵＶ
スペクトル１３５及び１３９を示した図である。
【０１０１】
図６６は、ピーク帰属部３７に係る帰属対象ピーク７３のＵＶスペクトル１３５と帰属
候補ピーク９５のＵＶスペクトル１３９の一致度について説明した図である。
【０１０２】
図６７は、ピーク帰属部３７に係る帰属対象ピーク７３と帰属候補ピーク９５のピーク
・パターンの一致度とＵＶスペクトルの一致度から算出する帰属候補ピークの一致度につ
いて説明した図である。
【０１０３】
図６８は、ピーク帰属部３７に係る対象ＦＰ４３における各ピークの基準群ＦＰ４５へ
の帰属を説明した図である。
【０１０４】
図６９は、ピーク帰属部３７に係る対象ＦＰ４３の各ピークが基準群ＦＰ４５に帰属さ
れた状況を示す対象ＦＰピーク特徴量４７を説明した図である。
【０１０５】
（基準ＦＰの選定）
前記基準ＦＰ選定部３３の機能を、図５〜図９を用いてさらに説明する。
【０１０６】
図５は、対象ＦＰ及び基準ＦＰのリテンション・タイムを示す図、図６は、対象ＦＰの
リテンション・タイム・出現パターンを示す図、図７は、基準ＦＰのリテンション・タイ
ム・出現パターンを示す図である。図８は、対象ＦＰと基準ＦＰのリテンション・タイム
・出現距離の一致数を示す図、図９は、対象ＦＰと基準ＦＰのリテンション・タイム・出
現パターンの一致度を示す図である。
【０１０７】
図５では、対象ＦＰ６１及び基準ＦＰ８３それぞれのリテンション・タイムを示す。図
６、図７では、対象ＦＰ６１及び基準ＦＰ８３それぞれのリテンション・タイムから全て
のリテンション・タイム間距離を算出し、それら距離を表形式にまとめたリテンション・
タイム・出現パターンを示す。図８では、これらの出現パターンからリテンション・タイ
ム・出現距離の一致数を算出し、それら一致数を表形式にまとめたリテンション・タイム
・出現距離の一致数を示す。図９では、この一致数をもとにリテンション・タイム・出現
パターンの一致度を算出し、それら一致度を表形式にまとめたリテンション・タイム・出
現パターンの一致度を示す。
【０１０８】
対象ＦＰ６１のピーク帰属処理において、対象ＦＰ６１とできるだけＦＰパターンが類
似した基準ＦＰで対象ＦＰ６１の各ピークを帰属する。この対象ＦＰ６１に類似した基準
ＦＰを複数の基準ＦＰから選定することが精度の高い帰属を行う上で重要なポイントであ
る。
【０１０９】
そこで、対象ＦＰ６１のＦＰパターンとの類似性を客観的かつ簡易的に評価する方法と
して、リテンション・タイム・出現パターンの一致度によりＦＰパターンの類似性を評価
する。
【０１１０】
例えば、対象ＦＰ６１及び基準ＦＰ８３のリテンション・タイムが図５のような場合、
対象ＦＰ６１及び基準ＦＰ８３それぞれのリテンション・タイム・出現パターンは、図６
、図７のようになる。図６、図７では、上段の対象ＦＰ６１及び基準ＦＰ８３に対し、下
段の図表のように、各セルの値がリテンション・タイム間距離で構成された表形式のパタ
ーンとして作成している。
【０１１１】
図６において、対象ＦＰ６１の各ピーク（６３、６５、６７、６９、７１、７３、７５
、７７、７９、８１）のリテンション・タイムは、（１０．２）、（１０．５）、（１０
．８）、（１１．１）、（１１．６）、（１２．１）、（１２．８）、（１３．１）、（
１３．６）、（１４．０）となっている。
【０１１２】
したがって、ピーク６３及びピーク６５間のリテンション・タイム間距離は、（１０．
５）−（１０．２）＝（０．３）となる。同様に、ピーク６３及びピーク６７間は、（０
．６）、ピーク６５及びピーク６７間は、（０．３）などとなる。以下、同様であり、図
６の下段図表の対象ＦＰ出現パターンとなる。
【０１１３】
図７において、基準ＦＰ８３の各ピーク（８５、８７、８９、９１、９３、９５、９７
、９９、１０１、１０３、１０５）のリテンション・タイムは、（１０．１）、（１０．
４）、（１０．７）、（１１．１）、（１１．７）、（１２．３）、（１２．７）、（１
３．１）、（１３．６）、（１４．１）、（１４．４）となっている。
【０１１４】
したがって、同様にリテンション・タイム間距離は、図７の下段図表の基準ＦＰ出現パ
ターンとなる。
【０１１５】
この図６、図７ようにパターン化した各ピークを総当たりで比較し一致数を求める。例
えば、図６下段図表の対象ＦＰ出現パターンの各セルの値と図７下段図表の基準ＦＰ出現
パターンの各セルの値とを比較し、図８のように一致数を得る。
【０１１６】
すなわち、対象ＦＰ６１と基準ＦＰ８３のリテンション・タイム・出現パターンの全て
のリテンション・タイム間距離を行単位で順番に総当たりで比較し、設定した範囲内で距
離が一致した数を算出した。
【０１１７】
例えば、図６、図７の対象及び基準ＦＰリテンション・タイム・出現パターンの１行を
比較すると、一致数は７個である。この７個の一致数が、図８の対象及び基準ＦＰリテン
ション・タイム・出現パターンの１行目に書き込まれる。図６、図7中のその他の行につ
いても同様であり、対象ＦＰリテンション・タイム・出現パターンの１行〜９行までと、
基準ＦＰリテンション・タイム・出現パターンの１行〜１０行までとを総当たりで比較し
、それぞれ一致数が得られる。
【０１１８】
図８に結果を示した。この図８において、丸で囲まれた左端の７の数値は、対象及び基
準ＦＰリテンション・タイム・出現パターンの１行目を比較した結果であり、その隣の７
の数値は、対象ＦＰリテンション・タイム・出現パターンの１行目と基準ＦＰリテンショ
ン・タイム・出現パターンの２行目とを比較した結果である。設定の範囲に限定はなく、
好ましくは０．０５分〜０．２分の範囲である。実施例１は、０．１分とした。
【０１１９】
リテンション・タイム・出現パターンの一致度をＲＰとすると、対象ＦＰ６１のｆ行目
のリテンション・タイム・出現パターンと基準ＦＰ８３のｇ行目のリテンション・タイム
・出現パターンの一致度（ＲＰ_ｆｇ）は、Tanimoto係数を用いて、
ＲＰ_ｆｇ ＝ ｛１−（ｍ／（ａ＋ｂ−ｍ））｝×（ａ−ｍ＋１）
として算出する。
【０１２０】
なお、式中ａは対象ＦＰ６１のピーク数（対象ＦＰピーク数）、ｂは基準ＦＰ８３のピ
ーク数（基準ＦＰピーク数）、ｍはリテンション・タイム・出現パターンの一致数（出現
距離の一致数）である（図８参照）。図８の一致数をもとに前記式により各リテンション
・タイム・出現パターンの一致度（ＲＰ）を算出した（図９参照）。
【０１２１】
これらＲＰの最小値であるＲＰ＿ｍｉｎを対象ＦＰ６１と基準ＦＰ８３とのリテンショ
ン・タイム・出現パターンの一致度とする。図９では、（０．５０）が対象ＦＰ６１の基
準ＦＰに対する一致度となる。
【０１２２】
この一致度を全ての基準ＦＰについて算出し、最も小さい一致度の基準ＦＰが選定され
、この基準ＦＰに対して対象ＦＰのピーク帰属を行わせる。
【０１２３】
基準ＦＰ選定部５は、対象ＦＰ６１及び基準ＦＰ８３を、ピーク高さ比でパターン化す
ることもできる。
【０１２４】
ピーク高さ比でパターン化した各ピークを総当たりで比較し、設定した範囲内で高さ比
が一致した数を算出する。この算出により図８と同様に一致数を得ることができる。
【０１２５】
なお、ピーク高さ比でパターン化する場合は、１行中に、同じような値が複数存在する
ケースがあり、これらを複数回カウントしないようにしなければならない。
【０１２６】
一致度は、Tanimoto係数を、「高さ比の一致数／（対象ＦＰピーク数＋基準ＦＰピーク
数−高さ比の一致数）」とし、（１−Tanimoto係数）が零に近いことで前記一致度を求め
ることができる。
【０１２７】
また、（１−Tanimoto係数）に、（対象ＦＰピーク数−出現パターン又は高さ比の一致
数＋１）の重み付けをし、「（１−Tanimoto係数）×（対象ＦＰピーク数−出現パターン
又は高さ比の一致数＋１」とし、重み付けにより対象ＦＰ６１のピーク（６３、６５、・
・・）がより多く一致している基準ＦＰを選ぶことができる。
【０１２８】
（ピーク・パターンによる特徴量化）
前記ピーク・パターン作成部３５の機能を、図１０〜図６７を用いてさらに説明する。
【０１２９】
図１０のように帰属対象ピーク７３を基準ＦＰ８３のいずれかのピークに帰属するとき
、何れのピークに帰属すべきか、ということになる。仮に、このピーク帰属を、ピーク、
リテンション・タイム又はＵＶスペクトルのいずれかの情報のみで行ったとすると、これ
ら３つの情報のいずれもが前記薬剤間誤差と分析誤差に起因する誤差を含んでいるため、
単独情報によるピーク帰属の精度には限界がある。
【０１３０】
また、図１３、図１４のように帰属対象ピーク７３と基準ＦＰ８３の各ピーク間でリテ
ンション・タイムのズレの許容幅を設定し、その許容幅内に存在する基準ＦＰ８３のピー
ク（以下、帰属候補ピーク）とＵＶスペクトル情報の２つの情報によるピーク帰属では、
すべての情報を総合して帰属先を判定しているため、上記単独情報によるピーク帰属に比
べ精度は向上する。
【０１３１】
しかし、３つの情報を使ったピーク帰属であったとしても、ＵＶスペクトルの特性とし
て、類似成分のＵＶスペクトルはほとんど同じとなってしまうため、帰属候補ピークに複
数の類似成分が含まれている場合は、結局ピーク情報のみでの帰属となってしまい、十分
な精度は得られない。そのため、より精度の高いピーク帰属を行うためには、これら３つ
の情報にプラスする情報が必要である。
【０１３２】
そこで、図１１、図１２のような周辺ピークの情報を含めたピーク・パターンを作成し
、このピーク・パターンの比較によってピークを帰属することにした。
【０１３３】
周辺のピークを含めたピーク・パターンとすると、これまでの３つの情報に周辺情報が
プラスされ、４つの情報によるピーク帰属が可能となり、より高い帰属精度が得られる。
【０１３４】
その結果として、一度の帰属処理で大量のピークを高精度かつ効率的に一斉帰属するこ
とができる。
【０１３５】
さらに、ピーク帰属で使用するデータを周辺情報を含めた４つの情報にすることで、既
存のピーク帰属の際に設定する制約条件（ピーク定義など）も必要なくなった。
【０１３６】
図１１では、帰属対象ピーク７３に対し、時間軸方向の双方に存在するピーク７１、７
５を含めたピーク・パターン１１５を作成した。
【０１３７】
図１２では、帰属対象ピーク７３に対し、時間軸方向の双方に存在するピーク６９、７
１、７５、７７を含めたピーク・パターン１２５を作成した。
【０１３８】
図１３、図１４では、帰属対象ピーク７３と基準ＦＰ８３の各ピーク間でリテンション
・タイムのズレの許容幅を設定し、その許容幅内に存在する基準ＦＰ８３のピークを帰属
対象ピーク７３に対応する候補ピーク（以下、帰属候補ピーク）とした。
【０１３９】
図１５では、帰属対象ピーク７３のピーク・パターン１１５と比較するピーク・パター
ンとして、帰属候補ピーク９３に対し、時間軸方向前後の双方に存在するピーク９１、９
５を含めたピーク・パターン１１７を作成した。
【０１４０】
図１６〜図１８では、帰属対象ピーク７３のピーク・パターン１１５と比較するピーク
・パターンとして、別な帰属候補ピーク９５、９７、９９に対し、それぞれ時間軸方向前
後の双方に存在するピークを含めたピーク・パターン１１９、１２１、１２３を作成した
。
【０１４１】
このピーク・パターンの比較をより高い精度で行わせるには、図１９〜図２２のように
、対象ＦＰならびに基準ＦＰの両方で周辺のピーク数を増加させたピーク・パターンを作
成することが肝要である。
【０１４２】
例えば、周辺のピーク４本を含めた計５本のピークによるピーク・パターンの比較にす
ると、より高い帰属精度が得られる。
【０１４３】
図１９では、帰属対象ピーク７３のピーク・パターン１２５と比較するピーク・パター
ンとして、帰属候補ピーク９３に対し、時間軸方向の双方に存在するピーク８９、９１、
９５、９７を含めたピーク・パターン１２７を作成した。
【０１４４】
図２０〜図２２では、帰属対象ピーク７３のピーク・パターン１２５と比較するピーク
・パターンとして、別な帰属候補ピーク９５、９７、９９に対し、それぞれ時間軸方向前
後の双方に存在するピークを含めたピーク・パターン１２９、１３１、１３３を作成した
。
【０１４５】
さらに、このピーク・パターンによる帰属をより高い精度で行わせるには、対象ＦＰと
基準ＦＰのピーク数に違いがある（つまり、どちらか一方に存在しないピークがある）場
合の対応が必要である。そのためには図２３〜図２５のように帰属対象ピークならびに帰
属候補ピークの両方でピーク・パターン構成ピークを網羅的に変化させたピーク・パター
ンを作成することが肝要である。
【０１４６】
具体的には、対象ＦＰの帰属対象ピークの周辺ピークの中から予めピーク・パターン構
成ピークの候補となるピーク（以下、ピーク・パターン構成候補ピーク）を設定し、この
ピーク・パターン構成候補ピークを順番にピーク・パターン構成ピークとしてピーク・パ
ターンを作成する。基準ＦＰの帰属候補ピークについても同様にピーク・パターン構成候
補ピークを設定し、このピーク・パターン構成候補ピークを順番にピーク・パターン構成
ピークとしてピーク・パターンを作成する。
【０１４７】
例えば、図２３のように帰属対象ピーク７３のピーク・パターン構成候補ピークとして
時間軸方向周辺４本（６９、７１、７５、７７）、帰属候補ピーク９３のピーク・パター
ン構成候補ピークとして時間軸方向周辺４本（８９、９１、９５、９７）とし、ピーク・
パターン構成ピークを任意の２本にそれぞれ設定する。この場合、図２４、図２５のよう
に帰属対象ピーク７３及び帰属候補ピーク９３それぞれでピーク・パターンが4Ｃ2（＝6
）パターン作成される。
【０１４８】
さらに、ピーク・パターン構成候補ピークを１０本とし、ピーク・パターン構成ピーク
を任意の２本に設定すると、帰属対象ピーク及び帰属候補ピークそれぞれで、10Ｃ2（＝
４５）パターンのピーク・パターンが作成される。ピーク・パターン構成ピークを任意の
４本に設定すると、帰属対象ピーク及び帰属候補ピークそれぞれで、10Ｃ4（＝２１０）
パターンのピーク・パターンが作成される。
【０１４９】
前記ピーク帰属部３７の機能を、図２６〜図６９を用いてさらに説明する。
【０１５０】
ピーク帰属部３７では、ピーク・パターン作成部３５で作成した帰属対象ピークと帰属
候補ピークの全ピーク・パターン間で、対応するピーク及びリテンション・タイムの差を
もとにピーク・パターンの一致度（以下、Ｐ＿Ｓｉｍ）を算出する。ピーク帰属部３７は
、Ｐ＿Ｓｉｍの最小値（以下、Ｐ＿Ｓｉｍ＿ｍｉｎ）を帰属対象ピークと帰属候補ピーク
のピーク・パターンの一致度とする。
【０１５１】
例えば、図２６〜図６１のように帰属対象ピーク７３及び帰属候補ピーク９３のそれぞ
れでピーク・パターン構成候補ピークを時間軸方向前後周辺４本とし、ピーク・パターン
構成ピークを任意の２本に設定する。この設定では、帰属対象ピーク及び帰属候補ピーク
それぞれで、4Ｃ2（＝６）パターンのピーク・パターンが作成される。したがって、帰属
対象ピーク７３及び帰属候補ピーク９３のＰ＿Ｓｉｍは、６パターン×６パターン（＝３
６）通りで算出され、これらＰ＿Ｓｉｍの最小値であるＰ＿Ｓｉｍ＿ｍｉｎを帰属対象ピ
ーク７３と帰属候補ピーク９３の一致度とする。
【０１５２】
ちなみに、帰属対象ピーク７３及び帰属候補ピーク９３のそれぞれでピーク・パターン
構成候補ピークを時間軸方向前後周辺１０本とし、ピーク・パターン構成ピークを任意の
２本に設定すると、帰属対象ピーク及び帰属候補ピークそれぞれで、10Ｃ2（＝４５）パ
ターンのピーク・パターンが作成される。したがって、帰属対象ピーク７３及び帰属候補
ピーク９３のＰ＿Ｓｉｍは、４５パターン×４５パターン（＝２０２５）通りで算出され
、これらＰ＿Ｓｉｍの最小値であるＰ＿Ｓｉｍ＿ｍｉｎを帰属対象ピーク７３と帰属候補
ピーク９３の一致度とする。また、ピーク・パターン構成ピークを任意の４本に設定する
と、帰属対象ピーク及び帰属候補ピークそれぞれで、10Ｃ4（＝２１０）パターンのピー
ク・パターンが作成される。したがって、帰属対象ピーク７３及び帰属候補ピーク９３の
Ｐ＿Ｓｉｍは、２１０パターン×２１０パターン（＝４４１００）通りで算出され、これ
らＰ＿Ｓｉｍの最小値であるＰ＿Ｓｉｍ＿ｍｉｎを帰属対象ピーク７３と帰属候補ピーク
９３の一致度とする。
【０１５３】
このＰ＿Ｓｉｍは、帰属対象ピーク７３の全ての帰属候補ピークについて同様に算出す
る。
【０１５４】
図６２、図６３で、３本のピークで構成されたピーク・パターンを比較するためのピー
ク・パターンの一致度の算出方法を説明する。この場合、帰属対象ピーク７３のピーク・
パターン１１５と帰属候補ピーク９５のピーク・パターン１１９を例にする。
【０１５５】
帰属対象ピーク７３のピーク・パターン１１５において、帰属対象ピーク７３のピーク
及びリテンション・タイムをｐ１及びｒ１、ピーク・パターン構成ピーク７１のピーク及
びリテンション・タイムをｄｎ１及びｃｎ１、ピーク・パターン構成ピーク７５のピーク
及びリテンション・タイムをｄｎ２及びｃｎ２とする。
【０１５６】
帰属候補ピーク９５のピーク・パターン１１９において、帰属候補ピーク９５のピーク
及びリテンション・タイムをｐ２及びｒ２、ピーク・パターン構成ピーク９３のピーク及
びリテンション・タイムをｆｎ１及びｅｎ１、ピーク・パターン構成ピーク９７のピーク
及びリテンション・タイムをｆｎ２及びｅｎ２とする。
【０１５７】
ピーク・パターンの一致度をＰ＿Ｓｉｍとすると、帰属対象ピーク７３と帰属候補ピー
ク９５の３本のピークで構成するピーク・パターンの一致度（Ｐ＿Ｓｉｍ（７３−９５）
）は、
Ｐ＿Ｓｉｍ（７３−９５） ＝ (|p1−p2|＋1)×(|(r1−(r2＋d)|＋1)
＋(|dn1−fn1|＋1)×(|(cn1−r1)−(en1−r2)|+1)
＋(|dn2−fn2|＋1)×(|(cn2−r1)−(en2−r2)|+1)
として算出する。
【０１５８】
なお、式中のｄは、リテンション・タイムのずれを補正する値である。
【０１５９】
図６４で、５本のピークで構成されたピーク・パターンを比較するためのピーク・パタ
ーンの一致度の算出方法を説明する。この場合、帰属対象ピーク７３のピーク・パターン
１２５と帰属候補ピーク９５のピーク・パターン１２９を例にする。
【０１６０】
帰属対象ピーク７３のピーク・パターン１２５において、帰属対象ピーク７３のピーク
及びリテンション・タイムをｐ１及びｒ１、ピーク・パターン構成ピーク６９、７１、７
５、７７のピーク及びリテンション・タイムをそれぞれｄｎ１及びｃｎ１、ｄｎ２及びｃ
ｎ２、ｄｎ３及びｃｎ３、ｄｎ４及びｃｎ４とする。
【０１６１】
帰属候補ピーク９５のピーク・パターン１２９において、帰属候補ピーク９５のピーク
及びリテンション・タイムをｐ２及びｒ２、ピーク・パターン構成ピーク９１、９３、９
７、９９のピーク及びリテンション・タイムをそれぞれｆｎ１及びｅｎ１、ｆｎ２及びｅ
ｎ２、ｆｎ３及びｅｎ３、ｆｎ４及びｅｎ４とする。
【０１６２】
帰属対象ピーク７３と帰属候補ピーク９５の５本のピークで構成するピーク・パターン
の一致度（Ｐ＿Ｓｉｍ（７３−９５））は、
Ｐ＿Ｓｉｍ（７３−９５） ＝ (|p1−p2|＋1)×(|(r1−(r2＋d)|＋1)
＋(|dn1−fn1|＋1)×(|(cn1−r1)−(en1−r2)|+1)
＋(|dn2−fn2|＋1)×(|(cn2−r1)−(en2−r2)|+1)
＋(|dn3−fn3|＋1)×(|(cn3−r1)−(en3−r2)|+1)
＋(|dn4−fn4|＋1)×(|(cn4−r1)−(en4−r2)|+1)
として算出する。
【０１６３】
なお、式中のｄは、リテンション・タイムのずれを補正する値である。
【０１６４】
ピーク帰属部３７では、図６７、図６８のように帰属対象ピークと帰属候補ピークでＵ
Ｖスペクトルの一致度を算出する。
【０１６５】
図６５は、帰属対象ピーク７３及び帰属候補ピーク９５のＵＶスペクトル（１３５と１
３９）図であり、図６６のようにこれら２つのＵＶスペクトルの一致度（ＵＶ＿Ｓｉｍ（
７３−９５））は、
ＵＶ＿Ｓｉｍ（７３−９５）＝ ＲＭＳＤ（１３５ ｖｓ １３９）
として算出する。
【０１６６】
ＲＭＳＤは、平均二乗偏差のことで、対応する２点の距離（ｄｉｓ）をそれぞれ二乗し
、その相加平均の平方根として定義される。つまり、
ＲＭＳＤ ＝ √｛ Σｄｉｓ² ／ ｎ ｝
として算出する。
【０１６７】
ｎは、ｄｉｓの数である。
【０１６８】
ここで、ＵＶスペクトルの波形は極大波長及び極小波長を含んでおり、この極大波長及び極小波長或いは何れかを比較することで一致度を算出することも可能である。しかし、吸収特性のない化合物あるいは吸収特性が類似した化合物等では、極大波長、極小波長は同じであるが、全体の波形がかなり異なる場合もあり、極大波長、極小波長の比較では、波形の一致度を算出できない恐れがある。
【０１６９】
これに対し、ＵＶスペクトルの波形によりＲＭＳＤを利用した場合には、波形全体の比較となるため、ＵＶスペクトルの波形の一致度の算出をより正確に行わせ、吸収特性のない化合物あるいは吸収特性が類似した化合物でも正確に識別できることになる。
【０１７０】
このＵＶスペクトルの一致度は、帰属対象ピーク７３の全ての帰属候補ピークについて
同様に算出する。
【０１７１】
さらに、ピーク帰属部３７では、図６７のように上記２つの一致度を統合した帰属候補
ピークの一致度を算出する。
【０１７２】
図６７のように帰属候補ピークの一致度（ＳＣＯＲＥ（７３−９５））はピーク・パタ
ーンとＵＶスペクトルのそれぞれの一致度を乗じて算出する。ピーク・パターン７３、９
５の一致度を示すスコアが、Ｐ＿Ｓｉｍ＿ｍｉｎ(７３−９５)であり、対応するＵＶ波形
データ１３５、１３９の一致度を示すスコアが、ＵＶ＿Ｓｉｍ（７３−９５）であるとす
る。このとき、帰属候補ピークの一致度 ＳＣＯＲＥ（７３−９５）は、
ＳＣＯＲＥ（７３−９５）
＝ Ｐ＿Ｓｉｍ＿ｍｉｎ(７３−９５)×ＵＶ＿Ｓｉｍ（７３−９５）
として算出する。
【０１７３】
この帰属候補ピークの一致度を帰属対象ピーク７３の全ての帰属候補ピークについて同
様に計算する。
【０１７４】
そして、全帰属候補ピーク間でこのＳＣＯＲＥを比較し、ＳＣＯＲＥが最小となる帰属
候補ピークを帰属対象ピーク７３の帰属ピークとして判定する。
【０１７５】
ピーク帰属部３７では、帰属対象ピークの帰属すべきピークを２つの観点を総合して判
定するため、正確なピーク帰属を実現することができる。
【０１７６】
また、対象ピーク特徴量作成部７では、対象ＦＰの基準ＦＰへの帰属結果をもとに、図
６８のように、対象ＦＰ４３の各ピークを基準群ＦＰ４５へ帰属する。
【０１７７】
対象ＦＰ４３の各ピークは、前記帰属処理により基準群ＦＰ４５を構成する基準ＦＰに
帰属される。この帰属結果をもとに最終的に基準群ＦＰ４５のピークに帰属する。
【０１７８】
尚、基準群ＦＰ４５は、正常品と定めた複数の基準ＦＰ全てを前記のように帰属処理し
作成したものであり、その各ピークは帰属されたピークの平均値（黒点）±標準偏差（縦
分割線）で表している。
【０１７９】
図６９は、対象ＦＰ４３を基準群ＦＰ４５に帰属した結果であり、この結果が対象ＦＰ
４３の対象ＦＰピーク特徴量４７である。
［ＦＰ領域分割特徴量作成の動作原理］
図７０〜図８６は、ＦＰ領域分割特徴量作成の動作原理を示し、図７０は、領域分割に
よる数量化を示す説明図、図７１は、リテンション・タイム等の変動との関係を示す説明
図、図７２は、領域の位置を変更して数量化する説明図、図７３は、ＦＰタイプ２のデー
タを示す図表、図７４は、ＦＰタイプ２のパターンを示す説明図、図７５は、縦・横分割
線での領域分割による領域ごとの特徴量化を示す説明図、図７６は、縦分割線（1本目）
の設定を示す説明図、図７７は、横分割線（1本目）の設定を示す説明図、図７８は、縦
・横分割線による領域分割を示す説明図、図７９は、特徴量化する領域の数を示す説明図
、図８０は、領域１の特定を示す説明図、図８１は、全ピークの高さ及び合計を示す図表
、図８２は、領域１のピーク高さの合計を示す説明図、図８３は、最初の1パターンによ
る全領域の特徴量を示す図表、図８４は、縦１本目の位置を順次変更してできた各領域で
の特徴量を示す図表、図８５は、横１本目の位置を順次変更してできた各領域での特徴量
を示す図表、図８６は、各縦・横分割線の位置を変更しない１通りでの特徴量を示す図表
である。
【０１８０】
前記対象ＦＰ領域分割特徴量作成部１１又は基準ＦＰ領域分割特徴量作成部２３は、前
記のように対象ＦＰタイプ２又は基準ＦＰタイプ２を分割した各領域に存在するピークの
存在率から対象ＦＰ領域分割特徴量又は基準ＦＰ領域分割特徴量を作成する。
【０１８１】
領域の分割は、例えば図７０のように行う。図７０では、例えば薬剤ＡのＦＰ５５につ
いて分割している。シグナル強度軸に平行な複数の縦分割線１４１と時間軸に平行な複
数の横分割線１４３とにより分割し、複数の領域である複数の格子１４５を作成する。
【０１８２】
複数の横分割線１４３は、本実施例では、シグナル強度が増大する方向へ等比間隔で
設定される。この設定により、ピークの密集する部分での領域分割を細分化し、ピークの
存在率をより正確に把握することが可能となる。但し、複数の横分割線１４３の本数を増
加するなどして等差間隔で設定することも可能である。
【０１８３】
この各格子１４５内に存在するピーク高さの割合で数量化し、特徴量とする。
【０１８４】
一方、図７１のように、分析条件の僅かなバラツキ等により、リテンション・タイムや
ピーク高さがＦＰ５５Ａ、５５Ｂのように変動する。この変動により各格子１４５内の値
が大きく変動する恐れがある。
【０１８５】
（基準ＦＰタイプ２の場合）
そこで、基準ＦＰタイプ２の場合には、図７２のように、各格子１４５の位置を変更（
シフト）し変更前後で数量化する。この操作により基準ＦＰ領域分割特徴量を正確に作成
することが可能となる。各縦・横分割線１４１、１４３を設定範囲内で平行移動させるよ
うに位置を変更設定することで前記各格子１４５の位置を変更する、
ここで、各格子１４５の位置を変更した数量化について、さらに説明する。
【０１８６】
図７３は、基準ＦＰタイプ２のデータｄ２０２、ｄ２０７、ｄ２０８を一例として示す
。このデータは、リテンション・タイム（ＲＴ）及びピーク高さ（Ｈｅｉｇｈｔ）の情報
のみの構成となっている。このデータは、前記基準ＦＰタイプ２作成部２１において、複
数の各基準ＦＰから特徴量化されたピークを除いて残ったピークとそのリテンション・タ
イムとで構成される基準ＦＰタイプ２に対応し、全ピークの各ＵＶスペクトルは除いてい
る。
【０１８７】
基準ＦＰタイプ２のデータｄ２０２、ｄ２０７、ｄ２０８のパターンは、図７４のよう
になっている。
【０１８８】
これらのＦＰパターンを縦・横分割線１４１、１４３により領域分割し、その領域ごと
に特徴量化する。
【０１８９】
（縦分割線（１本目）の設定）
縦分割線(１本目)の位置を設定するため、図７６のように、１本目のリテンション・タ
イム（ＲＴ）、振幅、刻みを指定する。
【０１９０】
この３つのパラメータをもとに、次の条件で縦１本目の位置が複数箇所設定される。
【０１９１】
縦分割線（１本目）＝ＲＴ−振幅＋（振幅×２／刻み）×ｉ
（ｉ＝０、１、２、・・・・、刻み−１）
例えば、ＲＴ＝１、振幅＝１、刻み回数＝１０ に指定すると、
縦分割線(１本目）= 0.0、 0.2、 0.4、 0.6、 0.8、 1.0、 1.2、 1.4、 1.6、 1.8
が設定される。
【０１９２】
（横分割線（１本目）の設定）
横(１本目)の位置を設定するため、図７７のように、１本目の高さ、振幅、刻みを指定
する。
【０１９３】
この3つのパラメータをもとに、次の条件で横１本目の位置が複数箇所設定される。
【０１９４】
例えば、高さ＝１、振幅＝０．５、刻み回数＝１０ に指定すると、
横分割線（１本目）＝高さ−振幅＋（振幅×２／刻み）×ｉ
（ｉ＝０、１、２、・・・・、刻み−１）
例えば、高さ＝１、振幅＝０．５、刻み回数＝１０ に指定すると、
横分割線(１本目）= 0.5、 0.6、 0.7、 0.8、 0.9、 1.0、 1.1、 1.2、 1.3、 1.4
が設定される。
【０１９５】
（縦・横分割線(１本目)の組み合わせ）
設定された縦・横分割線（１本目）の全ての組み合わせで、順次２本目以降の標本線を
設定し、領域分割する。
【０１９６】
前記の例では次のようになる。
【０１９７】
縦分割線(１本目）×横分割線(１本目）=
（0.0、 0.2、 0.4、 0.6、 0.8、 1.0、 1.2、 1.4、 1.6、 1.8）
×（0.5、 0.6、 0.7、 0.8、 0.9、 1.0、 1.1、 1.2、 1.3、 1.4）＝100通り
この１００通り全て組み合わせで順次２本目以降の分割線を設定し、領域を分割する。
【０１９８】
（縦・横分割線(２本目以降)の設定）
縦分割線２本目以降は指定した間隔（等差）で指定した本数になるまで設定する。
【０１９９】
縦分割線ｉ本目＝縦分割線（ｉ−１）本目＋間隔
（ｉ＝２、・・・・、指定した本数）
横分割線２本目以降は指定した間隔（等比）で指定した本数になるまで設定する。
【０２００】
横分割線ｉ本目＝横分割線（ｉ−１）本目＋間隔×２÷（ｉ−２）
（ｉ＝２、・・・・、指定した本数）
例えば、縦１本目＝０．０、縦間隔＝１０、縦本数＝７、横分割線1本目＝０．５、横
間隔＝１、横本数＝６ の場合、
縦分割線＝0、10、20、30、40、50、60
横分割線＝0.5、1.5、3.5、7.5、15.5、31.5
に設定される。
【０２０１】
（分割線による領域分割）
設定された縦と横の分割線を先程の例をもとＦＰ上に表記すると、図７８のようになる
。
【０２０２】
この縦と横の線で囲まれた領域ごとにＦＰを特徴量化する。
【０２０３】
領域は全部で３０個あるので、図７９のように３０特徴量が得られる。
【０２０４】
（領域ごとの特徴量化）
各領域は、次の式により特徴量化する。
【０２０５】
特徴量＝領域内ピーク高さ合計／全ピーク高さ合計
（特徴量化の方法）
以下に、上の式により、図８０に示すｄ２０２の領域１の特徴量を求める。
【０２０６】
まず、全ピークの高さ合計を算出すると、図８１のように、１５．５４５４７２となる
。
【０２０７】
次に、領域１のピーク高さ合計を、図８２のように算出する。
【０２０８】
従って、領域１の特徴量は、
特徴量 = 2 ／ 15.545472=0.128655
となる。
【０２０９】
（全領域の特徴量化）
上記特徴量化方法により、最初の１パターンによる全領域の特徴量を算出する。図８３
に算出結果を示す。
（縦分割線１本目を順次変更して特徴量化）
縦分割線１本目の位置を順次変更してできた各領域を上記方法で特徴量化する。図８４
に結果を示す。
【０２１０】
（横１本目を順次変更して特徴量化）
横分割線１本目の位置を変更するたびに縦１本目を１通り変更。できた各領域を上記方
法で特徴量化する。図８５に結果を示す。
【０２１１】
この処理により、縦・横分割線１本目が１０箇所ずつあった場合では、
１００行（１００通り）×３１列（ファイル名＋３０特徴量）
のデータとなる。
【０２１２】
（全基準データの特徴量化（基準ｔｙｐｅ２群ＦＰ））
これまでの処理を全基準データで実施する。例えば、基準データがｄ２０２、 ｄ２０
７、 ｄ２０８の３データであった場合は、
３００行（１００通り×３データ）×３１列（ファイル名＋３０特徴量）
となる。
【０２１３】
（対象ＦＰタイプ２の場合）
対象ＦＰタイプ２では、縦・横分割線（１本目）の組み合わせは、１通り(縦（ＲＴ）
＝１、横（高さ）＝１)となるので、この１通りでの特徴量を算出する。
［ＭＤ値］
図８７〜図９１は、前記のように、評価部２７に係る各種対象ＦＰとその評価値（ＭＤ
値）を示した図であり、前記のように各対象ＦＰを前記のように帰属処理することで、評
価部２７にて、上記ＭＴ法によりＭＤ値（ＭＤ値：０．２６、２．２０等）を求めること
ができる。
［多成分薬剤の評価方法］
図９２は、本発明実施例１のＦＰの特徴量作成方法を含む本発明実施例１の多成分薬剤
の評価方法を示す工程図である。
【０２１４】
図９２のように、多成分薬剤の評価方法は、ＦＰ作成工程１４８と、対象ＦＰピーク帰
属工程１４９と、対象ＦＰピーク特徴量作成工程１５１と、対象ＦＰタイプ２作成工程１
５３と、対象ＦＰ領域分割特徴量作成工程１５５と、対象ＦＰ特徴量統合工程１５７と、
基準ＦＰピーク帰属工程１５９と、基準ＦＰ帰属結果統合工程１６１と、基準ＦＰピーク
特徴量作成工程１６３と、基準ＦＰタイプ２作成工程１６５と、基準ＦＰ領域分割特徴量
作成工程１６７と、基準ＦＰ特徴量統合工程１６９と、評価工程１７１とを備えている。
【０２１５】
ＦＰ作成工程１４８は、対象ＦＰ作成工程１７３と、基準ＦＰ作成工程１７５とを備え
ている。
【０２１６】
対象ＦＰピーク帰属工程１４９は、基準ＦＰ選定工程１７７と、ピーク・パターン作成
工程１７９と、ピーク帰属工程１８１とを備えている。
【０２１７】
これら、ＦＰ作成工程１４８、対象ＦＰピーク帰属工程１４９、対象ＦＰピーク特徴量
作成工程１５１、対象ＦＰタイプ２作成工程１５３、対象ＦＰ領域分割特徴量作成工程１
５５、対象ＦＰ特徴量統合工程１５７、基準ＦＰピーク帰属工程１５９、基準ＦＰ帰属結
果統合工程１６１、基準ＦＰピーク特徴量作成工程１６３、基準ＦＰタイプ２作成工程１
６５、基準ＦＰ領域分割特徴量作成工程１６７、基準ＦＰ特徴量統合工程１６９、評価工
程１７１は、本実施例において前記多成分薬剤の評価装置１を用いて行われる。
【０２１８】
ＦＰ作成工程１４８は、図１のＦＰ作成部３の機能により行わせ、同様に対象ＦＰピー
ク帰属工程１４９、対象ＦＰピーク特徴量作成工程１５１、対象ＦＰタイプ２作成工程１
５３、対象ＦＰ領域分割特徴量作成工程１５５、対象ＦＰ特徴量統合工程１５７、基準Ｆ
Ｐピーク帰属工程１５９、基準ＦＰ帰属結果統合工程１６１、基準ＦＰピーク特徴量作成
工程１６３、基準ＦＰタイプ２作成工程１６５、基準ＦＰ領域分割特徴量作成工程１６７
、基準ＦＰ特徴量統合工程１６９、評価工程１７１を、対象ＦＰピーク帰属部５、対象Ｆ
Ｐピーク特徴量作成部７、対象ＦＰタイプ２作成部９、対象ＦＰ領域分割特徴量作成部１
１、対象ＦＰ特徴量統合部１３、基準ＦＰピーク帰属部１５、基準ＦＰ帰属結果統合部１
７、基準ＦＰピーク特徴量作成部１９、基準ＦＰタイプ２作成部２１、基準ＦＰ領域分割
特徴量作成部２３、基準ＦＰ特徴量統合部２５、評価部２７の各機能により行わせる。
【０２１９】
但し、各工程を、それぞれ別々のコンピュータで機能させることができ、例えば対象Ｆ
Ｐ作成工程１７３と、対象ＦＰピーク帰属工程１４９と、対象ＦＰピーク特徴量作成工程
１５１と、対象ＦＰタイプ２作成工程１５３と、対象ＦＰ領域分割特徴量作成工程１５５
と、対象ＦＰ特徴量統合工程１５７と、評価工程１７１とを一つのコンピュータで機能さ
せ、基準ＦＰ作成工程１７５と、基準ＦＰピーク帰属工程１５９と、基準ＦＰ帰属結果統
合工程１６１と、基準ＦＰピーク特徴量作成工程１６３と、基準ＦＰタイプ２作成工程１
６５と、基準ＦＰ領域分割特徴量作成工程１６７と、基準ＦＰ特徴量統合工程１６９とを
他のコンピュータで機能させることなどもできる。
【０２２０】
この場合、基準ＦＰ統合特徴量は、他のコンピュータで作成し、評価工程１７１に供給
されることになる。
【０２２１】
こうして、対象ＦＰタイプ２作成工程１５３は、時系列でピークが変化するパターンと
して対象ＦＰタイプ２を作成する。
対象ＦＰ領域分割特徴量作成工程１５５は、前記対象パターン・タイプ２を複数の領域に
分割し各領域に存在するピークの存在率から対象パターン領域分割特徴量を作成するパタ
ーン領域分割特徴量作成工程であるＦＰ領域分割特徴量作成工程を構成する。［多成分薬
剤の評価プログラム］
図９３〜図１０８は、多成分薬剤の評価プログラムに係るフローチャート、図１０９〜
図１１６は、基準データの作成に係るフローチャート、図１１７は、３Ｄクロマトのデー
タ例を示す図表、図１１８は、ピーク情報のデータ例を示す図表、図１１９は、ＦＰのデ
ータ例を示す図表、図１２０は、対象ＦＰの基準ＦＰへの帰属スコア計算結果（判定結果
ファイル）例を示す図表、図１２１は、対象ＦＰと基準ＦＰで対応するピークの照合過程
で作成する２つの中間ファイル（帰属候補ピークスコア表、帰属候補ピーク番号表）例を
示す図表、図１２２は、対象ＦＰと基準ＦＰで対応するピークを特定した結果である照合
結果ファイル例を示す図表、図１２３は、基準群ＦＰのデータ例を示す図表、図１２４は
、基準群ＦＰに帰属した対象ＦＰのピーク特徴量データのファイル例を示す図表、図１２
５は、対象及び基準ＦＰタイプ２のデータ例を示す図表、図１２６は、対象ＦＰ領域分割
特徴量ファイル例を示す図表、図１２７は、対象ＦＰ特徴量統合ファイル例を示す図表、
図１２８は、基準ｔｙｐｅ２群ＦＰ例を示す図表、図１２９は、基準群統合データ例を示
す図表である。
【０２２２】
図９３、図９４は、評価対象薬剤を評価するための処理全体のステップを示すフローチ
ャートであり、システム起動によりスタートし、ＦＰ作成部３のＦＰ作成機能と、対象Ｆ
Ｐピーク帰属部５の対象ＦＰピーク帰属機能と、対象ＦＰピーク特徴量作成部７の対象Ｆ
Ｐピーク特徴量作成機能と、対象ＦＰタイプ２作成部９の対象ＦＰタイプ２作成機能と、
対象ＦＰ領域分割特徴量作成部１１の対象ＦＰ領域分割特徴量作成機能と、対象ＦＰ特徴
量統合部１３の対象ＦＰ特徴量統合機能と、基準ＦＰピーク帰属部１５の基準ＦＰピーク
帰属機能と、基準ＦＰ帰属結果統合部１７の基準ＦＰ帰属結果統合機能と、基準ＦＰピー
ク特徴量作成部１９の基準ＦＰピーク特徴量作成機能と、基準ＦＰタイプ２作成部２１の
基準ＦＰタイプ２作成機能と、基準ＦＰ領域分割特徴量作成部２３の基準ＦＰ領域分割特
徴量作成機能と、基準ＦＰ特徴量統合部２５の基準ＦＰ特徴量統合機能と、評価部２７の
評価機能とをコンピュータに実現させる。
【０２２３】
ＦＰ作成機能は、ステップＳ１で実現される。対象ＦＰピーク帰属機能は、ステップＳ
２、Ｓ３、Ｓ４で実現される。対象ＦＰピーク特徴量作成機能は、ステップＳ５で実現さ
れる。対象ＦＰタイプ２作成機能は、ステップＳ６で実現される。対象ＦＰ領域分割特徴
量作成機能は、ステップＳ７で実現される。対象ＦＰ特徴量統合機能は、ステップＳ８で
実現される。評価機能は、ステップＳ９、Ｓ１０で実現される。
【０２２４】
ステップＳ１は、３Ｄクロマト及び特定の検出波長におけるピーク情報を入力データと
して「ＦＰ作成処理」が実行される。
【０２２５】
３Ｄクロマトは、評価対象薬剤をＨＰＬＣで分析することにより得られるデータであり
、図１１７の３Ｄクロマトのデータ例１８３で示すように、リテンション・タイム、検出
波長、ピーク（シグナル強度）の三次元情報で構成されたデータである。ピーク情報は、
同ＨＰＬＣ分析により得られる特定波長におけるクロマト・データをＨＰＬＣデータ解析
ツール（例えば、ＣｈｅｍＳｔａｔｉｏｎ等）で処理することで得られるデータであり、
図１１８のピーク情報例１８５で示すように、ピークとして検出された全ピークの極大値
及び面積値とその時点のリテンション・タイム等で構成されたデータである。
【０２２６】
ステップＳ１では、コンピュータのＦＰ作成部３の対象ＦＰ作成部２９（図１）が機能
し、３Ｄクロマト及びピーク情報から前記対象ＦＰ４３（図２）を作成し、そのデータを
ファイルとして出力する。この対象ＦＰ４３は、図１１９のＦＰのデータ例１８７で示す
ように、リテンション・タイムとピーク高さとピーク高さごとのＵＶスペクトルで構成さ
れたデータである。
【０２２７】
ステップＳ２は、ステップＳ１で出力した対象ＦＰ及び全基準ＦＰを入力として、「対
象ＦＰ帰属処理１」が実行される。
【０２２８】
ステップＳ２では、コンピュータの基準ＦＰ選定部３３が機能し、全基準ＦＰに対して
対象ＦＰ４３とのリテンション・タイム・出現パターンの一致度を算出し、対象ＦＰ４３
の帰属に適した基準ＦＰを選定する。
【０２２９】
基準ＦＰは、正常品と定めた薬剤の３Ｄクロマトとピーク情報から前記ステップＳ１と
同様の処理により作成されたＦＰである。なお、正常品は、安全性、有効性が確認された
薬剤（基準漢方薬）と定義され、製品ロットの異なる複数の薬剤が該当する。 基準ＦＰ
も、図１１９のＦＰのデータ例１８７と同様に構成されたデータである。
【０２３０】
ステップＳ３は、対象ＦＰ４３とステップＳ２で選定した基準ＦＰを入力とし、「対象
ＦＰ帰属処理２」が実行される。
【０２３１】
ステップＳ３では、コンピュータのピーク・パターン作成部３５（図１）及びピーク帰
属部３７（図１）が機能する。この機能により、対象ＦＰ４３とステップＳ２で選定した
基準ＦＰの全ピークで、図２３〜図６１のように網羅的にピーク・パターンを作成し、次
にそれらピーク・パターンの一致度（図６３または図６４のＰ＿Ｓｉｍ）を算出する。ま
た、対象ＦＰと基準ＦＰのピーク間でＵＶスペクトルの一致度（図６６のＵＶ＿Ｓｉｍ）
を算出する。さらにこれら２つの一致度から帰属候補ピークの一致度（図６７のＳＣＯＲ
Ｅ）を算出する。その算出結果は、図１２０の判定結果ファイル例１８９と同様なファイ
ルに出力される。
【０２３２】
ステップＳ４は、ステップＳ３で出力した判定結果ファイル１８９を入力とし、「対象
ＦＰ帰属処理３」が実行される。
【０２３３】
ステップＳ４では、コンピュータのピーク帰属部３７が機能し、対象ＦＰ４３と基準Ｆ
Ｐ間で、帰属候補ピークの一致度（ＳＣＯＲＥ）をもとに対象ＦＰの各ピークに対応する
基準ＦＰのピークを特定する。その結果は、図１２２の照合結果ファイル例１９５と同様
な照合結果ファイルに出力される。
【０２３４】
ステップＳ５は、ステップＳ４で出力した照合結果ファイルと基準群ＦＰ１９７とを入
力とし、「対象ＦＰ帰属処理４」が実行される。
【０２３５】
基準群ＦＰ１９７は、全基準ＦＰから前記ステップＳ２からステップＳ４と同様の処理
により作成された全基準ＦＰ間のピーク対応データである。
【０２３６】
ステップＳ５では、コンピュータの対象ＦＰピーク特徴量作成部７が機能し、対象ＦＰ
４３の照合結果ファイルをもとに、図６８、図６９のように、対象ＦＰ４３の各ピークを
基準群ＦＰ１９７のピークに帰属する。その結果は、図１２４のピーク・データ特徴量の
ファイル例１９９と同様なファイルに出力される。
【０２３７】
ステップＳ６は、ステップＳ５で出力したピーク・データ特徴量ファイルと対象ＦＰと
を入力とし、「ＦＰ_type２の作成」の処理が実行される。
【０２３８】
ステップＳ６では、コンピュータの対象ＦＰタイプ２作成部９が機能し、対象ＦＰピー
ク特徴量作成部７において特定されたピーク４７を元の対象ＦＰ４３から除いて残ったピ
ークとそのリテンション・タイムとで構成されるＦＰを対象ＦＰタイプ２（４９）として
作成する。その結果は、ＦＰタイプ２ファイル（図１２５のＦＰタイプ２ファイル例２０
１参照）に出力される。
【０２３９】
ステップＳ７では、「領域分割による対象ＦＰ_type２の特徴量化処理」が実行される
。この処理では、コンピュータの対象ＦＰ領域分割特徴量作成部１１が機能し、図７０の
領域分割により、対象ＦＰ領域分割特徴量が作成される。その結果は、対象ＦＰ領域分割
特徴量ファイル（図１２６の対象ＦＰ領域分割特徴量ファイル例２０３参照）に出力され
る。
【０２４０】
ステップＳ８では、「ピーク・データ特徴量と領域分割特徴量の統合」の処理が実行さ
れる。
【０２４１】
この処理では、コンピュータの対象ＦＰ特徴量統合部１３が機能し、対象ＦＰピーク特
徴量作成部７で作成された対象ＦＰピーク特徴量４７と対象ＦＰ領域分割特徴量作成部１
１で作成された対象ＦＰ領域分割特徴量５１とを統合して対象ＦＰ統合特徴量を作成する
。その結果は、対象ＦＰ特徴量統合ファイル（図１２７の対象ＦＰ特徴量統合ファイル例
２０５参照）に出力される。
【０２４２】
ステップＳ９では、コンピュータの評価部２７が機能し、ステップＳ８で出力した対象
ＦＰ統合特徴量と基準ＦＰ統合特徴量との同等性をＭＴ法により評価し、その評価結果を
図８７〜図９１のようなＭＤ値として出力する（図８７〜図９１）。
【０２４３】
ステップＳ１０は、ステップＳ９で出力したＭＤ値を入力として、「合否判定」が実行
される。
【０２４４】
ステップＳ１０では、コンピュータの評価部２７が機能し、ステップＳ９で出力したＭ
Ｄ値と予め設定した閾値（ＭＤ値の上限値）を比較し、合否を判定する（図２の評価結果
５３）。
［Ｓ１：ＦＰ作成処理（単一波長のみ利用）］
図９５は、図９３ステップＳ１「ＦＰ作成処理」の単一波長のピーク情報を利用した場
合のフローチャートである。
【０２４５】
図９５は、波長を単一波長、例えば２０３ｎｍとして評価対象のＦＰを作成するステッ
プの詳細である。この処理では、３Ｄクロマトと検出波長が２０３ｎｍにおけるピーク情
報から、２０３ｎｍで検出されたピークにおけるリテンション・タイムとピークならびに
それらピークのＵＶスペクトルで構成するＦＰを作成する。
【０２４６】
ステップＳ１０１では、「ピーク情報を読み込む」の処理が実行される。この処理では
、ＦＰの作成に必要な２つのデータのうちの1つ目としてピーク情報が読み込まれ、ステ
ップＳ１０２へ移行する。
【０２４７】
ステップＳ１０２では、「ピークのリテンション・タイム（Ｒ１）と対応するピーク・
データ（Ｐ１）を順番に取得」の処理が実行される。この処理では、ピーク情報から、ピ
ークのリテンション・タイム（Ｒ１）及びピーク・データ（Ｐ１）を１ピークずつ順番に
取得し、ステップＳ１０３へ移行する。
【０２４８】
ステップＳ１０３では、「３Ｄクロマトを読み込む」の処理が実行される。この処理で
は、ＦＰの作成に必要な２つのデータのうちの２つ目として３Ｄクロマトが読み込まれ、
ステップＳ１０４へ移行する。
【０２４９】
ステップＳ１０４では、「ピークのリテンション・タイム（Ｒ２）と対応するＵＶスペ
クトル（Ｕ１）を順番に取得」の処理が実行される。この処理では、３Ｄクロマトから、
リテンション・タイム（Ｒ２）及びＵＶスペクトル（Ｕ１）をＨＰＬＣ分析時のサンプリ
ング・レートごとに取得し、ステップＳ１０５へ移行する。
【０２５０】
ステップＳ１０５では、「｜Ｒ１−Ｒ２｜＜＝閾値？」の判断処理が実行される。この
処理では、ステップＳ１０２及びＳ１０４で読み込まれたＲ１とＲ２が閾値の範囲で対応
したものであるか否かが判断される。対応している（ＹＥＳ）場合は、２つのリテンショ
ン・タイムは同じであり、リテンション・タイムがＲ１のピークのＵＶスペクトルはＵ１
であると判断し、ステップＳ１０６へ移行する。対応していない（ＮＯ）場合は、２つの
リテンション・タイムは同じではなく、リテンション・タイムがＲ１のピークのＵＶスペ
クトルはＵ１ではないと判断し、３Ｄクロマトの次のデータとの比較のため、ステップＳ
１０４へ移行する。なお、この判断処理での閾値は、３Ｄクロマトにおける“サンプリン
グ・レート／２”とする。
【０２５１】
ステップＳ１０６では、「Ｕ１を最大値１で規格化」の処理が実行される。この処理で
は、Ｓ１０５でＲ１のＵＶスペクトルと判断したＵ１を最大値１で規格化し、ステップＳ
１０７へ移行する。
【０２５２】
ステップＳ１０７では、「Ｒ１とＰ１ならびに規格化したＵ１を出力（対象ＦＰ）」の
処理が実行される。この処理では、ピーク情報から取得したＲ１とＰ１ならびにＳ１０６
で規格化したＵ１を対象ＦＰに出力し、ステップＳ１０８へ移行する。
【０２５３】
ステップＳ１０８では、「全ピークの処理終了？」の判断処理が実行される。この処理
では、ピーク情報中の全てのピークに対して処理が行われたか否かが判断され、全ピーク
に対して処理が終了していなければ（ＮＯ）、未処理のピークを処理するため、ステップ
Ｓ１０２へ移行する。Ｓ１０２からＳ１０８までの処理は全ピークの処理が終了するまで
繰り返され、全ピークの処理が終了すると（ＹＥＳ）、ＦＰ作成処理を終了する。
［Ｓ１：ＦＰ作成処理（複数波長利用）］
図９６、図９７は、図９３のステップＳ１「ＦＰ作成処理」において、前記単一波長の
ピーク情報に代え、複数波長のピーク情報を利用した場合のフローチャートである。例え
ば２０３ｎｍを含めて、検出波長軸方向に複数（ｎ個）の波長を選択し、ＦＰを作成する
場合である。
【０２５４】
このＦＰ作成処理は、図９５のような単一波長では３Ｄクロマトで検出されている全ピ
ークを網羅できない場合に複数波長のピーク情報を利用し、３Ｄクロマトの全ピークを網
羅したＦＰを作成するためのものである。
【０２５５】
なお、図９６、図９７は、上記単一波長のみを利用したＦＰ作成処理で波長ごとのＦＰ
をｎ個作成後、それらＦＰから複数波長によるＦＰを作成するステップの詳細である。
【０２５６】
ステップＳ１１０では、「波長ごとにＦＰを作成」の処理が実行される。この処理では
、波長ごとに上記単一波長のみを利用したＦＰ作成処理が行われ、ｎ個のＦＰを作成し、
ステップＳ１１１へ移行する。
【０２５７】
ステップＳ１１１では、「ＦＰをピーク数（降順）でリスト化」の処理が実行される。
この処理では、ｎ個のＦＰをピーク数が多い順にリスト化し、ステップＳ１１２へ移行す
る。
【０２５８】
ステップＳ１１２では、ｎ個のＦＰを順番に処理するためのカウンタの初期化としてｎ
に１が代入（ｎ←１）され、ステップＳ１１３へ移行する。
【０２５９】
ステップＳ１１３では、「リストｎ番目のＦＰを読み込む」の処理が実行される。この
処理では、リストｎ番目のＦＰを読み込み、ステップＳ１１４へ移行する。
【０２６０】
ステップＳ１１４では、「全リテンション・タイム（Ｘ）を取得」の処理が実行される
。この処理では、Ｓ１１３で読み込んだＦＰのリテンション・タイム情報を全て取得し、
ステップＳ１１５へ移行する。
【０２６１】
ステップＳ１１５では、「ｎの更新（ｎ←ｎ＋１）」の処理が実行される。この処理で
は、処理を次のＦＰへ移行するため、ｎの更新としてｎにｎ＋１を代入し、ステップＳ１
１６へ移行する。
【０２６２】
ステップＳ１１６では、「リストｎ番目のＦＰを読み込む」の処理が実行される。この
処理では、リストｎ番目のＦＰを読み込み、ステップＳ１１７へ移行する。
【０２６３】
ステップＳ１１７では、「全リテンション・タイム（Ｙ）を取得」の処理が実行される
。この処理では、Ｓ１１６で読み込んだＦＰのリテンション・タイム情報を全て取得し、
ステップＳ１１８へ移行する。
【０２６４】
ステップＳ１１８では、「ＸとＹを重複なしで統合（Ｚ）」の処理が実行される。この
処理では、Ｓ１１４で取得したリテンション・タイム情報ＸとＳ１１７で取得したリテン
ション・タイム情報Ｙを重複なしで統合した後、Ｚに保存し、ステップＳ１１９へ移行す
る。
【０２６５】
ステップＳ１１９では、「Ｘの更新（Ｘ←Ｚ）」の処理が実行される。この処理では、
Ｘの更新としてＸにＳ１１８で保存したＺを代入し、ステップＳ１２０へ移行する。
【０２６６】
ステップＳ１２０では、「全ＦＰ処理終了？」の判断処理が実行される。この処理では
、Ｓ１１０で作成したｎ個のＦＰ全てが処理されたか否かが判断され、処理済み（ＹＥＳ
）の場合は、ステップＳ１２１へ移行する。未処理のＦＰがある（ＮＯ）の場合は、未処
理のＦＰに対してＳ１１５〜Ｓ１２０の処理を実行するため、Ｓ１１５へ移行する。全Ｆ
Ｐの処理が終了するまでＳ１１５〜Ｓ１２０の処理を繰り返す。
【０２６７】
ステップＳ１２１では、ｎ個のＦＰを再度順番に処理するためのカウンタの初期化とし
てｎに１が代入（ｎ←１）され、ステップＳ１２２へ移行する。
【０２６８】
ステップＳ１２２では、「リストｎ番目のＦＰを読み込む」の処理が実行される。この
処理では、リストｎ番目のＦＰを読み込み、ステップＳ１２３へ移行する。
【０２６９】
ステップＳ１２３では、「各ピークのリテンション・タイム（Ｒ１）、ピーク・データ
（Ｐ１）ならびにＵＶスペクトル（Ｕ１）を順番に取得」の処理が実行される。この処理
では、Ｓ１２２で読み込んだＦＰからリテンション・タイム（Ｒ１）、ピーク・データ（
Ｐ１）ならびにＵＶスペクトル（Ｕ１）を１ピークずつ順番に取得し、ステップＳ１２４
へ移行する。
【０２７０】
ステップＳ１２４では、「Ｘからリテンション・タイム（Ｒ２）を順番に取得」の処理
が実行される。この処理では、全ＦＰのリテンション・タイムが重複なしで保存されてい
るＸから1リテンション・タイム（Ｒ２）ずつ順番に取得し、ステップＳ１２５へ移行す
る。
【０２７１】
ステップＳ１２５では、「Ｒ１＝Ｒ２？」の判断処理が実行される。この処理では、Ｓ
１２３で取得したＲ１とＳ１２４で取得したＲ２が等しいか否かが判断され、等しい（Ｙ
ＥＳ）場合は、ステップＳ１２７へ移行する。等しくない（ＮＯ）場合は、ステップＳ１
２６へ移行する。
【０２７２】
ステップＳ１２６では、「Ｘの全リテンション・タイム比較終了？」の判断処理が実行
される。この処理では、Ｓ１２３で取得したＲ１に対して、Ｘの全リテンション・タイム
との比較が終了しているか否かが判断される。終了している（ＹＥＳ）場合は、リテンシ
ョン・タイムがＲ１のピークは処理済みであると判断し、次のピークへ処理を移行するた
め、ステップＳ１２３へ移行する。終了していない（ＮＯ）場合は、Ｘの次のリテンショ
ン。タイムに移行するため、ステップＳ１２４へ移行する。
【０２７３】
ステップＳ１２７では、「Ｒ１に（ｎ−１）×分析時間（Ｔ）を加算（Ｒ１←Ｒ１＋（
ｎ−１）×Ｔ）」の処理が実行される。この処理では、一番ピーク数の多いリスト１番目
のＦＰに存在するピークのリテンション・タイムはそのままで、リスト１番目のＦＰには
存在せず、リスト２番目のＦＰに存在するピークのリテンション・タイムはＲ１に分析時
間
（Ｔ）が加算され、さらに、リスト１〜ｎ−１番目のＦＰには存在せず、リストｎ番目の
ＦＰに存在するピークのリテンション・タイムはＲ１に（ｎ−１）×Ｔが加算され、ステ
ップＳ１２８へ移行する。
【０２７４】
ステップＳ１２８では、「Ｒ１、Ｐ１ならびにＵ１を出力（対象ＦＰ）」の処理が実行
される。この処理では、Ｓ１２７で処理したＲ１とＳ１２３で取得したＰ１ならびにＵ１
を対象ＦＰに出力し、ステップＳ１２９へ移行する。
【０２７５】
ステップＳ１２９では、「ＸからＲ２を削除」の処理が実行される。この処理では、リ
テンション・タイムがＲ１（＝Ｒ２）における処理がＳ１２７、Ｓ１２８で終了したため
、Ｘから処理済みリテンション・タイム（Ｒ２）を削除し、Ｓ１３０へ移行する。
【０２７６】
ステップＳ１３０では、「全ピーク処理終了？」の判断処理が実行される。この処理で
は、リストｎ番目のＦＰの全ピークに対して処理が終了しているか否かが判断され、処理
済み（ＹＥＳ）の場合は、リストｎ番目のＦＰにおけるＦＰ作成処理を終了し、ステップ
Ｓ１３１へ移行する。未処理のピークがある（ＮＯ）場合は、未処理のピークを処理する
ため、ステップＳ１２３へ移行する。全ピークの処理が終了するまでＳ１２３〜Ｓ１３０
の処理を繰り返す。
【０２７７】
ステップＳ１３１では、「ｎの更新（ｎ←ｎ＋１）」の処理が実行される。この処理で
は、処理を次のＦＰへ移行するため、ｎの更新としてｎにｎ＋１を代入し、ステップＳ１
３２へ移行する。
【０２７８】
ステップＳ１３２では、「全ＦＰ処理終了？」の判断処理が実行される。この処理では
、Ｓ１１０で作成したｎ個のＦＰ全てが処理されたか否かが判断され、処理済み（ＹＥＳ
）の場合は、ＦＰ作成処理を終了する。未処理のＦＰがある（ＮＯ）場合は、未処理のＦ
Ｐに対してＳ１２２〜Ｓ１３２の処理を実行するため、Ｓ１２２へ移行する。全ＦＰの処
理が終了するまでＳ１２２〜Ｓ１３２の処理を繰り返す。
［Ｓ２：対象ＦＰ帰属処理１］
図９８は、図９３ステップＳ２の「対象ＦＰ帰属処理１」の詳細を示すフローチャート
である。この処理は、帰属の前処理であり、正常品とされた複数の基準ＦＰから対象ＦＰ
４３の帰属に適した基準ＦＰを選定する。
【０２７９】
ステップＳ２０１では、「対象ＦＰを読み込む」の処理が実行される。この処理では、
帰属対象のＦＰを読み込み、ステップＳ２０２へ移行する。
【０２８０】
ステップＳ２０２では、「全リテンション・タイム（Ｒ１）を取得」の処理が実行され
る。この処理では、Ｓ２０１で読み込んだ対象ＦＰのリテンション・タイム情報を全て取
得し、ステップＳ２０３へ移行する。
【０２８１】
ステップＳ２０３では、「全基準ＦＰのファイル名をリスト化」の処理が実行される。
この処理では、後で全基準ＦＰを順番に処理するために予め全基準ＦＰのファイル名をリ
スト化し、ステップＳ２０４へ移行する。
【０２８２】
ステップＳ２０４では、全基準ＦＰを順番に処理するためのカウンタの初期値としてｎ
に１を代入（ｎ←１）し、ステップＳ２０５へ移行する。
【０２８３】
ステップＳ２０５では、「リストｎ番目の基準ＦＰ（基準ＦＰ_ｎ）を読み込む」の処理
が実行される。この処理では、Ｓ２０３でリスト化した全基準ＦＰのファイル名リストの
ｎ番目のＦＰを読み込み、ステップＳ２０６へ移行する。
【０２８４】
ステップＳ２０６では、「全リテンション・タイム（Ｒ２）を取得」の処理が実行され
る。この処理では、Ｓ２０５で読み込んだ基準ＦＰのリテンション・タイム情報を全て取
得し、ステップＳ２０７へ移行する。
【０２８５】
ステップＳ２０７では、「Ｒ１とＲ２のリテンション・タイム・出現パターンの一致度
を算出（ＲＰ_ｎ＿ｍｉｎ）」の処理が実行される。この処理では、Ｓ２０２で取得した対
象ＦＰのリテンション・タイムとＳ２０６で取得した基準ＦＰのリテンション・タイムか
らＲＰ_ｎ＿ｍｉｎを算出し、ステップＳ２０８へ移行する。なお、ＲＰ_ｎ＿ｍｉｎの詳細
な計算フローは、図１０３のサブルーチン１により別途説明する。
【０２８６】
ステップＳ２０８では、「ＲＰ_ｎ＿ｍｉｎの保存（ＲＰ_ａｌｌ＿ｍｉｎ」の処理が実行
される。この処理では、Ｓ２０７で算出したＲＰ_ｎ＿ｍｉｎをＲＰ_ａｌｌ＿ｍｉｎに保存
し、ステップＳ２０９へ移行する。
【０２８７】
ステップＳ２０９では、「ｎの更新（ｎ←ｎ＋１）」の処理が実行される。この処理で
は、処理を次のＦＰへ移行するためｎの更新としてｎにｎ＋１を代入し、ステップＳ２１
０へ移行する。
【０２８８】
ステップＳ２１０では、「全基準ＦＰ処理終了？」の判断処理が実行される。この処理
では、基準ＦＰ全てが処理されたか否かが判断され、処理済み（ＹＥＳ）の場合は、ステ
ップＳ２１１へ移行する。未処理の基準ＦＰがある（ＮＯ）場合は、未処理のＦＰに対し
てＳ２０５〜Ｓ２１０の処理を実行するため、Ｓ２０５へ移行する。全基準ＦＰの処理が
終了するまでＳ２０５〜Ｓ２１０の処理を繰り返す。
【０２８９】
ステップＳ２１１では、「ＲＰ_ａｌｌ＿ｍｉｎから一致度が最小となる基準ＦＰを選定
」の処理が実行される。この処理では、全基準ＦＰに対して算出したＲＰ_１＿ｍｉｎから
ＲＰ_ｎ＿ｍｉｎを比較し、対象ＦＰとのリテンション・タイム・出現パターンの一致度が
最小となる基準ＦＰを選定し、対象ＦＰ帰属処理１を終了する。
［Ｓ３：対象ＦＰ帰属処理２］
図９９は、図９３ステップＳ３の「対象ＦＰ帰属処理２」の詳細を示すフローチャート
である。この処理は、帰属の本処理であり、対象ＦＰ４３とステップＳ２で選定した基準
ＦＰとの間で、前記のようなピーク・パターン及びＵＶスペクトルの一致度から各帰属候
補ピークの一致度（ＳＣＯＲＥ）を算出する。
【０２９０】
ステップＳ３０１では、「対象ＦＰを読み込む」の処理が実行される。この処理では、
帰属対象のＦＰを読み込み、ステップＳ３０２へ移行する。
【０２９１】
ステップＳ３０２では、「帰属対象ピークのリテンション・タイム（Ｒ１）とピーク・
データ（Ｐ１）ならびにＵＶスペクトル（Ｕ１）を順番に取得」の処理が実行される。こ
の処理では、Ｓ３０１で読み込んだ対象ＦＰの各ピークを順番に帰属対象ピークとし、Ｒ
１とＰ１ならびにＵ１を取得し、ステップＳ３０３へ移行する。
【０２９２】
ステップＳ３０３では、「基準ＦＰを読み込む」の処理が実行される。この処理では、
図９８の［対象ＦＰ帰属処理１］で選定された基準ＦＰを読み込み、ステップＳ３０４へ
移行する。
【０２９３】
ステップＳ３０４では、「基準ＦＰのピークのリテンション・タイム（Ｒ２）とピーク
・データ（Ｐ２）ならびにＵＶスペクトル（Ｕ２）を順番に取得」の処理が実行される。
この処理では、Ｓ３０３で読み込んだ基準ＦＰからＲ２とＰ２ならびにＵ２を１ピークず
つ取得し、ステップＳ３０５へ移行する。
【０２９４】
ステップＳ３０５では、「｜Ｒ１−（Ｒ２＋ｄ）｜＜閾値？」の判断処理が実行される
。この処理では、ステップＳ３０２及びＳ３０４で読み込まれたＲ１とＲ２が閾値の範囲
内で対応したものであるか否かが判断される。対応している（ＹＥＳ）場合は、リテンシ
ョン・タイムがＲ２のピークはリテンション・タイムがＲ１のピークの帰属候補ピークで
あると判断し、帰属候補ピークの一致度（ＳＣＯＲＥ）を算出するためステップＳ３０６
へ移行する。対応していない（ＮＯ）場合は、リテンション・タイムがＲ２のピークとリ
テンション・タイムがＲ１のピークではリテンション・タイムが違いすぎるため、帰属候
補ピークにはならないと判断し、ステップＳ３０９へ移行する。なお、この判断処理での
ｄは、対象ＦＰと基準ＦＰのピークのリテンション・タイムを補正する値で、初期値は０
とし、処理を進める中で随時帰属したピーク間のリテンション・タイムの差を求め、その
値でｄを更新する。また、閾値は、帰属候補ピークとすべきかどうかを判断するためのリ
テンション・タイムの許容幅である。
【０２９５】
ステップＳ３０６では、「ＵＶスペクトルの一致度を算出（ＵＶ＿Ｓｉｍ）」の処理が
実行される。この処理では、Ｓ３０２で取得した帰属対象ピークのＵ１とＳ３０４で取得
した帰属候補ピークのＵ２から、ＵＶ＿Ｓｉｍを算出し、ステップＳ３０７へ移行する。
なお、ＵＶ＿Ｓｉｍの詳細な計算フローは、図８６のサブルーチン２に別途記載する。
【０２９６】
ステップＳ３０７では、「ピーク・パターンの一致度を算出（Ｐ＿Ｓｉｍ＿ｍｉｎ）」
の処理が実行される。この処理では、Ｓ３０２で取得した帰属対象ピークのＲ１及びＰ１
とＳ３０４で取得した帰属候補ピークのＲ２及びＰ２から、これらピークに対して網羅的
にピーク・パターンを作成する。これらピーク・パターンのＰ＿Ｓｉｍ＿ｍｉｎを算出し
、ステップＳ３０８へ移行する。なお、Ｐ＿Ｓｉｍ＿ｍｉｎの詳細な計算フローは、図８
７のサブルーチン３に別途記載する。
【０２９７】
ステップＳ３０８では、「帰属候補ピークの一致度を算出（ＳＣＯＲＥ）」の処理が実
行される。この処理では、Ｓ３０６で算出したＵＶ＿ＳｉｍとＳ３０７で算出したＰ＿Ｓ
ｉｍ＿ｍｉｎから、帰属対象ピークと帰属候補ピークのＳＣＯＲＥを、
ＳＣＯＲＥ ＝ ＵＶ＿Ｓｉｍ × Ｐ＿Ｓｉｍ＿ｍｉｎ
として算出し、ステップＳ３１０へ移行する。
【０２９８】
ステップＳ３０９では、「ＳＣＯＲＥに８８８８８８を代入（ＳＣＯＲＥ←８８８８８
８）」の処理が実行される。この処理では、帰属対象ピークの帰属候補ピークに該当しな
いピークのＳＣＯＲＥを８８８８８８とし、ステップＳ３１０へ移行する。
【０２９９】
ステップＳ３１０では、「ＳＣＯＲＥの保存（ＳＣＯＲＥ＿ａｌｌ）」の処理が実行さ
れる。この処理では、Ｓ３０８あるいはＳ３０９で得られたＳＣＯＲＥをＳＣＯＲＥ＿ａ
ｌｌに保存し、ステップＳ３１１へ移行する。
【０３００】
ステップＳ３１１では、「基準全ピークの処理終了？」の判断処理が実行される。この
処理では、基準ＦＰの全ピークが処理されたか否かが判断され、処理済み（ＹＥＳ）の場
合は、ステップＳ３１２へ移行する。未処理のピークがある（ＮＯ）場合は、未処理のピ
ークに対してＳ３０４〜Ｓ３１１の処理を実行するため、Ｓ３０４へ移行する。全ピーク
の処理が終了するまでＳ３０４〜Ｓ３１１の処理を繰り返す。
【０３０１】
ステップＳ３１２では、「判定結果ファイルにＳＣＯＲＥ＿ａｌｌを出力し、ＳＣＯＲ
Ｅ＿ａｌｌを初期化（空にする）」の処理が実行される。この処理では、判定結果ファイ
ルにＳＣＯＲＥ＿ａｌｌを出力後、ＳＣＯＲＥ＿ａｌｌを初期化（空にする）し、ステッ
プＳ３１３へ移行する。
【０３０２】
ステップＳ３１３では、「対象全ピークの処理終了？」の判断処理が実行される。この
処理では、対象ＦＰの全ピークが処理されたか否かが判断され、処理済み（ＹＥＳ）の場
合は、対象ＦＰ帰属処理２を終了する。未処理のピークがある（ＮＯ）場合は、未処理の
ピークに対してＳ３０２〜Ｓ３１３の処理を実行するため、Ｓ３０２へ移行する。全ピー
クの処理が終了するまでＳ３０２〜Ｓ３１３の処理を繰り返す。
【０３０３】
図１２０で出力した判定結果ファイル例１８９を示す。
［Ｓ４：対象ＦＰ帰属処理３］
図１００は、図９３ステップＳ４の「対象ＦＰ帰属処理３」の詳細を示すフローチャー
トである。この処理は、帰属の後処理であり、前記のように算出した帰属候補ピークの一
致度（ＳＣＯＲＥ）から対象ＦＰの各ピークに対応する基準ＦＰのピークを特定する。
【０３０４】
ステップＳ４０１では、「判定結果ファイルを読み込む」の処理が実行される。この処
理では、図８１の「対象ＦＰ帰属処理２」で作成した判定結果ファイルを読み込み、ステ
ップＳ４０２へ移行する。
【０３０５】
ステップＳ４０２では、「「ＳＣＯＲＥ＜閾値」の条件を満たしたデータで帰属候補ピ
ークスコア表を作成」の処理が実行される。この処理では、判定結果ファイルのＳＣＯＲ
Ｅをもとに帰属候補スコア表（図１２１上図の帰属候補スコア表１９１参照）を作成し、
ステップＳ４０３へ移行する。この帰属候補ピークスコア表は、基準ＦＰのピークごとに
、対象ＦＰ全ピークに対して算出されたＳＣＯＲＥから、閾値より小さいＳＣＯＲＥだけ
を昇順に並べた表である。ちなにみ、このＳＣＯＲＥは値が小さいほど帰属すべきピーク
の可能性が高い。なお、閾値は、帰属候補とすべきかどうかを判断するためのＳＣＯＲＥ
の上限値である。
【０３０６】
ステップＳ４０３では、「帰属候補ピーク番号表を作成」の処理が実行される。この処
理では、帰属候補ピークスコア表をもとに帰属候補ピーク番号表（図１２１下図の帰属候
補ピーク番号表１９３参照）を作成し、ステップＳ４０４へ移行する。この帰属候補ピー
ク番号表は、帰属候補ピークスコア表の各スコアをそのスコアに対応する対象ＦＰのピー
ク番号に置き換えた表である。このことから、この表は、基準ＦＰのピークごとに対応さ
せるべき対象ＦＰのピーク番号が順番に並んだ表となっている。
【０３０７】
ステップＳ４０４では、「帰属すべき対象ＦＰのピーク番号を取得」の処理が実行され
る。この処理では、Ｓ４０３で作成した帰属候補ピーク番号表から、基準ＦＰのピークご
とに一番上位に位置する対象ＦＰのピーク番号を取得し、ステップＳ４０５へ移行する。
【０３０８】
ステップＳ４０５では、「取得したピーク番号が降順（重複なし）に並んでいる？」の
判断処理が実行される。この処理では、Ｓ４０４で取得した対象ＦＰのピーク番号が重複
なしで降順に並んでいるか否かが判断される。並んでいる（ＹＥＳ）場合は、基準ＦＰの
各ピークに対応する対象ＦＰのピークが確定できたと判断し、ステップＳ４０８へ移行す
る。並んでいない（ＮＯ）場合は、問題のあった基準ＦＰのピークに帰属すべき対象ＦＰ
のピークを見直すため、ステップＳ４０６へ移行する。
【０３０９】
ステップＳ４０６では、「問題のあったピーク間でＳＣＯＲＥを比較し、帰属候補ピー
ク番号表を更新」の処理が実行される。この処理では、問題のあった対象ＦＰのピーク番
号に対応するＳＣＯＲＥを帰属候補スコア表で比較し、ＳＣＯＲＥが大きい方のピーク番
号を２番目に位置するピーク番号に置き換えた帰属候補ピーク番号表に更新し、ステップ
Ｓ４０７へ移行する。
【０３１０】
ステップＳ４０７では、「帰属候補ピークスコア表を更新」の処理が実行される。この
処理では、Ｓ４０６での帰属候補ピーク番号表の更新内容に沿って、帰属候補ピークスコ
ア表を更新し、ステップＳ４０４へ移行する。対象ＦＰのピーク番号に問題（重複あり、
降順に並んでいない）がなくなるまで、Ｓ４０４からＳ４０７の処理を繰り返す。
【０３１１】
ステップＳ４０８では、「帰属結果を保存（ＴＥＭＰ）」の処理が実行される。この処
理では、基準ＦＰの全ピークのピーク番号、リテンション・タイム、及びピークとこれら
ピークに対応するピークとして特定した対象ＦＰのピーク・データをＴＥＭＰに保存し、
ステップＳ４０９へ移行する。
【０３１２】
ステップＳ４０９では、「ＴＥＭＰに対象ＦＰの全ピークが入っている？」の判断処理
が実行される。この処理では、Ｓ４０８で保存したＴＥＭＰ中に対象ＦＰの全ピークのピ
ーク・データが入っているか否かが判断される。全て入っている（ＹＥＳ）場合は、対象
ＦＰの全てのピークで処理が終了したと判断し、Ｓ４１２へ移行する。入っていないピー
クがある（ＮＯ）場合は、入っていないピークのピーク・データをＴＥＭＰに追加するた
め、ステップＳ４１０へ移行する。
【０３１３】
ステップＳ４１０では、「ＴＥＭＰに入っていない対象ＦＰのピークのリテンション・
タイムを補正」の処理が実行される。この処理では、ＴＥＭＰに入っていない対象ＦＰの
ピーク（補正が必要な対象ＦＰのピーク）のリテンション・タイムは、
補正値 ＝ ｋ１＋（ｋ２−ｋ１）＊（ｔ０−ｔ１）／（ｔ２−ｔ１）
ｋ１：補正が必要な対象ＦＰのピーク近傍で帰属された２つの基準ＦＰ側のピ
ークのうちのリテンション・タイムが小さいピークのリテンション・タイム
ｋ２：補正が必要な対象ＦＰのピーク近傍で帰属された２つの基準ＦＰ側のピ
ークのうちのリテンション・タイムが大きいピークのリテンション・タイム
ｔ０：補正が必要な対象ＦＰのピークのリテンション・タイム
ｔ１：補正が必要な対象ＦＰのピーク近傍で帰属された２つの対象ＦＰ側のピ
ークのうちのリテンション・タイムが小さいピークのリテンション・タイム
ｔ２：補正が必要な対象ＦＰのピーク近傍で帰属された２つの対象ＦＰ側のピ
ークのうちのリテンション・タイムが大きいピークのリテンション・タイム
として基準ＦＰにおけるリテンション・タイムに補正し、ステップＳ４１１へ移行する。
。
【０３１４】
ステップＳ４１１では、「補正したリテンション・タイムとそのピークのピーク・デー
タをＴＥＭＰに追加し、ＴＥＭＰを更新」の処理が実行される。この処理では、Ｓ４１０
で補正したＴＥＭＰに入っていない対象ＦＰのピークのリテンション・タイムとＴＥＭＰ
中の基準ＦＰのリテンション・タイムとを比較し、ＴＥＭＰ中の妥当な位置にＴＥＭＰに
入っていない対象ＦＰのピークの補正したリテンション・タイムならびにピーク・データ
を追加し、ＴＥＭＰを更新し、Ｓ４０９へ移行する。対象ＦＰの全ピークが追加されるま
で、Ｓ４０９からＳ４１１の処理を繰り返す。
【０３１５】
ステップＳ４１２では、「ＴＥＭＰを照合結果ファイルに出力」の処理が実行される。
この処理では、基準ＦＰの全ピークと対象ＦＰの全ピークの対応関係を特定したＴＥＭＰ
を照合結果ファイルとして出力し、対象ＦＰ帰属処理３を終了する。
［Ｓ５：対象ＦＰ帰属処理４］
図１０１、図１０２は、図９３ステップＳ５の「対象ＦＰ帰属処理４」の詳細を示すフ
ローチャートである。この処理は、帰属の最終処理であり、図９３ステップＳ４で作成し
た照合結果ファイル（図１２２の照合結果ファイル例１９５参照）をもとに対象ＦＰの各
ピークを基準群ＦＰ（図１２３の基準群ＦＰのデータ例１９７参照）のピークに帰属する
。
【０３１６】
なお、基準群ＦＰ１９７は、前記のように全基準ＦＰ間でピークの対応関係を特定した
ＦＰであり、図１２３の基準群ＦＰのデータ例１９７のように、基準群ＦＰピーク番号と
基準群リテンション・タイムとピーク高さとで構成されたデータである。図２の基準群Ｆ
Ｐ４５で示すように各ピークは平均値（黒点）±標準偏差（縦分割線）で表すことができ
る。
【０３１７】
ステップＳ５０１では、「照合結果ファイルを読み込む」の処理が実行される。この処
理では、図１００のＳ４１２で出力した照合結果ファイルを読み込み、ステップＳ５０２
へ移行する。
【０３１８】
ステップＳ５０２では、「基準群ＦＰを読み込む」の処理が実行される。この処理では
、対象ＦＰの各ピークの最終的な帰属相手である基準群ＦＰ１９７を読み込み、ステップ
Ｓ５０３へ移行する。
【０３１９】
ステップＳ５０３では、「対象ＦＰと基準群ＦＰを統合し保存（ＴＥＭＰ）」の処理が
実行される。この処理では、照合結果ファイルと基準群ＦＰ１９７で共通に存在する基準
ＦＰのピーク・データをもとに２つのファイルを統合し、その結果をＴＥＭＰとして保存
し、ステップＳ５０４へ移行する。
【０３２０】
ステップＳ５０４では、「基準ＦＰに対応するピークがない対象ＦＰの全ピークのリテ
ンション・タイムを補正」の処理が実行される。この処理では、照合結果ファイルで基準
ＦＰに対応するピークがない対象ＦＰの全ピークのリテンション・タイムをＳ５０３で保
存したＴＥＭＰのリテンション・タイムに補正し、ステップＳ５０５へ移行する。なお、
リテンション・タイムの補正は、前記ステップＳ４の「対象ＦＰ帰属処理３」のステップ
Ｓ４１０と同様の方法で補正する。
【０３２１】
ステップＳ５０５では、「補正したリテンション・タイム（Ｒ１、Ｒ３）と対応するピ
ーク・データ（Ｐ１）を順番に取得」の処理が実行される。この処理では、Ｓ５０４で補
正したリテンション・タイムをＲ１ならびにＲ３として、対応するピークのピーク・デー
タをＰ１として順番に取得し、ステップＳ５０６へ移行する。
【０３２２】
ステップＳ５０６では、「ＴＥＭＰから対象ＦＰの帰属候補ピークのリテンション・タ
イム（Ｒ２）と対応するピーク・データ（Ｐ２）を順番に取得」の処理が実行される。こ
の処理では、Ｓ５０３で保存したＴＥＭＰから対象ＦＰのピークが帰属されていないリテ
ンション・タイムをＲ２として、対応するピーク・データをＰ２として順番に取得し、ス
テップＳ５０７へ移行する。
【０３２３】
ステップＳ５０７では、「｜Ｒ１−Ｒ２｜＜閾値１？」の判断処理が実行される。この
処理では、Ｓ５０５とＳ５０６で取得したＲ１とＲ２の差が閾値１より小さいか否かが判
断される。小さい（ＹＥＳ）場合は、対象ＦＰのリテンション・タイムがＲ１のピークと
基準ＦＰのリテンション・タイムがＲ２のピークが対応する可能性があると判断し、ステ
ップＳ５０８へ移行する。Ｒ１とＲ２の差が閾値１以上（ＮＯ）の場合は、対応する可能
性なしと判断し、ステップＳ５１２へ移行する。
【０３２４】
ステップＳ５０８では、「Ｒ１、Ｒ２に対応するＵＶスペクトルを取得（Ｕ１、Ｕ２）
」の処理が実行される。この処理では、Ｓ５０７で対応する可能性ありと判断されたリテ
ンション・タイムがＲ１とＲ２のピークに対応するＵＶスペクトルをそれぞれのＦＰから
取得し、ステップＳ５０９へ移行する。
【０３２５】
ステップＳ５０９では、「ＵＶスペクトルの一致度を算出（ＵＶ＿Ｓｉｍ）」の処理が
実行される。この処理では、Ｓ５０８で取得したＵＶスペクトルＵ１及びＵ２からステッ
プＳ３の「対象ＦＰ帰属処理２」のステップＳ３０６と同様な方法でＵＶ＿Ｓｉｍを算出
し、ステップＳ５１０へ移行する。なお、ＵＶ＿Ｓｉｍの詳細な計算フローは図１０４の
サブルーチン２で別途説明する。
【０３２６】
ステップＳ５１０では、「ＵＶ＿Ｓｉｍ＜閾値２？」の判断処理が実行される。この処
理では、Ｓ５０９で算出したＵＶ＿Ｓｉｍが閾値２より小さいか否かが判断される。小さ
い（ＹＥＳ）場合は、ＵＶスペクトルがＵ１のピークとＵ２のピークが対応していると判
断し、ステップＳ５１１へ移行する。ＵＶ＿Ｓｉｍが閾値２以上（ＮＯ）の場合は、対応
していないと判断し、ステップＳ５０７へ移行する。
【０３２７】
ステップＳ５１１では、「Ｒ３ ← Ｒ２、閾値２ ← ＵＶ＿Ｓｉｍ」の処理が実行
される。この処理では、Ｓ５１０で対応すると判断したリテンション・タイムがＲ３（つ
まり、Ｒ１）を対応相手のリテンション・タイムであるＲ２に更新したのち、閾値２をＵ
Ｖ＿Ｓｉｍの値に更新し、Ｓ５０７に移行する。
【０３２８】
ステップＳ５１２では、「全ての帰属候補ピークのリテンション・タイムを比較終了？
」の判断処理が実行される。この処理では、Ｒ１と全ての帰属候補ピークのリテンション
・タイムの比較が終了したか否かが判断され、終了している（ＹＥＳ）場合は、ステップ
Ｓ５１３へ移行する。終了していない（ＮＯ）場合は、ステップＳ５０７へ移行する。
【０３２９】
ステップＳ５１３では、「Ｒ１、Ｒ３とＰ１ならびに閾値２を保存（ＴＥＭＰ２）」の
処理が実行される。この処理では、Ｓ５１０で対応すると判断したリテンション・タイム
（Ｒ１）と対応相手のリテンション・タイム（Ｒ２）に更新したＲ３と対応するピーク（
Ｐ１）ならびに現時点の閾値２を保存（ＴＥＭＰ２）し、Ｓ５０７に移行する。
【０３３０】
ステップＳ５１４では、「全ての非対応ピークのリテンション・タイムを比較終了？」
の判断処理が実行される。この処理では、全ての非対応ピークのリテンション・タイムで
帰属候補ピークのリテンション・タイムとの比較が終了したか否かが判断される。終了し
ている（ＹＥＳ）場合は、全ての非対応ピークの帰属処理が終了したと判断し、ステップ
Ｓ５１６へ移行する。終了していない（ＮＯ）場合は、未処理の非対応ピークが残ってい
ると判断し、ステップＳ５１５へ移行する。
【０３３１】
ステップＳ５１５では、「閾値２ ← 初期値」の処理が実行される。この処理では、
Ｓ５１１でＵＶ＿Ｓｉｍに更新されている閾値２を初期値に戻し、ステップＳ５０５へ移
行する。
【０３３２】
ステップＳ５１６では、「ＴＥＭＰ２にＲ３の値が同じピークが存在する？」の判断処
理が実行される。この処理では、ＴＥＭＰ中の同じピークに複数の非対応ピークが帰属さ
れているか否かが判断される。同じピークに帰属された非対応ピークが存在する（ＹＥＳ
）場合は、ステップＳ５１７へ移行する。存在しない（ＮＯ）場合は、ステップＳ５１８
へ移行する。
【０３３３】
ステップＳ５１７では、「Ｒ３の値が同じピークの閾値２を比較し、値が大きいピーク
のＲ３を元の値（Ｒ１）に戻す」の処理が実行される。この処理では、ＴＥＭＰ２中のＲ
３の値が同じピークの閾値２を比較し、値が大きいピークのＲ３の値を元の値（つまり、
Ｒ１）に戻し、ステップＳ５１８へ移行する。
【０３３４】
ステップＳ５１８では、「ＴＥＭＰにＴＥＭＰ２のピークを追加（ＴＥＭＰのリテンシ
ョン・タイムとＲ３が一致したピークのみ）」の処理が実行される。この処理では、ＴＥ
ＭＰのリテンション・タイムとＲ３が一致したピークのみＴＥＭＰにＲ３に対応するピー
クを追加し、ステップＳ５１９へ移行する。Ｒ３がＴＥＭＰのリテンション・タイムと一
致しないピークは、基準群ＦＰに帰属相手となるピークが存在しないため、追加しない。
【０３３５】
ステップＳ５１９では、「ＴＥＭＰ中の対象ＦＰのピークを出力（ピーク特徴量ファイ
ル）」の処理が実行される。この処理では、基準群ＦＰ１９７に帰属された対象ＦＰのピ
ーク・データをピーク・データ特徴量ファイルとして出力し、対象ＦＰ帰属処理４を終了
する。
【０３３６】
図１２４に前記のように出力するピーク・データ特徴量のファイル例１９９を示す。
【０３３７】
[サブルーチン１]
図１０３は、図９８の「基準ＦＰ選定処理」における「サブルーチン１」の詳細を示す
フローチャートである。この処理は、ＦＰ間（例えば、対象ＦＰと基準ＦＰ）のリテンシ
ョン・タイム・出現パターンの一致度を計算する。
【０３３８】
ステップS１００１では、「ｘ←Ｒ１、ｙ←Ｒ２」の処理が実行される。この処理では
、図９８のＳ２０２とＳ２０６で取得したＲ１及びＲ２をそれぞれｘとｙに代入し、ステ
ップＳ１００２へ移行する。
【０３３９】
ステップＳ１００２では、「ｘ、ｙのデータ数を取得（ａ、ｂ）」の処理が実行される
。この処理では、ｘ、ｙのデータ数をそれぞれａ、ｂとして取得し、ステップＳ１００３
へ移行する。
【０３４０】
ステップＳ１００３では、ｘのリテンション・タイムを順番に呼び出すためのカウンタ
の初期値としてｉに１を代入（ｉ←１）し、ステップＳ１００４へ移行する。
【０３４１】
ステップＳ１００４では、「ｘｉ番目のリテンション・タイムからの全距離を取得（ｆ
）」の処理が実行される。この処理では、ｘｉ番目のリテンション・タイムとそれ以降の
全リテンション・タイム間距離をｆとして取得し、ステップＳ１００５へ移行する。
【０３４２】
ステップＳ１００５では、ｙのリテンション・タイムを順番に呼び出すためのカウンタ
の初期値としてｊに１を代入（ｊ←１）し、ステップＳ１００６へ移行する。
【０３４３】
ステップＳ１００６では、「ｙｊ番目のリテンション・タイムからの全距離を取得（ｇ
）」の処理が実行される。この処理では、ｙｊ番目のリテンション・タイムとそれ以降の
全リテンション・タイム間距離をｇとして取得し、ステップＳ１００７へ移行する。
【０３４４】
ステップＳ１００７では「”｜ｆの各リテンション・タイム間距離−ｇの各リテンショ
ン・タイム間距離｜＜閾値”の条件を満たしたデータ数を取得（ｍ）」の処理が実行され
る。この処理では、Ｓ１００４及びＳ１００６で取得したリテンション・タイム間距離ｆ
とｇを総当りで比較し、”｜ｆの各リテンション・タイム間距離−ｇの各リテンション・
タイム間距離｜＜閾値”の条件を満たしたデータ数をｍとして取得し、ステップＳ１００
８へ移行する。
【０３４５】
ステップＳ１００８では、「ｆとｇのリテンション・タイム・出現パターンの一致度を
算出（ＲＰ_ｆｇ）」の処理が実行される。この処理では、Ｓ１００２で取得したａ、ｂと
Ｓ１００７で取得したｍからＲＰ_ｆｇを、
ＲＰ_ｆｇ ＝ （１−（ｍ／（ａ＋ｂ−ｍ）））×（ａ−ｍ＋１）
として算出し、ステップＳ１００９へ移行する。
【０３４６】
ステップＳ１００９では、「ＲＰ_ｆｇを保存（ＲＰ＿ａｌｌ）」の処理が実行される。
この処理では、Ｓ１００８で算出した一致度をＲＰ＿ａｌｌに保存し、ステップＳ１０１
０へ移行する。
【０３４７】
ステップＳ１０１０では、「ｊの更新（ｊ←ｊ＋１）」の処理が実行される。この処理
では、ｙの処理を次のリテンション・タイムへ移行するためｊの更新としてｊにｊ＋１を
代入し、ステップＳ１０１１へ移行する。
【０３４８】
ステップＳ１０１１では、「ｙの全リテンション・タイムで処理終了？」の判断処理が
実行される。この処理では、ｙの全てのリテンション・タイムの処理が終了したか否かが
判断される。終了している（ＹＥＳ）場合は、ｙの全リテンション・タイムの処理が終了
したと判断し、ステップＳ１０１２へ移行する。終了していない（ＮＯ）場合は、ｙ中に
未処理のリテンション・タイムが残っていると判断し、ステップＳ１００６へ移行する。
つまり、Ｓ１００６〜Ｓ１０１１までの処理はｙの全てのリテンション・タイムが処理さ
れるまで繰り返す。
【０３４９】
ステップＳ１０１２では、「ｉの更新（ｉ←ｉ＋１）」の処理が実行される。この処理
では、ｘの処理を次のリテンション・タイムへ移行するためｉの更新としてｉにｉ＋１を
代入し、ステップＳ１０１３へ移行する。
【０３５０】
ステップＳ１０１３では、「ｘの全リテンション・タイムで処理終了？」の判断処理が
実行される。この処理では、ｘの全てのリテンション・タイムの処理が終了したか否かが
判断される。終了している（ＹＥＳ）場合は、ｘの全リテンション・タイムの処理が終了
したと判断し、ステップＳ１０１４へ移行する。終了していない（ＮＯ）場合は、ｘ中に
未処理のリテンション・タイムが残っていると判断し、ステップＳ１００４へ移行する。
つまり、Ｓ１００４〜Ｓ１０１３までの処理はｘの全てのリテンション・タイムが処理さ
れるまで繰り返す。
【０３５１】
ステップＳ１０１４では、「ＲＰ＿ａｌｌから最小値を取得（ＲＰ＿ｍｉｎ）」の処理
が実行される。この処理では、対象ＦＰと基準ＦＰとのリテンション・タイム・出現パタ
ーンの全組み合わせでのＲＰが保存されたＲＰ＿ａｌｌ中の最小値を、ＲＰ＿ｍｉｎとし
て取得し、そのＲＰ＿ｍｉｎを図９８のＳ２０７に渡し、リテンション・タイム・出現パ
ターンの一致度計算処理を終了する。
[サブルーチン２]
図１０４は、図９９の「対象ＦＰ帰属処理２」における「サブルーチン２」の詳細を示
すフローチャートである。この処理は、ＵＶスペクトルの一致度を計算する。
【０３５２】
ステップＳ２００１では、「ｘ←Ｕ１、ｙ←Ｕ２、ｚ←０」の処理が実行される。この
処理では、図９９のＳ３０２とＳ３０４で取得したＵＶスペクトルＵ１及びＵ２をそれぞ
れｘとｙに代入し、さらにＵＶスペクトル間の距離の二乗和（ｚ）の初期値として０を代
入し、ステップＳ２００２へ移行する。
【０３５３】
ステップＳ２００２では、「ｘデータ数を取得（ａ）」の処理が実行される。この処理
では、ｘのデータ数をａとして取得し、ステップＳ２００３へ移行する。
【０３５４】
ステップＳ２００３では、ＵＶスペクトルＵ１を構成する各検出波長における吸光度を
ｘから順番に呼び出すための初期値としてｉに１を代入し、ステップＳ２００４へ移行す
る。
【０３５５】
ステップＳ２００４では、「ｘｉ番目のデータを取得（ｂ）」の処理が実行される。こ
の処理では、ＵＶスペクトルＵ１を代入したｘのｉ番目の吸光度データをｂとして取得し
、ステップＳ２００５へ移行する。
【０３５６】
ステップＳ２００５では、「ｙｉ番目のデータを取得（ｃ）」の処理が実行される。こ
の処理では、ＵＶスペクトルＵ２を代入したｙのｉ番目の吸光度データをｃとして取得し
、ステップＳ２００６へ移行する。
【０３５７】
ステップＳ２００６では、「ＵＶスペクトル間距離（ｄ）とＵＶスペクトル間距離の二
乗和（ｚ）を算出」の処理が実行される。この処理では、ＵＶスペクトル間距離ｄとＵＶ
スペクトル間距離の二乗和ｚを、
ｄ ＝ ｂ−ｃ
ｚ ＝ ｚ＋ｄ²
として算出し、ステップＳ２００７へ移行する。
【０３５８】
ステップＳ２００７では、「ｉの更新（ｉ←ｉ＋１）」の処理が実行される。この処理
では、ｉの更新としてｉにｉ＋１を代入し、ステップＳ２００８へ移行する。
【０３５９】
ステップＳ２００８では、「ｘの全データで処理終了？」の判断処理が実行される。こ
の処理では、ｘとｙの全てのデータの処理が終了したか否かが判断される。終了している
（ＹＥＳ）場合は、ｘとｙの全データの処理が終了したと判断し、ステップＳ２００９へ
移行する。終了していない（ＮＯ）場合は、ｘとｙに未処理のデータが残っていると判断
し、ステップＳ２００４へ移行する。つまり、Ｓ２００４〜Ｓ２００８までの処理はｘと
ｙの全ての吸光度データが処理されるまで繰り返す。
【０３６０】
ステップＳ２００９では、「ｘとｙのＵＶスペクトルの一致度を計算（ＵＶ＿Ｓｉｍ）
」の処理が実行される。この処理では、ＵＶ＿Ｓｉｍを、前記ＵＶスペクトル間距離の二
乗和ｚとｘのデータ数ａから、
ＵＶ＿Ｓｉｍ ＝ √（ｚ／ａ）
として算出し、このＵＶ＿Ｓｉｍを図９９のステップＳ３０６に渡し、ＵＶスペクトルの
一致度計算処理を終了する。
[サブルーチン３]
図１０５は、図９９の「対象ＦＰ帰属処理２」における「サブルーチン３」の詳細を示
すフローチャートである。この処理は、ピーク・パターンの一致度を計算する。
【０３６１】
ステップＳ３００１では、「ピーク・パターン構成候補数（ｍ）とピーク・パターン構
成ピーク数（ｎ）を設定」の処理が実行される。この処理では、ピーク・パターンを網羅
的に作成するための設定として、ピーク・パターン構成候補数（ｍ）とピーク・パターン
構成ピーク数（ｎ）をそれぞれ設定し、ステップS３００２へ移行する。
【０３６２】
ステップＳ３００２では、「ｘ←対象ＦＰ名、ｒ１←Ｒ１、ｐ１←Ｐ１、ｙ←基準ＦＰ
名、ｒ２←Ｒ２、ｐ２←Ｐ２」の処理が実行される。この処理では、処理に必要とする対
象ＦＰ及び基準ＦＰのファイル名ならびに図９９のＳ３０２とＳ３０４で取得したリテン
ション・タイム及びピーク・データをそれぞれｘ、ｒ１、ｐ１とｙ、ｒ２、ｐ２に代入し
、ステップＳ３００３へ移行する。
【０３６３】
ステップＳ３００３では、「ｘの全リテンション・タイムを取得（ａ）」の処理が実行
される。この処理では、Ｓ３００２でｘに代入した名前のファイル（対象ＦＰ）を読み込
み、そのファイルの全リテンション・タイムをａとして取得し、ステップＳ３００４へ移
行する。
【０３６４】
ステップＳ３００４では、「ｙの全リテンション・タイムを取得（ｂ）」の処理が実行
される。この処理では、Ｓ３００２でｙに代入した名前のファイル（基準ＦＰ）を読み込
み、そのファイルの全リテンション・タイムをｂとして取得し、ステップＳ３００５へ移
行する。
【０３６５】
ステップＳ３００５では、「ａからｒ１のピーク・パターン構成候補ピークｍ本のリテ
ンション・タイム（ｃｍ）及びピーク・データ（ｄｍ）を取得」の処理が実行される。こ
の処理では、ａから帰属対象ピークのリテンション・タイムであるｒ１のピーク・パター
ン構成候補ピークｍ本のリテンション・タイムをｃｍ、ピーク・データをｄｍとしてそれ
ぞれ取得し、ステップＳ３００６へ移行する。なお、ピーク・パターン構成候補ピークｍ
本は、ｒ１とリテンション・タイムが近いｍ本である。
【０３６６】
ステップＳ３００６では、「ｂからｒ２のピーク・パターン構成候補ピークｍ本のリテ
ンション・タイム（ｅｍ）及びピーク・データ（ｆｍ）を取得」の処理が実行される。こ
の処理では、ｂから帰属候補ピークのリテンション・タイムであるｒ２のピーク・パター
ン構成候補ピークｍ本のリテンション・タイムをｅｍ、ピーク・データをｆｍとしてそれ
ぞれ取得し、ステップＳ３００７へ移行する。なお、ピーク・パターン構成候補ピークｍ
本は、ｒ２とリテンション・タイムが近いｍ本である。
【０３６７】
ステップＳ３００７では、「ｃｍ、ｄｍをリテンション・タイム順（昇順）に並べる」
の処理が実行される。この処理では、Ｓ３００５で取得したｃｍとｄｍをリテンション・
タイムが昇順になるように並び替え、ステップＳ３００８へ移行する。
【０３６８】
ステップＳ３００８では、「ｅｍ、ｆｍをリテンション・タイム順（昇順）に並べる」
の処理が実行される。この処理では、Ｓ３００６で取得したｅｍとｆｍをリテンション・
タイムが昇順になるように並び替え、ステップＳ３００９へ移行する。
【０３６９】
ステップＳ３００９では、「ｃｍ、ｄｍからピーク・パターン構成ピークｎ本のリテン
ション・タイム（ｃｎ）及びピーク・データ（ｄｎ）を順番に取得」の処理が実行される
。この処理では、ピーク・パターン構成候補ピークｍ本のｃｍ及びｄｍからピーク・パタ
ーン構成ピークｎ本のリテンション・タイムをｃｎ、ピーク・データをｄｎとして順番に
取得し、ステップＳ３０１０へ移行する。
【０３７０】
ステップＳ３０１０では、「ｅｍ、ｆｍからピーク・パターン構成ピークｎ本のリテン
ション・タイム（ｅｎ）及びピーク・データ（ｆｎ）を順番に取得」の処理が実行される
。この処理では、ピーク・パターン構成候補ピークｍ本のｅｍ及びｆｍからピーク・パタ
ーン構成ピークｎ本のリテンション・タイムをｅｎ、ピーク・データをｆｎとして順番に
取得し、ステップＳ３０１１へ移行する。
【０３７１】
ステップＳ３０１１では、「ピーク・パターンの一致度を計算（Ｐ＿Ｓｉｍ）」の処理
が実行される。この処理では、これまでに取得した帰属対象ピークのｒ１とｐ１及びその
ピーク・パターン構成ピークｎ本のｃｎとｄｎ、ならびに帰属候補ピークのｒ２とｐ２及
びそのピーク・パターン構成ピークｎ本のｅｎとｆｎから、ピーク・パターンの一致度（
Ｐ＿Ｓｉｍ）を、ｎ＝４の場合を例とすると、図６６中のように
Ｐ＿Ｓｉｍ＝ (|p1−p2|＋1)×(|(r1−(r2＋d)|＋1)
＋(|dn1−fn1|＋1)×(|(cn1−r1)−(en1−r2)|+1)
＋(|dn2−fn2|＋1)×(|(cn2−r1)−(en2−r2)|+1)
＋(|dn3−fn3|＋1)×(|(cn3−r1)−(en3−r2)|+1)
＋(|dn4−fn4|＋1)×(|(cn4−r1)−(en4−r2)|+1)
として算出し、ステップＳ３０１２へ移行する。
【０３７２】
ステップＳ３０１２では、「Ｐ＿Ｓｉｍを保存（Ｐ＿Ｓｉｍ＿ａｌｌ）」の処理が実行
される。この処理では、Ｓ３０１１で算出されたＰ＿Ｓｉｍを順次Ｐ＿Ｓｉｍ＿ａｌｌに
保存し、ステップＳ３０１３へ移行する。
【０３７３】
ステップＳ３０１３では、「ｅｍ中のｍ本からｎ本を取り出す全組み合わせ終了？」の
判断処理が実行される。この処理では、帰属候補ピークのピーク・パターン構成候補ピー
クｍ本からピーク・パターン構成ピークｎ本を取り出す全組み合わせで処理が終了したか
否かが判断される。終了した（ＹＥＳ）場合は、帰属候補ピークにおいて、網羅的なピー
ク・パターンの作成とそのパターンにおける一致度の計算が終了したと判断し、ステップ
Ｓ３０１４へ移行する。終了していない（ＮＯ）場合は、ｍ本からｎ本を取り出す組み合
わせが終了していないと判断し、ステップＳ３０１０へ移行する。つまり、Ｓ３０１０〜
Ｓ３０１３までの処理はｍ本からｎ本を取り出す全組み合わせで処理が終了するまで繰り
返す。
【０３７４】
ステップＳ３０１４では、「ｃｍ中のｍ本からｎ本を取り出す全組み合わせ終了？」の
判断処理が実行される。この処理では、帰属対象ピークのピーク・パターン構成候補ピー
クｍ本からピーク・パターン構成ピークｎ本を取り出す全組み合わせで処理が終了したか
否かが判断される。終了した（ＹＥＳ）場合は、帰属対象ピークにおいて、網羅的なピー
ク・パターンの作成とそのパターンにおける一致度の計算が終了したと判断し、ステップ
Ｓ３０１５へ移行する。終了していない（ＮＯ）場合は、ｍ本からｎ本を取り出す組み合
わせが終了していないと判断し、ステップＳ３００９へ移行する。つまり、Ｓ３００９〜
Ｓ３０１４までの処理はｍ本からｎ本を取り出す全組み合わせで処理が終了するまで繰り
返す。
【０３７５】
ステップＳ３０１５では、「Ｐ＿Ｓｉｍ＿ａｌｌから最小値を取得（Ｐ＿Ｓｉｍ＿ｍｉ
ｎ）」の処理が実行される。この処理では、Ｓ３０１２で保存したＰ＿Ｓｉｍ＿ａｌｌの
最小値をＰ＿Ｓｉｍ＿ｍｉｎとして取得し、このＰ＿Ｓｉｍ＿ｍｉｎを図９９のステップ
Ｓ３０７に渡し、ピーク・パターンの一致度計算処理を終了する。
［Ｓ６：対象ＦＰタイプ２の作成処理］
図１０６は、図９３ステップＳ６の「ＦＰ_type２作成」の詳細を示すフローチャート
である。
【０３７６】
ステップＳ６０１では、「対象ＦＰを読み込む」の処理が実行される。この処理では、
前記対象ＦＰ４３（図１１９のＦＰのデータ例１８７参照）のファイルを読み込み、ステ
ップＳ６０２へ移行する。
【０３７７】
ステップＳ６０２では、「ピーク・データ特徴量ファイルを読み込む」の処理が実行さ
れる。この処理では、前記対象ＦＰ４３に関し、ピーク・データ特徴量のファイル（図１
２４のピーク・データ特徴量のファイル例１９９参照）を読み込み、ステップＳ６０３へ
移行する。ピーク・データ特徴量ファイル例は、前記対象ＦＰピーク特徴量作成部７によ
り、基準群ＦＰ４５のピークに帰属された対象ＦＰ４３のピーク情報を含んでいる。
【０３７８】
ステップＳ６０３では、「対象ＦＰとピーク・データ特徴量ファイルを比較」の処理が
実行される。この処理では、対象ＦＰ４３のファイルとピーク・データ特徴量ファイルと
が比較される。この比較により、基準群ＦＰ４５のピークに帰属されなかった対象ＦＰ４
３の残存ピークが特定され、ステップＳ６０４へ移行する。
【０３７９】
ステップＳ６０４では、「対象ＦＰのみに存在するピークのリテンション・タイムとピ
ーク・データを出力」の処理が実行される。この処理では、対象ＦＰ４３の残存ピークの
リテンション・タイムとピーク・データとが、対象ＦＰタイプ２のデータ・ファイル（図
１２５の基準及び対象ＦＰタイプ２のデータ例２０１参照）に出力される。
［Ｓ７：領域分割による対象ＦＰ_type２の特徴量化処理］
図１０７は、図９４ステップＳ７の「領域分割による対象ＦＰ_type２の特徴量化処理
」の詳細を示すフローチャートである。
【０３８０】
ステップＳ７０１では、「ＦＰ空間の領域分割条件の設定」の処理が実行される。この
処理では、対象ＦＰタイプ２の領域を分割するため、縦・横の線（分割線）１本目の位置
がそれぞれ１個設定される。この設定により、ＦＰ空間に対し、例えば図７６、図７７の
ように、縦・横分割線（1本目）が分割線として設定される。但し、対象ＦＰタイプ２の
場合、領域の位置変更はないので、振幅は関係しない。ステップＳ７０１で縦・横分割線
（1本目）が設定されるとステップＳ７０２へ移行する。
【０３８１】
ステップＳ７０２では、「ＦＰ空間の領域分割パターンの作成」の処理が実行される。
この処理では、縦・横分割線１本目の全組み合わせで２本目以降の分割線の位置が設定さ
れ、分割パターン（１個）が作成される。この処理により、ＦＰ空間に対し、例えば図７
８のように、縦・横分割線による領域分割が行われる。領域分割が行われるとステップＳ
７０３へ移行する。
【０３８２】
ステップＳ７０３では、「対象ＦＰ_type２のファイルを読み込む」の処理が実行され
る。この処理により、対象ＦＰタイプ２のファイルが読み込まれ、ステップＳ７０４へ移
行する。
【０３８３】
ステップＳ７０４では、「ＦＰ空間全体の総ピーク・データを算出」の処理が実行され
る。この処理では、例えば、図７９のように分割された各格子１４５全てに存在するピー
クの高さ合計が算出され（図８１）、ステップＳ７０５へ移行する。
【０３８４】
ステップＳ７０５では、「ＦＰ空間を分割パターンで分割」の処理が実行される。、こ
の処理では、Ｓ７０２で設定された領域分割パターンによりＳ７０３で読み込まれた対象
ＦＰタイプ２を図７９のように領域分割し、ステップＳ７０６へ移行する。
【０３８５】
ステップＳ７０６では、「分割された領域内でのピーク・データの存在割合を算出」の
処理が実行される。この処理では、各格子１４５内でのピーク存在割合が前記特徴量＝領
域内ピーク高さ合計／全ピーク高さ合計として算出される。算出結果は、図８６のように
なる。算出が終了するとステップＳ７０７へ移行する。
【０３８６】
ステップＳ７０７では、「各領域の存在割合を特徴量として出力」の処理が実行される
。この処理では、１通りでの対象ＦＰ領域分割特徴量ファイル（図１２６で示す１通りで
の対象ＦＰ領域分割特徴量ファイル例２０３参照）が出力される。
［Ｓ８：ピーク・データ特徴量と領域分割特徴量の統合］
図１０８は、図９４ステップＳ８の「ピーク・データ特徴量と領域分割特徴量の統合」
の詳細を示すフローチャートである。
【０３８７】
ステップＳ８０１では、「ピーク・データ特徴量ファイルを読み込む」の処理が実行さ
れる。この処理により、図１２４で示すピーク・データ特徴量のファイル例１９９と同様
のファイルが読み込まれ、ステップＳ８０２へ移行する。
【０３８８】
ステップＳ８０２では、「領域分割特徴量ファイルを読み込む」の処理が実行される。
この処理により、図１２６で示す対象ＦＰ領域分割特徴量ファイル２０３が読み込まれ、
ステップＳ８０３へ移行する。 ステップＳ８０３では、「２つの特徴量データを横１行
のデータとして統合」の処理が実行される。この処理により、ピーク・データ特徴量のフ
ァイル（図１２４で示
すピーク・データ特徴量のファイル例１９９参照）及び対象ＦＰ領域分割特徴量ファイル
（図１２６で示す対象ＦＰ領域分割特徴量ファイル例２０３参照）が一行のである対象Ｆ
Ｐ特徴量統合ファイル（図１２７の対象ＦＰ特徴量統合ファイル例２０５参照）として統
合され、ステップＳ８０４へ移行する。
【０３８９】
ステップＳ８０４では、「統合したデータを出力」の処理が実行される。この処理では
、図１２７の対象ＦＰ特徴量統合ファイル２０５が出力される。
［基準ＦＰ帰属結果特徴量統合ファイルの作成］
対象ＦＰ特徴量統合データを基準ＦＰ特徴量統合データと比較するための基準ＦＰ特徴
量統合ファイルは、図１０９〜図１１６のように作成される。
【０３９０】
図１０９、図１１０は、基準ＦＰ特徴量統合ファイルを作成するためのフローチャート
であり、基準ＦＰ作成部３１のＦＰ作成機能と、基準ＦＰピーク帰属部１５の基準ＦＰピ
ーク帰属機能と、基準ＦＰ帰属結果統合部１７の基準ＦＰ帰属結果統合機能と、基準ＦＰ
ピーク特徴量作成部１９の基準ＦＰピーク特徴量作成機能と、基準ＦＰタイプ２作成部２
１の基準ＦＰタイプ２作成機能と、基準ＦＰ領域分割特徴量作成部２３の基準ＦＰ領域分
割特徴量作成機能と、基準ＦＰ特徴量統合部２５の基準ＦＰ特徴量統合機能とをコンピュ
ータに実現させる。
【０３９１】
基準ＦＰ作成機能は、ステップＳ１０００１で実現される。基準ＦＰピーク帰属機能は
、ステップＳ１０００２、Ｓ１０００３、Ｓ１０００４で実現される。基準ＦＰ帰属結果
統合機能は、ステップＳ１０００５で実現される。基準ＦＰピーク特徴量作成機能は、ス
テップＳ１０００６で実現される。基準ＦＰタイプ２作成機能は、ステップＳ１０００７
で実現される。基準ＦＰ領域分割特徴量作成機能は、ステップＳ１０００８で実現される
。基準ＦＰ特徴量統合機能は、ステップＳ１０００９で実現される。
【０３９２】
Ｓ１０００１〜Ｓ１０００４は、図９３、９４の対象ＦＰ特徴量統合ファイルの作成に
係るＳ１〜Ｓ４に対応し、Ｓ１００７〜Ｓ１０００９は、同Ｓ６〜Ｓ８に対応している。
【０３９３】
ステップＳ１０００１は、３Ｄクロマト及び特定の検出波長におけるピーク情報を入力
データとして「ＦＰ作成処理」が実行される。
【０３９４】
３Ｄクロマト及びピーク・データ共に、評価基準となる複数の評価基準薬剤（基準漢方
薬）ごとに備えられている。
【０３９５】
ステップＳ１０００１では、コンピュータのＦＰ作成部３の基準ＦＰ作成部３１（図１
）が機能し、３Ｄクロマト及びピーク情報から基準ＦＰが対象ＦＰ４３（図２）と同様に
作成され、基準ＦＰのデータがファイルとして出力される。
【０３９６】
ステップＳ１０００２は、ステップＳ１０００１で出力した全基準ＦＰを入力として、
「基準ＦＰ帰属処理１」が実行される。
【０３９７】
ステップＳ１０００２では、コンピュータの基準ＦＰピーク帰属部１５が機能し、全基
準ＦＰを対象とし、選択した組み合わせ且つ順番で帰属スコアを算出するために全基準Ｆ
Ｐから組み合わせを選択してステップＳ１０００３へ移行する。
【０３９８】
ステップＳ１０００３は、選択した基準ＦＰの組み合わせを入力とし、「基準ＦＰ帰属
処理２」が実行される。
【０３９９】
ステップＳ１０００３では、ステップＳ２で選定した基準ＦＰの組み合わせの全ピーク
で、図２３〜図６１のように網羅的にピーク・パターンを作成し、次にそれらピーク・パ
ターンの一致度（図６３または図６４のＰ＿Ｓｉｍ）を算出する。また、選定した基準Ｆ
Ｐの組み合わせのピーク間でＵＶスペクトルの一致度（図６６のＵＶ＿Ｓｉｍ）を算出す
る。さらにこれら２つの一致度から帰属候補ピークの一致度（図６７のＳＣＯＲＥ）を算
出する。その算出結果は、判定結果ファイル（図１２０の判定結果ファイル例１８９参照
）として出力される。
【０４００】
ステップＳ１０００４は、ステップＳ１０００３で出力した判定結果ファイルを入力と
し、「基準ＦＰ帰属処理３」が実行される。
【０４０１】
ステップＳ１０００４では、選定した基準ＦＰの組み合わせ間で、帰属候補ピークの一
致度（ＳＣＯＲＥ）をもとに選定した基準ＦＰの組み合わせ間で対応するピークを特定す
る。その結果は、基準ＦＰごとに基準ＦＰ帰属データとして出力される。
【０４０２】
ステップＳ１０００５は、ステップＳ１０００４で出力した全基準ＦＰ帰属データを入
力とし、「基準ＦＰ帰属結果統合処理」が実行される。
【０４０３】
ステップＳ１０００５では、コンピュータの基準ＦＰ帰属結果統合部１７が機能し、、
基準ＦＰピーク帰属部１５で特定した各基準ＦＰのピーク対応関係を参照し、全基準ＦＰ
帰属データを統合して基準ＦＰ対応表を作成し、ステップＳ１０００６へ移行する。
ステップＳ１０００６では、コンピュータの基準ＦＰピーク特徴量作成部１９が機能し、
基準ＦＰ帰属結果統合部１７で作成された基準ＦＰ対応表を基に全基準ＦＰによるピーク
特徴量（基準群ＦＰ）を作成する。この基準ＦＰピーク特徴量作成部１９での処理は、基
準ＦＰ対応表の各ピーク（列）で統計量（最大値、最小値、中央値、平均値、ｅｔｃ）を
算出し、その情報を基にピーク（列）を選定する。選定したピーク（列）は、基準群ＦＰ
（図１２３の基準群ＦＰ例１９７参照）として出力される。
【０４０４】
ステップＳ１０００７は、ステップＳ１０００６で出力した基準群ＦＰと全基準ＦＰと
を入力とし、「ＦＰ_type２の作成」の処理が実行される。
【０４０５】
ステップＳ１０００７では、コンピュータの基準ＦＰタイプ２作成部２１が対象ＦＰタ
イプ２作成部９と同様に機能し、図９３のステップＳ６と同様にして、複数の各基準ＦＰ
から前記特徴量化されたピークをそれぞれ除いて残ったピークとそのリテンション・タイ
ムとで構成されるＦＰを基準ＦＰタイプ２としてそれぞれ作成する（図１２５のＦＰタイ
プ２ファイル例２０１参照）。
【０４０６】
ステップＳ１０００８では、「基準ＦＰ_type２の特徴量化処理」が実行される。
この処理では、コンピュータの基準ＦＰ領域分割特徴量作成部２３が機能し、図７３〜図
８５の領域分割により、基準ＦＰ領域分割特徴量が作成される。その結果は、基準ｔｙｐ
ｅ２群ＦＰ（図１２８の基準ｔｙｐｅ２群ＦＰ例２０７参照）として出力される。
【０４０７】
ステップＳ１０００９では、「基準データ作成処理」の処理が実行される。この処理で
は、コンピュータの基準ＦＰ特徴量統合部２５が機能し、基準ＦＰピーク特徴量作成部１
９で作成された基準群ＦＰと基準ＦＰ領域分割特徴量作成部２３で作成された基準ｔｙｐ
ｅ２群ＦＰとを統合して全基準ＦＰの特徴量データを作成する。その結果は、基準群統合
データ（図１２９の基準群統合データ例２０９参照）として出力される。
［Ｓ１０００５：基準ＦＰ対応表の作成］
図１１１、図１１２は、図１１０ステップＳ１０００５の「基準ＦＰ帰属結果統合処理
（基準ＦＰ対応表の作成）」の詳細を示すフローチャートである。
【０４０８】
ステップＳ１０１０１では、「帰属順番１番の帰属データを統合データとして読み込む
」の処理が実行される。この処理では、Ｓ１０００４で１番目に帰属処理しピークの対応
関係を特定した基準ＦＰ帰属データを統合データとして読み込み、ステップＳ１０１０２
へ移行する。
【０４０９】
ステップＳ１０１０２では、「帰属データ２番目以降を順番に読み込む」の処理が実行
される。この処理では、まずＳ１０００４で２番目に帰属処理しピークの対応関係を特定
した基準ＦＰ帰属データを統合データのとして読み込み、ステップＳ１０１０３へ移行す
る。
【０４１０】
ステップＳ１０１０３では、「統合データと帰属データを共通するピーク・データで統
合」の処理が実行される。この処理では、統合データと帰属データで共通に存在する基準
ＦＰのピーク・データをもとに２つのファイルを統合し、その結果で統合データを更新し
、ステップＳ１０１０４へ移行する。
【０４１１】
ステップＳ１０１０４では、「帰属データ中の全てのピークを統合データに追加？」の
判断処理が実行される。この処理では、帰属データの全てのピークが統合データに追加さ
れたか否かが判断される。追加された（ＹＥＳ）場合は、ステップＳ１０１０５へ移行す
る。追加されていないピーク（欠落ピーク）がある（ＮＯ）場合は、この欠落ピークを統
合データに追加処理するため、ステップＳ１０１０７へ移行する。尚、欠落ピークの統合
データへの追加処理（Ｓ１０１０７−Ｓ１０１２０）は、前記Ｓ５（対象ＦＰ帰属処理４
）におけるステップＳ５０４−Ｓ５１７と同様の処理が行われる。
【０４１２】
ステップＳ１０１２１では、「統合データにＴＥＭＰ２のデータを追加（全てのリテン
ション・タイムとピーク）」の処理が実行される。この処理では、ＴＥＭＰ２の全てのリ
テンション・タイム（Ｒ３）とピーク（Ｐ１）を統合データの該当箇所に追加し、ステッ
プＳ１０１２２へ移行する。
【０４１３】
ステップＳ１０１２２では、「閾値２ ← 初期値、ＴＥＭＰ２内のデータを全て削除
」の処理が実行される。この処理では、ＵＶ＿Ｓｉｍに更新されている閾値２を初期値に
戻し、全欠落ピークのリテンション・タイムやピークなどのデータが入ったＴＥＭＰ２か
ら全てのデータを削除し、ステップＳ１０１０４へ戻る。
【０４１４】
ステップＳ１０１０４から移行するステップＳ１０１０５では、「全ての帰属データの
処理終了？」の判断処理が実行される。この処理では、全基準データの処理が終了したか
否かが判断される。処理が終了（ＹＥＳ）の場合は、全帰属データの統合結果である基準
ＦＰ対応表を出力するため、ステップＳ１０１０６へ移行する。全ての処理が終了してい
ない（ＮＯ）の場合は、ステップＳ１０１０２へ戻り、残りの帰属データを順次処理する
。
【０４１５】
ステップＳ１０１０６では、「統合データを出力（基準ＦＰ対応表）」の処理が実行さ
れる。この処理では、全帰属データを統合した結果を基準ＦＰ対応表として出力し、基準
ＦＰ対応表の作成処理を終了する。
［Ｓ１０００６：ピーク特徴量化処理］
図１１３は、図１０９ステップＳ１０００６の「ピーク特徴量化処理（基準群ＦＰの作
成）」の詳細を示すフローチャートである。
【０４１６】
ステップＳ１０２０１では、「基準ＦＰ対応表を読み込む」の処理が実行される。この
処理では、Ｓ１０００５で作成した基準ＦＰ対応表を読み込み、ステップＳ１０２０２へ
移行する。
【０４１７】
ステップＳ１０２０２では、「各ピーク（列）で統計量を算出」の処理が実行される。
この処理では、基準ＦＰ対応表の各ピーク(列)で統計量（最大値、最小値、中央値、平均
値、分散、標準偏差、存在数、存在率）を算出し、ステップＳ１０２０３へ移行する。
【０４１８】
ステップＳ１０２０３では、「算出した統計量を参考にピーク（列）を選定」の処理が
実行される。この処理では、Ｓ１０１０２で算出した統計量を参考にピークを選定し、ス
テップＳ１０２０４へ移行する。
【０４１９】
ステップＳ１０２０４では、「選定したピーク（列）を出力（基準群ＦＰ）」の処理が
実行される。この処理では、統計量によりピーク（列）の選定結果を基準群ＦＰとして出
力し、基準群ＦＰの作成処理を終了する。
【０４２０】
図１２３に前記のように出力する基準ＦＰ対応表例１９７を示す。
［Ｓ１０００７：基準ＦＰタイプ２の作成処理］
図１１４は、図１１０ステップＳ１０００７の「基準ＦＰ編集処理（基準ＦＰ_type２
の作成）」の詳細を示すフローチャートである。
【０４２１】
ステップＳ１０３０１では、「基準ＦＰを順番に読み込む」の処理が実行される。この
処理では、複数の基準ＦＰのファイル（図１１９のＦＰのデータ例１８７参照）を読み込
み、ステップＳ１０３０２へ移行する。
【０４２２】
ステップＳ１０３０２では、「基準群ＦＰを読み込む」の処理が実行される。この処理
では、基準群ＦＰのデータ・ファイル（図１２３の基準群ＦＰのデータ例１９７参照）が
読み込まれ、ステップＳ１０３０３へ移行する。
【０４２３】
ステップＳ１０３０３では、「基準群ＦＰから基準ＦＰのピーク・データ特徴量を取り
出す」の処理が実行される。この処理では、基準群ＦＰ４５のファイルから基準ＦＰの帰
属処理されたピーク・データ特徴量を取り出し、ステップＳ１０３０４へ移行する。
【０４２４】
ステップＳ１０３０４では、「基準ＦＰと取り出したピーク・データ特徴量ファイルを
比較」の処理が実行され、基準ＦＰがピーク・データ特徴量ファイルと比較され、ステッ
プＳ１０３０５へ移行する。
【０４２５】
ステップＳ１０３０５では、「基準ＦＰのみに存在するピークのリテンション・タイム
とピーク・データを出力」の処理が実行され、基準ＦＰからピーク・データ特徴量ファイ
ルのピークが除かれて、ステップＳ１０３０６へ移行する。
【０４２６】
ステップＳ１０３０６では、「全基準ＦＰで処理終了？」の判断処理が実行される。こ
の処理では、全基準ＦＰで処理終了したときは（ＹＥＳ）、Ｓ１０００７は終了し、全基
準ＦＰで処理終了していなければ（ＮＯ）、Ｓ１０３０１−Ｓ１０３０５が繰り返される
。したがって、複数の基準ＦＰが順に処理され、各基準ＦＰからピーク・データ特徴量フ
ァイルのピークが除かれて基準ＦＰタイプ２のファイル（図１２５に示す対象及び基準Ｆ
Ｐタイプ２のデータ例２０１参照）が作成される。
［Ｓ１０００８：領域分割による基準ＦＰ_type２の特徴量化処理］
図１１５は、図１１０ステップＳ１０００８の「領域分割による基準ＦＰ_type２の特
徴量化処理」の詳細を示すフローチャートである。
【０４２７】
ステップＳ１０４０１では、「ＦＰ空間の領域分割条件の設定」の処理が実行される。
この処理では、基準ＦＰタイプ２の領域を分割するため、縦・横の線（分割線）１本目の
位置がそれぞれ複数個設定される。この設定により、ＦＰ空間に対し、例えば図７６、図
７７のように、縦・横分割線（１本目）１４１、１４３が分割線として複数個設定される
。縦・横分割線（１本目）１４１、１４３が複数個設定されるとステップＳ１０４０２へ
移行する。
【０４２８】
ステップＳ１０４０２では、「ＦＰ空間の領域分割パターンの設定」の処理が実行され
る。この処理では、縦・横分割線１本目の全組み合わせで各２本目以降の分割線の位置が
設定され、分割パターン（ｍ×ｎ個）が作成される。この設定により、ＦＰ空間に対し、
例えば図７８のように、縦・横分割線１４１、１４３による領域分割のパターンが複数設
定される。領域分割が行われるとステップＳ１０４０３へ移行する。
【０４２９】
ステップＳ１０４０３では、「基準ＦＰ_type２のファイルを順番に読み込む」の処理
が実行される。この処理により、基準ＦＰタイプ２のファイルが読み込まれ、ステップＳ
１０４０４へ移行する。
【０４３０】
ステップＳ１０４０４では、「ＦＰ空間全体の総ピーク・データを算出」の処理が実行
される。この処理では、例えば、図７９のように分割された各格子１４５全てに存在する
ピークの高さ合計が算出され（図８１）、ステップＳ１０４０５へ移行する。
【０４３１】
ステップＳ１０４０５では、「ＦＰ空間を各分割パターンで順番に分割」の処理が実行
される。この処理では、Ｓ１０４０２で設定された複数の領域分割パターンでＦＰ空間が
順次分割され、ステップＳ１０４０６へ移行する。
【０４３２】
ステップＳ１０４０６では、「分割された領域内でのピーク・データの存在割合を算出
」の処理が実行される。この処理では、例えば、図７９のように分割された各格子１４５
全てに存在するピークの高さ合計が算出され（図８１）、図７９の各格子１４５内でのピ
ーク存在割合が前記特徴量＝領域内ピーク高さ合計／全ピーク高さ合計としてされる。算
出結果は、例えば図８３〜図８５のようになる。算出が終了するとステップＳ１０４０８
へ移行する。
【０４３３】
ステップＳ１０４０８では、「全分割パターンで分割終了」の処理が実行される。この
処理では、Ｓ１０４０２で設定された複数の全領域分割パターンでの特徴量処理が終了し
たか否かの判断が行われる。特徴量処理が終了すれば（ＹＥＳ）、ステップＳ１０４０９
へ移行し、特徴量処理が終了していなければ（ＮＯ）、ステップＳ１０４０５へ移行する
。Ｓ１０４０５〜Ｓ１０４０８は、全領域分割パターンでの特徴量処理が終了するまで繰
り返される。
【０４３４】
ステップＳ１０４０９では、「全基準ＦＰ_type２で処理終了？」の判断処理が実行さ
れる。この処理では、複数の基準ＦＰごとに作成された複数の基準ＦＰタイプ２の全てに
おいて特徴量処理が終了したか否かの判断が行われる。全基準ＦＰタイプ２が終了すれば
（ＹＥＳ）、Ｓ１０００８は終了し、全基準ＦＰタイプ２が終了していなければ（ＮＯ）
、ステップＳ１０４０３へ移行する。Ｓ１０４０３〜Ｓ１０４０９は、基準ＦＰタイプ２
での特徴量処理が終了するまで繰り返される。
【０４３５】
図１２８に、基準ｔｙｐｅ２群ＦＰ例２０７を示す。
［Ｓ１０００９：基準データ作成処理］
図１１６は、図１１０ステップＳ１０００９の「基準データの作成処理」の詳細を示す
フローチャートである。
【０４３６】
ステップＳ１０５０１では、「領域分割特徴量ファイルを読み込む」の処理が実行され
る。この処理により、基準ＦＰ領域分割特徴量ファイル（図１２８で示す基準ｔｙｐｅ２
群ＦＰ例２０７参照）が読み込まれ、ステップＳ１０５０２へ移行する。
【０４３７】
ステップＳ１０５０２では、「領域分割した際の分割パターン数を算出」の処理が実行
される。この処理により、領域分割の分割パターン数が算出される。分割パターン数は、
図７０〜図８０での説明のように、例えば１００通りと算出される。この算出後にステッ
プＳ１０５０３へ移行する。
【０４３８】
ステップＳ１０５０３では、「基準群ＦＰを読み込む」の処理が実行され、基準群ＦＰ
が読み込まれてステップＳ１０５０４へ移行する。
【０４３９】
ステップＳ１０５０４では、「基準群ＦＰの各行を分割パターン数分複製したファイル
（基準群ＦＰ２）を作成」の処理が実行される。この処理では、基準群ＦＰと領域分割特
徴量ファイルとを統合するために基準群ＦＰの行を分割パターン数に合わせて複製し、基
準群ＦＰ２を作成する。例えば、図１２３の基準群ＦＰのファイル例１９７を、図１２９
の基準群統合データ例２０９のピーク・データ特徴量（基準群ＦＰ２）と対応するように
複製する。この複製後に、ステップＳ１０５０５へ移行する。
【０４４０】
ステップＳ１０５０５では、「基準群ＦＰ２と領域分割特徴量ファイルを行ごとに統合
」の処理が実行される。この処理では、Ｓ１０５０４で複製された基準群ＦＰ２のデータ
と領域分割特徴量ファイルのデータとが行ごとに統合され、ステップＳ１０５０６へ移行
する。
【０４４１】
ステップＳ１０５０６では、「統合したデータを出力」の処理が実行される。この処理
では、統合結果による基準ＦＰ特徴量統合ファイル（図１２９の基準群統合データ例２０
９参照）が出力される。
［実施例１の効果］
本発明実施例１の多成分物質の評価方法では、前記ＦＰ作成工程１４８と、対象ＦＰピ
ーク帰属工程１４９と、対象ＦＰピーク特徴量作成工程１５１と、対象ＦＰタイプ２作成
工程１５３と、対象ＦＰ領域分割特徴量作成工程１５５と、対象ＦＰ特徴量統合工程１５
７と、基準ＦＰピーク帰属工程１５９と、基準ＦＰ帰属結果統合工程１６１と、基準ＦＰ
ピーク特徴量作成工程１６３と、基準ＦＰタイプ２作成工程１６５と、基準ＦＰ領域分割
特徴量作成工程１６７と、基準ＦＰ特徴量統合工程１６９と、評価工程１７１とを備えて
いる。
【０４４２】
前記ＦＰ作成工程１４８は、前記対象ＦＰ作成工程１７３と、基準ＦＰ作成工程１７５
とを備えている。
【０４４３】
前記対象ＦＰピーク帰属工程１４９は、前記基準ＦＰ選定工程１７７と、ピーク・パタ
ーン作成工程１７９と、ピーク帰属工程１８１とを備えている。
【０４４４】
評価対象となる多成分薬剤の３Ｄクロマト４１をこれら７工程（１７８、１４９、１５
１、１５３、１５５、１５７、１７１）で処理することで、評価対象薬剤の品質評価の精
度及び効率をより向上させることができる。
【０４４５】
特に、対象ＦＰ４３と複数の基準ＦＰとに基づき特徴量化された対象ＦＰピーク特徴量
を作成し、この特徴量化から漏れた対象ＦＰ４３の残存ピークとして対象ＦＰタイプ２を
作成し、この対象ＦＰタイプ２を複数の領域に分割し各領域に存在するピークの存在率か
ら対象ＦＰ領域分割特徴量を作成し、対象ＦＰピーク特徴量及び対象ＦＰ領域分割特徴量
を統合して対象ＦＰ統合特徴量を作成し、対象ＦＰ統合特徴量とこの対象ＦＰ統合特徴量
に対応し評価基準となる多成分物質の複数の基準ＦＰに基づく基準ＦＰ統合特徴量とを比
較評価するため、対象ＦＰピーク特徴量に含まれなかった対象ピークのピークをも含めて
評価することができ、評価対象薬剤の品質評価の精度を確実に向上させることができる。
【０４４６】
対象ＦＰ作成工程１７３で作成された対象ＦＰ４３は、３Ｄクロマト４１と同様に三次
元の情報（ピーク、リテンション・タイム及びＵＶスペクトル）で構成している。そのた
め、その薬剤特有の情報をそのまま継承したデータである。それにも係らず、データ容量
は約１／７０に圧縮されているため、３Ｄクロマト４１に比較して、処理すべき情報量を
大幅に減少させることができ処理速度を速めることができる。
【０４４７】
対象ＦＰ作成工程１７３は、検出波長の異なる複数のＦＰを合成したＦＰを作成する。
これにより、１つの波長で全ての成分を検出できない成分が組み合わさった多成分薬剤で
あっても、複数の検出波長のＦＰを合成することで、全ての成分を含んだ品質評価が可能
となる。
【０４４８】
対象ＦＰ作成工程１７３は、３Ｄクロマトで検出された全ピークを含んだＦＰを作成す
る。このため、多成分薬剤である漢方薬の品質評価に適している。
【０４４９】
基準ＦＰ選定工程１７７では、対象ＦＰの帰属に適した基準ＦＰをＦＰ間のリテンショ
ン・タイム・出現パターンを比較し、パターンの一致度が良い基準ＦＰ選定する。これに
より、ピーク帰属工程１８１において、パターンが類似したＦＰ間で帰属処理ができるた
め、精度の高い帰属が可能となる。
【０４５０】
ピーク・パターン作成工程１７９では、帰属対象ピーク及び帰属候補ピークそれぞれで
複数の周辺ピークを使って網羅的にピーク・パターンを作成する。これにより、対象ＦＰ
と基準ＦＰでＦＰ全体のパターンが多少異なっていても、ピーク帰属工程１８１で精度の
高い帰属が可能となる。
【０４５１】
ピーク帰属工程１８１では、ピーク・パターン作成工程１７９で作成されたピーク・パ
ターンの一致度に加え、帰属対象ピークと帰属候補ピークのＵＶスペクトルの一致度も加
味して、帰属すべきピークを特定している。そのため、精度の高い帰属が可能となる。
【０４５２】
ピーク帰属工程１８１では、対象ＦＰの全ピークを基準ＦＰのピークへ一斉帰属する。
そのため、効率のよい帰属処理が可能となる。
【０４５３】
評価工程１７１では、多次元データである多成分で構成されたＦＰをＭＴ法でＭＤ値と
して1次元に集約し、複数の評価対象ロットを簡単に比較評価する。このため、複数の成
分で構成される多成分系薬剤の評価に適している。
【０４５４】
対象ＦＰ領域分割特徴量作成工程１５５は、シグナル強度軸に平行な複数の縦分割線
１４１と時間軸に平行な複数の横分割線１４３とにより前記領域の分割を行う。
【０４５５】
このため、領域分割を簡素化し、処理速度を速めることができる。
【０４５６】
複数の横分割線１４３は、シグナル強度が増大する方向へ等比間隔で設定された。
【０４５７】
このため、ピーク密度の高い部分で領域を細分化することができ、領域分割によるピー
ク存在率の算出を効率的に行わせることができる。
【０４５８】
多成分物質の評価方法は、さらに前記基準ＦＰ作成工程１７５と、基準ＦＰピーク帰属
工程１５９と、基準ＦＰ帰属結果統合工程１６１と、基準ＦＰピーク特徴量作成工程１６
３と、基準ＦＰタイプ２作成工程１６５と、基準ＦＰ領域分割特徴量作成工程１６７と、
基準ＦＰ特徴量統合工程１６９とを備えた。
【０４５９】
このため、基準ＦＰピーク特徴量と基準ＦＰ領域分割特徴量とを統合した基準ＦＰ統合
特徴量を作成し、評価工程１７１です対象ＦＰ統合特徴量と比較することができ、評価対
象薬剤の品質評価の精度及び効率をより向上させることができる。
【０４６０】
基準ＦＰ領域分割特徴量作成工程１６７は、各領域の位置を変更し変更前後で基準ＦＰ
領域分割特徴量を作成することができる。
【０４６１】
このため、分析条件の僅かなバラツキ等により、リテンション・タイムやピーク高さが
変動し、単一パターンでは各格子１４５内の値が大きく変動する場合でも、この変動に係
わらず各格子１４５内のピークの存在量を捉えることができ、評価対象薬剤の品質評価の
精度及び効率をより向上させることができる。
【０４６２】
基準ＦＰ領域分割特徴量作成工程１６７は、シグナル強度軸に平行な複数の縦分割線
１４１と時間軸に平行な複数の横分割線１４３とにより前記領域の分割を行う。
【０４６３】
このため、領域分割を簡素化し、処理速度を速めることができる。
【０４６４】
複数の横分割線１４３は、シグナル強度が増大する方向へ等比間隔で設定された。
【０４６５】
このため、ピーク密度の高い部分で領域を細分化することができ、領域分割によるピー
ク存在率の算出を効率的に行わせることができる。
【０４６６】
基準ＦＰ領域分割特徴量作成工程１６７は、各縦・横分割線１４１、１４３を設定範囲
内で平行移動させるように位置を変更設定することで各領域１４５の位置を変更する。
【０４６７】
このため、各領域１４５の位置変更を簡単な処理で効率的に行わせることができる。
【０４６８】
本発明実施例の多成分薬剤の評価プログラムは、各機能をコンピュータに実現させ、評
価の精度及び効率をより向上させることができる。
【０４６９】
本発明実施例の多成分薬剤の評価装置は、各部３、５、７、９、１１、１３、１５、１
７、１９、２１、２３、２５、２７を作用させ、評価の精度及び効率をより向上させるこ
とができる。
［ピーク・パターンの一致度計算（Ｐ＿Ｓｉｍ）の変形例］
図６３、図６４、図１０５でのピーク・パターンの一致度計算（Ｐ＿Ｓｉｍ）は、ＦＰをピーク高さで作成した上記実施例の場合について適用し、比較対象のピーク高さの差に基づいて計算した。
【０４７０】
一方、本発明のパターン又はＦＰの特徴量作成方法、作成プログラム、及び作成装置おけるピークについては、上記のようにシグナル強度（高さ）の極大値を意味する場合と、シグナル強度の面積値（ピーク面積）を高さで表現したものを意味する場合の何れも含めることができる。
この場合、ＦＰをピーク面積で作成するときも、面積値を高さで表現してＦＰを作成するため、ＦＰとしては上記実施例のピーク高さで作成する場合と同様の表現となる。このため、ＦＰをピーク高さで作成した場合と同様の上記実施例の処理により同様に評価することができる。
【０４７１】
但し、ＦＰをピーク面積で作成したときは、比較対象のピーク値の差が大きくなるため、比に基づいた計算とし、取り扱いを容易とするのが適している。
【０４７２】
以下に、比に基づいて計算するピーク・パターンの一致度（Ｐ＿Ｓｉｍ）を、ｎ＝２、ｎ＝４の場合を例として示す。
ｎ＝２の場合
Ｐ＿Ｓｉｍ＝ (p1／p2^♯１)×(|(r1−(r2＋d)|＋1)
＋(dn1／fn1^♯１)×(|(cn1−r1)−(en1−r2)|+1)
＋(dn2／fn2^♯１)×(|(cn2−r1)−(en2−r2)|+1)
ｎ＝４の場合
Ｐ＿Ｓｉｍ＝ (p1／p2^♯１)×(|(r1−(r2＋d)|＋1)
＋(dn1／fn1^♯１)×(|(cn1−r1)−(en1−r2)|+1)
＋(dn2／fn2^♯１)×(|(cn2−r1)−(en2−r2)|+1)
＋(dn3／fn3^♯１)×(|(cn3−r1)−(en3−r2)|+1)
＋(dn4／fn4^♯１)×(|(cn4−r1)−(en4−r2)|+1)
ここに、^#1は、比較対象の2つの値の比（大きい値/小さい値） であることを示している。
【０４７３】
なお、ＦＰをピーク高さで作成した場合にも、比に基づいてピーク・パターンの一致度（Ｐ＿Ｓｉｍ）を計算することができ、ＦＰをピーク面積で作成した場合にも、上記ピーク高さの差同様に、ピーク面積値の差に基づいてピーク・パターンの一致度（Ｐ＿Ｓｉｍ）を得ることができる。
［サブルーチン２の変形例］
図１３０は、図１０４に代えて適用するサブルーチン２の変形例に係り、図９９の「対象ＦＰ帰属処理２」における「サブルーチン２」の変形例の詳細を示すフローチャートである。この変形例に係る処理により、ＵＶスペクトルの一致度を計算する。
【０４７４】
本サブルーチン２の変形例では、図１０４のサブルーチン２でのＲＭＳＤにＵＶパターンの移動平均の傾き情報（ＤＮＳ）を加味する処理ができるようにした。ＤＮＳは、後述の式で表わされ、ＵＶパターンにおける移動平均値の移動傾きを２つのパターンで比較した時の傾き符号（ ＋／− ）の不一致数として定義される。すなわち、ＤＮＳは、ＵＶパターンの極大、極小値の位置の一致具合を評価する値である。
このＤＮＳの情報を前記ＲＭＳＤに加味することで、ＵＶスペクトルの波形の一致度をより正確に算出することができる。
【０４７５】
図１３０の変形例に係るサブルーチン２において、ステップＳ２００１〜Ｓ２００８までは、図１０４のサブルーチン２とほぼ同一である。但し、ステップＳ２００１では、区間１←ｗ１、区間２←ｗ２の初期設定が追加して行われ、後述する移動平均、移動傾きの計算の区間に用いられる。
【０４７６】
本変形例のサブルーチン２では、ＤＮＳ加味のためにステップＳ２０１０〜Ｓ２０１３を追加し、ステップＳ２００９ＡでＤＮＳを加味した一致度の計算を可能とした。
【０４７７】
ステップＳ２０１０では、「ＤＮＳを加味する？」の判断処理が実行され、ＤＮＳを加味すると判断されたときは（ＹＥＳ）、ステップＳ２０１１へ移行し、ＤＮＳを加味しないと判断されたときは（ＮＯ）、ステップＳ２００９Ａへ移行する。ＤＮＳを加味するか否かの起因は、例えば初期の設定による。例えば、ＦＰを、ピーク面積で作成したときはＤＮＳ加味、ＦＰをピーク高さで作成したときは、ＤＮＳ非加味に設定する。
但し、ＦＰをピーク高さで作成した上記実施例の場合にも、ＤＮＳを加味する処理でＵＶパターン一致度を計算することができ、ＦＰをピーク面積で作成した場合にも、ＤＮＳを加味しない上記実施例の処理でＵＶパターン一致度を計算することができる。
【０４７８】
ステップＳ２０１１では、「区間１（ｗ１）でｘとｙの移動平均を計算」の処理が実行され、区間１（ｗ１）における移動平均が求められる。区間１（ｗ１）は、ＵＶデータの波長に関する区間であり、ステップＳ２００１の初期設定においてｗ１＝３であれば区間１（３）となり、３個の波長におけるＵＶ強度の平均が求められる。具体的には図１３１の図表で後述する。
ステップＳ２０１２では、「区間２（ｗ２）でｘとｙの移動傾きを計算」の処理が実行され、区間２（ｗ２）における移動傾きが求められる。区間２（ｗ２）は、ステップＳ２０１１で求めた移動平均に関する区間であり、ステップＳ２００１の初期設定においてｗ２＝３であれば区間２（３）となり、ステップＳ２０１１で計算した移動平均に基づき、３個の移動平均に渡る傾きの（±）が求められる。具体的には図１０１の図表で後述する。
【０４７９】
ステップＳ２０１３では、「ｘとｙの移動傾きの符号の不一致数を算出（ＤＮＳ）」の処理が実行され、ステップＳ２０１２で計算された移動傾きから傾きの（±）の一致数が計算される。移動傾きの＋は、図６６において右肩上がりを表し、移動傾きの−は、同右肩下がりを表す。
【０４８０】
ステップＳ２０１３からステップＳ２００９Ａへ移行すると、このステップＳ２００９Ａの処理においてＤＮＳを加味した一致度の計算が行われる。
【０４８１】
ステップＳ２００９Ａでは、「ｘとｙのＵＶスペクトルの一致度を計算（ＵＶ＿Ｓｉｍ）
」の処理が実行され、ＤＮＳを加味した一致度の計算では、ＵＶ＿Ｓｉｍを、前記ＵＶスペクトル間距離の二乗和ｚとｘのデータ数ａとＤＮＳとから、
ＵＶ＿Ｓｉｍ ＝ √（ｚ／ａ）×１．１^ＤＮＳ
として算出し、このＵＶ＿Ｓｉｍを図８１のステップＳ３０６に渡し、ＵＶスペクトルの
一致度計算処理を終了する。
【０４８２】
なお、ステップＳ２０１０からステップＳ２００９Ａへ移行した場合の処理は、図８６のステップＳ２００９と同一である。
【０４８３】
図１３１は、移動平均及び移動傾きの計算例を示す図表である。
【０４８４】
図１３１の上段は、ＵＶデータ例、中段は、移動平均の計算例、下段は、移動傾きの計算例を示す。ＵＶデータ例は、具体的な数値に代え、ＵＶ強度をａ１〜ａ７で表記している。例えば、２２０ｎｍのＵＶ強度がａ１、２２１ｎｍのＵＶ強度がａ２等となる。移動平均の計算例及び移動傾きの計算例も具体的数値に代え、ＵＶ強度ａ１〜ａ７を使用している。
【０４８５】
移動平均は、区間１（ｗ１＝３）を例とし、ステップＳ２０１２（図１３０）において、区間（ａ１、ａ２、ａ３）、区間（ａ２、ａ３、ａ４）、・・・毎に計算した値としてｍ１、ｍ２、・・・が計算される。移動傾きも区間２（３）を例とし、ステップＳ２０１３（図１３０）において、区間（ｍ１、ｍ２、ｍ３）、区間（ｍ２、ｍ３、ｍ４）、・・・毎に計算した値としてｓ１、・・・が計算される。例えば、移動平均の差ｍ３−ｍ１が移動傾きとなり、その（±）を取り出す。
【０４８６】
こうして、ＦＰをピーク面積で作成したときは、基準群ＦＰへの帰属処理ならびに基準ＦＰ帰属結果統合処理において、ＤＮＳを加味した処理でＵＶパターン一致度を計算することができる。この計算により、図６６で示す対応する２点の距離（ｄｉｓ）がピーク高さで作成したＦＰに比較して大きくなっても、取り扱いを容易とし、ＵＶパターン一致度を正確に計算することができる。
［その他］
本発明実施例のパターン又はＦＰの特徴量作成方法、作成プログラム、及び作成装置では、ＦＰをピーク面積で作成するときは、シグナル強度軸を面積値軸、シグナル強度を面積値として同様に適用することができる。
【０４８７】
本発明実施例は、多成分薬剤として漢方薬の評価について適用したが、その他の多成分
物質の評価にも適用することができる。
【０４８８】
本実施例では、対象ＦＰタイプ２又は基準ＦＰタイプ２について領域分割特徴量を作成
したが、対象ＦＰ、基準ＦＰについて領域分割特徴量を作成することも可能である。
【０４８９】
また、時系列でピークが変化するパターンを複数の領域に分割し各領域に存在するピー
クの存在率又は存在量からパターン領域分割特徴量を作成する対象パターン領域分割特徴
量作成工程を備えるものであれば、広く適用することができる。
【０４９０】
上記実施例のＦＰは、３Ｄクロマト上での全ピークを対象としたが、細かいデータ、例えば３Ｄクロマト上でピーク面積が５％未満のピークを除いてＦＰを作成することもできる。
【０４９１】
上記実施例のＦＰは、ピーク高さに基づいて作成し、図８７〜図９１の評価を得たが、ピーク面積に基づいてＦＰを作成した場合についても、ピーク高さに基づいて作成した上記実施例同様の手順によりＭＴ法によりＭＤ値を求め、図８７〜図９１と同様に評価を得ることができる。
【０４９２】
クロマトは、３Ｄクロマトに限らず、ＦＰとしてＵＶスペクトルを除いたピークとそのリテンション・タイムとで構成されたものを用いることもできる。この場合、ＵＶスペクトルの一致度を除き、上記実施例と同様に行わせることができる。
【符号の説明】
【０４９３】
１ 多成分薬剤の評価装置（パターンの特徴量作成装置、ＦＰの特徴量作成装置）
３ ＦＰ作成部
５ 対象ＦＰピーク帰属部
７ 対象ＦＰピーク特徴量作成部
９ 対象ＦＰタイプ２作成部
１１ 対象ＦＰ領域分割特徴量作成部
１３ 対象ＦＰ特徴量統合部
１５ 基準ＦＰピーク帰属部
１７ 基準ＦＰ帰属結果統合部
１９ 基準ＦＰピーク特徴量作成部
２１ 基準ＦＰタイプ２作成部
２３ 基準ＦＰ領域分割特徴量作成部
２５ 基準ＦＰ特徴量統合部
２７ 評価部
３１ 基準ＦＰ作成部
３３ 基準ＦＰ選定部
３５ ピーク・パターン作成部
３７ ピーク帰属部
３９ 漢方薬
４１ ３Ｄクロマト
４２ 対象ＦＰに含まれるピークのＵＶスペクトル
４３ 対象ＦＰ
４５ 基準群ＦＰ
４７ 基準群ＦＰに帰属した対象ＦＰ
４９ 対象ＦＰタイプ２（ＦＰ）
５１ 対象ＦＰ領域分割特徴量
５３ 対象ＦＰの評価結果
５５ 薬剤ＡのＦＰ
５７ 薬剤ＢのＦＰ
５９ 薬剤ＣのＦＰ
６１ 対象ＦＰ（リテンション・タイム１０．０−１４．５分）
６３、６５、６７、６９、７１、７３、７５、７７、７９、８１ 対象ＦＰ（リテンショ
ン・タイム１０．０−１４．５分）中の各ピーク
８３ 基準ＦＰ（リテンション・タイム１０．０−１４．５分）
８５、８７、８９、９１、９３、９５、９７、９９、１０１、１０３、１０５ 基準ＦＰ
（リテンション・タイム１０．０−１４．５分）中の各ピーク
１０７ 対象ＦＰリテンション・タイム・出現パターン
１０９ 基準ＦＰリテンション・タイム・出現パターン
１１１ リテンション・タイム・出現距離の一致数
１１３ リテンション・タイム・出現パターンの一致度
１１５ 対象ＦＰ帰属対象ピークのピーク・パターン（３本）
１１７、１１９、１２１、１２３ 基準ＦＰ帰属候補ピークのピーク・パターン（３本）
１２５ 対象ＦＰ帰属対象ピークのピーク・パターン（５本）
１２７、１２９、１３１、１３３ 基準ＦＰ帰属候補ピークのピーク・パターン（５本）
１３５ 帰属対象ピークのＵＶスペクトル
１３９ 帰属候補ピークのＵＶスペクトル
１４１ 縦領域分割線
１４３ 横領域分割線
１４５ 縦・横領域分割線により分割した各領域（格子）
１４７ 各領域をピーク高さ
１４８ ＦＰ作成工程
１４９ 対象ＦＰピーク帰属工程
１５１ 対象ＦＰピーク特徴量作成工程
１５３ 対象ＦＰタイプ２作成工程
１５５ 対象ＦＰ領域分割特徴量作成工程（パターン領域分割特徴量作成工程
、ＦＰ領域分割特徴量作成工程）
１５７ 対象ＦＰ特徴量統合工程
１５９ 基準ＦＰピーク帰属工程
１６１ 基準ＦＰ帰属結果統合工程
１６３ 基準ＦＰピーク特徴量作成工程
１６５ 基準ＦＰタイプ２作成工程
１６７ 基準ＦＰ領域分割特徴量作成工程（パターン領域分割特徴量作成工程
、ＦＰ領域分割特徴量作成工程）
１６９ 基準ＦＰ特徴量統合工程
１７１ 評価工程
１７３ 対象ＦＰ作成工程
１７５ 基準ＦＰ作成工程
１７７ 基準ＦＰ選定工程
１７９ ピーク・パターン作成工程
１８１ ピーク帰属工程
１８３ ３Ｄクロマト・データ例
１８５ ピーク情報データ例
１８７ ＦＰデータ例
１８９ 判定結果ファイル例
１９１ 帰属候補ピークスコア表例
１９３ 帰属候補ピーク番号表例
１９５ 照合結果ファイル例
１９７ 基準群ＦＰデータ例
１９９ 対象ＦＰピーク特徴量ファイル例
２０１ ＦＰタイプ２データ例
２０３ 対象ＦＰ領域分割特徴量ファイル例
２０５ 対象ＦＰ統合特徴量ファイル例
２０７ 基準ｔｙｐｅ２群ＦＰデータ例
２０９ 基準群統合データ例

【特許請求の範囲】
【請求項１】
時系列でピークが変化するパターンを複数の領域に分割し各領域に存在するピークの存
在率又は存在量からパターン領域分割特徴量を作成するパターン領域分割特徴量作成工程、
を備えたことを特徴とするパターンの特徴量作成方法。
【請求項２】
多成分物質のクロマトから検出されたピークとそのリテンション・タイムとで構成され
るＦＰを複数の領域に分割し各領域に存在するピークの存在率又は存在量からＦＰ領域分
割特徴量を作成するＦＰ領域分割特徴量作成工程、
を備えたことを特徴とするＦＰの特徴量作成方法。
【請求項３】
請求項２記載のＦＰの特徴量作成方法であって、
前記多成分物質は、多成分薬剤である、
ことを特徴とするＦＰの特徴量作成方法。
【請求項４】
請求項３記載のＦＰの特徴量作成方法であって、
前記多成分薬剤は、生薬、生薬の組合せ、それらの抽出物、漢方薬の何れかである、
ことを特徴とするＦＰの特徴量作成方法。
【請求項５】
請求項２〜４の何れかに記載のＦＰの特徴量作成方法であって、
前記ＦＰ領域分割特徴量作成工程は、シグナル強度軸又は面積値軸に平行な複数の縦分割線と時間軸に平行な複数の横分割線とにより前記領域の分割を行う、
ことを特徴とするＦＰの特徴量作成方法。
【請求項６】
請求項５記載のＦＰの特徴量作成方法であって、
前記複数の横分割線は、シグナル強度又は面積値が増大する方向へ等比間隔で設定された、
ことを特徴とするＦＰの特徴量作成方法。
【請求項７】
請求項２〜６の何れかに記載のＦＰの特徴量作成方法であって、
前記ＦＰ領域分割特徴量作成工程は、前記各領域の位置を変更し変更前後で前記ＦＰ領
域分割特徴量を作成する、
ことを特徴とするＦＰの特徴量作成方法。
【請求項８】
請求項７記載のＦＰの特徴量作成方法であって、
前記ＦＰ領域分割特徴量作成工程は、前記各縦・横分割線を設定範囲内で平行移動させ
るように位置を変更設定することで前記各領域の位置を変更する、
ことを特徴とするＦＰの特徴量作成方法。
【請求項９】
時系列でピークが変化するパターンを複数の領域に分割し各領域に存在するピークの存
在率又は存在量からパターン領域分割特徴量を作成するパターン領域分割特徴量作成機能
、
をコンピュータに実現させることを特徴とするパターンの特徴量作成プログラム。
【請求項１０】
多成分物質のクロマトから検出されたピークとそのリテンション・タイムとで構成され
るＦＰを複数の領域に分割し各領域に存在するピークの存在率又は存在量からＦＰ領域分
割特徴量を作成するＦＰ領域分割特徴量作成機能、
を備えたことを特徴とするＦＰの特徴量作成プログラム。
【請求項１１】
請求項１０記載のＦＰの特徴量作成プログラムであって、
前記多成分物質は、多成分薬剤である、
ことを特徴とするＦＰの特徴量作成プログラム。
【請求項１２】
請求項１１記載のＦＰの特徴量作成プログラムであって、
前記多成分薬剤は、生薬、生薬の組合せ、それらの抽出物、漢方薬の何れかである、
ことを特徴とするＦＰの特徴量作成プログラム。
【請求項１３】
請求項１０〜１２の何れかに記載のＦＰの特徴量作成プログラムであって、
前記ＦＰ領域分割特徴量作成機能は、シグナル強度軸又は面積値軸に平行な複数の縦分割線と時間軸に平行な複数の横分割線とにより前記領域の分割を行う、
ことを特徴とするＦＰの特徴量作成プログラム。
【請求項１４】
請求項１３記載のＦＰの特徴量作成プログラムであって、
前記複数の横分割線は、シグナル強度又は面積値が増大する方向へ等比間隔で設定された、
ことを特徴とするＦＰの特徴量作成プログラム。
【請求項１５】
請求項１０〜１４の何れかに記載のＦＰの特徴量作成プログラムであって、
前記ＦＰ領域分割特徴量作成機能は、前記各領域の位置を変更し変更前後で前記ＦＰ領
域分割特徴量を作成する、
ことを特徴とするＦＰの特徴量作成プログラム。
【請求項１６】
請求項１５記載のＦＰの特徴量作成プログラムであって、
前記ＦＰ領域分割特徴量作成機能は、前記各縦・横分割線を設定範囲内で平行移動させ
るように位置を変更設定することで前記各領域の位置を変更する、
ことを特徴とするＦＰの特徴量作成プログラム。
【請求項１７】
時系列でピークが変化するパターンを複数の領域に分割し各領域に存在するピークの存
在率又は存在量からパターン領域分割特徴量を作成するパターン領域分割特徴量作成部、
を備えたことを特徴とするパターンの特徴量作成装置。
【請求項１８】
多成分物質のクロマトから検出されたピークとそのリテンション・タイムとで構成され
るＦＰを複数の領域に分割し各領域に存在するピークの存在率又は存在量からＦＰ領域分
割特徴量を作成するＦＰ領域分割特徴量作成部、
を備えたことを特徴とするＦＰの特徴量作成装置。
【請求項１９】
請求項１８記載のＦＰの特徴量作成装置であって、
前記多成分物質は、多成分薬剤である、
ことを特徴とするＦＰの特徴量作成装置。
【請求項２０】
請求項１９記載のＦＰの特徴量作成装置であって、
前記多成分薬剤は、生薬、生薬の組合せ、それらの抽出物、漢方薬の何れかである、
ことを特徴とするＦＰの特徴量作成装置。
【請求項２１】
請求項1８〜２０の何れかに記載のＦＰの特徴量作成装置であって、
前記ＦＰ領域分割特徴量作成部は、シグナル強度軸又は面積値軸に平行な複数の縦分割線と時間軸に平行な複数の横分割線とにより前記領域の分割を行う、
ことを特徴とするＦＰの特徴量作成装置。
【請求項２２】
請求項２１記載のＦＰの特徴量作成装置であって、
前記複数の横分割線は、シグナル強度又は面積値が増大する方向へ等比間隔で設定された、
ことを特徴とするＦＰの特徴量作成装置。
【請求項２３】
請求項１８〜２２の何れかに記載のＦＰの特徴量作成装置であって、
前記ＦＰ領域分割特徴量作成部は、前記各領域の位置を変更し変更前後で前記ＦＰ領域
分割特徴量を作成する、
ことを特徴とするＦＰの特徴量作成装置。
【請求項２４】
請求項２３記載のＦＰの特徴量作成装置であって、
前記ＦＰ領域分割特徴量作成部は、前記各縦・横分割線を設定範囲内で平行移動させる
ように位置を変更設定することで前記各領域の位置を変更する、
ことを特徴とするＦＰの特徴量作成装置。

【図１】

【図４】

【図５】

【図６】

【図７】

【図８】

【図９】

【図１０】

【図１１】

【図１２】

【図１３】

【図１４】

【図１５】

【図１６】

【図１７】

【図１８】

【図１９】

【図２０】

【図２１】

【図２２】

【図２３】

【図２４】

【図２５】

【図２６】

【図２７】

【図２８】

【図２９】

【図３０】

【図３１】

【図３２】

【図３３】

【図３４】

【図３５】

【図３６】

【図３７】

【図３８】

【図３９】

【図４０】

【図４１】

【図４２】

【図４３】

【図４４】

【図４５】

【図４６】

【図４７】

【図４８】

【図４９】

【図５０】

【図５１】

【図５２】

【図５３】

【図５４】

【図５５】

【図５６】

【図５７】

【図５８】

【図５９】

【図６０】

【図６１】

【図６２】

【図６３】

【図６４】

【図６５】

【図６６】

【図６７】

【図７０】

【図７１】

【図７２】

【図７３】

【図７４】

【図７５】

【図７６】

【図７７】

【図７８】

【図７９】

【図８０】

【図８１】

【図８２】

【図８３】

【図８４】

【図８５】

【図８６】

【図８７】

【図８８】

【図８９】

【図９０】

【図９１】

【図９２】

【図９３】

【図９４】

【図９５】

【図９６】

【図９７】

【図９８】

【図９９】

【図１００】

【図１０１】

【図１０２】

【図１０３】

【図１０４】

【図１０５】

【図１０６】

【図１０７】

【図１０８】

【図１０９】

【図１１０】

【図１１１】

【図１１２】

【図１１３】

【図１１４】

【図１１５】

【図１１６】

【図１１７】

【図１１８】

【図１１９】

【図１２０】

【図１２１】

【図１２２】

【図１２３】

【図１２４】

【図１２５】

【図１２６】

【図１２７】

【図１２８】

【図１２９】

【図１３０】

【図１３１】

【図２】

【図３】

【図６８】

【図６９】

【公開番号】特開２０１３−１１５９８（Ｐ２０１３−１１５９８Ａ）
【公開日】平成２５年１月１７日（２０１３．１．１７）
【国際特許分類】
【出願番号】特願２０１２−１２５３８２（Ｐ２０１２−１２５３８２）
【出願日】平成２４年５月３１日（２０１２．５．３１）
【公序良俗違反の表示】
（特許庁注：以下のものは登録商標）
１．ＷＩＮＤＯＷＳ
【出願人】（５０５４２０５４１）アングルトライ株式会社 (6)
【出願人】（０００００３６６５）株式会社ツムラ (43)