JP5915539B2 - パターンの評価方法、多成分物質の評価方法、評価プログラム、及び評価装置 - Google Patents

パターンの評価方法、多成分物質の評価方法、評価プログラム、及び評価装置 Download PDF

Info

Publication number
JP5915539B2
JP5915539B2 JP2012549177A JP2012549177A JP5915539B2 JP 5915539 B2 JP5915539 B2 JP 5915539B2 JP 2012549177 A JP2012549177 A JP 2012549177A JP 2012549177 A JP2012549177 A JP 2012549177A JP 5915539 B2 JP5915539 B2 JP 5915539B2
Authority
JP
Japan
Prior art keywords
peak
target
pattern
peaks
attribution
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2012549177A
Other languages
English (en)
Other versions
JPWO2012164946A1 (ja
Inventor
芳和 森
芳和 森
桂一 野田
桂一 野田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tsumura and Co
Original Assignee
Tsumura and Co
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tsumura and Co filed Critical Tsumura and Co
Publication of JPWO2012164946A1 publication Critical patent/JPWO2012164946A1/ja
Application granted granted Critical
Publication of JP5915539B2 publication Critical patent/JP5915539B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02—Column chromatography
    • G01N30/86—Signal analysis
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/20—Identification of molecular entities, parts thereof or of chemical compositions
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02—Column chromatography
    • G01N30/86—Signal analysis
    • G01N30/8675—Evaluation, i.e. decoding of the signal into analytical information
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02—Column chromatography
    • G01N30/86—Signal analysis
    • G01N30/8675—Evaluation, i.e. decoding of the signal into analytical information
    • G01N30/8686—Fingerprinting, e.g. without prior knowledge of the sample components
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02—Column chromatography
    • G01N30/88—Integrated analysis systems specially adapted therefor, not covered by a single one of the groups G01N30/04 - G01N30/86
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02—Column chromatography
    • G01N30/86—Signal analysis
    • G01N30/8603—Signal analysis with integration or differentiation
    • G01N2030/862—Other mathematical operations for data preprocessing
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02—Column chromatography
    • G01N30/62—Detectors specially adapted therefor
    • G01N30/74—Optical detectors
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02—Column chromatography
    • G01N30/86—Signal analysis
    • G01N30/8624—Detection of slopes or peaks; baseline correction
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02—Column chromatography
    • G01N30/86—Signal analysis
    • G01N30/8624—Detection of slopes or peaks; baseline correction
    • G01N30/8631—Peaks
    • G01N30/8634—Peak quality criteria
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02—Column chromatography
    • G01N30/86—Signal analysis
    • G01N30/8665—Signal analysis for calibrating the measuring apparatus
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02—Column chromatography
    • G01N30/86—Signal analysis
    • G01N30/8693—Models, e.g. prediction of retention times, method development and validation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/08—Feature extraction
    • G06F2218/10—Feature extraction by analysing the shape of a waveform, e.g. extracting parameters relating to peaks
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/12—Classification; Matching
    • G06F2218/14—Classification; Matching by matching peak patterns

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Library & Information Science (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Computing Systems (AREA)
  • Theoretical Computer Science (AREA)
  • Food Science & Technology (AREA)
  • Medicinal Chemistry (AREA)
  • Information Retrieval, Db Structures And Fs Structures Therefor (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
  • Investigating Or Analysing Biological Materials (AREA)

Description

本発明は、パターンの評価方法、多成分物質、例えば多成分薬剤である漢方薬の品質評価を行うための多成分物質の評価方法、評価プログラム、及び評価装置に関する。
多成分物質として、例えば多成分で構成される薬剤（以下、多成分薬剤）である漢方薬などの天然物由来の薬剤がある。これら薬剤の定量的、定性的プロフィールは、使用する原料生薬の、地質学的要因、生態学的要因、収集時期、収集場所、収集年代、生育期の天候等が原因で変化する。
そのため、これら多成分薬剤などについては、その安全性および有効性を担保するための品質として一定の基準を規定し、その基準に基づいて、国家の監督機関、化学的組織、製造業者等が品質評価を行っている。
しかしながら、多成分薬剤の品質等の判定基準は、多成分薬剤中のある特徴的な１ないしは数成分を選択して、その含有量などに基づいて設定することが一般的であった。
例えば、非特許文献１では、多成分薬剤において有効成分の同定ができていない場合、定量分析が可能、水に溶けやすい、熱水中で分解しない、他の成分と化学反応をしない等の物性を持った複数の成分を選択し、化学分析により得られるそれら成分の含有量を評価の基準としている。
また、多成分薬剤にクロマトグラフィーを適用し、保持時間ごとに紫外可視吸収スペクトルを得て、その中の一部の成分情報から評価の基準を設定することも知られている。
例えば、特許文献１では、ＨＰＬＣクロマトグラム・データ（以下、クロマト）中の一部のピークを選択し、バーコード化することによって多成分薬剤を評価している。
しかしながら、これらの方法は、評価の対象が「特定成分の含有量」又は「特定成分のクロマト・ピーク」に限定されており、多成分薬剤が含有する成分の一部を評価対象にしているに過ぎない。このため、多成分薬剤については、評価の対象外となっている成分が多数存在することから、多成分薬剤の評価方法として精度的に不十分である。
多成分薬剤の品質を正確に評価するには、全ピーク情報或いは数％の細かい情報を除外した全ピーク情報に近いピーク情報を網羅したパターンの評価が必要であり、そのため、多成分薬剤間で全ピーク或いはこれに近い個数のピークを対応させる必要がある。
ところが、複数のピークを高精度で効率的に対応させることが困難であり、このことが多成分薬剤の高精度で効率的な評価の妨げとなっていた。
さらに説明すると、同じ製品名の多成分薬剤であっても生薬は天然物であるがために構成成分が僅かに異なる場合がある。このため、同じ品質の薬剤であっても構成成分の含量比率が異なったり、ある薬剤に存在する成分が他の薬剤には存在しないこと（以下、薬剤間誤差）がある。さらに、クロマトにおけるピーク強度やピークの溶出時間には厳密な再現性がない（以下、分析誤差）などの要因もある。これらにより、多成分薬剤間の全ピーク或いはこれに近い個数のピークについて同じ成分に由来するピークに対応させること（以下、ピーク帰属）ができないことが原因となって高精度で効率的な評価の妨げとなっていた。
特開２００２−２１４２１５号公報
月刊薬事 ｖｏｌ.２８, Ｎｏ.３, ６７-７１ (１９８６)
解決しようとする問題点は、既存の評価方法ではピーク情報を網羅したパターンの評価を高精度で効率的に評価することに限界があった点である。
本発明は、評価の精度及び効率を向上させることを可能とするために、時系列でピークが変化する評価対象の対象パターンを取得するパターン取得工程と、前記対象パターンの帰属対象ピークと前記対象パターンに対応し評価基準となる基準パターンにおいて前記帰属対象ピークの出現時間との差が設定した範囲内の帰属候補ピークそれぞれに対し時間軸方向前後位置の少なくとも一方に存在する周辺ｎ本のピークを含めたｎ＋１本のピークによりピーク・パターンを作成するピーク・パターン作成工程と、前記帰属対象ピーク及び帰属候補ピークのピーク・パターンを比較し対応する前記帰属対象ピークと帰属候補ピークを特定するピーク帰属工程と、このピーク帰属工程で特定されて帰属したピークと評価基準となる複数の基準パターンのピークとを比較評価する評価工程とを備えたことをパターンの評価方法の特徴とする。
本発明は、評価対象となる多成分物質のクロマトから検出されたピークとそのリテンション・タイムとで構成される対象ＦＰを作成するＦＰ作成工程と、前記対象ＦＰの帰属対象ピークと前記対象ＦＰに対応し評価基準となる多成分物質のクロマトから検出されたピークとそのリテンション・タイムとで構成される基準ＦＰにおいて前記帰属対象ピークのリテンション・タイムとの差が設定した範囲内の帰属候補ピークそれぞれに対し時間軸方向前後位置の少なくとも一方に存在する周辺ｎ本のピークを含めたｎ＋１本のピークによりピーク・パターンを作成するピーク・パターン作成工程と、前記帰属対象ピーク及び帰属候補ピークのピーク・パターンを比較し対応する前記帰属対象ピークと帰属候補ピークを特定するピーク帰属工程と、このピーク帰属工程で特定されて帰属したピークと評価基準となる複数の基準ＦＰのピークとを比較評価する評価工程とを備えたことを多成分物質の評価方法の特徴とする。
本発明は、時系列でピークが変化する評価対象の対象パターンを取得するパターン取得機能と、前記対象パターンの帰属対象ピークと前記対象パターンに対応し評価基準となる基準パターンにおいて前記帰属対象ピークの出現時間との差が設定した範囲内の帰属候補ピークそれぞれに対し時間軸方向前後位置の少なくとも一方に存在する周辺ｎ本のピークを含めたｎ＋１本のピークによりピーク・パターンを作成するピーク・パターン作成機能と、前記帰属対象ピーク及び帰属候補ピークのピーク・パターンを比較し対応する前記帰属対象ピークと帰属候補ピークを特定するピーク帰属機能と、このピーク帰属機能で特定されて帰属したピークと評価基準となる複数の基準パターンのピークとを比較評価する評価機能とをコンピュータに実現させることをパターンの評価プログラムの特徴とする。
本発明は 評価対象となる多成分物質のクロマトから検出されたピークとそのリテンション・タイムとで構成される対象ＦＰを作成するＦＰ作成機能と、前記対象ＦＰの帰属対象ピークと前記対象ＦＰに対応し評価基準となる多成分物質のクロマトから検出されたピークとそのリテンション・タイムとで構成される基準ＦＰにおいて前記帰属対象ピークのリテンション・タイムとの差が設定した範囲内の帰属候補ピークそれぞれに対し時間軸方向前後位置の少なくとも一方に存在する周辺ｎ本のピークを含めたｎ＋１本のピークによりピーク・パターンを作成するピーク・パターン作成機能と、前記帰属対象ピーク及び帰属候補ピークのピーク・パターンを比較し対応する前記帰属対象ピークと帰属候補ピークを特定するピーク帰属機能と、このピーク帰属機能で特定されて帰属したピークと評価基準となる複数の基準ＦＰのピークとを比較評価する評価機能とを備えたことを多成分物質の評価プログラムの特徴とする。
本発明は、時系列でピークが変化する評価対象の対象パターンを取得するパターン取得部と、前記対象パターンの帰属対象ピークと前記対象パターンに対応し評価基準となる基準パターンにおいて前記帰属対象ピークの出現時間との差が設定した範囲内の帰属候補ピークそれぞれに対し時間軸方向前後位置の少なくとも一方に存在する周辺ｎ本のピークを含めたｎ＋１本のピークによりピーク・パターンを作成するピーク・パターン作成部と、前記帰属対象ピーク及び帰属候補ピークのピーク・パターンを比較し対応する前記帰属対象ピークと帰属候補ピークを特定するピーク帰属部と、このピーク帰属部で特定されて帰属したピークと評価基準となる複数の基準パターンのピークとを比較評価する評価部とを備えたことをパターンの評価装置の特徴とする。
評価対象となる多成分物質のクロマトから検出されたピークとそのリテンション・タイムとで構成される対象ＦＰを作成するＦＰ作成部と、前記対象ＦＰの帰属対象ピークと前記対象ＦＰに対応し評価基準となる多成分物質のクロマトから検出されたピークとそのリテンション・タイムとで構成される基準ＦＰにおいて前記帰属対象ピークのリテンション・タイムとの差が設定した範囲内の帰属候補ピークそれぞれに対し時間軸方向前後位置の少なくとも一方に存在する周辺ｎ本のピークを含めたｎ＋１本のピークによりピーク・パターンを作成するピーク・パターン作成部と、前記帰属対象ピーク及び帰属候補ピークのピーク・パターンを比較し対応する前記帰属対象ピークと帰属候補ピークを特定するピーク帰属部と、このピーク帰属部で特定されて帰属したピークと評価基準となる複数の基準ＦＰのピークとを比較評価する評価部とを備えたことを多成分物質の評価装置の特徴とする。
本発明のパターンの評価方法は、上記構成であるから、ピーク帰属をピーク・パターンのパターン比較で行わせることができ、対象パターンの各ピークを基準パターンの各ピークに対して精度よく、効率よく帰属させることが可能となり、評価の精度及び効率をより向上させることができる。
本発明の多成分物質の評価方法も同様にして、多成分物質の評価の精度及び効率をより向上させることができる。
本発明のパターン及び多成分物質の評価プログラムは、上記構成であるから、各機能をコンピュータに実現させ、評価の精度及び効率をより向上させることができる。
本発明のパターン及び多成分物質の評価装置は、上記構成であるから、各部を作用させ、評価の精度及び効率をより向上させることができる。
多成分薬剤の評価装置のブロック図である。（実施例１） 多成分薬剤の評価手順を示すブロック図である。（実施例１） 三次元クロマトグラム・データ（以下、３Ｄクロマト）から作成したＦＰの説明図である。（実施例１） 薬剤ごとのＦＰを示し、（Ａ）は薬剤Ａ、（Ｂ）は薬剤Ｂ、（Ｃ）は薬剤Ｃを示すグラフである。（実施例１） 対象ＦＰ及び基準ＦＰのリテンション・タイムを示す図である。（実施例１） 対象ＦＰのリテンション・タイム・出現パターンを示す図である。（実施例１） 基準ＦＰのリテンション・タイム・出現パターンを示す図である。（実施例１） 対象ＦＰと基準ＦＰのリテンション・タイム・出現距離の一致数を示す図である。（実施例１） 対象ＦＰと基準ＦＰのリテンション・タイム・出現パターンの一致度を示す図である。（実施例１） 対象ＦＰの帰属対象ピークを示す図である。（実施例１） 帰属対象ピークを含めたピーク３本によるピーク・パターン図である。（実施例１） 帰属対象ピークを含めたピーク５本によるピーク・パターン図である。（実施例１） 帰属対象ピークの許容幅を示す図である。（実施例１） 帰属対象ピークに対する基準ＦＰの帰属候補ピークを示す図である。（実施例１） 帰属対象ピークと帰属候補ピークのピーク３本によるピーク・パターン図である。（実施例１） 帰属対象ピークと別な帰属候補ピークのピーク３本によるピーク・パターン図である。（実施例１） 帰属対象ピークと別な帰属候補ピークのピーク３本によるピーク・パターン図である。（実施例１）） 帰属対象ピークと別な帰属候補ピークのピーク３本によるピーク・パターン図である。（実施例１） 帰属対象ピークと帰属候補ピークのピーク５本によるピーク・パターン図である。（実施例１） 帰属対象ピークと別な帰属候補ピークのピーク５本によるピーク・パターン図である。（実施例１） 帰属対象ピークと別な帰属候補ピークのピーク５本によるピーク・パターン図である。（実施例１）） 帰属対象ピークと別な帰属候補ピークのピーク５本によるピーク・パターン図である。（実施例１） 帰属対象ピークと帰属候補ピークのピーク・パターン構成候補ピークを示す図である。（実施例１） ピーク・パターン構成候補ピークを4本とした場合の帰属対象ピークの全ピーク・パターン数を示す図である。（実施例１） ピーク・パターン構成候補ピークを4本とした場合の帰属候補ピークの全ピーク・パターン数を示す図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークのピーク・パターンに対する帰属候補ピークのピーク・パターンの網羅的比較の説明図である。（実施例１） 帰属対象ピークと帰属候補ピークのピーク３本によるピーク・パターンの一致度の算出方法を示す図である。（実施例１） 帰属対象ピークと帰属候補ピークのピーク３本によるピーク・パターンの一致度の算出方法を示す図である。（実施例１） 帰属対象ピークと帰属候補ピークのピーク５本によるピーク・パターンの一致度の算出方法を示す図である。（実施例１） 帰属対象ピーク及び帰属候補ピークのＵＶスペクトルを示す図である。（実施例１） 帰属対象ピークと帰属候補ピークのＵＶスペクトルの一致度の説明図である。（実施例１） ピーク・パターンとＵＶスペクトルの両方の比較による帰属候補ピークの一致度の説明図である。（実施例１） 対象ＦＰの基準群ＦＰへの帰属を示す説明図である。（実施例１） 対象ＦＰが基準群ＦＰへ帰属された状況を示す図である。（実施例１） 各種対象ＦＰとその評価値（ＭＤ値）を示す図である。（実施例１） 各種対象ＦＰとその評価値（ＭＤ値）を示す図である。（実施例１） 各種対象ＦＰとその評価値（ＭＤ値）を示す図である。（実施例１） 各種対象ＦＰとその評価値（ＭＤ値）を示す図である。（実施例１） 各種対象ＦＰとその評価値（ＭＤ値）を示す図である。（実施例１） 多成分薬剤の評価方法を示す工程図である。（実施例１） 多成分薬剤の評価フロー図である。（実施例１） 単波長によるＦＰ作成機能におけるデータ処理フローチャートである。（実施例１） 複数波長によるＦＰ作成機能におけるデータ処理フローチャートである。（実施例１） 複数波長によるＦＰ作成機能におけるデータ処理フローチャートである。（実施例１） ピーク帰属処理１（基準ＦＰの選定）におけるデータ処理フローチャートである。（実施例１） ピーク帰属処理２（帰属スコアの算出）におけるデータ処理フローチャートである。（実施例１） ピーク帰属処理３（対応ピークの特定）におけるデータ処理フローチャートである。（実施例１） ピーク帰属処理４（基準群ＦＰへの帰属）におけるデータ処理フローチャートである。（実施例１） ピーク帰属処理４（基準群ＦＰへの帰属）におけるデータ処理フローチャートである。（実施例１） ピーク帰属処理１（基準ＦＰの選定）におけるリテンション・タイム・出現パターンの一致度計算処理のフローチャートである。（実施例１） ピーク帰属処理２（帰属スコアの算出）におけるＵＶスペクトルの一致度計算処理のフローチャートである。（実施例１） ピーク帰属処理２（帰属スコアの算出）におけるピーク・パターンの一致度計算処理のフローチャートである。（実施例１） 基準ＦＰ特徴量ファイルを作成するためのフローチャートである。（実施例１） 「基準ＦＰ帰属結果統合処理（基準ＦＰ対応表の作成）」の詳細を示すフローチャートである。（実施例１） 「基準ＦＰ帰属結果統合処理（基準ＦＰ対応表の作成）」の詳細を示すフローチャートである。（実施例１） 「ピーク特徴量化処理（基準群ＦＰの作成）」の詳細を示すフローチャートである。（実施例１） ３Ｄクロマトのデータ例を示す図表である。（実施例１） ピーク情報のデータ例を示す図表である。（実施例１） ＦＰのデータ例を示す図表である。（実施例１） 対象ＦＰの基準ＦＰへの帰属スコア計算結果（判定結果）ファイル例を示す図表である。（実施例１） 対象ＦＰと基準ＦＰで対応するピークの照合過程を示す図表である。（実施例１） 照合結果ファイル例を示す図表である。（実施例１） 基準群ＦＰのデータ例を示す図表である。（実施例１） 対象ＦＰピーク特徴量ファイル例を示す図表である。（実施例１） 図８６に代えて適用するサブルーチン２の変形例の詳細を示すフローチャートである。（実施例１） 移動平均及び移動傾きの計算例を示す図表である。（実施例１）
評価の精度及び効率を向上させることを可能にするという目的は、対象パターンの各ピークの基準群パターンの各ピークへの高精度かつ効率的なピーク一斉帰属により得られた結果を特徴量として比較評価することにより実現した。
高精度で効率的なピーク一斉帰属は、対象パターンと基準パターンのピークを時間軸方向前後位置の少なくとも一方に存在するピークでパターン化したピーク・パターンを比較することにより実現した。
本発明の実施例１は、パターンの評価方法、評価プログラム、評価装置を、多成分物質、例えば多成分薬剤を評価する多成分薬剤の評価方法、評価プログラム、評価装置として適用したものである。
多成分薬剤は、複数の有効化学成分を含有する薬剤と定義され、限定はされないが、生薬、生薬の組合せ、それらの抽出物、漢方薬等が含まれる。また剤形も特に限定されず、例えば、第１５改正日本薬局方の製剤総則で規定されている液剤、エキス剤、カプセル剤、顆粒剤、丸剤、懸濁剤・乳剤、散剤、酒精剤、錠剤、浸剤・煎剤、チンキ剤、トローチ剤、芳香水剤、流エキス剤等が含まれる。多成分物質としては、薬剤以外のものも含まれる。
漢方薬の具体例は、医療用漢方製剤１４８処方「使用上の注意」の業界統一と自主改訂、一般用漢方処方の手引き(１９７８年)に記載されている。
多成分薬剤の評価では、評価対象薬剤が正常品と定めた複数の薬剤と同等であるかどうかを評価するため、まず、評価対象薬剤の三次元クロマトグラム・データ（以下、３Ｄクロマト）からその薬剤特有の情報を抽出した対象ＦＰを作成する。
次に、対象ＦＰの各ピークを、全基準ＦＰをピーク帰属処理し作成した全基準ＦＰのピーク対応データ（以下、基準群ＦＰ）に帰属し、ピーク特徴量を得る。
次に、基準群ＦＰのピークと帰属した対象ＦＰのピーク（以下、対象ＦＰ帰属ピーク）の同等性をＭＴ法で評価する。最後に、得られた評価値（以下、ＭＤ値）と予め設定しておいた判定値（ＭＤ値の上限値）を比較し、評価対象薬剤が正常品と同等であるかどうかを判定する。
なお、３Ｄクロマトは、評価対象となる多成分物質である多成分薬剤のＨＰＬＣクロマトグラム・データ（以下、クロマト）であり、ＵＶスペクトルを含んでいる。
ＦＰとは、特定の波長で検出されたピークにおけるシグナル強度（高さ）の極大値又は面積値（以下、ピーク）とそのピークの出現時間（以下、リテンション・タイム）、とで構成するフィンガー・プリント・データである。
対象ＦＰは、評価対象である漢方薬の三次元クロマトグラム・データである３Ｄクロマトから、特定の検出波長における複数のピークとそのリテンション・タイムならびにＵＶスペクトルを抽出したものである。
基準ＦＰは、対象ＦＰに対応しており、正常品と定めた多成分物質である多成分薬剤としての漢方薬のＦＰである。
［多成分薬剤の評価装置］
図１は、多成分薬剤の評価装置のブロック図、図２は、多成分薬剤の評価手順を示すブロック図、図３は、３Ｄクロマトから作成したＦＰの説明図、図４（Ａ）は、薬剤Ａ、（Ｂ）は、薬剤Ｂ、（Ｃ）は、薬剤ＣのＦＰである。
図1、図２のように、パターンの評価装置としての多成分薬剤の評価装置１は、パターン取得部としてのＦＰ作成部３と、基準ＦＰ選定部５と、ピーク・パターン作成部７と、ピーク帰属部９と、評価部１１とを備えている。この多成分薬剤の評価装置１は、コンピュータで構成され、図示はしないが、ＣＰＵ、ＲＯＭ、ＲＡＭなどを備えている。多成分薬剤の評価装置１は、コンピュータにインストールされた多成分物質の評価プログラムとして多成分薬剤の評価プログラムを実現させ、多成分薬剤の評価を行わせることができる。但し、多成分薬剤の評価プログラムは、これを記録した多成分薬剤の評価プログラム記録媒体を用い、コンピュータで構成された多成分薬剤の評価装置１にこれを読み取らせることで、多成分薬剤の評価を実現させることもできる。
本実施例において、ＦＰ作成部３と、基準ＦＰ選定部５と、ピーク・パターン作成部７と、ピーク帰属部９と、評価部１１とは、一つのコンピュータで構成している。但し、ＦＰ作成部３と、基準ＦＰ選定部５と、ピーク・パターン作成部７と、ピーク帰属部９と、評価部１１とを、それぞれ別々のコンピュータで構成し、或いはＦＰ作成部３及び基準ＦＰ選定部５とピーク・パターン作成部７及びピーク帰属部９と評価部１１とを別々のコンピュータで構成することなどもできる。
ＦＰ作成部３は、時系列でピークが変化する評価対象の対象パターンを取得するものである。具体的にＦＰ作成部３は、例えば、対象ＦＰ１７（以下単に、「ＦＰ１７」ということもある。）を対象パターンとして作成し取得する機能部である。このＦＰ１７は、３Ｄクロマト１５と同様に三次元の情報（ピーク、リテンション・タイム及びＵＶスペクトル）で構成している。
このため、ＦＰ１７は、その薬剤特有の情報をそのまま継承したデータである。それにも係らず、データ容量は約１／７０に圧縮されているため、３Dクロマト１５に比較して、処理すべき情報量を大幅に減少させることができ処理速度を速めることができる。
３Ｄクロマト１５は、漢方薬１３（図２参照）に高速液体クロマトグラフィー（ＨＰＬＣ）を適用した結果である。この３Ｄクロマト１５は、各成分の移動速度として現れ、それを特定時間における移動距離として表し、或いはカラム末端から時系列に現れる様をチャートに表したものである。ＨＰＬＣにおいては、時間軸に対する検出器応答をプロットしたもので、ピークの出現時間を保持時間（リテンション・タイム）と呼んでいる。
検出器としては、特に限定されないが、光学的性質を利用した吸光度検出器(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（＝６）パターンのピーク・パターンが作成される。したがって、帰属対象ピーク４５及び帰属候補ピーク６５のＰ＿Ｓｉｍは、６パターン×６パターン（＝３６）通りで算出され、これらＰ＿Ｓｉｍの最小値であるＰ＿Ｓｉｍ＿ｍｉｎを帰属対象ピーク４５と帰属候補ピーク６５の一致度とする。
ちなみに、帰属対象ピーク４５及び帰属候補ピーク６５のそれぞれでピーク・パターン構成候補ピークを時間軸方向前後周辺１０本とし、ピーク・パターン構成ピークを任意の２本に設定すると、帰属対象ピーク及び帰属候補ピークそれぞれで、4Ｃ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)
として算出する。
なお、式中のｄは、リテンション・タイムのずれを補正する値である。
ピーク帰属部９では、図６５、図６６のように帰属対象ピークと帰属候補ピークでＵＶスペクトルの一致度を算出する。
図６５は、帰属対象ピーク４５及び帰属候補ピーク６７のＵＶスペクトル（１０７と１１１）図であり、図６６のようにこれら２つのＵＶスペクトルの一致度（ＵＶ＿Ｓｉｍ（４５−６７））は、
ＵＶ＿Ｓｉｍ（４５−６７）＝ ＲＭＳＤ（１０７ ｖｓ １１１）
として算出する。
ＲＭＳＤは、平均二乗偏差のことで、対応する２点の距離（ｄｉｓ）をそれぞれ二乗し、その相加平均の平方根として定義される。つまり、
ＲＭＳＤ ＝ √｛ Σｄｉｓ2 ／ ｎ ｝
として算出する。
ｎは、ｄｉｓの数である。
ここで、ＵＶスペクトルの波形は極大波長及び極小波長を含んでおり、この極大波長及び極小波長或いは何れかを比較することで一致度を算出することも可能である。しかし、吸収特性のない化合物あるいは吸収特性が類似した化合物等では、極大波長、極小波長は同じであるが、全体の波形がかなり異なる場合もあり、極大波長、極小波長の比較では、波形の一致度を算出できない恐れがある。
これに対し、ＵＶスペクトルの波形によりＲＭＳＤを利用した場合には、波形全体の比較となるため、ＵＶスペクトルの波形の一致度の算出をより正確に行わせ、吸収特性のない化合物あるいは吸収特性が類似した化合物でも正確に識別できることになる。
このＵＶスペクトルの一致度は、帰属対象ピーク４５の全ての帰属候補ピークについて同様に算出する。
さらに、ピーク帰属部９では、図６７のように上記２つの一致度を統合した帰属候補ピークの一致度を算出する。
図６７のように帰属候補ピークの一致度（ＳＣＯＲＥ（４５−６７））はピーク・パターンとＵＶスペクトルのそれぞれの一致度を乗じて算出する。ピーク・パターン４５、６７の一致度を示すスコアが、Ｐ＿Ｓｉｍ＿ｍｉｎ(４５−６７)であり、対応するＵＶスペクトル１０７、１１１の一致度を示すスコアが、ＵＶ＿Ｓｉｍ（４５−６７）であるとする。このとき、帰属候補ピークの一致度 ＳＣＯＲＥ（４５−６７）は、
ＳＣＯＲＥ（４５−６７）
＝ Ｐ＿Ｓｉｍ＿ｍｉｎ(４５−６７)×ＵＶ＿Ｓｉｍ（４５−６７）
として算出する。
この帰属候補ピークの一致度を帰属対象ピーク４５の全ての帰属候補ピークについて同様に計算する。
そして、全帰属候補ピーク間でこのＳＣＯＲＥを比較し、ＳＣＯＲＥが最小となる帰属候補ピークを帰属対象ピーク４５の帰属ピークとして判定する。
ピーク帰属部９では、帰属対象ピークの帰属すべきピークを２つの観点を総合して判定するため、正確なピーク帰属を実現することができる。
また、ピーク帰属部９では、対象ＦＰの基準ＦＰへの帰属結果をもとに、図６８のように、対象ＦＰ１７の各ピークを基準群ＦＰ１９へ帰属する。
対象ＦＰ１７の各ピークは、前記帰属処理により基準群ＦＰを構成する基準ＦＰに帰属される。この帰属結果をもとに最終的に基準群ＦＰ１９のピークに帰属する。
尚、基準群ＦＰ１９は、正常品と定めた複数の基準ＦＰ全てを前記のように帰属処理し作成したものであり、その各ピークは帰属されたピークの平均値（黒点）±標準偏差（縦線）で表している。
図６９は、対象ＦＰ１７を基準群ＦＰ１９に帰属した結果であり、この結果が対象ＦＰ１７の帰属処理の最終結果である。
［ＭＤ値］
この結果から、前記のようにＭＴ法によりＭＤ値（ＭＤ値：０．２５、２．９９等）を求めることができる（図７０〜図７４参照）。
［多成分薬剤の評価方法］
図７５は、本発明実施例１の多成分薬剤の評価方法を示す工程図である。
図７５のように、パターンの評価方法としての多成分薬剤の評価方法は、パターン取得工程としてのＦＰ作成工程１１３と、基準パターン選定工程としての基準ＦＰ選定工程１１５と、ピーク・パターン作成工程１１７と、ピーク帰属工程１１９と、評価工程１２１とを備えている。これら、ＦＰ作成工程１１３、基準ＦＰ選定工程１１５、ピーク・パターン作成工程１１７、ピーク帰属工程１１９、評価工程１２１は、本実施例において前記多成分薬剤の評価装置１を用いて行われ、ＦＰ作成工程１１３は、ＦＰ作成部３の機能により行わせ、同様に基準ＦＰ選定工程１１５、ピーク・パターン作成工程１１７、ピーク帰属工程１１９、評価工程１２１を、基準ＦＰ選定部５、ピーク・パターン特定部７、ピーク帰属部９、評価部１１の機能により行わせることができる。
［多成分薬剤の評価プログラム］
図７６〜図９１は、多成分薬剤の評価プログラムに係るフローチャート、図９２は、３Ｄクロマトのデータ例を示す図表、図９３は、ピーク情報データ例を示す図表、図９４は、ＦＰデータ例を示す図表、図９５は、ステップＳ３で作成した判定結果ファイル例を示す図表、図９６は、対象ＦＰと基準ＦＰで対応するピークを特定する過程で作成する２つの中間ファイル（帰属候補ピークスコア表、帰属候補ピーク番号表）例を示す図表、図９７は、対象ＦＰと基準ＦＰで対応するピークを特定した結果である照合結果ファイル例を示す図表、図９８は、基準群ＦＰデータ例を示す図表、図９９は、対象ＦＰ帰属ピークのデータである対象ＦＰのピーク・データ特徴量ファイル例を示す図表である。
図７６は、評価対象薬剤を評価するための処理全体のステップを示すフローチャートであり、システム起動によりスタートし、ＦＰ作成部３のＦＰ作成機能と、基準ＦＰ選定部５の基準ＦＰ選定機能と、ピーク・パターン作成部７のピーク・パターン作成機能と、ピーク帰属部９のピーク帰属機能と、評価部１１の評価機能とをコンピュータに実現させる。
ＦＰ作成機能は、ステップＳ１で実現され、基準ＦＰ選定機能は、ステップＳ２で実現され、ピーク・パターン作成機能は、ステップＳ３で実現され、ピーク帰属機能は、ステップＳ３〜Ｓ５で実現され、評価機能は、ステップＳ６、Ｓ７で実現される。
ステップＳ１は、３Ｄクロマト及び特定の検出波長におけるピーク情報を入力データとして「ＦＰ作成処理」が実行される。
３Ｄクロマトは、評価対象薬剤をＨＰＬＣで分析することにより得られるデータであり、図９２の３Ｄクロマトのデータ例１２３で示すように、リテンション・タイム、検出波長、ピーク（シグナル強度）の三次元情報で構成されたデータである。ピーク情報は、同ＨＰＬＣ分析により得られる特定波長におけるクロマト・データをＨＰＬＣデータ解析ツール（例えば、ＣｈｅｍＳｔａｔｉｏｎ等）で処理することで得られるデータであり、図９３のピーク情報データ例１２５のように、ピークとして検出された全ピークの極大値及び面積値とその時点のリテンション・タイム等で構成されたデータである。
ステップＳ１では、コンピュータのＦＰ作成部３（図１）が機能し、３Ｄクロマト及びピーク情報から前記対象ＦＰ１７（図２）を作成し、そのデータをファイルとして出力する。この対象ＦＰ１７は、図９４のＦＰデータ例１２７のように、リテンション・タイムとピーク高さとピーク高さごとのＵＶスペクトルで構成されたデータである。
ステップＳ２は、ステップＳ１で出力した対象ＦＰ及び全基準ＦＰを入力として、「対象ＦＰ帰属処理１」が実行される。
ステップＳ２では、コンピュータの基準ＦＰ選定部５が機能し、全基準ＦＰに対して対象ＦＰ１７とのリテンション・タイム・出現パターンの一致度を算出し、対象ＦＰ１７の帰属に適した基準ＦＰを選定する。
基準ＦＰは、正常品と定めた薬剤の３Ｄクロマトとピーク情報から前記ステップＳ１と同様の処理により作成されたＦＰである。なお、正常品は、安全性、有効性が確認された薬剤と定義され、製品ロットの異なる複数の薬剤が該当する。 基準ＦＰも、図９４のＦＰデータ例１２７と同様に構成されたデータである。
ステップＳ３は、対象ＦＰ１７とステップＳ２で選定した基準ＦＰを入力とし、「対象ＦＰ帰属処理２」が実行される。
ステップＳ３では、コンピュータのピーク・パターン作成部７（図１）及びピーク帰属部９（図１）が機能する。この機能により、対象ＦＰ１７とステップＳ２で選定した基準ＦＰの全ピークで、図２３〜図６１のように網羅的にピーク・パターンを作成し、次にそれらピーク・パターンの一致度（図６３または図６４のＰ＿Ｓｉｍ）を算出する。また、対象ＦＰと基準ＦＰのピーク間でＵＶスペクトルの一致度（図６６のＵＶ＿Ｓｉｍ）を算出する。さらにこれら２つの一致度から帰属候補ピークの一致度（図６７のＳＣＯＲＥ）を算出し、その算出結果をファイル（判定結果ファイル）に出力する。
ステップＳ４は、ステップＳ３で出力した判定結果ファイルを入力とし、「対象ＦＰ帰属処理３」が実行される。
ステップＳ４では、コンピュータのピーク帰属部７が機能し、対象ＦＰ１７と基準ＦＰ間で、帰属候補ピークの一致度（ＳＣＯＲＥ）をもとに対象ＦＰの各ピークに対応する基準ＦＰのピークを特定し、その結果をファイル（照合結果ファイル）に出力する。
ステップＳ５は、ステップＳ４で出力した照合結果ファイルと基準群ＦＰを入力とし、「対象ＦＰ帰属処理４」が実行される。
基準群ＦＰは、全基準ＦＰから前記ステップＳ２からステップＳ４と同様の処理により作成された全基準ＦＰ間のピーク対応データである。
ステップＳ５では、コンピュータのピーク帰属部７が機能し、対象ＦＰ１７の照合結果ファイルをもとに、図６８、図６９のように、対象ＦＰ１７の各ピークを基準群ＦＰのピークに帰属し、その結果をファイル（ピーク・データ特徴量ファイル）に出力する。
ステップＳ６は、ステップＳ５で出力したピーク・データ特徴量ファイルと基準群ＦＰを入力とし、「ＦＰ評価処理」が実行される。
ステップＳ６では、コンピュータの評価部１１が機能し、ステップＳ５で出力したピーク・データ特徴量データと基準群ＦＰとの同等性をＭＴ法により評価し、その評価結果をＭＤ値をとして出力する（図７０〜図７４）。
ステップＳ７は、ステップＳ６で出力したＭＤ値を入力として、「合否判定」が実行される。
ステップＳ７では、コンピュータの評価部１１が機能し、ステップＳ６で出力したＭＤ値と予め設定した閾値（ＭＤ値の上限値）を比較し、合否を判定する（図２のグラフ２３）。
［Ｓ１：ＦＰ作成処理（単一波長のみ利用）］
図７７は、図７６ステップＳ１「ＦＰ作成処理」の単一波長のピーク情報を利用した場合のフローチャートである。
図７７は、波長を単一波長、例えば２０３ｎｍとして評価対象のＦＰを作成するステップの詳細である。この処理では、３Ｄクロマトと検出波長が２０３ｎｍにおけるピーク情報から、２０３ｎｍで検出されたピークにおけるリテンション・タイムとピークならびにそれらピークのＵＶスペクトルで構成するＦＰを作成する。
ステップＳ１０１では、「ピーク情報を読み込む」の処理が実行される。この処理では、ＦＰの作成に必要な２つのデータのうちの1つ目としてピーク情報が読み込まれ、ステップＳ１０２へ移行する。
ステップＳ１０２では、「ピークのリテンション・タイム（Ｒ１）と対応するピーク・データ（Ｐ１）を順番に取得」の処理が実行される。この処理では、ピーク情報から、ピークのリテンション・タイム（Ｒ１）及びピーク・データ（Ｐ１）を１ピークずつ順番に取得し、ステップＳ１０３へ移行する。
ステップＳ１０３では、「３Ｄクロマトを読み込む」の処理が実行される。この処理では、ＦＰの作成に必要な２つのデータのうちの２つ目として３Ｄクロマトが読み込まれ、ステップＳ１０４へ移行する。
ステップＳ１０４では、「ピークのリテンション・タイム（Ｒ２）と対応するＵＶスペクトル（Ｕ１）を順番に取得」の処理が実行される。この処理では、３Ｄクロマトから、リテンション・タイム（Ｒ２）及びＵＶスペクトル（Ｕ１）をＨＰＬＣ分析時のサンプリング・レート／２ごとに取得し、ステップＳ１０５へ移行する。
ステップＳ１０５では、「｜Ｒ１−Ｒ２｜＜＝閾値？」の判断処理が実行される。この処理では、ステップＳ１０２及びＳ１０４で読み込まれたＲ１とＲ２が閾値の範囲で対応したものであるか否かが判断される。対応している（ＹＥＳ）場合は、２つのリテンション・タイムは同じであり、リテンション・タイムがＲ１のピークのＵＶスペクトルはＵ１
であると判断し、ステップＳ１０６へ移行する。対応していない（ＮＯ）場合は、２つのリテンション・タイムは同じではなく、リテンション・タイムがＲ１のピークのＵＶスペクトルはＵ１ではないと判断し、３Ｄクロマトの次のデータとの比較のため、ステップＳ１０４へ移行する。なお、この判断処理での閾値は、３Ｄクロマトにおける“サンプリング・レート”とする。
ステップＳ１０６では、「Ｕ１を最大値１で規格化」の処理が実行される。この処理では、Ｓ１０５でＲ１のＵＶスペクトルと判断したＵ１を最大値１で規格化し、ステップＳ１０７へ移行する。
ステップＳ１０７では、「Ｒ１とＰ１ならびに規格化したＵ１を出力（対象ＦＰ）」の処理が実行される。この処理では、ピーク情報から取得したＲ１とＰ１ならびにＳ１０６で規格化したＵ１を対象ＦＰに出力し、ステップＳ１０８へ移行する。
ステップＳ１０８では、「全ピークの処理終了？」の判断処理が実行される。この処理では、ピーク情報中の全てのピークに対して処理が行われたか否かが判断され、全ピークに対して処理が終了していなければ（ＮＯ）、未処理のピークを処理するため、ステップＳ１０２へ移行する。Ｓ１０２からＳ１０８までの処理は全ピークの処理が終了するまで繰り返され、全ピークの処理が終了すると（ＹＥＳ）、ＦＰ作成処理を終了する。
［Ｓ１：ＦＰ作成処理（複数波長利用）］
図７８、図７９は、図７６ステップＳ１「ＦＰ作成処理」において、前記単一波長のピーク情報に代え、複数波長のピーク情報を利用した場合のフローチャートである。例えば２０３ｎｍを含めて、検出波長軸方向に複数（ｎ個）の波長を選択し、ＦＰを作成する場合である。
このＦＰ作成処理は、図７７のような単一波長では３Ｄクロマトで検出されている全ピークを網羅できない場合に複数波長のピーク情報を利用し、３Ｄクロマトの全ピークを網羅したＦＰを作成するためのものである。
なお、図７８、図７９は、上記単一波長のみを利用したＦＰ作成処理で波長ごとのＦＰをｎ個作成後、それらＦＰから複数波長によるＦＰを作成するステップの詳細である。
ステップＳ１１０では、「波長ごとにＦＰを作成」の処理が実行される。この処理では、波長ごとに上記単一波長のみを利用したＦＰ作成処理が行われ、ｎ個のＦＰを作成し、ステップＳ１１１へ移行する。
ステップＳ１１１では、「ＦＰをピーク数（降順）でリスト化」の処理が実行される。この処理では、ｎ個のＦＰをピーク数が多い順にリスト化し、ステップＳ１１２へ移行する。
ステップＳ１１２では、ｎ個のＦＰを順番に処理するためのカウンタの初期化としてｎに１が代入（ｎ←１）され、ステップＳ１１３へ移行する。
ステップＳ１１３では、「リストｎ番目のＦＰを読み込む」の処理が実行される。この処理では、リストｎ番目のＦＰを読み込み、ステップＳ１１４へ移行する。
ステップＳ１１４では、「全リテンション・タイム（Ｘ）を取得」の処理が実行される。この処理では、Ｓ１１３で読み込んだＦＰのリテンション・タイム情報を全て取得し、ステップＳ１１５へ移行する。
ステップＳ１１５では、「ｎの更新（ｎ←ｎ＋１）」の処理が実行される。この処理では、処理を次のＦＰへ移行するため、ｎの更新としてｎにｎ＋１を代入し、ステップＳ１１６へ移行する。
ステップＳ１１６では、「リストｎ番目のＦＰを読み込む」の処理が実行される。この処理では、リストｎ番目のＦＰを読み込み、ステップＳ１１７へ移行する。
ステップＳ１１７では、「全リテンション・タイム（Ｙ）を取得」の処理が実行される。この処理では、Ｓ１１６で読み込んだＦＰのリテンション・タイム情報を全て取得し、ステップＳ１１８へ移行する。
ステップＳ１１８では、「ＸとＹを重複なしで統合（Ｚ）」の処理が実行される。この処理では、Ｓ１１４で取得したリテンション・タイム情報ＸとＳ１１７で取得したリテンション・タイム情報Ｙを重複なしで統合した後、Ｚに保存し、ステップＳ１１９へ移行する。
ステップＳ１１９では、「Ｘの更新（Ｘ←Ｚ）」の処理が実行される。この処理では、Ｘの更新としてＸにＳ１１８で保存したＺを代入し、ステップＳ１２０へ移行する。
ステップＳ１２０では、「全ＦＰ処理終了？」の判断処理が実行される。この処理では、Ｓ１１０で作成したｎ個のＦＰ全てが処理されたか否かが判断され、処理済み（ＹＥＳ）の場合は、ステップＳ１２１へ移行する。未処理のＦＰがある（ＮＯ）の場合は、未処理のＦＰに対してＳ１１５〜Ｓ１２０の処理を実行するため、Ｓ１１５へ移行する。全ＦＰの処理が終了するまでＳ１１５〜Ｓ１２０の処理を繰り返す。
ステップＳ１２１では、ｎ個のＦＰを再度順番に処理するためのカウンタの初期化としてｎに１が代入（ｎ←１）され、ステップＳ１２２へ移行する。
ステップＳ１２２では、「リストｎ番目のＦＰを読み込む」の処理が実行される。この処理では、リストｎ番目のＦＰを読み込み、ステップＳ１２３へ移行する。
ステップＳ１２３では、「各ピークのリテンション・タイム（Ｒ１）、ピーク・データ（Ｐ１）ならびにＵＶスペクトル（Ｕ１）を順番に取得」の処理が実行される。この処理では、Ｓ１２２で読み込んだＦＰからリテンション・タイム（Ｒ１）、ピーク・データ（Ｐ１）ならびにＵＶスペクトル（Ｕ１）を１ピークずつ順番に取得し、ステップＳ１２４へ移行する。
ステップＳ１２４では、「Ｘからリテンション・タイム（Ｒ２）を順番に取得」の処理が実行される。この処理では、全ＦＰのリテンション・タイムが重複なしで保存されているＸから1リテンション・タイム（Ｒ２）ずつ順番に取得し、ステップＳ１２５へ移行する。
ステップＳ１２５では、「Ｒ１＝Ｒ２？」の判断処理が実行される。この処理では、Ｓ１２３で取得したＲ１とＳ１２４で取得したＲ２が等しいか否かが判断され、等しい（ＹＥＳ）場合は、ステップＳ１２７へ移行する。等しくない（ＮＯ）場合は、ステップＳ１２６へ移行する。
ステップＳ１２６では、「Ｘの全リテンション・タイム比較終了？」の判断処理が実行される。この処理では、Ｓ１２３で取得したＲ１に対して、Ｘの全リテンション・タイムとの比較が終了しているか否かが判断される。終了している（ＹＥＳ）場合は、リテンション・タイムがＲ１のピークは処理済みであると判断し、次のピークへ処理を移行するため、ステップＳ１２３へ移行する。終了していない（ＮＯ）場合は、Ｘの次のリテンション。タイムに移行するため、ステップＳ１２４へ移行する。
ステップＳ１２７では、「Ｒ１に（ｎ−１）×分析時間（Ｔ）を加算（Ｒ１←Ｒ１＋（ｎ−１）×Ｔ）」の処理が実行される。この処理では、一番ピーク数の多いリスト１番目のＦＰに存在するピークのリテンション・タイムはそのままで、リスト１番目のＦＰには存在せず、リスト２番目のＦＰに存在するピークのリテンション・タイムはＲ１に分析時間（Ｔ）が加算され、さらに、リスト１〜ｎ−１番目のＦＰには存在せず、リストｎ番目のＦＰに存在するピークのリテンション・タイムはＲ１に（ｎ−１）×Ｔが加算され、ステップＳ１２８へ移行する。
ステップＳ１２８では、「Ｒ１、Ｐ１ならびにＵ１を出力（対象ＦＰ）」の処理が実行される。この処理では、Ｓ１２７で処理したＲ１とＳ１２３で取得したＰ１ならびにＵ１を対象ＦＰに出力し、ステップＳ１２９へ移行する。
ステップＳ１２９では、「ＸからＲ２を削除」の処理が実行される。この処理では、リテンション・タイムがＲ１（＝Ｒ２）における処理がＳ１２７、Ｓ１２８で終了したため、Ｘから処理済みリテンション・タイム（Ｒ２）を削除し、Ｓ１３０へ移行する。
ステップＳ１３０では、「全ピーク処理終了？」の判断処理が実行される。この処理では、リストｎ番目のＦＰの全ピークに対して処理が終了しているか否かが判断され、処理済み（ＹＥＳ）の場合は、リストｎ番目のＦＰにおけるＦＰ作成処理を終了し、ステップＳ１３１へ移行する。未処理のピークがある（ＮＯ）場合は、未処理のピークを処理するため、ステップＳ１２３へ移行する。全ピークの処理が終了するまでＳ１２３〜Ｓ１３０の処理を繰り返す。
ステップＳ１３１では、「ｎの更新（ｎ←ｎ＋１）」の処理が実行される。この処理では、処理を次のＦＰへ移行するため、ｎの更新としてｎにｎ＋１を代入し、ステップＳ１３２へ移行する。
ステップＳ１３２では、「全ＦＰ処理終了？」の判断処理が実行される。この処理では、Ｓ１１０で作成したｎ個のＦＰ全てが処理されたか否かが判断され、処理済み（ＹＥＳ）の場合は、ＦＰ作成処理を終了する。未処理のＦＰがある（ＮＯ）場合は、未処理のＦＰに対してＳ１２２〜Ｓ１３２の処理を実行するため、Ｓ１２２へ移行する。全ＦＰの処理が終了するまでＳ１２２〜Ｓ１３２の処理を繰り返す。
［Ｓ２：対象ＦＰ帰属処理１］
図８０は、図７６ステップＳ２の「対象ＦＰ帰属処理１」の詳細を示すフローチャートである。この処理は、帰属の前処理であり、正常品とされた複数の基準ＦＰから対象ＦＰ１７の帰属に適した基準ＦＰを選定する。
ステップＳ２０１では、「対象ＦＰを読み込む」の処理が実行される。この処理では、帰属対象のＦＰを読み込み、ステップＳ２０２へ移行する。
ステップＳ２０２では、「全リテンション・タイム（Ｒ１）を取得」の処理が実行される。この処理では、Ｓ２０１で読み込んだ対象ＦＰのリテンション・タイム情報を全て取得し、ステップＳ２０３へ移行する。
ステップＳ２０３では、「全基準ＦＰのファイル名をリスト化」の処理が実行される。この処理では、後で全基準ＦＰを順番に処理するために予め全基準ＦＰのファイル名をリスト化し、ステップＳ２０４へ移行する。
ステップＳ２０４では、全基準ＦＰを順番に処理するためのカウンタの初期値としてｎに１を代入（ｎ←１）し、ステップＳ２０５へ移行する。
ステップＳ２０５では、「リストｎ番目の基準ＦＰ（基準ＦＰｎ）を読み込む」の処理が実行される。この処理では、Ｓ２０３でリスト化した全基準ＦＰのファイル名リストのｎ番目のＦＰを読み込み、ステップＳ２０６へ移行する。
ステップＳ２０６では、「全リテンション・タイム（Ｒ２）を取得」の処理が実行される。この処理では、Ｓ２０５で読み込んだ基準ＦＰのリテンション・タイム情報を全て取得し、ステップＳ２０７へ移行する。
ステップＳ２０７では、「Ｒ１とＲ２のリテンション・タイム・出現パターンの一致度を算出（ＲＰｎ＿ｍｉｎ）」の処理が実行される。この処理では、Ｓ２０２で取得した対象ＦＰのリテンション・タイムとＳ２０６で取得した基準ＦＰのリテンション・タイムからＲＰｎ＿ｍｉｎを算出し、ステップＳ２０８へ移行する。なお、ＲＰｎ＿ｍｉｎの詳細な計算フローは、図８５のサブルーチン１により別途説明する。
ステップＳ２０８では、「ＲＰｎ＿ｍｉｎの保存（ＲＰａｌｌ＿ｍｉｎ」の処理が実行される。この処理では、Ｓ２０７で算出したＲＰｎ＿ｍｉｎをＲＰａｌｌ＿ｍｉｎに保存し、ステップＳ２０９へ移行する。
ステップＳ２０９では、「ｎの更新（ｎ←ｎ＋１）」の処理が実行される。この処理では、処理を次のＦＰへ移行するためｎの更新としてｎにｎ＋１を代入し、ステップＳ２１０へ移行する。
ステップＳ２１０では、「全基準ＦＰ処理終了？」の判断処理が実行される。この処理では、基準ＦＰ全てが処理されたか否かが判断され、処理済み（ＹＥＳ）の場合は、ステップＳ２１１へ移行する。未処理の基準ＦＰがある（ＮＯ）場合は、未処理のＦＰに対してＳ２０５〜Ｓ２１０の処理を実行するため、Ｓ２０５へ移行する。全基準ＦＰの処理が終了するまでＳ２０５〜Ｓ２１０の処理を繰り返す。
ステップＳ２１１では、「ＲＰａｌｌ＿ｍｉｎから一致度が最小となる基準ＦＰを選定」の処理が実行される。この処理では、全基準ＦＰに対して算出したＲＰ１＿ｍｉｎからＲＰｎ＿ｍｉｎを比較し、対象ＦＰとのリテンション・タイム・出現パターンの一致度が最小となる基準ＦＰを選定し、対象ＦＰ帰属処理１を終了する。
［Ｓ３：対象ＦＰ帰属処理２］
図８１は、図７６ステップＳ３の「対象ＦＰ帰属処理２」の詳細を示すフローチャートである。この処理は、帰属の本処理であり、対象ＦＰ１７とステップＳ２で選定した基準ＦＰとの間で、前記のようなピーク・パターン及びＵＶスペクトルの一致度から各帰属候補ピークの一致度（ＳＣＯＲＥ）を算出する。
ステップＳ３０１では、「対象ＦＰを読み込む」の処理が実行される。この処理では、帰属対象のＦＰを読み込み、ステップＳ３０２へ移行する。
ステップＳ３０２では、「帰属対象ピークのリテンション・タイム（Ｒ１）とピーク・データ（Ｐ１）ならびにＵＶスペクトル（Ｕ１）を順番に取得」の処理が実行される。この処理では、Ｓ３０１で読み込んだ対象ＦＰの各ピークを順番に帰属対象ピークとし、Ｒ１とＰ１ならびにＵ１を取得し、ステップＳ３０３へ移行する。
ステップＳ３０３では、「基準ＦＰを読み込む」の処理が実行される。この処理では、図８０の［対象ＦＰ帰属処理１］で選定された基準ＦＰを読み込み、ステップＳ３０４へ移行する。
ステップＳ３０４では、「基準ＦＰのピークのリテンション・タイム（Ｒ２）とピーク・データ（Ｐ２）ならびにＵＶスペクトル（Ｕ２）を順番に取得」の処理が実行される。この処理では、Ｓ３０３で読み込んだ基準ＦＰからＲ２とＰ２ならびにＵ２を１ピークずつ取得し、ステップＳ３０５へ移行する。
ステップＳ３０５では、「｜Ｒ１−（Ｒ２＋ｄ）｜＜閾値？」の判断処理が実行される。この処理では、ステップＳ３０２及びＳ３０４で読み込まれたＲ１とＲ２が閾値の範囲内で対応したものであるか否かが判断される。対応している（ＹＥＳ）場合は、リテンション・タイムがＲ２のピークはリテンション・タイムがＲ１のピークの帰属候補ピークであると判断し、帰属候補ピークの一致度（ＳＣＯＲＥ）を算出するためステップＳ３０６へ移行する。対応していない（ＮＯ）場合は、リテンション・タイムがＲ２のピークとリテンション・タイムがＲ１のピークではリテンション・タイムが違いすぎるため、帰属候補ピークにはならないと判断し、ステップＳ３０９へ移行する。なお、この判断処理でのｄは、対象ＦＰと基準ＦＰのピークのリテンション・タイムを補正する値で、初期値は０とし、処理を進める中で随時帰属したピーク間のリテンション・タイムの差を求め、その値でｄを更新する。また、閾値は、帰属候補ピークとすべきかどうかを判断するためのリテンション・タイムの許容幅である。
ステップＳ３０６では、「ＵＶスペクトルの一致度を算出（ＵＶ＿Ｓｉｍ）」の処理が実行される。この処理では、Ｓ３０２で取得した帰属対象ピークのＵ１とＳ３０４で取得した帰属候補ピークのＵ２から、ＵＶ＿Ｓｉｍを算出し、ステップＳ３０７へ移行する。なお、ＵＶ＿Ｓｉｍの詳細な計算フローは、図８６のサブルーチン２に別途記載する。
ステップＳ３０７では、「ピーク・パターンの一致度を算出（Ｐ＿Ｓｉｍ＿ｍｉｎ）」の処理が実行される。この処理では、Ｓ３０２で取得した帰属対象ピークのＲ１及びＰ１とＳ３０４で取得した帰属候補ピークのＲ２及びＰ２から、これらピークに対して網羅的にピーク・パターンを作成する。これらピーク・パターンのＰ＿Ｓｉｍ＿ｍｉｎを算出し、ステップＳ３０８へ移行する。なお、Ｐ＿Ｓｉｍ＿ｍｉｎの詳細な計算フローは、図８７のサブルーチン３に別途記載する。
ステップＳ３０８では、「帰属候補ピークの一致度を算出（ＳＣＯＲＥ）」の処理が実行される。この処理では、Ｓ３０６で算出したＵＶ＿ＳｉｍとＳ３０７で算出したＰ＿Ｓｉｍ＿ｍｉｎから、帰属対象ピークと帰属候補ピークのＳＣＯＲＥを、
ＳＣＯＲＥ ＝ ＵＶ＿Ｓｉｍ × Ｐ＿Ｓｉｍ＿ｍｉｎ
として算出し、ステップＳ３１０へ移行する。
ステップＳ３０９では、「ＳＣＯＲＥに８８８８８８を代入（ＳＣＯＲＥ←８８８８８８）」の処理が実行される。この処理では、帰属対象ピークの帰属候補ピークに該当しないピークのＳＣＯＲＥを８８８８８８とし、ステップＳ３１０へ移行する。
ステップＳ３１０では、「ＳＣＯＲＥの保存（ＳＣＯＲＥ＿ａｌｌ）」の処理が実行される。この処理では、Ｓ３０８あるいはＳ３０９で得られたＳＣＯＲＥをＳＣＯＲＥ＿ａｌｌに保存し、ステップＳ３１１へ移行する。
ステップＳ３１１では、「基準全ピークの処理終了？」の判断処理が実行される。この処理では、基準ＦＰの全ピークが処理されたか否かが判断され、処理済み（ＹＥＳ）の場合は、ステップＳ３１２へ移行する。未処理のピークがある（ＮＯ）場合は、未処理のピークに対してＳ３０４〜Ｓ３１１の処理を実行するため、Ｓ３０４へ移行する。全ピークの処理が終了するまでＳ３０４〜Ｓ３１１の処理を繰り返す。
ステップＳ３１２では、「判定結果ファイルにＳＣＯＲＥ＿ａｌｌを出力し、ＳＣＯＲＥ＿ａｌｌを初期化（空にする）」の処理が実行される。この処理では、判定結果ファイルにＳＣＯＲＥ＿ａｌｌを出力後、ＳＣＯＲＥ＿ａｌｌを初期化（空にする）し、ステップＳ３１３へ移行する。
ステップＳ３１３では、「対象全ピークの処理終了？」の判断処理が実行される。この処理では、対象ＦＰの全ピークが処理されたか否かが判断され、処理済み（ＹＥＳ）の場合は、対象ＦＰ帰属処理２を終了する。未処理のピークがある（ＮＯ）場合は、未処理のピークに対してＳ３０２〜Ｓ３１３の処理を実行するため、Ｓ３０２へ移行する。全ピークの処理が終了するまでＳ３０２〜Ｓ３１３の処理を繰り返す。
図９５で出力した判定結果ファイル例１２９を示す。
［Ｓ４：対象ＦＰ帰属処理３］
図８２は、図７６ステップＳ４の「対象ＦＰ帰属処理３」の詳細を示すフローチャートである。この処理は、帰属の後処理であり、前記のように算出した帰属候補ピークの一致度（ＳＣＯＲＥ）から対象ＦＰの各ピークに対応する基準ＦＰのピークを特定する。
ステップＳ４０１では、「判定結果ファイルを読み込む」の処理が実行される。この処理では、図８１の「対象ＦＰ帰属処理２」で作成した判定結果ファイルを読み込み、ステップＳ４０２へ移行する。
ステップＳ４０２では、「「ＳＣＯＲＥ＜閾値」の条件を満たしたデータで帰属候補ピークスコア表を作成」の処理が実行される。この処理では、判定結果ファイルのＳＣＯＲＥをもとに図９６（上図）の帰属候補スコア表１３１を作成し、ステップＳ４０３へ移行する。この帰属候補ピークスコア表は、基準ＦＰのピークごとに、対象ＦＰ全ピークに対して算出されたＳＣＯＲＥから、閾値より小さいＳＣＯＲＥだけを昇順に並べた表である。ちなにみ、このＳＣＯＲＥは値が小さいほど帰属すべきピークの可能性が高い。なお、閾値は、帰属候補とすべきかどうかを判断するためのＳＣＯＲＥの上限値である。
ステップＳ４０３では、「帰属候補ピーク番号表を作成」の処理が実行される。この処理では、帰属候補ピークスコア表をもとに図９６（下図）の帰属候補ピーク番号表１３３を作成し、ステップＳ４０４へ移行する。この帰属候補ピーク番号表は、帰属候補ピークスコア表の各スコアをそのスコアに対応する対象ＦＰのピーク番号に置き換えた表である。このことから、この表は、基準ＦＰのピークごとに対応させるべき対象ＦＰのピーク番号が順番に並んだ表となっている。
ステップＳ４０４では、「帰属すべき対象ＦＰのピーク番号を取得」の処理が実行される。この処理では、Ｓ４０３で作成した帰属候補ピーク番号表から、基準ＦＰのピークごとに一番上位に位置する対象ＦＰのピーク番号を取得し、ステップＳ４０５へ移行する。
ステップＳ４０５では、「取得したピーク番号が降順（重複なし）に並んでいる？」の判断処理が実行される。この処理では、Ｓ４０４で取得した対象ＦＰのピーク番号が重複なしで降順に並んでいるか否かが判断される。並んでいる（ＹＥＳ）場合は、基準ＦＰの各ピークに対応する対象ＦＰのピークが確定できたと判断し、ステップＳ４０８へ移行す
る。並んでいない（ＮＯ）場合は、問題のあった基準ＦＰのピークに帰属すべき対象ＦＰのピークを見直すため、ステップＳ４０６へ移行する。
ステップＳ４０６では、「問題のあったピーク間でＳＣＯＲＥを比較し、帰属候補ピーク番号表を更新」の処理が実行される。この処理では、問題のあった対象ＦＰのピーク番号に対応するＳＣＯＲＥを帰属候補スコア表で比較し、ＳＣＯＲＥが大きい方のピーク番号を２番目に位置するピーク番号に置き換えた帰属候補ピーク番号表に更新し、ステップＳ４０７へ移行する。
ステップＳ４０７では、「帰属候補ピークスコア表を更新」の処理が実行される。この処理では、Ｓ４０６での帰属候補ピーク番号表の更新内容に沿って、帰属候補ピークスコア表を更新し、ステップＳ４０４へ移行する。対象ＦＰのピーク番号に問題（重複あり、降順に並んでいない）がなくなるまで、Ｓ４０４からＳ４０７の処理を繰り返す。
ステップＳ４０８では、「帰属結果を保存（ＴＥＭＰ）」の処理が実行される。この処理では、基準ＦＰの全ピークのピーク番号、リテンション・タイム、及びピークとこれらピークに対応するピークとして特定した対象ＦＰのピーク・データをＴＥＭＰに保存し、ステップＳ４０９へ移行する。
ステップＳ４０９では、「ＴＥＭＰに対象ＦＰの全ピークが入っている？」の判断処理が実行される。この処理では、Ｓ４０８で保存したＴＥＭＰ中に対象ＦＰの全ピークのピーク・データが入っているか否かが判断される。全て入っている（ＹＥＳ）場合は、対象ＦＰの全てのピークで処理が終了したと判断し、Ｓ４１２へ移行する。入っていないピークがある（ＮＯ）場合は、入っていないピークのピーク・データをＴＥＭＰに追加するため、ステップＳ４１０へ移行する。ステップＳ４１０では、「ＴＥＭＰに入っていない対象ＦＰのピークのリテンション・タイムを補正」の処理が実行される。この処理では、ＴＥＭＰに入っていない対象ＦＰのピーク（補正が必要な対象ＦＰのピーク）のリテンション・タイムは、
補正値 ＝ ｋ１＋（ｋ２−ｋ１）＊（ｔ０−ｔ１）／（ｔ２−ｔ１）
ｋ１：補正が必要な対象ＦＰのピーク近傍で帰属された２つの基準ＦＰ側のピークのうちのリテンション・タイムが小さいピークのリテンション・タイム
ｋ２：補正が必要な対象ＦＰのピーク近傍で帰属された２つの基準ＦＰ側のピークのうちのリテンション・タイムが大きいピークのリテンション・タイム
ｔ０：補正が必要な対象ＦＰのピークのリテンション・タイム
ｔ１：補正が必要な対象ＦＰのピーク近傍で帰属された２つの対象ＦＰ側のピークのうちのリテンション・タイムが小さいピークのリテンション・タイム
ｔ２：補正が必要な対象ＦＰのピーク近傍で帰属された２つの対象ＦＰ側のピークのうちのリテンション・タイムが大きいピークのリテンション・タイム
として基準ＦＰにおけるリテンション・タイムに補正し、ステップＳ４１１へ移行する。
ステップＳ４１１では、「補正したリテンション・タイムとそのピークのピーク・データをＴＥＭＰに追加し、ＴＥＭＰを更新」の処理が実行される。この処理では、Ｓ４１０で補正したＴＥＭＰに入っていない対象ＦＰのピークのリテンション・タイムとＴＥＭＰ中の基準ＦＰのリテンション・タイムとを比較し、ＴＥＭＰ中の妥当な位置にＴＥＭＰに入っていない対象ＦＰのピークの補正したリテンション・タイムならびにピーク・データを追加し、ＴＥＭＰを更新し、Ｓ４０９へ移行する。対象ＦＰの全ピークが追加されるまで、Ｓ４０９からＳ４１１の処理を繰り返す。
ステップＳ４１２では、「ＴＥＭＰを照合結果ファイルに出力」の処理が実行される。この処理では、基準ＦＰの全ピークと対象ＦＰの全ピークの対応関係を特定したＴＥＭＰを照合結果ファイルとして出力し、対象ＦＰ帰属処理３を終了する。
図９７に前記のように出力した照合結果ファイル例１３５を示す。
［Ｓ５：対象ＦＰ帰属処理４］
図８３、図８４は、図７６ステップＳ５の「対象ＦＰ帰属処理４」の詳細を示すフローチャートである。この処理は、帰属の最終処理であり、図７６ステップＳ４で作成した照合結果ファイルをもとに対象ＦＰの各ピークを基準群ＦＰのピークに帰属する。なお、基準群ＦＰは、前記のように全基準ＦＰ間でピークの対応関係を特定したＦＰであり、図９８の基準群ＦＰデータ例１３７のように、基準群ＦＰピーク番号と基準群リテンション・タイムとピーク高さとで構成されたデータである。図２の基準群ＦＰ１９で示すように各ピークは平均値（黒点）±標準偏差（縦線）で表すことができる。
ステップＳ５０１では、「照合結果ファイルを読み込む」の処理が実行される。この処理では、図８２のＳ４１２で出力した照合結果ファイルを読み込み、ステップＳ５０２へ移行する。
ステップＳ５０２では、「基準群ＦＰを読み込む」の処理が実行される。この処理では、対象ＦＰの各ピークの最終的な帰属相手である基準群ＦＰを読み込み、ステップＳ５０３へ移行する。
ステップＳ５０３では、「対象ＦＰと基準群ＦＰを統合し保存（ＴＥＭＰ）」の処理が実行される。この処理では、照合結果ファイルと基準群ＦＰで共通に存在する基準ＦＰのピーク・データをもとに２つのファイルを統合し、その結果をＴＥＭＰとして保存し、ステップＳ５０４へ移行する。
ステップＳ５０４では、「基準ＦＰに対応するピークがない対象ＦＰのピークのリテンション・タイムを補正」の処理が実行される。この処理では、照合結果ファイルで基準ＦＰに対応するピークがない対象ＦＰの全ピークのリテンション・タイムをＳ５０３で保存したＴＥＭＰのリテンション・タイムに補正し、ステップＳ５０５へ移行する。なお、リテンション・タイムの補正は、前記ステップＳ４の「対象ＦＰ帰属処理３」のステップＳ４１０と同様の方法で補正する。
ステップＳ５０５では、「補正したリテンション・タイム（Ｒ１、Ｒ３）と対応するピーク・データ（Ｐ１）を順番に取得」の処理が実行される。この処理では、Ｓ５０４で補正したリテンション・タイムをＲ１ならびにＲ３として、対応するピークのピーク・データをＰ１として順番に取得し、ステップＳ５０６へ移行する。
ステップＳ５０６では、「ＴＥＭＰから対象ＦＰの帰属候補ピークのリテンション・タイム（Ｒ２）と対応するピーク・データ（Ｐ２）を順番に取得」の処理が実行される。この処理では、Ｓ５０３で保存したＴＥＭＰから対象ＦＰのピークが帰属されていないリテンション・タイムをＲ２として、対応するピーク・データをＰ２として順番に取得し、ステップＳ５０７へ移行する。
ステップＳ５０７では、「｜Ｒ１−Ｒ２｜＜閾値１？」の判断処理が実行される。この処理では、Ｓ５０５とＳ５０６で取得したＲ１とＲ２の差が閾値１より小さいか否かが判断される。小さい（ＹＥＳ）場合は、対象ＦＰのリテンション・タイムがＲ１のピークと基準ＦＰのリテンション・タイムがＲ２のピークが対応する可能性があると判断し、ステップＳ５０８へ移行する。Ｒ１とＲ２の差が閾値１以上（ＮＯ）の場合は、対応する可能性なしと判断し、ステップＳ５１２へ移行する。
ステップＳ５０８では、「Ｒ１、Ｒ２に対応するＵＶスペクトルを取得（Ｕ１、Ｕ２）」の処理が実行される。この処理では、Ｓ５０７で対応する可能性ありと判断されたリテンション・タイムがＲ１とＲ２のピークに対応するＵＶスペクトルをそれぞれのＦＰから取得し、ステップＳ５０９へ移行する。
ステップＳ５０９では、「ＵＶスペクトルの一致度を算出（ＵＶ＿Ｓｉｍ）」の処理が実行される。この処理では、Ｓ５０８で取得したＵＶスペクトルＵ１及びＵ２からステップＳ３の「対象ＦＰ帰属処理２」のステップＳ３０６と同様な方法でＵＶ＿Ｓｉｍを算出し、ステップＳ５１０へ移行する。なお、ＵＶ＿Ｓｉｍの詳細な計算フローは図８６のサブルーチン２で別途説明する。
ステップＳ５１０では、「ＵＶ＿Ｓｉｍ＜閾値２？」の判断処理が実行される。この処理では、Ｓ５０９で算出したＵＶ＿Ｓｉｍが閾値２より小さいか否かが判断される。小さい（ＹＥＳ）場合は、ＵＶスペクトルがＵ１のピークとＵ２のピークが対応していると判断し、ステップＳ５１１へ移行する。ＵＶ＿Ｓｉｍが閾値２以上（ＮＯ）の場合は、対応していないと判断し、ステップＳ５０７へ移行する。
ステップＳ５１１では、「Ｒ３ ← Ｒ２、閾値２ ← ＵＶ＿Ｓｉｍ」の処理が実行される。この処理では、Ｓ５１０で対応すると判断したリテンション・タイムがＲ３（つまり、Ｒ１）を対応相手のリテンション・タイムであるＲ２に更新したのち、閾値２をＵＶ＿Ｓｉｍの値に更新し、Ｓ５０７に移行する。
ステップＳ５１２では、「全ての帰属候補ピークのリテンション・タイムを比較終了？」の判断処理が実行される。この処理では、Ｒ１と全ての帰属候補ピークのリテンション・タイムの比較が終了したか否かが判断され、終了している（ＹＥＳ）場合は、ステップＳ５１３へ移行する。終了していない（ＮＯ）場合は、ステップＳ５０７へ移行する。
ステップＳ５１３では、「Ｒ１、Ｒ３とＰ１ならびに閾値２を保存（ＴＥＭＰ２）」の処理が実行される。この処理では、Ｓ５１０で対応すると判断したリテンション・タイム（Ｒ１）と対応相手のリテンション・タイム（Ｒ２）に更新したＲ３と対応するピーク（Ｐ１）ならびに現時点の閾値２を保存（ＴＥＭＰ２）し、Ｓ５０７に移行する。
ステップＳ５１４では、「全ての非対応ピークのリテンション・タイムを比較終了？」の判断処理が実行される。この処理では、全ての非対応ピークのリテンション・タイムで帰属候補ピークのリテンション・タイムとの比較が終了したか否かが判断される。終了している（ＹＥＳ）場合は、全ての非対応ピークの帰属処理が終了したと判断し、ステップＳ５１６へ移行する。終了していない（ＮＯ）場合は、未処理の非対応ピークが残っていると判断し、ステップＳ５１５へ移行する。
ステップＳ５１５では、「閾値２ ← 初期値」の処理が実行される。この処理では、Ｓ５１１でＵＶ＿Ｓｉｍに更新されている閾値２を初期値に戻し、ステップＳ５０５へ移行する。
ステップＳ５１６では、「ＴＥＭＰ２にＲ３の値が同じピークが存在する？」の判断処理が実行される。この処理では、ＴＥＭＰ中の同じピークに複数の非対応ピークが帰属されているか否かが判断される。同じピークに帰属された非対応ピークが存在する（ＹＥＳ）場合は、ステップＳ５１７へ移行する。存在しない（ＮＯ）場合は、ステップＳ５１８へ移行する。
ステップＳ５１７では、「Ｒ３の値が同じピークの閾値２を比較し、値が大きいピークのＲ３を元の値（Ｒ１）に戻す」の処理が実行される。この処理では、ＴＥＭＰ２中のＲ３の値が同じピークの閾値２を比較し、値が大きいピークのＲ３の値を元の値（つまり、Ｒ１）に戻し、ステップＳ５１８へ移行する。
ステップＳ５１８では、「ＴＥＭＰにＴＥＭＰ２のピークを追加（ＴＥＭＰのリテンション・タイムとＲ３が一致したピークのみ）」の処理が実行される。この処理では、ＴＥＭＰのリテンション・タイムとＲ３が一致したピークのみＴＥＭＰにＲ３に対応するピークを追加し、ステップＳ５１９へ移行する。Ｒ３がＴＥＭＰのリテンション・タイムと一致しないピークは、基準群ＦＰに帰属相手となるピークが存在しないため、追加しない。
ステップＳ５１９では、「ＴＥＭＰ中の対象ＦＰのピークを出力（ピーク特徴量ファイル）」の処理が実行される。この処理では、基準群ＦＰ１３７に帰属された対象ＦＰのピーク・データをピーク・データ特徴量ファイルとして出力し、対象ＦＰ帰属処理４を終了する。
図９９に前記のように出力するピーク・データ特徴量のファイル例１３９を示す。
[サブルーチン１]
図８５は、図８０の「基準ＦＰ選定処理」における「サブルーチン１」の詳細を示すフローチャートである。この処理は、ＦＰ間（例えば、対象ＦＰと基準ＦＰ）のリテンション・タイム・出現パターンの一致度を計算する。
ステップS１００１では、「ｘ←Ｒ１、ｙ←Ｒ２」の処理が実行される。この処理では、図８０のＳ２０２とＳ２０６で取得したＲ１及びＲ２をそれぞれｘとｙに代入し、ステップＳ１００２へ移行する。
ステップＳ１００２では、「ｘ、ｙのデータ数を取得（ａ、ｂ）」の処理が実行される。この処理では、ｘ、ｙのデータ数をそれぞれａ、ｂとして取得し、ステップＳ１００３へ移行する。
ステップＳ１００３では、ｘのリテンション・タイムを順番に呼び出すためのカウンタの初期値としてｉに１を代入（ｉ←１）し、ステップＳ１００４へ移行する。
ステップＳ１００４では、「ｘｉ番目のリテンション・タイムからの全距離を取得（ｆ）」の処理が実行される。この処理では、ｘｉ番目のリテンション・タイムとそれ以降の全リテンション・タイム間距離をｆとして取得し、ステップＳ１００５へ移行する。
ステップＳ１００５では、ｙのリテンション・タイムを順番に呼び出すためのカウンタの初期値としてｊに１を代入（ｊ←１）し、ステップＳ１００６へ移行する。
ステップＳ１００６では、「ｙｊ番目のリテンション・タイムからの全距離を取得（ｇ）」の処理が実行される。この処理では、ｙｊ番目のリテンション・タイムとそれ以降の全リテンション・タイム間距離をｇとして取得し、ステップＳ１００７へ移行する。
ステップＳ１００７では「”｜ｆの各リテンション・タイム間距離−ｇの各リテンション・タイム間距離｜＜閾値”の条件を満たしたデータ数を取得（ｍ）」の処理が実行される。この処理では、Ｓ１００４及びＳ１００６で取得したリテンション・タイム間距離ｆとｇを総当りで比較し、”｜ｆの各リテンション・タイム間距離−ｇの各リテンション・タイム間距離｜＜閾値”の条件を満たしたデータ数をｍとして取得し、ステップＳ１００８へ移行する。
ステップＳ１００８では、「ｆとｇのリテンション・タイム・出現パターンの一致度を算出（ＲＰｆｇ）」の処理が実行される。この処理では、Ｓ１００２で取得したａ、ｂとＳ１００７で取得したｍからＲＰｆｇを、
ＲＰｆｇ ＝ （１−（ｍ／（ａ＋ｂ−ｍ）））×（ａ−ｍ＋１）
として算出し、ステップＳ１００９へ移行する。
ステップＳ１００９では、「ＲＰｆｇを保存（ＲＰ＿ａｌｌ）」の処理が実行される。この処理では、Ｓ１００８で算出した一致度をＲＰ＿ａｌｌに保存し、ステップＳ１０１０へ移行する。
ステップＳ１０１０では、「ｊの更新（ｊ←ｊ＋１）」の処理が実行される。この処理では、ｙの処理を次のリテンション・タイムへ移行するためｊの更新としてｊにｊ＋１を代入し、ステップＳ１０１１へ移行する。
ステップＳ１０１１では、「ｙの全リテンション・タイムで処理終了？」の判断処理が実行される。この処理では、ｙの全てのリテンション・タイムの処理が終了したか否かが判断される。終了している（ＹＥＳ）場合は、ｙの全リテンション・タイムの処理が終了したと判断し、ステップＳ１０１２へ移行する。終了していない（ＮＯ）場合は、ｙ中に未処理のリテンション・タイムが残っていると判断し、ステップＳ１００６へ移行する。つまり、Ｓ１００６〜Ｓ１０１１までの処理はｙの全てのリテンション・タイムが処理されるまで繰り返す。
ステップＳ１０１２では、「ｉの更新（ｉ←ｉ＋１）」の処理が実行される。この処理では、ｘの処理を次のリテンション・タイムへ移行するためｉの更新としてｉにｉ＋１を代入し、ステップＳ１０１３へ移行する。
ステップＳ１０１３では、「ｘの全リテンション・タイムで処理終了？」の判断処理が実行される。この処理では、ｘの全てのリテンション・タイムの処理が終了したか否かが判断される。終了している（ＹＥＳ）場合は、ｘの全リテンション・タイムの処理が終了したと判断し、ステップＳ１０１４へ移行する。終了していない（ＮＯ）場合は、ｘ中に未処理のリテンション・タイムが残っていると判断し、ステップＳ１００４へ移行する。つまり、Ｓ１００４〜Ｓ１０１３までの処理はｘの全てのリテンション・タイムが処理されるまで繰り返す。
ステップＳ１０１４では、「ＲＰ＿ａｌｌから最小値を取得（ＲＰ＿ｍｉｎ）」の処理が実行される。この処理では、対象ＦＰと基準ＦＰとのリテンション・タイム・出現パターンの全組み合わせでのＲＰが保存されたＲＰ＿ａｌｌ中の最小値を、ＲＰ＿ｍｉｎとして取得し、そのＲＰ＿ｍｉｎを図８０のＳ２０７に渡し、リテンション・タイム・出現パターンの一致度計算処理を終了する。
[サブルーチン２]
図８６は、図８１の「対象ＦＰ帰属処理２」における「サブルーチン２」の詳細を示すフローチャートである。この処理は、ＵＶスペクトルの一致度を計算する。
ステップＳ２００１では、「ｘ←Ｕ１、ｙ←Ｕ２、ｚ←０」の処理が実行される。この処理では、図８１のＳ３０２とＳ３０４で取得したＵＶスペクトルＵ１及びＵ２をそれぞれｘとｙに代入し、さらにＵＶスペクトル間の距離の二乗和（ｚ）の初期値として０を代入し、ステップＳ２００２へ移行する。
ステップＳ２００２では、「ｘデータ数を取得（ａ）」の処理が実行される。この処理では、ｘのデータ数をａとして取得し、ステップＳ２００３へ移行する。
ステップＳ２００３では、「ｉ←１」の処理が実行される。この処理では、ＵＶスペクトルＵ１及びＵ２を構成する各検出波長における吸光度をｘ及びｙから順番に呼び出すための初期値としてｉに１を代入し、ステップＳ２００４へ移行する。
ステップＳ２００４では、「ｘｉ番目のデータを取得（ｂ）」の処理が実行される。この処理では、ＵＶスペクトルＵ１を代入したｘのｉ番目の吸光度データをｂとして取得し、ステップＳ２００５へ移行する。
ステップＳ２００５では、「ｙｉ番目のデータを取得（ｃ）」の処理が実行される。この処理では、ＵＶスペクトルＵ２を代入したｙのｉ番目の吸光度データをｃとして取得し、ステップＳ２００６へ移行する。
ステップＳ２００６では、「ＵＶスペクトル間距離（ｄ）とＵＶスペクトル間距離の二乗和（ｚ）を算出」の処理が実行される。この処理では、ＵＶスペクトル間距離ｄとＵＶスペクトル間距離の二乗和ｚを、
ｄ ＝ ｂ−ｃ
ｚ ＝ ｚ＋ｄ2
として算出し、ステップＳ２００７へ移行する。
ステップＳ２００７では、「ｉの更新（ｉ←ｉ＋１）」の処理が実行される。この処理では、ｉの更新としてｉにｉ＋１を代入し、ステップＳ２００８へ移行する。
ステップＳ２００８では、「ｘの全データで処理終了？」の判断処理が実行される。この処理では、ｘとｙの全てのデータの処理が終了したか否かが判断される。終了している（ＹＥＳ）場合は、ｘとｙの全データの処理が終了したと判断し、ステップＳ２００９へ移行する。終了していない（ＮＯ）場合は、ｘとｙに未処理のデータが残っていると判断し、ステップＳ２００４へ移行する。つまり、Ｓ２００４〜Ｓ２００８までの処理はｘとｙの全ての吸光度データが処理されるまで繰り返す。
ステップＳ２００９では、「ｘとｙのＵＶスペクトルの一致度を計算（ＵＶ＿Ｓｉｍ）」の処理が実行される。この処理では、ＵＶ＿Ｓｉｍを、前記ＵＶスペクトル間距離の二乗和ｚとｘのデータ数ａから、
ＵＶ＿Ｓｉｍ ＝ √（ｚ／ａ）
として算出し、このＵＶ＿Ｓｉｍを図８１のステップＳ３０６に渡し、ＵＶスペクトルの一致度計算処理を終了する。
[サブルーチン３]
図８７は、図８１の「対象ＦＰ帰属処理２」における「サブルーチン３」の詳細を示すフローチャートである。この処理は、ピーク・パターンの一致度を計算する。
ステップＳ３００１では、「ピーク・パターン構成候補数（ｍ）とピーク・パターン構成ピーク数（ｎ）を設定」の処理が実行される。この処理では、ピーク・パターンを網羅的に作成するための設定として、ピーク・パターン構成候補数（ｍ）とピーク・パターン構成ピーク数（ｎ）をそれぞれ設定し、ステップ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)
として算出し、ステップＳ３０１２へ移行する。
ステップＳ３０１２では、「Ｐ＿Ｓｉｍを保存（Ｐ＿Ｓｉｍ＿ａｌｌ）」の処理が実行される。この処理では、Ｓ３０１１で算出されたＰ＿Ｓｉｍを順次Ｐ＿Ｓｉｍ＿ａｌｌに保存し、ステップＳ３０１３へ移行する。
ステップＳ３０１３では、「ｅｍ中のｍ本からｎ本を取り出す全組み合わせ終了？」の判断処理が実行される。この処理では、帰属候補ピークのピーク・パターン構成候補ピークｍ本からピーク・パターン構成ピークｎ本を取り出す全組み合わせで処理が終了したか否かが判断される。終了した（ＹＥＳ）場合は、帰属候補ピークにおいて、網羅的なピーク・パターンの作成とそのパターンにおける一致度の計算が終了したと判断し、ステップＳ３０１４へ移行する。終了していない（ＮＯ）場合は、ｍ本からｎ本を取り出す組み合わせが終了していないと判断し、ステップＳ３０１０へ移行する。つまり、Ｓ３０１０〜Ｓ３０１３までの処理はｍ本からｎ本を取り出す全組み合わせで処理が終了するまで繰り返す。
ステップＳ３０１４では、「ｃｍ中のｍ本からｎ本を取り出す全組み合わせ終了？」の判断処理が実行される。この処理では、帰属対象ピークのピーク・パターン構成候補ピークｍ本からピーク・パターン構成ピークｎ本を取り出す全組み合わせで処理が終了したか否かが判断される。終了した（ＹＥＳ）場合は、帰属対象ピークにおいて、網羅的なピーク・パターンの作成とそのパターンにおける一致度の計算が終了したと判断し、ステップＳ３０１５へ移行する。終了していない（ＮＯ）場合は、ｍ本からｎ本を取り出す組み合わせが終了していないと判断し、ステップＳ３００９へ移行する。つまり、Ｓ３００９〜Ｓ３０１４までの処理はｍ本からｎ本を取り出す全組み合わせで処理が終了するまで繰り返す。
ステップＳ３０１５では、「Ｐ＿Ｓｉｍ＿ａｌｌから最小値を取得（Ｐ＿Ｓｉｍ＿ｍｉｎ）」の処理が実行される。この処理では、Ｓ３０１２で保存したＰ＿Ｓｉｍ＿ａｌｌの最小値をＰ＿Ｓｉｍ＿ｍｉｎとして取得し、このＰ＿Ｓｉｍ＿ｍｉｎを図８１のステップＳ３０７に渡し、ピーク・パターンの一致度計算処理を終了する。
［基準群ＦＰの作成］
対象ＦＰ特徴量データを基準ＦＰ特徴量データと比較するための基準ＦＰ特徴量ファイルは、図８８〜図９１のように作成される。
図８８は、基準ＦＰ特徴量統合ファイルを作成するためのフローチャートであり、基準ＦＰ作成部のＦＰ作成機能と、基準ＦＰピーク帰属部の基準ＦＰピーク帰属機能と、基準ＦＰ帰属結果統合部の基準ＦＰ帰属結果統合機能と、基準ＦＰピーク特徴量作成部の基準ＦＰピーク特徴量作成機能とをコンピュータに実現させる。
基準ＦＰ作成機能は、ステップＳ１０００１で実現される。基準ＦＰピーク帰属機能は、ステップＳ１０００２、Ｓ１０００３、Ｓ１０００４で実現される。基準ＦＰ帰属結果統合機能は、ステップＳ１０００５で実現される。基準ＦＰピーク特徴量作成機能は、ステップＳ１０００６で実現される。
Ｓ１０００１〜Ｓ１０００４は、図７６の対象ＦＰ特徴量統合ファイルの作成に係るＳ１〜Ｓ４に対応している。
ステップＳ１０００１は、３Ｄクロマト及び特定の検出波長におけるピーク情報を入力データとして「ＦＰ作成処理」が実行される。
３Ｄクロマト及びピーク・データ共に、評価基準となる複数の評価基準薬剤（基準漢方薬）ごとに備えられている。
ステップＳ１０００１では、コンピュータの基準ＦＰ作成部が機能し、３Ｄクロマト及びピーク情報から基準ＦＰが対象ＦＰ１７（図２）と同様に作成され、基準ＦＰのデータがファイルとして出力される。
ステップＳ１０００２は、ステップＳ１０００１で出力した全基準ＦＰを入力として、「基準ＦＰ帰属処理１」が実行される。
ステップＳ１０００２では、コンピュータの基準ＦＰピーク帰属部が機能し、全基準ＦＰを対象とし、選択した組み合わせ且つ順番で帰属スコアを算出するために全基準ＦＰから組み合わせを選択してステップＳ１０００３へ移行する。
ステップＳ１０００３は、選択した基準ＦＰの組み合わせを入力とし、「基準ＦＰ帰属処理２」が実行される。
ステップＳ１０００３では、ステップＳ２で選定した基準ＦＰの組み合わせの全ピークで、図２３〜図６１のように網羅的にピーク・パターンを作成し、次にそれらピーク・パターンの一致度（図６３または図６４のＰ＿Ｓｉｍ）を算出する。また、選定した基準ＦＰの組み合わせのピーク間でＵＶスペクトルの一致度（図６６のＵＶ＿Ｓｉｍ）を算出する。さらにこれら２つの一致度から帰属候補ピークの一致度（図６７のＳＣＯＲＥ）を算出する。その算出結果は、判定結果ファイル（図９５の判定結果ファイル例１２９参照）として出力される。
ステップＳ１０００４は、ステップＳ１０００３で出力した判定結果ファイルを入力とし、「基準ＦＰ帰属処理３」が実行される。
ステップＳ１０００４では、選定した基準ＦＰの組み合わせ間で、帰属候補ピークの一致度（ＳＣＯＲＥ）をもとに選定した基準ＦＰの組み合わせ間で対応するピークを特定する。その結果は、基準ＦＰごとに基準ＦＰ帰属データとして出力される。
ステップＳ１０００５は、ステップＳ１０００４で出力した全基準ＦＰ帰属データを入力とし、「基準ＦＰ帰属結果統合処理」が実行される。
ステップＳ１０００５では、コンピュータの基準ＦＰ帰属結果統合部が機能し、基準ＦＰピーク帰属部で特定した各基準ＦＰのピーク対応関係を参照し、全基準ＦＰ帰属データを統合して基準ＦＰ対応表を作成し、ステップＳ１０００６へ移行する。
ステップＳ１０００６では、コンピュータの基準ＦＰピーク特徴量作成部が機能し、基準ＦＰ帰属結果統合部で作成された基準ＦＰ対応表を基に全基準ＦＰによるピーク特徴量（基準群ＦＰ）を作成する。この基準ＦＰピーク特徴量作成部での処理は、基準ＦＰ対応表の各ピーク（列）で統計量（最大値、最小値、中央値、平均値、ｅｔｃ）を算出し、その情報を基にピーク（列）を選定する。選定したピーク（列）は、基準群ＦＰ（図９８の基準群ＦＰ例１３７参照）として出力される。
［Ｓ１０００５：基準ＦＰ対応表の作成］
図８９、９０は、ステップＳ１０００５の「基準ＦＰ帰属結果統合処理（基準ＦＰ対応表の作成）」の詳細を示すフローチャートである。
ステップＳ１０１０１では、「帰属順番１番の帰属データを統合データとして読み込む」の処理が実行される。この処理では、Ｓ１０００４で１番目に帰属処理しピークの対応関係を特定した基準ＦＰ帰属データを統合データとして読み込み、ステップＳ１０１０２へ移行する。
ステップＳ１０１０２では、「帰属データ２番目以降を順番に読み込む」の処理が実行される。この処理では、まずＳ１０００４で２番目に帰属処理しピークの対応関係を特定した基準ＦＰ帰属データを統合データのとして読み込み、ステップＳ１０１０３へ移行する。
ステップＳ１０１０３では、「統合データと帰属データを共通するピーク・データで統合」の処理が実行される。この処理では、統合データと帰属データで共通に存在する基準ＦＰのピーク・データをもとに２つのファイルを統合し、その結果で統合データを更新し、ステップＳ１０１０４へ移行する。
ステップＳ１０１０４では、「帰属データ中の全てのピークを統合データに追加？」の判断処理が実行される。この処理では、帰属データの全てのピークが統合データに追加されたか否かが判断される。追加された（ＹＥＳ）場合は、ステップＳ１０１０５へ移行する。追加されていないピーク（欠落ピーク）がある（ＮＯ）場合は、この欠落ピークを統合データに追加処理するため、ステップＳ１０１０７へ移行する。尚、欠落ピークの統合データへの追加処理（Ｓ１０１０７−Ｓ１０１２０）は、前記Ｓ５（対象ＦＰ帰属処理４）におけるステップＳ５０４−Ｓ５１７と同様の処理が行われる。
ステップＳ１０１２１では、「統合データにＴＥＭＰ２のデータを追加（全てのリテンション・タイムとピーク）」の処理が実行される。この処理では、ＴＥＭＰ２の全てのリテンション・タイム（Ｒ３）とピーク（Ｐ１）を統合データの該当箇所に追加し、ステップＳ１０１２２へ移行する。
ステップＳ１０１２２では、「閾値２ ← 初期値、ＴＥＭＰ２内のデータを全て削除」の処理が実行される。この処理では、ＵＶ＿Ｓｉｍに更新されている閾値２を初期値に戻し、全欠落ピークのリテンション・タイムやピークなどのデータが入ったＴＥＭＰ２から全てのデータを削除し、ステップＳ１０１０４へ戻る。
ステップＳ１０１０４から移行するステップＳ１０１０５では、「全ての帰属データの処理終了？」の判断処理が実行される。この処理では、全基準データの処理が終了したか否かが判断される。処理が終了（ＹＥＳ）の場合は、全帰属データの統合結果である基準ＦＰ対応表を出力するため、ステップＳ１０１０６へ移行する。全ての処理が終了していない（ＮＯ）の場合は、ステップＳ１０１０２へ戻り、残りの帰属データを順次処理する。
ステップＳ１０１０６では、「統合データを出力（基準ＦＰ対応表）」の処理が実行される。この処理では、全帰属データを統合した結果を基準ＦＰ対応表として出力し、基準ＦＰ対応表の作成処理を終了する。
［Ｓ１０００６：ピーク特徴量化処理］
図９１は、図８８ステップＳ１０００６の「ピーク特徴量化処理（基準群ＦＰの作成）」の詳細を示すフローチャートである。
ステップＳ１０２０１では、「基準ＦＰ対応表を読み込む」の処理が実行される。この処理では、Ｓ１０００５で作成した基準ＦＰ対応表を読み込み、ステップＳ１０２０２へ移行する。
ステップＳ１０２０２では、「各ピーク（列）で統計量を算出」の処理が実行される。この処理では、基準ＦＰ対応表の各ピーク(列)で統計量（最大値、最小値、中央値、平均値、分散、標準偏差、存在数、存在率）を算出し、ステップＳ１０２０３へ移行する。
ステップＳ１０２０３では、「算出した統計量を参考にピーク（列）を選定」の処理が実行される。この処理では、Ｓ１０１０２で算出した統計量を参考にピークを選定し、ステップＳ１０２０４へ移行する。
ステップＳ１０２０４では、「選定したピーク（列）を出力（基準群ＦＰ）」の処理が実行される。この処理では、統計量によりピーク（列）の選定結果を基準群ＦＰとして出力し、基準群ＦＰの作成処理を終了する。
図９８に前記のように出力する基準ＦＰ対応表例１３７を示す。
［実施例１の効果］
本発明実施例１では、評価対象となる多成分薬剤の３Ｄクロマト１５から特定の波長、例えば２０３ｎｍで検出されたピークとそのピークのリテンション・タイム、ならびにそのピークのＵＶスペクトルで構成する対象ＦＰ１７を作成するＦＰ作成工程１１３と、対象ＦＰ１７のピーク帰属に適した基準ＦＰを複数の基準ＦＰから選定する基準ＦＰ選定工程１１５と、対象ＦＰと選定された基準ＦＰのピークそれぞれに対し、例えば時間軸方向前後の少なくとも一方に存在する２本のピークを含めた３本のピークで構成されるピーク・パターンを作成するピーク・パターン作成工程１１７と、前記ピーク・パターン及びピークのＵＶスペクトルを比較し対応するピークを特定するピーク帰属工程１１９と、帰属したピークと複数の基準ＦＰのピークを例えばＭＴ法により比較評価する評価工程１２１とを備えた。
評価対象となる多成分薬剤の３Ｄクロマト１５をこれら５工程（１１３、１１５、１１７、１１９、１２１）で処理することで、評価対象薬剤の品質評価の精度及び効率をより向上させることができる。
ＦＰ作成工程１１３で作成された対象ＦＰ１７は、３Ｄクロマト１５と同様に三次元の情報（ピーク、リテンション・タイム及びＵＶスペクトル）で構成している。そのため、その薬剤特有の情報をそのまま継承したデータである。それにも係らず、データ容量は約１／７０に圧縮されているため、３Ｄクロマト１５に比較して、処理すべき情報量を大幅に減少させることができ処理速度を速めることができる。
ＦＰ作成工程１１３は、検出波長の異なる複数のＦＰを合成したＦＰを作成する。これにより、１つの波長で全ての成分を検出できない成分が組み合わさった多成分薬剤であっても、複数の検出波長のＦＰを合成することで、全ての成分を含んだ品質評価が可能となる。
ＦＰ作成工程１１３は、３Ｄクロマトで検出された全ピークを含んだＦＰを作成する。このため、多成分薬剤である漢方薬の品質評価に適している。
基準ＦＰ選定工程１１５では、対象ＦＰの帰属に適した基準ＦＰをＦＰ間のリテンション・タイム・出現パターンを比較し、パターンの一致度が良い基準ＦＰ選定する。これにより、ピーク帰属処理工程１１９において、パターンが類似したＦＰ間で帰属処理ができるため、精度の高い帰属が可能となる。
ピーク・パターン作成工程１１７では、帰属対象ピーク及び帰属候補ピークそれぞれで複数の周辺ピークを使って網羅的にピーク・パターンを作成する。これにより、対象ＦＰと基準ＦＰでＦＰ全体のパターンが多少異なっていても、ピーク帰属処理工程１１９で精度の高い帰属が可能となる。
ピーク帰属工程１１９では、ピーク・パターン作成工程１１７で作成されたピーク・パターンの一致度に加え、帰属対象ピークと帰属候補ピークのＵＶスペクトルの一致度も加味して、帰属すべきピークを特定している。そのため、精度の高い帰属が可能となる。
ピーク帰属工程１１９では、対象ＦＰの全ピークを基準ＦＰのピークへ一斉帰属する。そのため、効率のよい帰属処理が可能となる。
評価工程１２１では、多次元データである多成分で構成されたＦＰをＭＴ法でＭＤ値として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つの値の比（大きい値/小さい値） であることを示している。
なお、ＦＰをピーク高さで作成した場合にも、比に基づいてピーク・パターンの一致度（Ｐ＿Ｓｉｍ）を計算することができ、ＦＰをピーク面積で作成した場合にも、上記ピーク高さの差同様に、ピーク面積値の差に基づいてピーク・パターンの一致度（Ｐ＿Ｓｉｍ）を得ることができる。
［サブルーチン２の変形例］
図１００は、図８６に代えて適用するサブルーチン２の変形例に係り、図８１の「対象ＦＰ帰属処理２」における「サブルーチン２」の変形例の詳細を示すフローチャートである。この変形例に係る処理により、ＵＶスペクトルの一致度を計算する。
本サブルーチン２の変形例では、図８６のサブルーチン２でのＲＭＳＤにＵＶパターンの移動平均の傾き情報（ＤＮＳ）を加味する処理ができるようにした。ＤＮＳは、後述の式で表わされ、ＵＶパターンにおける移動平均値の移動傾きを２つのパターンで比較した時の傾き符号（ ＋／− ）の不一致数として定義される。すなわち、ＤＮＳは、ＵＶパターンの極大、極小値の位置の一致具合を評価する値である。
このＤＮＳの情報を前記ＲＭＳＤに加味することで、ＵＶスペクトルの波形の一致度をより正確に算出することができる。
図１００の変形例に係るサブルーチン２において、ステップＳ２００１〜Ｓ２００８までは、図８６のサブルーチン２とほぼ同一である。但し、ステップＳ２００１では、区間１←ｗ１、区間２←ｗ２の初期設定が追加して行われ、後述する移動平均、移動傾きの計算の区間に用いられる。
本変形例のサブルーチン２では、ＤＮＳ加味のためにステップＳ２０１０〜Ｓ２０１３を追加し、ステップＳ２００９ＡでＤＮＳを加味した一致度の計算を可能とした。
ステップＳ２０１０では、「ＤＮＳを加味する？」の判断処理が実行され、ＤＮＳを加味すると判断されたときは（ＹＥＳ）、ステップＳ２０１１へ移行し、ＤＮＳを加味しないと判断されたときは（ＮＯ）、ステップＳ２００９Ａへ移行する。ＤＮＳを加味するか否かの起因は、例えば初期の設定による。例えば、ＦＰを、ピーク面積で作成したときはＤＮＳ加味、ＦＰをピーク高さで作成したときは、ＤＮＳ非加味に設定する。
但し、ＦＰをピーク高さで作成した上記実施例の場合にも、ＤＮＳを加味する処理でＵＶパターン一致度を計算することができ、ＦＰをピーク面積で作成した場合にも、ＤＮＳを加味しない上記実施例の処理でＵＶパターン一致度を計算することができる。
ステップＳ２０１１では、「区間１（ｗ１）でｘとｙの移動平均を計算」の処理が実行され、区間１（ｗ１）における移動平均が求められる。区間１（ｗ１）は、ＵＶデータの波長に関する区間であり、ステップＳ２００１の初期設定においてｗ１＝３であれば区間１（３）となり、３個の波長におけるＵＶ強度の平均が求められる。具体的には図１０１の図表で後述する。
ステップＳ２０１２では、「区間２（ｗ２）でｘとｙの移動傾きを計算」の処理が実行され、区間２（ｗ２）における移動傾きが求められる。区間２（ｗ２）は、ステップＳ２０１１で求めた移動平均に関する区間であり、ステップＳ２００１の初期設定においてｗ２＝３であれば区間２（３）となり、ステップＳ２０１１で計算した移動平均に基づき、３個の移動平均に渡る傾きの（±）が求められる。具体的には図１０１の図表で後述する。
ステップＳ２０１３では、「ｘとｙの移動傾きの符号の不一致数を算出（ＤＮＳ）」の処理が実行され、ステップＳ２０１２で計算された移動傾きから傾きの（±）の一致数が計算される。移動傾きの＋は、図６６において右肩上がりを表し、移動傾きの−は、同右肩下がりを表す。
ステップＳ２０１３からステップＳ２００９Ａへ移行すると、このステップＳ２００９Ａの処理においてＤＮＳを加味した一致度の計算が行われる。
ステップＳ２００９Ａでは、「ｘとｙのＵＶスペクトルの一致度を計算（ＵＶ＿Ｓｉｍ）」の処理が実行され、ＤＮＳを加味した一致度の計算では、ＵＶ＿Ｓｉｍを、前記ＵＶスペクトル間距離の二乗和ｚとｘのデータ数ａとＤＮＳとから、
ＵＶ＿Ｓｉｍ ＝ √（ｚ／ａ）×１．１ＤＮＳ
として算出し、このＵＶ＿Ｓｉｍを図８１のステップＳ３０６に渡し、ＵＶスペクトルの一致度計算処理を終了する。
なお、ステップＳ２０１０からステップＳ２００９Ａへ移行した場合の処理は、図８６のステップＳ２００９と同一である。
図１０１は、移動平均及び移動傾きの計算例を示す図表である。
図１０１の上段は、ＵＶデータ例、中段は、移動平均の計算例、下段は、移動傾きの計算例を示す。ＵＶデータ例は、具体的な数値に代え、ＵＶ強度をａ１〜ａ７で表記している。例えば、２２０ｎｍのＵＶ強度がａ１、２２１ｎｍのＵＶ強度がａ２等となる。移動平均の計算例及び移動傾きの計算例も具体的数値に代え、ＵＶ強度ａ１〜ａ７を使用している。
移動平均は、区間１（ｗ１＝３）を例とし、ステップＳ２０１２（図１００）において、区間（ａ１、ａ２、ａ３）、区間（ａ２、ａ３、ａ４）、・・・毎に計算した値としてｍ１、ｍ２、・・・が計算される。移動傾きも区間２（３）を例とし、ステップＳ２０１３（図１００）において、区間（ｍ１、ｍ２、ｍ３）、区間（ｍ２、ｍ３、ｍ４）、・・・毎に計算した値としてｓ１、・・・が計算される。例えば、移動平均の差ｍ３−ｍ１が移動傾きとなり、その（±）を取り出す。
こうして、ＦＰをピーク面積で作成したときは、基準群ＦＰへの帰属処理ならびに基準ＦＰ帰属結果統合処理において、ＤＮＳを加味した処理でＵＶパターン一致度を計算することができる。この計算により、図６６で示す対応する２点の距離（ｄｉｓ）がピーク高さで作成したＦＰに比較して大きくなっても、取り扱いを容易とし、ＵＶパターン一致度を正確に計算することができる。
［その他］
本発明実施例は、多成分薬剤として漢方薬の評価について適用したが、その他の多成分物質の評価にも適用することができる。
本発明のパターンの評価方法は、時系列でピークが変化する評価対象の対象パターンを取得するパターン取得工程と、前記対象パターンとこの対象パターンに対応し評価基準となる基準パターンのピークそれぞれに対し時間軸方向前後位置の少なくとも一方に存在するｎ本のピークを含めたｎ＋１本のピークによりそれぞれのピーク・パターンを作成するピーク・パターン作成工程と、前記それぞれのピーク・パターンを比較し対応するピークを特定するピーク帰属工程と、このピーク帰属工程で特定されて帰属したピークと評価基準となる複数の基準パターンのピークとを比較評価する評価工程とを備えて、対象パターンと基準パターンとの同等性の評価等に広く適用することができる。
上記実施例のＦＰは、３Ｄクロマト上での全ピークを対象としたが、細かいデータ、例えば３Ｄクロマト上でピーク面積が５％未満のピークを除いてＦＰを作成することもできる。
上記実施例のＦＰは、ピーク高さに基づいて作成し、図７０〜図７４の評価を得たが、ピーク面積に基づいてＦＰを作成した場合についても、ピーク高さに基づいて作成した上記実施例同様の手順によりＭＴ法によりＭＤ値を求め、図７０〜図７４と同様に評価を得ることができる。
クロマトは、３Ｄクロマトに限らず、ＦＰとしてＵＶスペクトルを除いたピークとそのリテンション・タイムとで構成されたものを用いることもできる。この場合、ＵＶスペクトルの一致度を除き、上記実施例と同様に行わせることができる。
１ 多成分薬剤の評価装置（パターンの評価装置）
３ ＦＰ作成部（パターン取得部）
５ 基準ＦＰ選定部（基準パターン選定部）
７ ピーク・パターン作成部
９ ピーク帰属部
１１ 評価部
１３ 漢方薬
１５ ３Ｄクロマト
１７ 対象ＦＰ（対象パターン）
１９ 基準群ＦＰ
２１ 基準群ＦＰに帰属した対象ＦＰ
２３ 対象ＦＰの評価結果
２５ 対象ＦＰに含まれるピークのＵＶスペクトル
２７ 薬剤ＡのＦＰ
２９ 薬剤ＢのＦＰ
３１ 薬剤ＣのＦＰ
３３ 対象ＦＰ（リテンション・タイム１０．０−１４．５分）
３５、３７、３９、４１、４３、４５、４７、４９、５１、５３ 対象ＦＰ（リテンション・タイム１０．０−１４．５分）中の各ピーク
５５ 基準ＦＰ（リテンション・タイム１０．０−１４．５分）
５７、５９、６１、６３、６５、６７、６９、７１、７３、７５、７７ 基準ＦＰ（リテンション・タイム１０．０−１４．５分）中の各ピーク
７９ 対象ＦＰリテンション・タイム・出現パターン
８１ 基準ＦＰリテンション・タイム・出現パターン
８３ リテンション・タイム・出現距離の一致数
８５ リテンション・タイム・出現パターンの一致度
８７ 対象ＦＰ帰属対象ピークのピーク・パターン（３本）
８９、９１、９３、９５ 基準ＦＰ帰属候補ピークのピーク・パターン（３本）
９７ 対象ＦＰ帰属対象ピークのピーク・パターン（５本）
９９、１０１、１０３、１０５ 基準ＦＰ帰属候補ピークのピーク・パターン（５本）
１０７ 帰属対象ピークのＵＶスペクトル
１１１ 帰属候補ピークのＵＶスペクトル
１１３ ＦＰ作成工程（パターン取得工程）
１１５ 基準ＦＰ選定工程（基準パターン選定工程）
１１７ ピーク・パターン作成工程
１１９ ピーク帰属工程
１２１ 評価工程
１２３ ３Ｄクロマト・データ例
１２５ ピーク情報データ例
１２７ ＦＰデータ例
１２９ 判定結果ファイル例
１３１ 帰属候補ピークスコア表
１３３ 帰属候補ピーク番号表
１３５ 照合結果ファイル例
１３７ 基準群ＦＰデータ例
１３９ 対象ＦＰピーク特徴量ファイル例

Claims (51)

  1. 時系列でピークが変化する評価対象の対象パターンを取得するパターン取得工程と、
    前記対象パターンの帰属対象ピークと前記対象パターンに対応し評価基準となる基準パターンにおいて前記帰属対象ピークの出現時間との差が設定した範囲内の帰属候補ピークそれぞれに対し時間軸方向前後位置の少なくとも一方に存在する周辺ｎ本のピークを含めたｎ＋１本のピークによりピーク・パターンを作成するピーク・パターン作成工程と、
    前記帰属対象ピーク及び帰属候補ピークのピーク・パターンを比較し対応する前記帰属対象ピークと帰属候補ピークを特定するピーク帰属工程と、
    このピーク帰属工程で特定されて帰属したピークと評価基準となる複数の基準パターンのピークとを比較評価する評価工程と、
    を備えたことを特徴とするパターンの評価方法。
  2. 評価対象となる多成分物質のクロマトから検出されたピークとそのリテンション・タイムとで構成される対象ＦＰを作成するＦＰ作成工程と、
    前記対象ＦＰの帰属対象ピークと前記対象ＦＰに対応し評価基準となる多成分物質のクロマトから検出されたピークとそのリテンション・タイムとで構成される基準ＦＰにおいて前記帰属対象ピークのリテンション・タイムとの差が設定した範囲内の帰属候補ピークそれぞれに対し時間軸方向前後位置の少なくとも一方に存在する周辺ｎ本のピークを含めたｎ＋１本のピークによりピーク・パターンを作成するピーク・パターン作成工程と、
    前記帰属対象ピーク及び帰属候補ピークのピーク・パターンを比較し対応する前記帰属対象ピークと帰属候補ピークを特定するピーク帰属工程と、
    このピーク帰属工程で特定されて帰属したピークと評価基準となる複数の基準ＦＰのピークとを比較評価する評価工程と、
    を備えたことを特徴とする多成分物質の評価方法。
  3. 請求項２記載の多成分物質の評価方法であって、
    前記ＦＰ作成工程は、評価対象となる多成分物質の３Ｄクロマトから検出されたピークとその出現時間とそのＵＶスペクトルとで構成される対象ＦＰを作成し、
    前記ピーク・パターン作成工程は、前記対象ＦＰとこの対象ＦＰに対応し評価基準となる多成分物質の３Ｄクロマトから検出されたピークとそのリテンション・タイムとそのＵＶスペクトルとで構成される基準ＦＰのピークそれぞれに対し時間軸方向前後位置の少なくとも一方に存在するｎ本のピークを含めたｎ＋１本のピークによりそれぞれのピーク・パターンを作成する、
    ことを特徴とする多成分物質の評価方法
  4. 請求項２又は３記載の多成分物質の評価方法であって、
    前記ＦＰ作成工程は、前記評価対象のクロマト又は３Ｄクロマトから検出され検出波長の異なる複数のＦＰを合成して前記対象ＦＰを作成する、
    ことを特徴とする多成分物質の評価方法。
  5. 請求項２〜４の何れかに記載の多成分物質の評価方法であって、
    前記ＦＰ作成工程は、前記評価対象のクロマト又は３Ｄクロマトの全ピークを抽出して前記対象ＦＰを作成する、
    ことを特徴とする多成分物質の評価方法。
  6. 請求項２〜５の何れかに記載の多成分物質の評価方法であって、
    前記ピーク・パターン作成工程は、前記対象ＦＰの帰属対象となるピークならびに前記基準ＦＰの帰属候補となるピークの両方でピーク・パターンを構成するピークを変化させて網羅的にピーク・パターンを作成する、
    ことを特徴とする多成分物質の評価方法。
  7. 請求項２〜６の何れかに記載の多成分物質の評価方法であって、
    前記ピーク・パターン作成工程は、前記対象ＦＰ及び基準ＦＰの全てのピークに対し前記ピーク・パターンを作成する、
    ことを特徴とする多成分物質の評価方法。
  8. 請求項２〜７の何れかに記載の多成分物質の評価方法であって、
    前記ピーク帰属工程は、前記ピーク・パターンを構成する全てのピーク及びその出現時間から算出する前記それぞれのピーク・パターンの一致度に基づき前記特定を行う、
    ことを特徴とする多成分物質の評価方法。
  9. 請求項３記載の多成分物質の評価方法であって、
    前記ピーク帰属工程は、前記対象ＦＰ及び基準ＦＰのピークのＵＶスペクトルの比較を考慮して前記特定を行う、
    ことを特徴とする多成分物質の評価方法。
  10. 請求項９記載の多成分物質の評価方法であって、
    前記ピーク帰属工程は、前記ピーク・パターンを構成する全てのピーク及びその出現時間から算出する前記それぞれのピーク・パターンの一致度と前記ＵＶスペクトルの一致度とを総合して算出するピークの一致度に基づき前記特定を行う、
    ことを特徴とする多成分物質の評価方法。
  11. 請求項２〜１０の何れかに記載の多成分物質の評価方法であって、
    前記ピーク帰属工程は、前記対象ＦＰ及び基準ＦＰの全てのピークに対し前記特定を行う、
    ことを特徴とする多成分物質の評価方法。
  12. 請求項２〜１１の何れかに記載の多成分物質の評価方法であって、
    前記ピーク帰属工程は、前記対象ＦＰを基準ＦＰへ帰属したピークを前記複数の基準ＦＰに基づく基準群ＦＰの各ピークに対応させて前記対象ＦＰのピークを特定する、
    ことを特徴とする多成分物質の評価方法。
  13. 請求項２〜１２の何れかに記載の多成分物質の評価方法であって、
    前記評価工程は、前記基準群ＦＰに帰属された対象ＦＰのピークを特徴量としてＭＴ法によりＭＤ値を求め、評価基準となる複数のＭＤ値に対して比較評価する、
    ことを特徴とする多成分物質の評価方法。
  14. 請求項２〜１３の何れかに記載の多成分物質の評価方法であって、
    前記ピーク・パターン作成工程で用いる基準ＦＰを、前記複数の基準ＦＰから選定する基準ＦＰ選定工程を備えた、
    ことを特徴とする多成分物質の評価方法。
  15. 請求項１４記載の多成分物質の評価方法であって、
    前記基準ＦＰ選定工程は、前記複数の全基準ＦＰと前記対象ＦＰとのピークの出現時間及び出現パターンを比較し前記対象ＦＰの帰属相手として適した基準ＦＰを選定する、
    ことを特徴とする多成分物質の評価方法。
  16. 請求項２〜１５の何れかに記載の多成分物質の評価方法であって、
    前記多成分物質は、多成分薬剤である、
    ことを特徴とする多成分物質の評価方法。
  17. 請求項１６記載の多成分物質の評価方法であって、
    前記多成分薬剤は、生薬、生薬の組合せ、それらの抽出物、漢方薬の何れかである、
    ことを特徴とする多成分物質の評価方法。
  18. 時系列でピークが変化する評価対象の対象パターンを取得するパターン取得機能と、
    前記対象パターンの帰属対象ピークと前記対象パターンに対応し評価基準となる基準パターンにおいて前記帰属対象ピークの出現時間との差が設定した範囲内の帰属候補ピークそれぞれに対し時間軸方向前後位置の少なくとも一方に存在する周辺ｎ本のピークを含めたｎ＋１本のピークによりピーク・パターンを作成するピーク・パターン作成機能と、
    前記帰属対象ピーク及び帰属候補ピークのピーク・パターンを比較し対応する前記帰属対象ピークと帰属候補ピークを特定するピーク帰属機能と、
    このピーク帰属機能で特定されて帰属したピークと評価基準となる複数の基準パターンのピークとを比較評価する評価機能と、
    をコンピュータに実現させることを特徴とするパターンの評価プログラム。
  19. 評価対象となる多成分物質のクロマトから検出されたピークとそのリテンション・タイムとで構成される対象ＦＰを作成するＦＰ作成機能と、
    前記対象ＦＰの帰属対象ピークと前記対象ＦＰに対応し評価基準となる多成分物質のクロマトから検出されたピークとそのリテンション・タイムとで構成される基準ＦＰにおいて前記帰属対象ピークのリテンション・タイムとの差が設定した範囲内の帰属候補ピークそれぞれに対し時間軸方向前後位置の少なくとも一方に存在する周辺ｎ本のピークを含めたｎ＋１本のピークによりピーク・パターンを作成するピーク・パターン作成機能と、
    前記帰属対象ピーク及び帰属候補ピークのピーク・パターンを比較し対応する前記帰属対象ピークと帰属候補ピークを特定するピーク帰属機能と、
    このピーク帰属機能で特定されて帰属したピークと評価基準となる複数の基準ＦＰのピークとを比較評価する評価機能と、
    を備えたことを特徴とする多成分物質の評価プログラム。
  20. 請求項１９記載の多成分物質の評価プログラムであって、
    前記ＦＰ作成機能は、評価対象となる多成分物質の３Ｄクロマトから検出されたピークとその出現時間とそのＵＶスペクトルとで構成される対象ＦＰを作成し、
    前記ピーク・パターン作成機能は、前記対象ＦＰとこの対象ＦＰに対応し評価基準となる多成分物質の３Ｄクロマトから検出されたピークとそのリテンション・タイムとそのＵＶスペクトルとで構成される基準ＦＰのピークそれぞれに対し時間軸方向前後位置の少なくとも一方に存在するｎ本のピークを含めたｎ＋１本のピークによりそれぞれのピーク・パターンを作成する、
    ことを特徴とする多成分物質の評価プログラム
  21. 請求項１９又は２０記載の多成分物質の評価プログラムであって、
    前記ＦＰ作成機能は、前記評価対象のクロマト又は３Ｄクロマトから検出され検出波長の異なる複数のＦＰを合成して前記対象ＦＰを作成する、
    ことを特徴とする多成分物質の評価プログラム。
  22. 請求項１９〜２１の何れかに記載の多成分物質の評価プログラムであって、
    前記ＦＰ作成機能は、前記評価対象のクロマト又は３Ｄクロマトの全ピークを抽出して前記対象ＦＰを作成する、
    ことを特徴とする多成分物質の評価プログラム。
  23. 請求項１９〜２２の何れかに記載の多成分物質の評価プログラムであって、
    前記ピーク・パターン作成機能は、前記対象ＦＰの帰属対象となるピークならびに前記
    基準ＦＰの帰属候補となるピークの両方でピーク・パターンを構成するピークを変化させて網羅的にピーク・パターンを作成する、
    ことを特徴とする多成分物質の評価プログラム。
  24. 請求項１９〜２３の何れかに記載の多成分物質の評価プログラムであって、
    前記ピーク・パターン作成機能は、前記対象ＦＰ及び基準ＦＰの全てのピークに対し前記ピーク・パターンを作成する、
    ことを特徴とする多成分物質の評価プログラム。
  25. 請求項１９〜２４の何れかに記載の多成分物質の評価プログラムであって、
    前記ピーク帰属機能は、前記ピーク・パターンを構成する全てのピーク及びその出現時間から算出する前記それぞれのピーク・パターンの一致度に基づき前記特定を行う、
    ことを特徴とする多成分物質の評価プログラム。
  26. 請求項２０記載の多成分物質の評価プログラムであって、
    前記ピーク帰属機能は、前記対象ＦＰ及び基準ＦＰのピークのＵＶスペクトルの比較を考慮して前記特定を行う、
    ことを特徴とする多成分物質の評価プログラム。
  27. 請求項２６記載の多成分物質の評価プログラムであって、
    前記ピーク帰属機能は、前記ピーク・パターンを構成する全てのピーク及びその出現時間から算出する前記それぞれのピーク・パターンの一致度と前記ＵＶスペクトルの一致度とを総合して算出するピークの一致度に基づき前記特定を行う、
    ことを特徴とする多成分物質の評価プログラム。
  28. 請求項１９〜２７の何れかに記載の多成分物質の評価プログラムであって、
    前記ピーク帰属機能は、前記対象ＦＰ及び基準ＦＰの全てのピークに対し前記特定を行う、
    ことを特徴とする多成分物質の評価プログラム。
  29. 請求項１９〜２８の何れかに記載の多成分物質の評価プログラムであって、
    前記ピーク帰属機能は、前記対象ＦＰを基準ＦＰへ帰属したピークを前記複数の基準ＦＰに基づく基準群ＦＰの各ピークに対応させて前記対象ＦＰのピークを特定する、
    ことを特徴とする多成分物質の評価プログラム。
  30. 請求項１９〜２９の何れかに記載の多成分物質の評価プログラムであって、
    前記評価機能は、前記基準群ＦＰに帰属された対象ＦＰのピークを特徴量としてＭＴ法によりＭＤ値を求め、評価基準となる複数のＭＤ値に対して比較評価する、
    ことを特徴とする多成分物質の評価プログラム。
  31. 請求項１９〜３０の何れかに記載の多成分物質の評価プログラムであって、
    前記ピーク・パターン作成機能で用いる基準ＦＰを、前記複数の基準ＦＰから選定する基準ＦＰ選定機能を備えた、
    ことを特徴とする多成分物質の評価プログラム。
  32. 請求項３１記載の多成分物質の評価プログラムであって、
    前記基準ＦＰ選定機能は、前記複数の全基準ＦＰと前記対象ＦＰとのピークの出現時間及び出現パターンを比較し前記対象ＦＰの帰属相手として適した基準ＦＰを選定する、
    ことを特徴とする多成分物質の評価プログラム。
  33. 請求項１９〜３２の何れかに記載の多成分物質の評価プログラムであって、
    前記多成分物質は、多成分薬剤である、
    ことを特徴とする多成分物質の評価プログラム。
  34. 請求項３３記載の多成分物質の評価プログラムであって、
    前記多成分薬剤は、生薬、生薬の組合せ、それらの抽出物、漢方薬の何れかである、
    ことを特徴とする多成分物質の評価プログラム。
  35. 時系列でピークが変化する評価対象の対象パターンを取得するパターン取得部と、
    前記対象パターンの帰属対象ピークと前記対象パターンに対応し評価基準となる基準パターンにおいて前記帰属対象ピークの出現時間との差が設定した範囲内の帰属候補ピークそれぞれに対し時間軸方向前後位置の少なくとも一方に存在する周辺ｎ本のピークを含めたｎ＋１本のピークによりピーク・パターンを作成するピーク・パターン作成部と、
    前記帰属対象ピーク及び帰属候補ピークのピーク・パターンを比較し対応する前記帰属対象ピークと帰属候補ピークを特定するピーク帰属部と、
    このピーク帰属部で特定されて帰属したピークと評価基準となる複数の基準パターンのピークとを比較評価する評価部と、
    を備えたことを特徴とするパターンの評価装置。
  36. 評価対象となる多成分物質のクロマトから検出されたピークとそのリテンション・タイムとで構成される対象ＦＰを作成するＦＰ作成部と、
    前記対象ＦＰの帰属対象ピークと前記対象ＦＰに対応し評価基準となる多成分物質のクロマトから検出されたピークとそのリテンション・タイムとで構成される基準ＦＰにおいて前記帰属対象ピークのリテンション・タイムとの差が設定した範囲内の帰属候補ピークそれぞれに対し時間軸方向前後位置の少なくとも一方に存在する周辺ｎ本のピークを含めたｎ＋１本のピークによりピーク・パターンを作成するピーク・パターン作成部と、
    前記帰属対象ピーク及び帰属候補ピークのピーク・パターンを比較し対応する前記帰属対象ピークと帰属候補ピークを特定するピーク帰属部と、
    このピーク帰属部で特定されて帰属したピークと評価基準となる複数の基準ＦＰのピークとを比較評価する評価部と、
    を備えたことを特徴とする多成分物質の評価装置。
  37. 請求項３６記載の多成分物質の評価装置であって、
    前記ＦＰ作成部は、評価対象となる多成分物質の３Ｄクロマトから検出されたピークとその出現時間とそのＵＶスペクトルとで構成される対象ＦＰを作成し、
    前記ピーク・パターン作成部は、前記対象ＦＰとこの対象ＦＰに対応し評価基準となる多成分物質の３Ｄクロマトから検出されたピークとそのリテンション・タイムとそのＵＶスペクトルとで構成される基準ＦＰのピークそれぞれに対し時間軸方向前後位置の少なくとも一方に存在するｎ本のピークを含めたｎ＋１本のピークによりそれぞれのピーク・パターンを作成する、
    ことを特徴とする多成分物質の評価装置
  38. 請求項３６又は３７記載の多成分物質の評価装置であって、
    前記ＦＰ作成部は、前記評価対象のクロマト又は３Ｄクロマトから検出され検出波長の異なる複数のＦＰを合成して前記対象ＦＰを作成する、
    ことを特徴とする多成分物質の評価装置。
  39. 請求項３６〜３８の何れかに記載の多成分物質の評価装置であって、
    前記ＦＰ作成部は、前記評価対象のクロマト又は３Ｄクロマトの全ピークを抽出して前記対象ＦＰを作成する、
    ことを特徴とする多成分物質の評価装置。
  40. 請求項３６〜３９の何れかに記載の多成分物質の評価装置であって、
    前記ピーク・パターン作成部は、前記対象ＦＰの帰属対象となるピークならびに前記基
    準ＦＰの帰属候補となるピークの両方でピーク・パターンを構成するピークを変化させて網羅的にピーク・パターンを作成する、
    ことを特徴とする多成分物質の評価装置。
  41. 請求項３６〜４０の何れかに記載の多成分物質の評価装置であって、
    前記ピーク・パターン作成部は、前記対象ＦＰ及び基準ＦＰの全てのピークに対し前記ピーク・パターンを作成する、
    ことを特徴とする多成分物質の評価装置。
  42. 請求項３６〜４１の何れかに記載の多成分物質の評価装置であって、
    前記ピーク帰属部は、前記ピーク・パターンを構成する全てのピーク及びその出現時間から算出する前記それぞれのピーク・パターンの一致度に基づき前記特定を行う、
    ことを特徴とする多成分物質の評価装置。
  43. 請求項３７記載の多成分物質の評価装置であって、
    前記ピーク帰属部は、前記対象ＦＰ及び基準ＦＰのピークのＵＶスペクトルの比較を考慮して前記特定を行う、
    ことを特徴とする多成分物質の評価装置。
  44. 請求項４３記載の多成分物質の評価装置であって、
    前記ピーク帰属部は、前記ピーク・パターンを構成する全てのピーク及びその出現時間から算出する前記それぞれのピーク・パターンの一致度と前記ＵＶスペクトルの一致度とを総合して算出するピークの一致度に基づき前記特定を行う、
    ことを特徴とする多成分物質の評価装置。
  45. 請求項３６〜４２の何れかに記載の多成分物質の評価装置であって、
    前記ピーク帰属部は、前記対象ＦＰ及び基準ＦＰの全てのピークに対し前記特定を行う、
    ことを特徴とする多成分物質の評価装置。
  46. 請求項３６〜４５の何れかに記載の多成分物質の評価装置であって、
    前記ピーク帰属部は、前記対象ＦＰを基準ＦＰへ帰属したピークを前記複数の基準ＦＰに基づく基準群ＦＰの各ピークに対応させて前記対象ＦＰのピークを特定する、
    ことを特徴とする多成分物質の評価装置。
  47. 請求項３６〜４６の何れかに記載の多成分物質の評価装置であって、
    前記評価部は、前記基準群ＦＰに帰属された対象ＦＰのピークを特徴量としてＭＴ法によりＭＤ値を求め、評価基準となる複数のＭＤ値に対して比較評価する、
    ことを特徴とする多成分物質の評価装置。
  48. 請求項３６〜４７の何れかに記載の多成分物質の評価装置であって、
    前記ピーク・パターン作成部で用いる基準ＦＰを、前記複数の基準ＦＰから選定する基準ＦＰ選定部を備えた、
    ことを特徴とする多成分物質の評価装置。
  49. 請求項４８記載の多成分物質の評価装置であって、
    前記基準ＦＰ選定部は、前記複数の全基準ＦＰと前記対象ＦＰとのピークの出現時間及び出現パターンを比較し前記対象ＦＰの帰属相手として適した基準ＦＰを選定する、
    ことを特徴とする多成分物質の評価装置。
  50. 請求項３６〜４９の何れかに記載の多成分物質の評価装置であって、
    前記多成分物質は、多成分薬剤である、
    ことを特徴とする多成分物質の評価装置。
  51. 請求項４９記載の多成分物質の評価装置であって、
    前記多成分薬剤は、生薬、生薬の組合せ、それらの抽出物、漢方薬の何れかである、
    ことを特徴とする多成分物質の評価装置。
JP2012549177A 2011-06-01 2012-05-31 パターンの評価方法、多成分物質の評価方法、評価プログラム、及び評価装置 Active JP5915539B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2011123844 2011-06-01
JP2011123844 2011-06-01
PCT/JP2012/003602 WO2012164946A1 (ja) 2011-06-01 2012-05-31 パターンの評価方法、多成分物質の評価方法、評価プログラム、及び評価装置

Publications (2)

Publication Number Publication Date
JPWO2012164946A1 JPWO2012164946A1 (ja) 2015-02-23
JP5915539B2 true JP5915539B2 (ja) 2016-05-11

Family

ID=47258818

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2012549177A Active JP5915539B2 (ja) 2011-06-01 2012-05-31 パターンの評価方法、多成分物質の評価方法、評価プログラム、及び評価装置

Country Status (8)

Country Link
US (1) US20140149051A1 (ja)
EP (1) EP2717049B1 (ja)
JP (1) JP5915539B2 (ja)
KR (1) KR101492895B1 (ja)
CN (1) CN102959396B (ja)
HK (1) HK1181117A1 (ja)
TW (1) TWI528033B (ja)
WO (1) WO2012164946A1 (ja)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2954015C (en) 2014-07-08 2022-08-09 Hollister Incorporated Trans anal irrigation platform with bed module
US11788995B2 (en) 2016-08-10 2023-10-17 Mls Acq, Inc. System and method for background removal in spectrometry system
CN113237985B (zh) * 2021-03-24 2023-10-20 迈克医疗电子有限公司 血红蛋白中组分类型的确定方法、装置和计算机设备

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2658344B2 (ja) * 1989-01-26 1997-09-30 株式会社島津製作所 クロマトグラフ用データ処理装置
WO1997039347A1 (en) * 1996-04-12 1997-10-23 Waters Investments Limited Chromatographic pattern analysis system employing chromatographic variability characterization
JP3541577B2 (ja) * 1996-08-30 2004-07-14 株式会社島津製作所 Ｐｄａ検出器用データ処理装置
JP3094921B2 (ja) * 1996-09-26 2000-10-03 株式会社島津製作所 クロマトグラフ用データ処理装置
JP4886933B2 (ja) * 2001-01-12 2012-02-29 カウンセル オブ サイエンティフィック アンド インダストリアル リサーチ クロマトグラフフィンガープリントならびに単一の医薬および処方物の標準化のための新規な方法
US6691053B2 (en) * 2001-11-30 2004-02-10 Agilent Technologies, Inc. Method of matching retention times among multiple chromatographic system
US20090149335A1 (en) * 2002-02-22 2009-06-11 Biolife Solutions Inc. Method and use of microarray technology and proteogenomic analysis to predict efficacy of human and xenographic cell, tissue and organ transplant
US7660678B2 (en) * 2004-02-05 2010-02-09 Medpro Holdings, Llc On-site method of providing analysis of potency and purity of pharmaceutical compounds
JP4719694B2 (ja) * 2004-02-13 2011-07-06 ウオーターズ・テクノロジーズ・コーポレイシヨン 化学物質を追跡し、定量化するためのシステムおよび方法
KR20070038494A (ko) * 2004-06-07 2007-04-10 가부시키가이샤 쓰무라 다성분계 약제의 평가방법
CN1670529B (zh) * 2005-03-28 2010-04-21 中山大学 一种复方血栓通制剂hplc指纹图谱的构建方法
US20090012723A1 (en) * 2005-06-09 2009-01-08 Chemlmage Corporation Adaptive Method for Outlier Detection and Spectral Library Augmentation
JP2007315941A (ja) * 2006-05-26 2007-12-06 Univ Of Miyazaki 植物品種判定装置、植物品種判定方法及び植物品種判定用プログラム
EP2221617A1 (en) * 2009-02-20 2010-08-25 Nederlandse Organisatie voor toegepast -natuurwetenschappelijk onderzoek TNO A method, a system and a computer program for data processing for an automatic extraction of respective peaks in a chromatographic spectrum

Also Published As

Publication number Publication date
KR101492895B1 (ko) 2015-02-12
EP2717049B1 (en) 2024-04-17
KR20130037681A (ko) 2013-04-16
JPWO2012164946A1 (ja) 2015-02-23
EP2717049A4 (en) 2014-11-19
TW201303291A (zh) 2013-01-16
EP2717049A1 (en) 2014-04-09
TWI528033B (zh) 2016-04-01
US20140149051A1 (en) 2014-05-29
CN102959396A (zh) 2013-03-06
WO2012164946A1 (ja) 2012-12-06
CN102959396B (zh) 2015-12-02
HK1181117A1 (zh) 2013-11-01

Similar Documents

Publication Publication Date Title
JP5915540B2 (ja) ピーク帰属方法、帰属プログラム、及び帰属装置
JP5915539B2 (ja) パターンの評価方法、多成分物質の評価方法、評価プログラム、及び評価装置
JP5912880B2 (ja) パターン又はｆｐの特徴量作成方法、作成プログラム、及び作成装置
JP5954180B2 (ja) パターンの評価方法、多成分物質の評価方法、評価プログラム、及び評価装置
JP5910506B2 (ja) 集合データの類似性評価方法、類似性評価プログラム、及び類似性評価装置
JP5895847B2 (ja) Ｆｐ作成方法、作成プログラム、作成装置、及びｆｐ

Legal Events

Date Code Title Description
A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20150818

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20151015

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20160308

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20160321

R150 Certificate of patent or registration of utility model

Ref document number: 5915539

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250