WO2012164956A1

WO2012164956A1 - パターン又はｆｐの特徴量作成方法、作成プログラム、及び作成装置

Info

Publication number: WO2012164956A1
Application number: PCT/JP2012/003618
Authority: WO
Inventors: 手島　昌一; 芳和森; 桂一野田
Original assignee: 株式会社ツムラ
Priority date: 2011-06-01
Filing date: 2012-05-31
Publication date: 2012-12-06
Also published as: EP2717045B1; EP2717045A1; CN102959395A; JP5912880B2; US10605792B2; HK1181119A1; KR101442117B1; US20170007502A1; KR20130029406A; TW201312109A; CN102959395B; EP2717045A4; JP2013011598A; TWI561818B; US20130204539A1

Abstract

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

Description

パターン又はＦＰの特徴量作成方法、作成プログラム、及び作成装置

　本発明は、パターンの特徴量作成方法、多成分物質、例えば多成分薬剤である漢方薬の品質評価を行うための多成分物質のＦＰの特徴量作成方法、作成プログラム、及び作成装置に関する。

　多成分物質として、例えば多成分で構成される薬剤（以下、多成分薬剤）である漢方薬などの天然物由来の薬剤がある。これら薬剤の定量的、定性的プロフィールは、使用する原料生薬の、地質学的要因、生態学的要因、収集時期、収集場所、収集年代、生育期の天候等が原因で変化する。

　そのため、これら多成分薬剤などについては、その安全性および有効性を担保するための品質として一定の基準を規定し、その基準に基づいて、国家の監督機関、化学的組織、製造業者等が品質評価を行っている。

　しかしながら、多成分薬剤の品質等の判定基準は、多成分薬剤中のある特徴的な１ないしは数成分を選択して、その含有量などに基づいて設定することが一般的であった。

　例えば、非特許文献１では、多成分薬剤において有効成分の同定ができていない場合、定量分析が可能、水に溶けやすい、熱水中で分解しない、他の成分と化学反応をしない等の物性を持った複数の成分を選択し、化学分析により得られるそれら成分の含有量を評価の基準としている。

　また、多成分薬剤にクロマトグラフィーを適用し、保持時間ごとに紫外可視吸収スペクトルを得て、その中の一部の成分情報から評価の基準を設定することも知られている。

　例えば、特許文献１では、ＨＰＬＣクロマトグラム・データ中の一部のピークを選択し、バーコード化することによって多成分薬剤を評価している。

　しかしながら、これらの方法は、評価の対象が「特定成分の含有量」又は「特定成分のクロマト・ピーク」に限定されており、多成分薬剤が含有する成分の一部が評価対象になっているに過ぎない。このため、多成分薬剤については、評価の対象外となっている成分が多数存在することから、多成分薬剤の評価方法として精度的に不十分である。

　多成分薬剤の品質を正確に評価するには、全ピーク情報或いは数％の細かい情報を除外した全ピーク情報に近いピーク情報を網羅した評価が必要であり、そのため、多成分薬剤間で全ピーク或いはこれに近い個数のピークを対応させる必要がある。

　ところが、複数のピークを高精度で効率的に対応させることが困難であり、このことが多成分薬剤の高精度で効率的な評価の妨げとなっていた。

　さらに説明すると、同じ製品名の多成分薬剤であっても生薬は天然物であるがために構成成分が僅かに異なる場合がある。このため、同じ品質の薬剤であっても構成成分の含量比率が異なったり、ある薬剤に存在する成分が他の薬剤には存在しないこと（以下、薬剤間誤差）がある。さらに、クロマトにおけるピーク強度やピークの溶出時間には厳密な再現性がない（以下、分析誤差）などの要因もある。これらにより、多成分薬剤間の全ピーク或いはこれに近い個数のピークについて同じ成分に由来するピークに対応させること（以下、ピーク帰属）ができないことが原因となって高精度で効率的な評価の妨げとなっていた。

特開２００２－２１４２１５号公報

月刊薬事ｖｏｌ．２８，Ｎｏ．３，６７-７１ (１９８６)

　解決しようとする問題点は、既存の評価方法では多成分物質の品質等を高精度で効率的に評価することに限界があった点である。

　本発明は、評価の精度及び効率を向上させることに寄与可能とするために、時系列でピークが変化するパターンを複数の領域に分割し各領域に存在するピークの存在率又は存在量からパターン領域分割特徴量を作成するパターン領域分割特徴量作成工程を備えたことをパターンの特徴量作成方法の特徴とする。

　本発明は、多成分物質のクロマトから検出されたピークとそのリテンション・タイムとで構成されるＦＰを複数の領域に分割し各領域に存在するピークの存在率又は存在量からＦＰ領域分割特徴量を作成するＦＰ領域分割特徴量作成工程を備えたことをＦＰの特徴量作成方法の特徴とする。

　本発明は、時系列でピークが変化するパターンを複数の領域に分割し各領域に存在するピークの存在率又は存在量からパターン領域分割特徴量を作成するパターン領域分割特徴量作成機能をコンピュータに実現させることをパターンの特徴量作成プログラムの特徴とする。

　本発明は、多成分物質のクロマトから検出されたピークとそのリテンション・タイムとで構成されるＦＰを複数の領域に分割し各領域に存在するピークの存在率又は存在量からＦＰ領域分割特徴量を作成するＦＰ領域分割特徴量作成機能を備えたことをＦＰの特徴量作成プログラムの特徴とする。

　本発明は、時系列でピークが変化するパターンを複数の領域に分割し各領域に存在するピークの存在率又は存在量からパターン領域分割特徴量を作成するパターン領域分割特徴量作成部を備えたことをパターンの特徴量作成装置の特徴とする。本発明は、多成分物質のクロマトから検出されたピークとそのリテンション・タイムとで構成されるＦＰを複数の領域に分割し各領域に存在するピークの存在率又は存在量からＦＰ領域分割特徴量を作成するＦＰ領域分割特徴量作成部を備えたことをＦＰの特徴量作成装置の特徴とする。

　本発明のパターン又はＦＰの特徴量作成方法は、上記構成であるから、パターン又はＦＰの特徴量を領域分割により簡単に得ることが可能となる。したがって、例えば細かいピークをも捉えて特徴量を作成することができる。

　本発明のパターン又はＦＰの特徴量作成プログラムは、上記構成であるから、各機能をコンピュータに実現させ、パターン又はＦＰの特徴量を簡単に得ることが可能となる。本発明のパターン又はＦＰの特徴量作成装置は、上記構成であるから、各部を作用させ、パターン又はＦＰの特徴量を簡単に得ることが可能となる。

多成分薬剤の評価装置のブロック図である。（実施例１）多成分薬剤の評価手順を示すブロック図である。（実施例１）三次元クロマトグラム・データ（以下、３Ｄクロマト）から作成したＦＰの説明図である。（実施例１）（Ａ）は、薬剤Ａ、（Ｂ）は、薬剤Ｂ、（Ｃ）は、薬剤ＣのＦＰである。（実施例１）対象ＦＰ及び基準ＦＰのリテンション・タイムを示す図である。（実施例１）対象ＦＰのリテンション・タイム・出現パターンを示す図である。（実施例１）基準ＦＰのリテンション・タイム・出現パターンを示す図である。（実施例１）対象ＦＰと基準ＦＰのリテンション・タイム・出現距離の一致数を示す図である。（実施例１）対象ＦＰと基準ＦＰのリテンション・タイム・出現パターンの一致度を示す図である。（実施例１）対象ＦＰの帰属対象ピークを示す図である。（実施例１）帰属対象ピークを含めたピーク３本によるピーク・パターン図である。（実施例１）帰属対象ピークを含めたピーク５本によるピーク・パターン図である。（実施例１）帰属対象ピークの許容幅を示す図である。（実施例１）帰属対象ピークに対する基準ＦＰの帰属候補ピークを示す図である。（実施例１）帰属対象ピークと帰属候補ピークのピーク３本によるピーク・パターン図である。（実施例１）帰属対象ピークと別な帰属候補ピークのピーク３本によるピーク・パターン図である。（実施例１）帰属対象ピークと別な帰属候補ピークのピーク３本によるピーク・パターン図である。（実施例１））帰属対象ピークと別な帰属候補ピークのピーク３本によるピーク・パターン図である。（実施例１）帰属対象ピークと帰属候補ピークのピーク５本によるピーク・パターン図である。（実施例１）帰属対象ピークと別な帰属候補ピークのピーク５本によるピーク・パターン図である。（実施例１）帰属対象ピークと別な帰属候補ピークのピーク５本によるピーク・パターン図である。（実施例１））帰属対象ピークと別な帰属候補ピークのピーク５本によるピーク・パターン図である。（実施例１）帰属対象ピークと帰属候補ピークのピーク・パターン構成候補ピークを示す図である。（実施例１）ピーク・パターン構成候補ピークを４本とした場合の帰属対象ピークの全ピーク・パターン数を示す図である。（実施例１）ピーク・パターン構成候補ピークを４本とした場合の帰属候補ピークの全ピーク・パターン数を示す図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１）帰属対象ピークと帰属候補ピークのピーク３本によるピーク・パターンの一致度の算出方法を示す図である。（実施例１）帰属対象ピークと帰属候補ピークのピーク３本によるピーク・パターンの一致度の算出方法を示す図である。（実施例１）帰属対象ピークと帰属候補ピークのピーク５本によるピーク・パターンの一致度の算出方法を示す図である。（実施例１）帰属対象ピーク及び帰属候補ピークのＵＶスペクトルを示す図である。（実施例１）帰属対象ピークと帰属候補ピークのＵＶスペクトルの一致度の説明図である。（実施例１）ピーク・パターンとＵＶスペクトルの両方の比較による帰属候補ピークの一致度計算の説明図である。（実施例１）対象ＦＰの基準群ＦＰへの帰属を示す説明図である。（実施例１）対象ＦＰが基準群ＦＰへ帰属された状況を示す図である。（実施例１）領域分割による数量化を示す説明図である。（実施例１）リテンション・タイム等の変動との関係を示す説明図である。（実施例１）領域の位置を変更して数量化する説明図である。（実施例１）ＦＰタイプ２のデータを示す図表である。（実施例１）ＦＰタイプ２のパターンを示す説明図である。（実施例１）縦・横分割線での領域分割による領域ごとの特徴量化を示す説明図である。（実施例１）縦分割線（1本目）の設定を示す説明図である。（実施例１）横分割線（1本目）の設定を示す説明図である。（実施例１）縦・横分割線による領域分割を示す説明図である。（実施例１）特徴量化する領域の数を示す説明図である。（実施例１）領域１の特定を示す説明図である。（実施例１）全ピークの高さ及び合計を示す図表である。（実施例１）領域１のピーク高さの合計を示す説明図である。（実施例１）全領域の特徴量を示す図表である。（実施例１）縦１本目の位置を順次変更してできた各領域での特徴量を示す図表である。（実施例１）横１本目の位置を順次変更してできた各領域での特徴量を示す図表である。（実施例１）各縦・横分割線の位置を変更しない１通りでの特徴量を示す図表である。（実施例１）各種対象ＦＰとその評価値（ＭＤ値）を示す図である。（実施例１）各種対象ＦＰとその評価値（ＭＤ値）を示す図である。（実施例１）各種対象ＦＰとその評価値（ＭＤ値）を示す図である。（実施例１）各種対象ＦＰとその評価値（ＭＤ値）を示す図である。（実施例１）各種対象ＦＰとその評価値（ＭＤ値）を示す図である。（実施例１）多成分薬剤の評価方法を示す工程図である。（実施例１）多成分薬剤の品質評価フロー図である。（実施例１）多成分薬剤の品質評価フロー図である。（実施例１）単波長によるＦＰ作成機能におけるデータ処理フローチャートである。（実施例１）複数波長によるＦＰ作成機能におけるデータ処理フローチャートである。（実施例１）複数波長によるＦＰ作成機能におけるデータ処理フローチャートである。（実施例１）ピーク帰属処理１（基準ＦＰの選定）におけるデータ処理フローチャートである。（実施例１）ピーク帰属処理２（帰属スコアの算出）におけるデータ処理フローチャートである。（実施例１）ピーク帰属処理３（対応ピークの特定）におけるデータ処理フローチャートである。（実施例１）ピーク帰属処理４（基準群ＦＰへの帰属）におけるデータ処理フローチャートである。（実施例１）ピーク帰属処理４（基準群ＦＰへの帰属）におけるデータ処理フローチャートである。（実施例１）ピーク帰属処理１（基準ＦＰの選定）におけるリテンション・タイム・出現パターンの一致度計算処理のフローチャートである。（実施例１）ピーク帰属処理２（帰属スコアの算出）におけるＵＶスペクトルの一致度計算処理のフローチャートである。（実施例１）ピーク帰属処理２（帰属スコアの算出）におけるピーク・パターンの一致度計算処理のフローチャートである。（実施例１）「ＦＰ_type２作成」の詳細を示すフローチャートである。（実施例１）「領域分割による対象ＦＰ_type２の特徴量化処理」の詳細を示すフローチャートである。（実施例１）「対象ＦＰのピーク特徴量と領域分割特徴量の統合」の詳細を示すフローチャートである。（実施例１）基準ＦＰ特徴量統合ファイルを作成するためのフローチャートである。（実施例１）基準ＦＰ特徴量統合ファイルを作成するためのフローチャートである。（実施例１）「基準ＦＰ帰属結果統合処理（基準ＦＰ対応表の作成）」の詳細を示すフローチャートである。（実施例１）「基準ＦＰ帰属結果統合処理（基準ＦＰ対応表の作成）」の詳細を示すフローチャートである。（実施例１）「ピーク特徴量化処理（基準群ＦＰの作成）」の詳細を示すフローチャートである。（実施例１）「基準ＦＰ_type２の作成処理」の詳細を示すフローチャートである。（実施例１）「領域分割による基準ＦＰの特徴量化処理」の詳細を示すフローチャートである。（実施例１）基準ＦＰの特徴量統合処理に係るフローチャートである。（実施例１）３Ｄクロマトのデータ例を示す図表である。（実施例１）ピーク情報のデータ例を示す図表である。（実施例１）ＦＰのデータ例を示す図表である。（実施例１）対象ＦＰの基準ＦＰへの帰属スコア計算結果（判定結果ファイル例）を示す図表である。（実施例１）対象ＦＰと基準ＦＰで対応するピークの照合過程を示す図表である。（実施例１）対象ＦＰと基準ＦＰで対応するピークを特定した結果（照合結果ファイル）例を示す図表である。（実施例１）基準群ＦＰのデータ例を示す図表である。（実施例１）対象ＦＰピーク特徴量ファイル例を示す図表である。（実施例１）対象及び基準ＦＰタイプ２のデータ例を示す図表である。（実施例１）対象ＦＰ領域分割特徴量ファイル例を示す図表である。（実施例１）対象ＦＰ統合特徴量ファイル例を示す図表である。（実施例１）基準ｔｙｐｅ２群ＦＰ例を示す図表である。（実施例１）基準群統合データ例を示す図表である。（実施例１）図１０４に代えて適用するサブルーチン２の変形例の詳細を示すフローチャートである。（実施例１）移動平均及び移動傾きの計算例を示す図表である。（実施例１）

　評価の精度及び効率を向上させることに寄与可能にするという目的を、多成分物質のクロマトから検出されたピークとそのリテンション・タイムとで構成されるＦＰを複数の領域に分割し各領域に存在するピークの存在率又は存在量からＦＰ領域分割特徴量を作成することにより実現した。

　本発明の実施例１は、多成分物質、例えば多成分薬剤を評価する多成分薬剤の評価方法、評価プログラム、パターンであるＦＰの特徴量作成方法、作成プログラム、作成装置である。

　多成分薬剤は、複数の有効化学成分を含有する薬剤と定義され、限定はされないが、生薬、生薬の組合せ、それらの抽出物、漢方薬等が含まれる。また剤形も特に限定されず、例えば、第１５改正日本薬局方の製剤総則で規定されている液剤、エキス剤、カプセル剤、顆粒剤、丸剤、懸濁剤・乳剤、散剤、酒精剤、錠剤、浸剤・煎剤、チンキ剤、トローチ剤、芳香水剤、流エキス剤等が含まれる。多成分物質としては、薬剤以外のものも含まれる。

　漢方薬の具体例は、医療用漢方製剤１４８処方「使用上の注意」の業界統一と自主改訂、一般用漢方処方の手引き(１９７８年)に記載されている。

　多成分薬剤の評価では、評価対象薬剤が正常品と定めた複数の薬剤と同等であるかどうかを評価するため、まず、評価対象薬剤の三次元クロマトグラム・データ（以下、３Ｄクロマト）からその薬剤特有の情報を抽出した対象ＦＰを作成する。

　次に、対象ＦＰの各ピークを、全基準ＦＰをピーク帰属処理し作成した全基準ＦＰのピーク対応データ（以下、基準群ＦＰ）に帰属し、ピーク特徴量を得る。

　さらに、対象ＦＰから帰属されたピークを除き、残ったピークでＦＰタイプ２を作成し、このＦＰタイプ２を領域分割し領域分割特徴量を得る。

　これら２つの特徴量を統合し、対象ＦＰ統合特徴量を得る。

　この対象ＦＰ統合特徴量と全基準ＦＰから得た基準ＦＰ統合特徴量により、基準群ＦＰと対象ＦＰの同等性をＭＴ法で評価する。最後に、得られた評価値（以下、ＭＤ値）と予め設定しておいた判定値（ＭＤ値の上限値）を比較し、評価対象薬剤が正常品と同等であるかどうかを判定する。
［多成分薬剤の評価装置］
　図１は、多成分薬剤の評価装置のブロック図、図２は、多成分薬剤の評価手順を示すブロック図、図３は、３Ｄクロマトから作成したＦＰの説明図、図４（Ａ）は、薬剤Ａ、（Ｂ）は、薬剤Ｂ、（Ｃ）は、薬剤ＣのＦＰである。

　図1のように、多成分薬剤の評価装置１は、ＦＰ作成部３と、対象ＦＰピーク帰属部５と、対象ＦＰピーク特徴量作成部７と、対象ＦＰタイプ２作成部９と、対象ＦＰ領域分割特徴量作成部１１と、対象ＦＰ特徴量統合部１３と、基準ＦＰピーク帰属部１５と、基準ＦＰ帰属結果統合部１７と、基準ＦＰピーク特徴量作成部１９と、基準ＦＰタイプ２作成部２１と、基準ＦＰ領域分割特徴量作成部２３と、基準ＦＰ特徴量統合部２５と、評価部２７とを備えている。多成分薬剤の評価装置１は、パターンであるＦＰの特徴量作成装置を含む。

　ＦＰ作成部３は、対象ＦＰ作成部２９と、基準ＦＰ作成部３１とを備えている。

　対象ＦＰピーク帰属部５は、基準ＦＰ選定部３３と、ピーク・パターン作成部３５と、ピーク帰属部３７とを備えている。

　この多成分薬剤の評価装置１は、例えば、一つのコンピュータで構成され、図示はしないが、ＣＰＵ、ＲＯＭ、ＲＡＭなどを備えている。多成分薬剤の評価装置１は、コンピュータにインストールされたパターンの特徴量作成プログラムとしてＦＰの特徴量作成プログラムを実現させ、ＦＰの特徴量を得ることができる。但し、ＦＰの特徴量作成プログラムは、これを記録したＦＰの特徴量作成プログラム記録媒体を用い、コンピュータで構成された多成分薬剤の評価装置１にこれを読み取らせることで、ＦＰの特徴量を得ることを実現させることもできる。

　多成分薬剤の評価装置１は、各部をそれぞれ別々のコンピュータで構成することができ、例えば対象ＦＰピーク帰属部５と、対象ＦＰピーク特徴量作成部７と、対象ＦＰタイプ２作成部９と、対象ＦＰ領域分割特徴量作成部１１と、対象ＦＰ特徴量統合部１３と、評価部２７とを一つのコンピュータで構成し、基準ＦＰ作成部３１と、基準ＦＰピーク帰属部１５と、基準ＦＰ帰属結果統合部１７と、基準ＦＰピーク特徴量作成部１９と、基準ＦＰタイプ２作成部２１と、基準ＦＰ領域分割特徴量作成部２３と、基準ＦＰ特徴量統合部２５を他のコンピュータで構成することなどもできる。

　この場合、基準ＦＰ統合特徴量は、他のコンピュータで作成し、評価装置１の評価部２７に入力されることになる。

　そして、対象ＦＰ作成部２９と、対象ＦＰピーク帰属部５と、対象ＦＰピーク特徴量作成部７と、対象ＦＰタイプ２作成部９と、対象ＦＰ領域分割特徴量作成部１１と、対象ＦＰ特徴量統合部１３とにより対象ＦＰ統合特徴量が作成され、基準ＦＰ作成部３１と、基準ＦＰピーク帰属部１５と、基準ＦＰ帰属結果統合部１７と、基準ＦＰピーク特徴量作成部１９と、基準ＦＰタイプ２作成部２１と、基準ＦＰ領域分割特徴量作成部２３と、基準ＦＰ特徴量統合部２５とにより基準ＦＰ統合特徴量が作成され、これらが比較評価されて対象ＦＰ４３と基準群ＦＰ４５との同等性を評価する。

　ＦＰ作成部３の対象ＦＰ作成部２９は、例えば、図２、図３のように、漢方薬３９のクロマトとして三次元クロマトグラム・データである３Ｄクロマト４１から、特定の検出波長における複数のピークとそのリテンション・タイムならびにＵＶスペクトルを抽出した対象ＦＰ４３（以下単に、「ＦＰ４３」ということもある。）を作成する機能部である。

　このＦＰ４３は、３Ｄクロマト４１と同様に三次元の情報（ピーク、リテンション・タイム及びＵＶスペクトル）で構成している。

　このため、ＦＰ４３は、その薬剤特有の情報をそのまま継承したデータである。それにも係らず、データ容量は約１／７０に圧縮されているため、３Ｄクロマト４１に比較して、処理すべき情報量を大幅に減少させることができ処理速度を速めることができる。

　３Ｄクロマト４１は、漢方薬３９に高速液体クロマトグラフィー（ＨＰＬＣ）を適用した結果である。この３Ｄクロマト４１は、各成分の移動速度として現れ、それを特定時間における移動距離として表し、或いはカラム末端から時系列に現れる様をチャートに表したものである。ＨＰＬＣにおいては、時間軸に対する検出器応答をプロットしたもので、ピークの出現時間を保持時間（リテンション・タイム）と呼んでいる。

　検出器としては、特に限定されないが、光学的性質を利用した吸光度検出器(Absorbance 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次元に集約し、複数の評価対象ロットを簡単に比較評価する。このため、複数の成分で構成される多成分系薬剤の評価に適している。

　対象ＦＰ領域分割特徴量作成工程１５５は、シグナル強度軸に平行な複数の縦分割線１４１と時間軸に平行な複数の横分割線１４３とにより前記領域の分割を行う。

　このため、領域分割を簡素化し、処理速度を速めることができる。

　複数の横分割線１４３は、シグナル強度が増大する方向へ等比間隔で設定された。

　このため、ピーク密度の高い部分で領域を細分化することができ、領域分割によるピーク存在率の算出を効率的に行わせることができる。

　多成分物質の評価方法は、さらに前記基準ＦＰ作成工程１７５と、基準ＦＰピーク帰属工程１５９と、基準ＦＰ帰属結果統合工程１６１と、基準ＦＰピーク特徴量作成工程１６３と、基準ＦＰタイプ２作成工程１６５と、基準ＦＰ領域分割特徴量作成工程１６７と、基準ＦＰ特徴量統合工程１６９とを備えた。

　このため、基準ＦＰピーク特徴量と基準ＦＰ領域分割特徴量とを統合した基準ＦＰ統合特徴量を作成し、評価工程１７１です対象ＦＰ統合特徴量と比較することができ、評価対象薬剤の品質評価の精度及び効率をより向上させることができる。

　基準ＦＰ領域分割特徴量作成工程１６７は、各領域の位置を変更し変更前後で基準ＦＰ領域分割特徴量を作成することができる。

　このため、分析条件の僅かなバラツキ等により、リテンション・タイムやピーク高さが変動し、単一パターンでは各格子１４５内の値が大きく変動する場合でも、この変動に係わらず各格子１４５内のピークの存在量を捉えることができ、評価対象薬剤の品質評価の精度及び効率をより向上させることができる。

　基準ＦＰ領域分割特徴量作成工程１６７は、シグナル強度軸に平行な複数の縦分割線１４１と時間軸に平行な複数の横分割線１４３とにより前記領域の分割を行う。

　基準ＦＰ領域分割特徴量作成工程１６７は、各縦・横分割線１４１、１４３を設定範囲内で平行移動させるように位置を変更設定することで各領域１４５の位置を変更する。

　このため、各領域１４５の位置変更を簡単な処理で効率的に行わせることができる。

　本発明実施例の多成分薬剤の評価プログラムは、各機能をコンピュータに実現させ、評価の精度及び効率をより向上させることができる。

　本発明実施例の多成分薬剤の評価装置は、各部３、５、７、９、１１、１３、１５、１７、１９、２１、２３、２５、２７を作用させ、評価の精度及び効率をより向上させることができる。
［ピーク・パターンの一致度計算（Ｐ＿Ｓｉｍ）の変形例］
　図６３、図６４、図１０５でのピーク・パターンの一致度計算（Ｐ＿Ｓｉｍ）は、ＦＰをピーク高さで作成した上記実施例の場合について適用し、比較対象のピーク高さの差に基づいて計算した。

　一方、本発明のパターン又はＦＰの特徴量作成方法、作成プログラム、及び作成装置おけるピークについては、上記のようにシグナル強度（高さ）の極大値を意味する場合と、シグナル強度の面積値（ピーク面積）を高さで表現したものを意味する場合の何れも含めることができる。

　この場合、ＦＰをピーク面積で作成するときも、面積値を高さで表現してＦＰを作成するため、ＦＰとしては上記実施例のピーク高さで作成する場合と同様の表現となる。このため、ＦＰをピーク高さで作成した場合と同様の上記実施例の処理により同様に評価することができる。

　但し、ＦＰをピーク面積で作成したときは、比較対象のピーク値の差が大きくなるため、比に基づいた計算とし、取り扱いを容易とするのが適している。

　なお、ＦＰをピーク高さで作成した場合にも、比に基づいてピーク・パターンの一致度（Ｐ＿Ｓｉｍ）を計算することができ、ＦＰをピーク面積で作成した場合にも、上記ピーク高さの差同様に、ピーク面積値の差に基づいてピーク・パターンの一致度（Ｐ＿Ｓｉｍ）を得ることができる。
［サブルーチン２の変形例］
　図１３０は、図１０４に代えて適用するサブルーチン２の変形例に係り、図９９の「対象ＦＰ帰属処理２」における「サブルーチン２」の変形例の詳細を示すフローチャートである。この変形例に係る処理により、ＵＶスペクトルの一致度を計算する。

　本サブルーチン２の変形例では、図１０４のサブルーチン２でのＲＭＳＤにＵＶパターンの移動平均の傾き情報（ＤＮＳ）を加味する処理ができるようにした。ＤＮＳは、後述の式で表わされ、ＵＶパターンにおける移動平均値の移動傾きを２つのパターンで比較した時の傾き符号（＋／－）の不一致数として定義される。すなわち、ＤＮＳは、ＵＶパターンの極大、極小値の位置の一致具合を評価する値である。

　このＤＮＳの情報を前記ＲＭＳＤに加味することで、ＵＶスペクトルの波形の一致度をより正確に算出することができる。

　図１３０の変形例に係るサブルーチン２において、ステップＳ２００１～Ｓ２００８までは、図１０４のサブルーチン２とほぼ同一である。但し、ステップＳ２００１では、区間１←ｗ１、区間２←ｗ２の初期設定が追加して行われ、後述する移動平均、移動傾きの計算の区間に用いられる。

　本変形例のサブルーチン２では、ＤＮＳ加味のためにステップＳ２０１０～Ｓ２０１３を追加し、ステップＳ２００９ＡでＤＮＳを加味した一致度の計算を可能とした。

　ステップＳ２０１０では、「ＤＮＳを加味する？」の判断処理が実行され、ＤＮＳを加味すると判断されたときは（ＹＥＳ）、ステップＳ２０１１へ移行し、ＤＮＳを加味しないと判断されたときは（ＮＯ）、ステップＳ２００９Ａへ移行する。ＤＮＳを加味するか否かの起因は、例えば初期の設定による。例えば、ＦＰを、ピーク面積で作成したときはＤＮＳ加味、ＦＰをピーク高さで作成したときは、ＤＮＳ非加味に設定する。

　但し、ＦＰをピーク高さで作成した上記実施例の場合にも、ＤＮＳを加味する処理でＵＶパターン一致度を計算することができ、ＦＰをピーク面積で作成した場合にも、ＤＮＳを加味しない上記実施例の処理でＵＶパターン一致度を計算することができる。

　ステップＳ２０１１では、「区間１（ｗ１）でｘとｙの移動平均を計算」の処理が実行され、区間１（ｗ１）における移動平均が求められる。区間１（ｗ１）は、ＵＶデータの波長に関する区間であり、ステップＳ２００１の初期設定においてｗ１＝３であれば区間１（３）となり、３個の波長におけるＵＶ強度の平均が求められる。具体的には図１３１の図表で後述する。

　ステップＳ２０１２では、「区間２（ｗ２）でｘとｙの移動傾きを計算」の処理が実行され、区間２（ｗ２）における移動傾きが求められる。区間２（ｗ２）は、ステップＳ２０１１で求めた移動平均に関する区間であり、ステップＳ２００１の初期設定においてｗ２＝３であれば区間２（３）となり、ステップＳ２０１１で計算した移動平均に基づき、３個の移動平均に渡る傾きの（±）が求められる。具体的には図１０１の図表で後述する。

　ステップＳ２０１３では、「ｘとｙの移動傾きの符号の不一致数を算出（ＤＮＳ）」の処理が実行され、ステップＳ２０１２で計算された移動傾きから傾きの（±）の一致数が計算される。移動傾きの＋は、図６６において右肩上がりを表し、移動傾きの－は、同右肩下がりを表す。

　ステップＳ２０１３からステップＳ２００９Ａへ移行すると、このステップＳ２００９Ａの処理においてＤＮＳを加味した一致度の計算が行われる。

　ステップＳ２００９Ａでは、「ｘとｙのＵＶスペクトルの一致度を計算（ＵＶ＿Ｓｉｍ）」の処理が実行され、ＤＮＳを加味した一致度の計算では、ＵＶ＿Ｓｉｍを、前記ＵＶスペクトル間距離の二乗和ｚとｘのデータ数ａとＤＮＳとから、
　ＵＶ＿Ｓｉｍ　＝　√（ｚ／ａ）×１．１^ＤＮＳ
として算出し、このＵＶ＿Ｓｉｍを図８１のステップＳ３０６に渡し、ＵＶスペクトルの一致度計算処理を終了する。

　なお、ステップＳ２０１０からステップＳ２００９Ａへ移行した場合の処理は、図８６のステップＳ２００９と同一である。

　図１３１は、移動平均及び移動傾きの計算例を示す図表である。

　図１３１の上段は、ＵＶデータ例、中段は、移動平均の計算例、下段は、移動傾きの計算例を示す。ＵＶデータ例は、具体的な数値に代え、ＵＶ強度をａ１～ａ７で表記している。例えば、２２０ｎｍのＵＶ強度がａ１、２２１ｎｍのＵＶ強度がａ２等となる。移動平均の計算例及び移動傾きの計算例も具体的数値に代え、ＵＶ強度ａ１～ａ７を使用している。

　移動平均は、区間１（ｗ１＝３）を例とし、ステップＳ２０１２（図１３０）において、区間（ａ１、ａ２、ａ３）、区間（ａ２、ａ３、ａ４）、・・・毎に計算した値としてｍ１、ｍ２、・・・が計算される。移動傾きも区間２（３）を例とし、ステップＳ２０１３（図１３０）において、区間（ｍ１、ｍ２、ｍ３）、区間（ｍ２、ｍ３、ｍ４）、・・・毎に計算した値としてｓ１、・・・が計算される。例えば、移動平均の差ｍ３－ｍ１が移動傾きとなり、その（±）を取り出す。

　こうして、ＦＰをピーク面積で作成したときは、基準群ＦＰへの帰属処理ならびに基準ＦＰ帰属結果統合処理において、ＤＮＳを加味した処理でＵＶパターン一致度を計算することができる。この計算により、図６６で示す対応する２点の距離（ｄｉｓ）がピーク高さで作成したＦＰに比較して大きくなっても、取り扱いを容易とし、ＵＶパターン一致度を正確に計算することができる。
［その他］
　本発明実施例のパターン又はＦＰの特徴量作成方法、作成プログラム、及び作成装置では、ＦＰをピーク面積で作成するときは、シグナル強度軸を面積値軸、シグナル強度を面積値として同様に適用することができる。

　本発明実施例は、多成分薬剤として漢方薬の評価について適用したが、その他の多成分物質の評価にも適用することができる。

　本実施例では、対象ＦＰタイプ２又は基準ＦＰタイプ２について領域分割特徴量を作成したが、対象ＦＰ、基準ＦＰについて領域分割特徴量を作成することも可能である。

　また、時系列でピークが変化するパターンを複数の領域に分割し各領域に存在するピークの存在率又は存在量からパターン領域分割特徴量を作成する対象パターン領域分割特徴量作成工程を備えるものであれば、広く適用することができる。

　上記実施例のＦＰは、３Ｄクロマト上での全ピークを対象としたが、細かいデータ、例えば３Ｄクロマト上でピーク面積が５％未満のピークを除いてＦＰを作成することもできる。

　上記実施例のＦＰは、ピーク高さに基づいて作成し、図８７～図９１の評価を得たが、ピーク面積に基づいてＦＰを作成した場合についても、ピーク高さに基づいて作成した上記実施例同様の手順によりＭＴ法によりＭＤ値を求め、図８７～図９１と同様に評価を得ることができる。

　クロマトは、３Ｄクロマトに限らず、ＦＰとしてＵＶスペクトルを除いたピークとそのリテンション・タイムとで構成されたものを用いることもできる。この場合、ＵＶスペクトルの一致度を除き、上記実施例と同様に行わせることができる。

１　多成分薬剤の評価装置（パターンの特徴量作成装置、ＦＰの特徴量作成装置）
３　ＦＰ作成部
５　対象ＦＰピーク帰属部
７　対象ＦＰピーク特徴量作成部
９　対象ＦＰタイプ２作成部
１１　対象ＦＰ領域分割特徴量作成部
１３　対象ＦＰ特徴量統合部
１５　基準ＦＰピーク帰属部
１７　基準ＦＰ帰属結果統合部
１９　基準ＦＰピーク特徴量作成部
２１　基準ＦＰタイプ２作成部
２３　基準ＦＰ領域分割特徴量作成部
２５　基準ＦＰ特徴量統合部
２７　評価部
３１　基準ＦＰ作成部
３３　基準ＦＰ選定部
３５　ピーク・パターン作成部
３７　ピーク帰属部
３９　漢方薬
４１　３Ｄクロマト
４２　対象ＦＰに含まれるピークのＵＶスペクトル
４３　対象ＦＰ
４５　基準群ＦＰ
４７　基準群ＦＰに帰属した対象ＦＰ
４９　対象ＦＰタイプ２（ＦＰ）
５１　対象ＦＰ領域分割特徴量
５３　対象ＦＰの評価結果
５５　薬剤ＡのＦＰ
５７　薬剤ＢのＦＰ
５９　薬剤ＣのＦＰ
６１　対象ＦＰ（リテンション・タイム１０．０－１４．５分）
６３、６５、６７、６９、７１、７３、７５、７７、７９、８１　対象ＦＰ（リテンション・タイム１０．０－１４．５分）中の各ピーク
８３　基準ＦＰ（リテンション・タイム１０．０－１４．５分）
８５、８７、８９、９１、９３、９５、９７、９９、１０１、１０３、１０５　基準ＦＰ（リテンション・タイム１０．０－１４．５分）中の各ピーク
１０７　対象ＦＰリテンション・タイム・出現パターン
１０９　基準ＦＰリテンション・タイム・出現パターン
１１１　リテンション・タイム・出現距離の一致数
１１３　リテンション・タイム・出現パターンの一致度
１１５　対象ＦＰ帰属対象ピークのピーク・パターン（３本）
１１７、１１９、１２１、１２３　基準ＦＰ帰属候補ピークのピーク・パターン（３本）１２５　対象ＦＰ帰属対象ピークのピーク・パターン（５本）
１２７、１２９、１３１、１３３　基準ＦＰ帰属候補ピークのピーク・パターン（５本）１３５　帰属対象ピークのＵＶスペクトル
１３９　帰属候補ピークのＵＶスペクトル
１４１　縦領域分割線
１４３　横領域分割線
１４５　縦・横領域分割線により分割した各領域（格子）
１４７　各領域をピーク高さ
１４８　ＦＰ作成工程
１４９　対象ＦＰピーク帰属工程
１５１　対象ＦＰピーク特徴量作成工程
１５３　対象ＦＰタイプ２作成工程
１５５　対象ＦＰ領域分割特徴量作成工程（パターン領域分割特徴量作成工程、ＦＰ領域分割特徴量作成工程）
１５７　対象ＦＰ特徴量統合工程
１５９　基準ＦＰピーク帰属工程
１６１　基準ＦＰ帰属結果統合工程
１６３　基準ＦＰピーク特徴量作成工程
１６５　基準ＦＰタイプ２作成工程
１６７　基準ＦＰ領域分割特徴量作成工程（パターン領域分割特徴量作成工程、ＦＰ領域分割特徴量作成工程）
１６９　基準ＦＰ特徴量統合工程
１７１　評価工程
１７３　対象ＦＰ作成工程
１７５　基準ＦＰ作成工程
１７７　基準ＦＰ選定工程
１７９　ピーク・パターン作成工程
１８１　ピーク帰属工程
１８３　３Ｄクロマト・データ例
１８５　ピーク情報データ例
１８７　ＦＰデータ例
１８９　判定結果ファイル例
１９１　帰属候補ピークスコア表例
１９３　帰属候補ピーク番号表例
１９５　照合結果ファイル例
１９７　基準群ＦＰデータ例
１９９　対象ＦＰピーク特徴量ファイル例
２０１　ＦＰタイプ２データ例
２０３　対象ＦＰ領域分割特徴量ファイル例
２０５　対象ＦＰ統合特徴量ファイル例
２０７　基準ｔｙｐｅ２群ＦＰデータ例
２０９　基準群統合データ例

Claims

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