KR101412655B1 - Fp작성방법, 작성프로그램, 작성장치, 및 fp - Google Patents
Fp작성방법, 작성프로그램, 작성장치, 및 fp Download PDFInfo
- Publication number
- KR101412655B1 KR101412655B1 KR1020127032269A KR20127032269A KR101412655B1 KR 101412655 B1 KR101412655 B1 KR 101412655B1 KR 1020127032269 A KR1020127032269 A KR 1020127032269A KR 20127032269 A KR20127032269 A KR 20127032269A KR 101412655 B1 KR101412655 B1 KR 101412655B1
- Authority
- KR
- South Korea
- Prior art keywords
- peak
- retention time
- data
- peaks
- attribution
- Prior art date
Links
Images
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
- 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/04—Preparation or injection of sample to be analysed
-
- 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/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
- 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/26—Conditioning of the fluid carrier; Flow patterns
- G01N30/38—Flow patterns
- G01N30/46—Flow patterns using more than one column
- G01N30/461—Flow patterns using more than one column with serial coupling of separation columns
- G01N30/463—Flow patterns using more than one column with serial coupling of separation columns for multidimensional chromatography
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
Landscapes
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Pathology (AREA)
- Immunology (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- Engineering & Computer Science (AREA)
- Library & Information Science (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Theoretical Computer Science (AREA)
- Computing Systems (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Crystallography & Structural Chemistry (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)
- Investigating Or Analysing Biological Materials (AREA)
Abstract
평가대상의 평가의 정도 및 효율을 향상시키는 것에 기여가능하게 한다.
평가 대상이 되는 다성분 약제의 3D크로마토로부터 특정의 파장에서 검출된 피크의 피크와 그 피크의 리텐션타임, 그 피크의 UV스펙트럼으로 구성되는 대상FP를 작성하는 FP작성부(3)를 구비하여, 이 FP작성부(3)를 포함하는 다성분 약제의 평가장치는 대상FP의 피크귀속에 적합한 기준FP를 복수의 기준FP로부터 선정하는 기준FP선정부(5)와, 대상FP로 선정된 기준FP의 피크 각각에 대하여, 예를 들면 시간축방향 전후의 적어도 일방에 존재하는 2개의 피크를 포함한 3개의 피크로 구성되는 피크패턴을 작성하는 피크패턴작성부(7)와, 피크패턴 및 피크의 UV스펙트럼을 비교하여 대응하는 피크를 특정하는 피크귀속부(9)와, 귀속한 피크와 복수의 기준FP의 피크를 예를 들면 MT법에 의해 비교평가하는 평가부(11)를 구비한 것을 특징으로 한다.
평가 대상이 되는 다성분 약제의 3D크로마토로부터 특정의 파장에서 검출된 피크의 피크와 그 피크의 리텐션타임, 그 피크의 UV스펙트럼으로 구성되는 대상FP를 작성하는 FP작성부(3)를 구비하여, 이 FP작성부(3)를 포함하는 다성분 약제의 평가장치는 대상FP의 피크귀속에 적합한 기준FP를 복수의 기준FP로부터 선정하는 기준FP선정부(5)와, 대상FP로 선정된 기준FP의 피크 각각에 대하여, 예를 들면 시간축방향 전후의 적어도 일방에 존재하는 2개의 피크를 포함한 3개의 피크로 구성되는 피크패턴을 작성하는 피크패턴작성부(7)와, 피크패턴 및 피크의 UV스펙트럼을 비교하여 대응하는 피크를 특정하는 피크귀속부(9)와, 귀속한 피크와 복수의 기준FP의 피크를 예를 들면 MT법에 의해 비교평가하는 평가부(11)를 구비한 것을 특징으로 한다.
Description
본 발명은 평가대상, 예를 들면 다성분 약제(藥劑)인 한방약의 품질평가를 행하기 위해 사용하는 FP의 FP작성방법, 작성프로그램, 작성장치 및 FP에 관한 것이다.
평가대상으로서, 예를 들면, 다성분으로 구성되는 약제(이하, 다성분 약제)인 한방약 등의 천연물 유래(由來)의 약제가 있다. 이들 약제의 정량적, 정성적 프로필은 사용하는 원료 생약의 지질학적 요인, 생태학적 요인, 수집시기, 수집장소, 수집연대, 생육기의 기후 등이 원인으로 변화한다.
그 때문에, 이들 다성분 약제 등에 대해서는, 그 안전성 및 유효성을 담보하기 위한 품질로서 일정의 기준을 규정하고, 그 기준에 기초하여, 국가의 감독기관, 화학적 조직, 제조업자 등이 품질평가를 행하고 있다.
그러나, 다성분 약제의 품질 등의 판정기준은 다성분 약제 중에 있는 특징적인 1 내지 수 성분을 선택하여, 그 함유량 등에 기초하여 설정하는 것이 일반적이었다.
예를 들면, 비특허문헌 1에서는 다성분 약제에서 유효성분의 동정(同定)이 되어 있지 않은 경우, 정량분석이 가능하고, 물에 용해하기 쉬우며, 열수(熱水) 중에서 분해하지 않으며, 다른 성분과 화학반응을 하지 않는 등의 물성을 가진 복수의 성분을 선택하고, 화학분석에 의해 얻어지는 그들 성분의 함유량을 평가의 기준으로 하고 있다.
또한, 다성분 약제에 크로마토그래피를 적용하여, 보지(保持)시간마다 자외가시(紫外可視) 흡수 스펙트럼을 얻어서, 그 중의 일부의 성분 정보로부터 평가의 기준을 설정하는 것도 알려져 있다.
예를 들면, 특허문헌 1에서는 HPLC 크로마토그램 데이타(이하, 크로마토) 중의 일부의 피크를 선택하여 바코드화하는 것에 의해 다성분 약제를 평가하고 있다.
그러나, 이들의 방법은 평가의 대상이 '특정 성분의 함유량' 또는 '특정 성분의 크로마토 피크'에 한정되어 있어서, 다성분 약제가 함유하는 성분의 일부를 평가대상으로 하고 있는 것에 지나지 않는다. 이 때문에, 다성분 약제에 대해서는 평가의 대상 외로 되어 있는 성분이 다수 존재하기 때문에, 다성분 약제의 평가방법으로서 정도(精度)적으로 불충분하다.
다성분 약제의 품질을 정확히 평가하기 위해서는, 전(全) 피크정보 혹은 수 %의 미세한 정보를 제외한 전 피크정보에 가까운 피크정보를 망라한 패턴의 평가가 필요하고, 이를 위해, 다성분 약제 사이에 전 피크 혹은 이에 가까운 개수의 피크를 대응시킬 필요가 있다.
그런데, 복수의 피크를 고정도로 효율적으로 대응시키는 것이 곤란하며, 이 것이 다성분 약제의 고정도의 효율적인 평가에 방해가 되고 있다.
보다 상세히 설명하면, 동일 약품명의 다성분 약제라고 해도, 원료가 되는 생약이 천연물이기 때문에 구성성분이 약간 상이한 경우가 있다. 이 때문에, 같은 품질의 약제라고 해도 구성 성분의 함량비율이 상이하다던가, 어느 약제에 존재하는 성분이 다른 약제에는 존재하지 않는 것(이하, 약제간 오차)이 있다. 더욱이, 크로마토에 있어서의 피크강도나 피크의 용출시간에는 엄밀한 재현성이 없는(이하, 분석오차) 등의 요인도 있다. 이들에 의해, 다성분 약제의 전 피크 혹은 이들에 가까운 개수의 피크에 대해서 동일 성분에 유래하는 피크에 대응시키는 것(이하, 피크귀속)이 가능하지 않은 것이 원인이 되어 고정도의 효율적인 평가에 방해가 되고 있다.
월간약사(月刊藥事) vol. 28, No. 3, 67 - 71 (1986)
해결하려는 문제점은 기존의 평가방법에서는 평가대상의 품질을 고정도로 효율적으로 평가하는 것에 한계가 있었다는 점이다.
본 발명은 평가의 정도 및 효율을 향상시키는 것에 기여하기 위해, 평가대상의 크로마토로부터 검출된 피크와 그 리텐션타임으로 구성하는 FP를 작성하는 FP작성공정을 구비한 FP작성방법에 있어서, 상기 FP작성공정은 상기 크로마토를 리텐션타임, 검출파장, 및 피크를 데이타로 하는 3D크로마토로 하고, 이 3D크로마토로부터 특정파장에서 검출된 피크와 그 리텐션타임 및 그 피크의 UV스펙트럼으로 상기 FP를 작성하는 것을 FP작성방법의 특징으로 한다.
또한, 평가대상의 3D크로마토로부터 검출된 피크와 그 리텐션타임으로 구성하는 FP를 작성하는 FP작성기능을 컴퓨터에 실현시키는 평가대상의 FP작성프로그램에 있어서, 상기 FP작성기능은 상기 크로마토를 리텐션타임, 검출파장, 및 피크를 데이타로 하는 3D크로마토로 하고, 이 3D크로마토로부터 특정파장에서 검출된 피크와 그 리텐션타임 및 그 피크의 UV스펙트럼으로 상기 FP를 구성하는 것을 FP작성프로그램의 특징으로 한다.
또한, 평가대상의 3D크로마토로부터 검출된 피크와 그 리텐션타임으로 구성하는 FP를 작성하는 FP작성부를 구비한 평가대상의 FP작성장치에 있어서, 상기 FP작성부는 상기 크로마토를 리텐션타임, 검출파장, 및 피크를 데이타로 하는 3D크로마토로 하고, 이 3D크로마토로부터 특정파장에서 검출된 피크와 그 리텐션타임 및 그 피크의 UV스펙트럼으로 구성한 것을 FP를 구성하는 것을 FP작성장치의 특징으로 한다.
또한, 상기 FP작성방법에 따라서 작성된 FP에 있어서, 상기 3D크로마토로부터 특정파장에서 검출된 피크와 그 리텐션타임 및 그 피크의 UV스펙트럼으로 구성한 것을 FP의 특징으로 한다.
본 발명의 평가대상의 FP작성방법은 상기 구성이기 때문에, 평가대상의 삼차원 크로마토그램 데이타(이하, 3D크로마토)로부터 특정의 파장에서 검출된 피크에 있어서의 시그널강도(높이)의 극대치 또는 면적치(이하, 피크)와 그 피크의 출현시간(이하, 리텐션타임), 및 그 피크의 자외가시흡수스펙트럼(이하, UV스텍트럼)으로 구성하는 데이타(이하, 핑거프린트데이타 : FP)를 작성할 수 있다.
이 FP는 3D크로마토와 마찬가지로 삼차원의 정보(피크, 리텐션타임 및 UV스펙트럼)로 구성되어 있다.
이 때문에, FP는 그 약제 특유의 정보를 그대로 계승한 데이타이다. 그럼에도 불구하고, 데이타용량은 대폭으로 압축되어 있기 때문에, 3D크로마토에 비교하여, 처리해야할 정보량을 대폭으로 감소시킬 수 있어서 처리속도를 빠르게 할 수 있다.
결과적으로 대상FP의 각 피크를 기준FP의 각 피크에 대하여 정도 좋고, 효율 좋게 귀속시키는 것이 가능하게 되어, 평가의 정도 및 효율을 보다 향상시키는 것에 기여할 수 있다.
본 발명의 FP작성프로그램은 상기 구성이기 때문에, 각 기능을 컴퓨터에 실현시켜서, 삼차원정보의 FP를 작성하여, 평가의 정도 및 효율을 보다 향상시키는 것에 기여할 수 있다.
본 발명의 FP작성장치는 상기 구성이기 때문에, 각 부를 작용시켜서, 삼차원정보의 FP를 작성하여, 평가의 정도 및 효율을 보다 향상시키는 것에 기여할 수 있다.
본 발명의 FP는 상기 구성이기 때문에, 평가의 정도 및 효율을 보다 향상시키는 것에 기여할 수 있다.
도 1은 다성분 약제의 평가장치의 블록도이다. (실시예 1)
도 2는 다성분 약제의 평가수순을 도시하는 블록도이다. (실시예 1)
도 3은 삼차원 크로마토그램 데이타(이하, 3D크로마토)로부터 작성한 FP의 설명도이다. (실시예 1)
도 4는 약제마다의 FP를 도시하며, (A)는 약제A, (B)는 약제B, (C)는 약제C를 도시하는 그래프이다. (실시예 1)
도 5는 대상FP 및 기준FP의 리텐션타임을 도시하는 도면이다. (실시예 1)
도 6은 대상FP의 리텐션타임출현패턴을 도시하는 도면이다. (실시예 1)
도 7은 기준FP의 리텐션타임출현패턴을 도시하는 도면이다. (실시예 1)
도 8은 대상FP와 기준FP의 리텐션타임출현거리의 일치수를 도시하는 도면이다. (실시예 1)
도 9는 대상FP와 기준FP의 리텐션타임출현패턴의 일치도를 도시하는 도면이다. (실시예 1)
도 10은 대상FP의 귀속대상피크를 도시하는 도면이다. (실시예 1)
도 11은 귀속대상피크를 포함한 피크 3개에 따른 피크패턴도이다. (실시예 1)
도 12는 귀속대상피크를 포함한 피크 5개에 따른 피크패턴도이다. (실시예 1)
도 13은 귀속대상피크의 허용폭을 도시하는 도면이다. (실시예 1)
도 14는 귀속대상피크에 대한 기준FP의 귀속후보피크를 도시하는 도면이다. (실시예 1)
도 15는 귀속대상피크와 귀속후보피크의 피크 3개에 따른 피크패턴도이다. (실시예 1)
도 16은 귀속대상피크와 별도의 귀속후보피크의 피크 3개에 따른 피크패턴도이다. (실시예 1)
도 17은 귀속대상피크와 별도의 귀속후보피크의 피크 3개에 따른 피크패턴도이다. (실시예 1)
도 18은 귀속대상피크와 별도의 귀속후보피크의 피크 3개에 따른 피크패턴도이다. (실시예 1)
도 19는 귀속대상피크와 귀속후보피크의 피크 5개에 따른 피크패턴도이다. (실시예 1)
도 20은 귀속대상피크와 별도의 귀속후보피크의 피크 5개에 따른 피크패턴도이다. (실시예 1)
도 21은 귀속대상피크와 별도의 귀속후보피크의 피크 5개에 따른 피크패턴도이다. (실시예 1)
도 22는 귀속대상피크와 별도의 귀속후보피크의 피크 5개에 따른 피크패턴도이다. (실시예 1)
도 23은 귀속대상피크와 귀속후보피크의 피크패턴구성후보피크를 도시하는 도면이다. (실시예 1)
도 24는 피크패턴구성후보피크를 4개로 한 경우의 귀속대상피크의 전 피크패턴수를 도시하는 도면이다. (실시예 1)
도 25는 피크패턴구성후보피크를 4개로 한 경우의 귀속후보피크의 전 피크패턴수를 도시하는 도면이다. (실시예 1)
도 26은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 27은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 28은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 29는 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 30은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 31은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 32는 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 33은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 34는 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 35는 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 36은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 37은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 38은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 39는 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 40은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 41은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 42는 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 43은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 44는 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 45는 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 46은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 47은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 48은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 49는 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 50은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 51은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 52는 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 53은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 54는 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 55는 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 56은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 57은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 58은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 59는 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 60은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 61은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 62는 귀속대상피크와 귀속후보피크의 피크 3개에 따른 피크패턴의 일치도의 산출방법을 도시하는 도면이다. (실시예 1)
도 63은 귀속대상피크와 귀속후보피크의 피크 3개에 따른 피크패턴의 일치도의 산출방법을 도시하는 도면이다. (실시예 1)
도 64는 귀속대상피크와 귀속후보피크의 피크 5개에 따른 피크패턴의 일치도의 산출방법을 도시하는 도면이다. (실시예 1)
도 65는 귀속대상피크와 귀속후보피크의 UV스펙트럼을 도시하는 도면이다. (실시예 1)
도 66은 귀속대상피크와 귀속후보피크의 UV스펙트럼의 일치도의 설명도이다. (실시예 1)
도 67은 피크패턴과 UV스펙트럼의 양방의 비교에 따른 귀속후보피크의 일치도의 설명도이다. (실시예 1)
도 68은 대상FP의 기준군FP로의 귀속을 도시하는 설명도이다. (실시예 1)
도 69는 대상FP가 기준군FP로 귀속된 상황을 도시하는 도면이다. (실시예 1)
도 70은 각종 대상FP와 그 평가치(MD값)를 도시하는 도면이다. (실시예 1)
도 71은 각종 대상FP와 그 평가치(MD값)를 도시하는 도면이다. (실시예 1)
도 72는 각종 대상FP와 그 평가치(MD값)를 도시하는 도면이다. (실시예 1)
도 73은 각종 대상FP와 그 평가치(MD값)을 도시하는 도면이다. (실시예 1)
도 74는 각종 대상FP와 그 평가치(MD값)를 도시하는 도면이다. (실시예 1)
도 75는 다성분 약제의 평가방법을 도시하는 공정도이다. (실시예 1)
도 76은 다성분 약제의 평가플로우도이다. (실시예 1)
도 77은 단파장에 따른 FP작성기능에 있어서의 데이타처리 플로우차트이다. (실시예 1)
도 78은 복수 파장에 따른 FP작성기능에 있어서의 데이타처리 플로우차트이다. (실시예 1)
도 79는 복수 파장에 따른 FP작성기능에 있어서의 데이타처리 플로우차트이다. (실시예 1)
도 80은 피크귀속처리1(기준FP의 선정)에 있어서의 데이타처리 플로우차트이다. (실시예 1)
도 81은 피크귀속처리2(귀속스코어의 산출)에 있어서의 데이타처리 플로우차트이다. (실시예 1)
도 82는 피크귀속처리3(대응피크의 특정)에 있어서의 데이타처리 플로우차트이다. (실시예 1)
도 83은 피크귀속처리4(기준군FP로의 귀속)에 있어서의 데이타처리 플로우차트이다. (실시예 1)
도 84는 피크귀속처리4(기준군FP로의 귀속)에 있어서의 데이타처리 플로우차트이다. (실시예 1)
도 85는 피크귀속처리1(기준FP의 선정)에 있어서의 리텐션타임출현패턴의 일치도계산처리의 플로우차트이다. (실시예 1)
도 86은 피크귀속처리2(귀속스코어의 산출)에 있어서의 UV스펙트럼의 일치도계산처리의 플로우차트이다. (실시예 1)
도 87은 피크귀속처리2(귀속스코어의 산출)에 있어서의 피크패턴의 일치도계산처리의 플로우차트이다. (실시예 1)
도 88은 기준FP특징량파일을 작성하기 위한 플로우차트이다. (실시예 1)
도 89는 '기준FP귀속결과통합처리(기준FP대응표의 작성)'의 상세를 도시하는 플로우차트이다. (실시예 1)
도 90는 '기준FP귀속결과통합처리(기준FP대응표의 작성)'의 상세를 도시하는 플로우차트이다. (실시예 1)
도 91은 '피크특징량화처리(기준군FP의 작성)'의 상세를 도시하는 플로우차트이다. (실시예 1)
도 92는 3D크로마토의 데이타 예를 도시하는 도표이다. (실시예 1)
도 93은 피크정보의 데이타 예를 도시하는 도표이다. (실시예 1)
도 94는 FP의 데이타 예를 도시하는 도표이다. (실시예 1)
도 95는 대상FP의 기준FP로의 귀속스코어계산결과(판정결과) 파일 예를 도시하는 도표이다. (실시예 1)
도 96은 대상FP와 기준FP에서 대응하는 피크의 조합(照合)과정을 도시하는 도표이다. (실시예 1)
도 97은 조합결과 파일 예를 도시하는 도표이다. (실시예 1)
도 98은 기준군FP의 데이타 예를 도시하는 도표이다. (실시예 1)
도 99는 대상FP피크특징량파일 예를 도시하는 도표이다. (실시예 1)
도 100은 도 86 대신에 적용하는 섭루틴2의 변형례의 상세를 도시하는 플로우차트이다. (실시예 1)
도 101은 이동평균 및 이동경사의 계산 예를 도시하는 도표이다. (실시예 1)
도 2는 다성분 약제의 평가수순을 도시하는 블록도이다. (실시예 1)
도 3은 삼차원 크로마토그램 데이타(이하, 3D크로마토)로부터 작성한 FP의 설명도이다. (실시예 1)
도 4는 약제마다의 FP를 도시하며, (A)는 약제A, (B)는 약제B, (C)는 약제C를 도시하는 그래프이다. (실시예 1)
도 5는 대상FP 및 기준FP의 리텐션타임을 도시하는 도면이다. (실시예 1)
도 6은 대상FP의 리텐션타임출현패턴을 도시하는 도면이다. (실시예 1)
도 7은 기준FP의 리텐션타임출현패턴을 도시하는 도면이다. (실시예 1)
도 8은 대상FP와 기준FP의 리텐션타임출현거리의 일치수를 도시하는 도면이다. (실시예 1)
도 9는 대상FP와 기준FP의 리텐션타임출현패턴의 일치도를 도시하는 도면이다. (실시예 1)
도 10은 대상FP의 귀속대상피크를 도시하는 도면이다. (실시예 1)
도 11은 귀속대상피크를 포함한 피크 3개에 따른 피크패턴도이다. (실시예 1)
도 12는 귀속대상피크를 포함한 피크 5개에 따른 피크패턴도이다. (실시예 1)
도 13은 귀속대상피크의 허용폭을 도시하는 도면이다. (실시예 1)
도 14는 귀속대상피크에 대한 기준FP의 귀속후보피크를 도시하는 도면이다. (실시예 1)
도 15는 귀속대상피크와 귀속후보피크의 피크 3개에 따른 피크패턴도이다. (실시예 1)
도 16은 귀속대상피크와 별도의 귀속후보피크의 피크 3개에 따른 피크패턴도이다. (실시예 1)
도 17은 귀속대상피크와 별도의 귀속후보피크의 피크 3개에 따른 피크패턴도이다. (실시예 1)
도 18은 귀속대상피크와 별도의 귀속후보피크의 피크 3개에 따른 피크패턴도이다. (실시예 1)
도 19는 귀속대상피크와 귀속후보피크의 피크 5개에 따른 피크패턴도이다. (실시예 1)
도 20은 귀속대상피크와 별도의 귀속후보피크의 피크 5개에 따른 피크패턴도이다. (실시예 1)
도 21은 귀속대상피크와 별도의 귀속후보피크의 피크 5개에 따른 피크패턴도이다. (실시예 1)
도 22는 귀속대상피크와 별도의 귀속후보피크의 피크 5개에 따른 피크패턴도이다. (실시예 1)
도 23은 귀속대상피크와 귀속후보피크의 피크패턴구성후보피크를 도시하는 도면이다. (실시예 1)
도 24는 피크패턴구성후보피크를 4개로 한 경우의 귀속대상피크의 전 피크패턴수를 도시하는 도면이다. (실시예 1)
도 25는 피크패턴구성후보피크를 4개로 한 경우의 귀속후보피크의 전 피크패턴수를 도시하는 도면이다. (실시예 1)
도 26은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 27은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 28은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 29는 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 30은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 31은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 32는 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 33은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 34는 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 35는 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 36은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 37은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 38은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 39는 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 40은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 41은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 42는 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 43은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 44는 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 45는 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 46은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 47은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 48은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 49는 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 50은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 51은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 52는 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 53은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 54는 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 55는 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 56은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 57은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 58은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 59는 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 60은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 61은 귀속대상피크의 피크패턴에 대한 귀속후보피크의 피크패턴의 망라적 비교의 설명도이다. (실시예 1)
도 62는 귀속대상피크와 귀속후보피크의 피크 3개에 따른 피크패턴의 일치도의 산출방법을 도시하는 도면이다. (실시예 1)
도 63은 귀속대상피크와 귀속후보피크의 피크 3개에 따른 피크패턴의 일치도의 산출방법을 도시하는 도면이다. (실시예 1)
도 64는 귀속대상피크와 귀속후보피크의 피크 5개에 따른 피크패턴의 일치도의 산출방법을 도시하는 도면이다. (실시예 1)
도 65는 귀속대상피크와 귀속후보피크의 UV스펙트럼을 도시하는 도면이다. (실시예 1)
도 66은 귀속대상피크와 귀속후보피크의 UV스펙트럼의 일치도의 설명도이다. (실시예 1)
도 67은 피크패턴과 UV스펙트럼의 양방의 비교에 따른 귀속후보피크의 일치도의 설명도이다. (실시예 1)
도 68은 대상FP의 기준군FP로의 귀속을 도시하는 설명도이다. (실시예 1)
도 69는 대상FP가 기준군FP로 귀속된 상황을 도시하는 도면이다. (실시예 1)
도 70은 각종 대상FP와 그 평가치(MD값)를 도시하는 도면이다. (실시예 1)
도 71은 각종 대상FP와 그 평가치(MD값)를 도시하는 도면이다. (실시예 1)
도 72는 각종 대상FP와 그 평가치(MD값)를 도시하는 도면이다. (실시예 1)
도 73은 각종 대상FP와 그 평가치(MD값)을 도시하는 도면이다. (실시예 1)
도 74는 각종 대상FP와 그 평가치(MD값)를 도시하는 도면이다. (실시예 1)
도 75는 다성분 약제의 평가방법을 도시하는 공정도이다. (실시예 1)
도 76은 다성분 약제의 평가플로우도이다. (실시예 1)
도 77은 단파장에 따른 FP작성기능에 있어서의 데이타처리 플로우차트이다. (실시예 1)
도 78은 복수 파장에 따른 FP작성기능에 있어서의 데이타처리 플로우차트이다. (실시예 1)
도 79는 복수 파장에 따른 FP작성기능에 있어서의 데이타처리 플로우차트이다. (실시예 1)
도 80은 피크귀속처리1(기준FP의 선정)에 있어서의 데이타처리 플로우차트이다. (실시예 1)
도 81은 피크귀속처리2(귀속스코어의 산출)에 있어서의 데이타처리 플로우차트이다. (실시예 1)
도 82는 피크귀속처리3(대응피크의 특정)에 있어서의 데이타처리 플로우차트이다. (실시예 1)
도 83은 피크귀속처리4(기준군FP로의 귀속)에 있어서의 데이타처리 플로우차트이다. (실시예 1)
도 84는 피크귀속처리4(기준군FP로의 귀속)에 있어서의 데이타처리 플로우차트이다. (실시예 1)
도 85는 피크귀속처리1(기준FP의 선정)에 있어서의 리텐션타임출현패턴의 일치도계산처리의 플로우차트이다. (실시예 1)
도 86은 피크귀속처리2(귀속스코어의 산출)에 있어서의 UV스펙트럼의 일치도계산처리의 플로우차트이다. (실시예 1)
도 87은 피크귀속처리2(귀속스코어의 산출)에 있어서의 피크패턴의 일치도계산처리의 플로우차트이다. (실시예 1)
도 88은 기준FP특징량파일을 작성하기 위한 플로우차트이다. (실시예 1)
도 89는 '기준FP귀속결과통합처리(기준FP대응표의 작성)'의 상세를 도시하는 플로우차트이다. (실시예 1)
도 90는 '기준FP귀속결과통합처리(기준FP대응표의 작성)'의 상세를 도시하는 플로우차트이다. (실시예 1)
도 91은 '피크특징량화처리(기준군FP의 작성)'의 상세를 도시하는 플로우차트이다. (실시예 1)
도 92는 3D크로마토의 데이타 예를 도시하는 도표이다. (실시예 1)
도 93은 피크정보의 데이타 예를 도시하는 도표이다. (실시예 1)
도 94는 FP의 데이타 예를 도시하는 도표이다. (실시예 1)
도 95는 대상FP의 기준FP로의 귀속스코어계산결과(판정결과) 파일 예를 도시하는 도표이다. (실시예 1)
도 96은 대상FP와 기준FP에서 대응하는 피크의 조합(照合)과정을 도시하는 도표이다. (실시예 1)
도 97은 조합결과 파일 예를 도시하는 도표이다. (실시예 1)
도 98은 기준군FP의 데이타 예를 도시하는 도표이다. (실시예 1)
도 99는 대상FP피크특징량파일 예를 도시하는 도표이다. (실시예 1)
도 100은 도 86 대신에 적용하는 섭루틴2의 변형례의 상세를 도시하는 플로우차트이다. (실시예 1)
도 101은 이동평균 및 이동경사의 계산 예를 도시하는 도표이다. (실시예 1)
평가의 정도 및 효율을 향상시키는 것에 기여한다는 목적은, 삼차원의 정보(피크, 리텐션타임 및 UV스펙트럼)로 작성한 FP에 의해 실현했다.
실시예 1
본 발명의 실시예 1은 다성분 물질 등의 평가대상, 예를 들면 다성분 약제의 FP작성방법, 작성프로그램, 작성장치 및 FP이다.
다성분 약제는 복수의 유효화학성분을 함유하는 약제로 정의되며, 한정되지 않지만, 생약, 생약의 조합, 그들의 추출물, 한방약 등이 포함된다. 또한, 제형(劑形)도 특히 한정되지 않으며, 예를 들어, 제 15 개정 일본약국방(藥局方)의 제제(製劑)총칙에서 한정되어 있는 액제(液劑), 엑기스제, 캅셀제, 과립제, 환제, 현탁제·유제, 산제(散劑), 주정제(酒精劑), 정제(錠劑), 침제(浸劑)·전제(煎劑), 틴크제, 트로치제, 방향수제, 유(流)엑기스제 등이 포함된다. 다성분 물질로서는 약제 이외의 것도 포함된다.
한방제의 구체 예는 의료용 한방제제(製劑) 148 처방 '사용상의 주의'의 업계통일과 자주개정(自主改訂), 일반용 한방처방의 안내(1978년)에 기재되어 있다.
다성분 약제의 평가에는 평가대상 약제가 정상품으로 정해진 복수의 약제와 동등한가 어떤가를 평가하기 위해, 먼저 평가대상 약제의 3차원 크로마토그램 데이타(이하, 3D크로마토)로부터 그 약제 특유의 정보를 추출한 대상FP를 작성한다.
이어서, 대상FP의 각 피크를 전(全) 기준FP를 피크귀속처리하여 작성한 전 기준FP의 피크 대응 데이타(이하, 기준군FP)에 귀속하여, 피크 특정량을 얻는다.
이어서, 기준군FP의 피크와 귀속한 대상FP의 피크(이하, 대상FP귀속피크)의 동등성을 MT법으로 평가한다. 마지막으로, 얻어진 평가치(이하, MD값)와 미리 설정해 둔 판정치(MD값의 상한치)를 비교하여, 평가대상 약제가 정상품과 동등한가 어떤가를 판정한다.
또한, 3D크로마토는 평가대상이 되는 다성분 물질인 다성분 약제의 HPLC 크로노마토그램 데이타(이하, 크로마토)이며, UV스펙트럼을 포함한다.
FP란 특정의 파장에서 검출된 피크에 있어서의 시그널강도(높이)의 극대치 또는 면적치(이하, 피크)와 그 피크의 출현시간(이하, 리텐션타임)으로 구성하는 핑거프린트 데이타이다.
대상FP는 평가대상인 한방약의 삼차원 크로마토그램 데이타인 3D크로마토로부터 특정의 검출파장에 있어서의 복수의 피크와 그 리텐션타임 및 UV스펙트럼을 추출한 것이다.
기준FP는 대상FP에 대응하며, 정상품으로 정한 다성분 물질인 다성분 약제로서의 한방약의 FP이다.
[다성분 약제의 평가장치]
도 1은 다성분 약제의 평가장치의 블록도, 도 2는 다성분 약제의 평가수순을 도시하는 블록도, 도 3은 3D크로마토로부터 작성한 FP의 설명도, 도 4 (A)는 약제A, (B)는 약제B, (C)는 약제C의 FP이다.
도 1, 도 2와 같이, 패턴의 평가장치로서의 다성분 약제의 평가장치(1)는 패턴취득부로서의 FP작성부(3)와, 기준FP선정부(5)와, 피크패턴작성부(7)와, 피크귀속부(9)와, 평가부(11)를 구비하고 있다. 이 다성분 약제의 평가장치(1)는 컴퓨터로 구성되며, 도시하지는 않았지만, CPU, ROM, RAM 등을 구비하고 있다.
본 실시예에서, FP작성부(3)와, 기준FP선정부(5)와, 피크패턴작성부(7)와, 피크귀속부(9)와, 평가부(11)는 하나의 컴퓨터로 구성되어 있다. 단, FP작성부(3)와, 기준FP선정부(5)와, 피크패턴작성부(7)와, 피크귀속부(9)와, 평가부(11)를 각각 별도의 컴퓨터로 구성하거나, 혹은 FP작성부(3) 및 기준FP선정부(5)와, 피크패턴작성부(7) 및 피크귀속부(9)와, 평가부(11)를 별도의 컴퓨터로 구성하는 것 등도 가능하다.
FP작성부(3)는 다성분 약제의 평가장치(1)의 일부로서 구성된 FP작성장치에 구비되어, 컴퓨터에 인스톨된 FP작성프로그램에 의해, 예를 들면, 한방약(13, 도 2 참조)의 크로마토로서 도 3과 같은 삼차원 크로마토그램 데이타인 3D크로마토(15)로부터 특정의 검출파장에 있어서의 복수의 피크와 그 리텐션타임 및 UV스펙트럼을 추출한 대상FP(17)(이하, 단지 'FP(17')로 표기하는 경우도 있음)를 작성하는 기능부이다. 또한, FP작성프로그램은 이것을 기록한 FP작성프로그램 기록매체를 이용하여, 컴퓨터로 구성된 FP작성부(3)에 이것을 판독시키는 것으로 FP작성을 실현시킬 수도 있다.
이 FP(17)는 3D크로마토(15)와 마찬가지로 삼차원의 정보(피크, 리텐션타임 및 UV스펙트럼)로 구성되어 있다.
이 때문에, FP(17)는 그 약제 특유의 정보를 그대로 계승한 데이타이다. 그럼에도 불구하고, 데이타 용량은 약 1/70으로 압축되어 있기 때문에, 3D크로마토(15)에 비교하여, 처리해야할 정보량을 대폭으로 감소시킬 수 있어서 처리속도를 빠르게 할 수 있다.
3D크로마토(15)는 한방약(13)(도 2 참조)에 고속액체 크로마토그래피(HPLC)를 적용한 결과이다. 이 3D크로마토(15)는 각 성분의 이동속도로 나타나며, 그것을 특정시간에 있어서의 이동거리로 표시하고, 혹은 컬럼 말단으로부터 시계열에서 나타나는 상태를 차트로 표시한 것이다. HPLC에 있어서는, 시간축에 대한 검출기응답을 플롯한 것이고, 피크의 출현시간을 보지(保持)시간(리텐션타임)이라고 부르고 있다.
검출기로는, 특히 한정되지 않지만, 광학적 성질을 이용한 흡광도검출기(Absorbance Detector)가 사용되고, 피크는 자외선(UV)의 검출파장에 따른 시그널강도로 삼차원적으로 얻어진 것이다. 광학적 성질을 이용한 것으로는 투과율검출기(Transmittance Detector)를 이용할 수도 있다.
검출파장에 한정은 없으나, 바람직하게는 150nm ~ 900nm의 범위이고, 특히 바람직하게는 200nm ~ 400nm의 자외가시흡(吸)영역, 더욱이 바람직하게는 200nm ~ 300nm로부터 선택된 복수의 파장이다.
그리고, 3D크로마토(15)는 적어도 한방약의 번호(로트 번호), 리텐션타임, 검출파장, 및 피크를 데이타로 가지는 것이다.
또한, 3D크로마토(15)는 시판의 장치에 의해서도 얻어질 수 있으며, 이러한 시판의 장치로는 Agilent 1100 시스템 등을 들 수 있다. 또한, 크로노마토그래피는 HPLC에 한정되는 것은 아니며, 다양한 것을 채용할 수가 있다.
3D크로마토(15)는, 도 2, 도 3과 같이 x축을 리텐션타임, y축을 검출파장, z축을 시그널강도로 하여 표시한다.
FP(17)는 적어도 한방약의 번호(로트 번호), 리텐션타임, 특정의 파장에 있어서의 피크, 및 UV스펙트럼을 데이타로 가지는 것이다.
FP(17)는, 도 2, 도 3과 같이 x축을 리텐션타임, y축을 특정의 검출파장에 있어서의 피크로 한 이차원에서 표시하지만, 도 3과 같이 1 피크로 표현하는 UV스펙트럼(25)과 같은 UV스펙트럼 정보를 피크마다에 가진 데이타이다. FP(17)를 작성하는 특정의 검출파장은 특히 한정되지 않으며 다양하게 선택할 수 있다. 단, FP(17)에는 3D크로마토 중의 모든 피크를 포함시키는 것이 정보를 계승한다는 점에서 중요하다. 이 때문에, 본 실시예 1에서는 검출파장을 3D크로마토 중의 모든 피크를 포함하고 있는 203nm로 했다.
한편, 단독의 파장에서는 모든 피크를 포함시킬 수 없는 것이 있다. 그러한 경우는 검출파장을 복수로 하고, 후술하는 바와 같이 복수의 파장을 조합하여 모든 피크를 포함한 FP를 작성한다.
본 실시예 1에서는 피크를 시스널강도(피크높이)의 극대치로 했지만, 피크로서 면적치를 채용할 수도 있다. 또한, FP에 UV스펙트럼을 포함시키지 않고, x축을 리텐션타임, y축을 특정의 검출파장에 있어서의 피크로 한 이차원의 표시정보만으로 하는 것도 가능하다. 이 경우는, FP를 한방약의 번호(로트 번호), 리텐션타임을 데이타로 가지는 크로마토로서의 2D크로마토로부터 작성할 수도 있다.
도 4 (A)는 약제A, (B)는 약제B, (C)는 약제C의 FP이다.
기준FP선정부(5)는 피크패턴작성부(7)에서 이용하는 기준FP를 복수의 기준FP로부터 선정하는 기능부이다. 이 기준FP선정부(5)는 대상FP의 피크귀속에 적합한 다성분 물질의 FP를 복수의 기준FP로부터 선정한다. 즉, 대상FP의 각 피크의 피크귀속을 높은 정도로 행하기 위해, 도 5 ~ 도 9와 같이 대상FP와 기준FP 사이에서 피크의 리텐션타임출현패턴의 일치도를 산출하고, 이 일치도가 최소로 되는 기준FP를 전 기준FP로부터 선정한다. 구체적으로는 후술한다. 피크패턴작성부(7)는, 도 10 ~ 도 12와 같이 대상FP(33)에 있어서의 귀속의 대상이 되는 피크(이하, 귀속대상피크)에 대하여, 시간축방향 전후의 적어도 일방에 존재하는 n개의 피크를 포함한 합계 n+1개의 피크로 구성되는 피크패턴을 귀속대상피크의 피크패턴으로 작성하는 기능부이다. n은 자연수이다. 구체적으로는 후술한다.
또한, 도 11에서는 시간축방향 전후의 적어도 일방에 존재하는 2개의 피크를 포함한 합계 3개의 피크로 구성되는 피크패턴을, 도 12에서는 시간축방향 전후의 적어도 일방에 존재하는 4개의 피크를 포함한 합계 5개의 피크로 구성되는 피크패턴을 도시한다.
또한, 피크패턴작성부(7)는, 도 13 ~ 도 22(후술)와 같이 기준FP(55)에 있어서 귀속대상피크의 리텐션타임과의 차가 설정한 범위(허용폭) 내의 모든 피크(이하, 귀속후보피크)에 대하여, 시간축방향 전후의 적어도 일방에 존재하는 n개의 피크를 포함한 합계 n+1개의 피크로 구성되는 피크패턴을 귀속후보피크의 피크패턴으로 작성하는 기능부이다. 또한, 도 15 ~ 도 18(후술)에서는 시간축방향 전후의 적어도 일방에 존재하는 2개의 피크를 포함한 합계 3개의 피크로 구성되는 피크패턴을 도시한다. 도 19 ~ 도 22(후술)에서는 시간축방향 전후의 적어도 일방에 존재하는 4개의 피크를 포함한 합계 5개의 피크로 구성되는 피크패턴을 도시한다.
허용폭에 한정은 없으나, 정도와 효율의 면에서 0.5분 ~ 2분이 바람직하다. 실시예 1에서는 1분으로 했다.
또한, 피크패턴작성부(7)에서는 대상FP(33)와 기준FP(55)의 피크수에 차이가 있는(즉, 어느 쪽인가 일방에 존재하지 않는 피크가 있는) 경우에 대하여도 유연하게 대응할 수 있도록 한다. 이를 위해, 도 23 ~ 도 61(후술)과 같이, 귀속대상피크 및 귀속후보피크의 양방에서 피크패턴을 구성하는 피크(이하, 피크패턴구성피크)를 변화시켜서 망라적(網羅的)으로 피크패턴을 작성한다. 또한, 도 23 ~ 도 61에서는 시간축방향 전후의 적어도 일방에 존재하는 2개의 피크를 포함한 합계 3개의 피크로 구성되는 피크패턴의 경우에 대해서 도시한다.
피크귀속부(9)는 대상FP, 기준FP 각각의 피크패턴을 비교하여 대응하는 피크를 특정하는 기능부이다. 실시예에서는 귀속대상피크의 피크패턴과 귀속후보피크의 피크패턴의 일치도 및 UV스펙트럼의 일치도를 산출하여 대응하는 피크를 특정한다. 구체적으로는 후술한다.
또한, 이 2개의 일치도를 통합한 귀속후보피크의 일치도를 산출하고, 이 일치도에 기초하여, 대상FP(33)의 각 피크를 기준FP(55)의 각 피크에 귀속하는 기능부이다.
또한, 대상FP(33)와 기준FP(55)의 귀속결과에 기초하여, 최종적으로 도 68, 도 69(후술)와 같이 대상FP의 각 피크를 기준군FP의 각 피크에 귀속하는 기능부이다.
피크귀속부(9)에서 피크패턴의 일치도는, 도 62 ~ 도 64(후술)와 같이 귀속대상피크와 귀속후보피크의 피크패턴 사이의 대응하는 피크 및 리텐션타임의 차를 바탕으로 산출한다. 또한, UV스펙트럼의 일치도는, 도 65, 도 66(후술)과 같이 귀속대상피크(45)의 UV스펙트럼(107)과 귀속후보피크(67)의 UV스펙트럼(111)의 각 파장에 있어서의 흡광도(吸光度)의 차를 기초로 산출한다. 또한, 도 67(후술)과 같이 이들 2개의 일치도를 곱하여 귀속후보피크(67)의 일치도를 산출한다.
평가부(11)는 피크귀속부(9)에서 특정되어 귀속한 피크와 상기 복수의 기준FP의 피크를 비교평가하는 기능부이다. 실시예에서는, MT법으로 대상FP귀속피크(21)와 기준군FP(19)의 동등성을 평가하는 기능부이다.
MT법은 품질공학에서 일반적으로 알려져 있는 계산법을 의미한다. 예를 들면, '품질공학의 수리(數理)' 일본규격협회발행(2000) 제 136 - 138 면, 품질공학응용강좌 '화학·약학·생물학의 기술개발' 일본규격협회편(1999) 제 454 - 456 면 및 품질공학 11 (5), 78 - 84 (2003), 입문 MT시스템(2008)에 기재가 있다.
또한, 일반에 시판되고 있는 MT법 프로그램 소프트도 사용할 수 있다. 시판의 MT법 프로그램 소프트로서는 앵클드라이(주)의 ATMTS; (재) 일본규격협회의 TM-ANOVA; (주) 오켄의 MT법 for windows 등을 들 수 있다.
평가부(11)는 대상FP(17) 중, 한방약의 로트 번호와, 리텐션타임 또는 UV검출파장의 일방에 대하여 MT법에서의 변수축을 할당하고, 피크를 MT법에서의 특징량으로 한다.
변수축의 할당에는 특히 한정은 없지만, MT법에서의 소위 항목축에 리텐션타임을 할당하고, 소위 번호열축에 다성분계 약제의 번호를 할당하고, MT법에서의 소위 특징량에 피크를 할당하는 것이 바람직하다.
여기서, 상기 항목축과 번호열축은 이하와 같이 정의된다. 즉, MT법에 있어서는, 데이타 세트 Xij에 대해서, 평균치mj와 표준편차σj를 구하고, Xij를 규준화(規準化)한 값인 xij = (Xij - mj)/σj로부터, i와 j의 상관계수 r을 구하여, 단위공간이나 마하라노비스의 거리를 얻는다. 이때, '평균치mj와 표준편차σj는 항목축의 값마다, 번호열축의 값을 변화시켜서 구한다.'와 같이 항목축과 번호열축이 정의된다.
축이 할당된 데이타와 특징량으로부터, MT법을 이용하여, 기준점과 단위량(이하, '단위공간'으로 약기하는 경우도 있다.)를 얻는다. 여기서, 기준점, 단위량 및 단위공간은 상기 MT법의 문헌의 기재에 따라서 정의된다.
MT법에 의해, 평가해야하는 약제의 단위공간과의 상위의 정도를 표시하는 값으로서 MD값을 얻는다. 여기서 MD값은 MT법의 문헌의 설명과 마찬가지로 정의되고, 또한, MD값은 문헌에 기재된 방법으로 구해진다.
이와 같이 하여 얻어진 MD값을 이용하여, 평가해야할 약제는 정상품으로서 정해진 복수의 약제로부터의 상위(相違)의 정도를 판정하여 평가할 수가 있다.
예를 들면, 도 70 ~ 도 74의 각 대상FP를 상기와 같이 귀속처리하는 것으로, 상기 MT법에 의해 MD값(MD값 : 0.25, 2.99 등)을 구할 수 있다.
이 MD값을 정상품의 MD값에 대하여 평가하는 경우, 정상품으로 정해진 복수의 약제에서 마찬가지로 MD값을 구한다. 이 정상품의 MD값으로부터 역치를 설정하고, 도 2의 평가부(11)의 평가결과(23)에 도시하는 바와 같이 평가대상약제의 MD값을 플로트하여, 정상품인지 비정상품인지의 판정을 행할 수 있다. 도 2의 평가부(11)의 평가결과(23)에서는, 예를 들면 MD값 10 이하를 정상품으로 하고 있다.
또한, 평가부(11)는 대상FP귀속피크(21)와 기준군FP(19)의 동등성을 비교평가할 수 있다면 충분하기 때문에, MT법 이외의 패턴인식수법 등을 적용하는 것도 가능하다.
[피크패턴처리의 동작원리]
도 5 ~ 도 67은 상기 기준FP선정부(5), 피크패턴작성부(7), 피크귀속부(9) 및 평가부(11)의 동작원리를 설명하는 것이다.
도 5 ~ 도 9는 기준FP선정부(5)에 관한 대상FP와 기준FP의 리텐션타임출현패턴의 일치도에 대하여 설명한 도면이다. 도 5는 대상FP 및 기준FP의 리텐션타임을 도시하는 도면, 도 6은 대상FP의 리텐션타임을 도시하는 도면, 도 7은 기준FP의 리텐션타임출현패턴을 도시하는 도면, 도 8은 대상FP와 기준FP의 리텐션타임출현거리의 일치수를 도시하는 도면, 도 9는 대상FP와 기준FP의 리텐션타임출현패턴의 일치도를 도시하는 도면이다.
도 5에서는, 대상FP(33) 및 기준FP(55) 각각의 리텐션타임을 도시한다. 도 6, 도 7에서는 대상FP(33) 및 기준FP(55) 각각의 리텐션타임으로부터 모든 리텐션타임간 거리를 산출하고, 그들 거리를 표형식으로 정리한 리텐션타임출현패턴을 도시한다. 도 8에서는, 이들 출현패턴으로부터 리텐션타임출현거리의 일치수를 산출하고, 이들 일치수를 표형식으로 정리한 리텐션타임출현거리의 일치수를 도시한다. 도 9에서는, 이 일치수를 기초로 하여 리텐션타임출현패턴의 일치도를 산출하고, 이들 일치도를 표형식으로 정리한 리텐션타임출현패턴의 일치도를 도시한다.
도 10 ~ 도 12는 피크패턴작성부(7)에 관련한 귀속대상피크와 그 주변피크에서 작성한 피크패턴에 대하여 설명한 도면이다. 도 10은 기준FP의 귀속대상피크를 도시하는 도면, 도 11은 주변피크 2개를 포함한 피크 3개로 작성한 피크패턴에 대하여, 도 12는 주변피크 4개를 포함한 피크 5개로 작성한 피크패턴에 대하여 설명한 도면이다.
도 13 , 도 14는 피크패턴작성부(7)에 관련한 귀속대상피크와 귀속후보피크의 관계에 대해 설명하며, 도 13은 귀속대상피크의 허용폭을 도시하는 도면, 도 14는 귀속대상피크에 대한 기준FP의 귀속후보피크를 도시하는 도면이다.
도 15 ~ 도 18은 피크패턴작성부(7)에 관련한 피크 3개로 작성한 귀속대상피크와 귀속후보피크의 피크패턴의 예이다. 도 15는 귀속대상피크와 귀속후보피크의 피크 3개에 따른 피크패턴도, 도 16은 귀속대상피크와 별도의 귀속후보피크의 피크 3개에 따른 피크패턴도, 도 17은 귀속대상피크와 별도의 귀속후보피크의 피크 3개에 따른 피크패턴도, 도 18은 귀속대상피크와 별도의 귀속후보피크의 피크 3개에 따른 피크패턴도이다.
도 19 ~ 도 22는 피크패턴작성부(7)에 관련한 피크 5개로 작성한 귀속대상피크와 귀속후보피크의 피크패턴도이다.
도 23 ~ 도 61은 피크패턴작성부(7)에 관련한 귀속대상피크 및 귀속후보피크의 피크패턴을 망라적으로 작성하고, 비교하는 망라적 비교의 원리를 설명하기 위한 도면이다.
도 62, 63은 피크귀속부(9)에 관련한 피크 3개로 작성한 피크패턴의 일치도의 산출방법에 대하여 설명하는 도면이다.
도 64는 피크귀속부(9)에 관련한 피크 5개로 작성한 피크패턴의 일치도의 산출방법에 대하여 설명하는 도면이다.
도 65는 피크귀속부(9)에 관련한 귀속대상피크(45) 및 귀속후보피크(67)의 UV스펙트럼(107 및 111)을 도시하는 도면이다.
도 66은 피크귀속부(9)에 관련한 귀속대상피크(45)의 UV스펙트럼(107)과 귀속후보피크(67)의 UV스펙트럼(111)의 일치도에 대하여 설명한 도면이다.
도 67은 피크귀속부(9)에 관련한 귀속대상피크(45)와 귀속후보피크(67)의 피크패턴의 일치도와 UV스펙트럼의 일치도로부터 산출하는 귀속후보피크의 일치도에 대해 설명한 도면이다.
도 68은 피크귀속부(9)에 관련한 대상FP(17)에 있어서의 각 피크의 기준군FP(19)로의 귀속을 설명한 도면이다.
도 69는 피크귀속부(9)에 관련한 대상FP(17)의 각 피크가 기준군FP(19)에 귀속된 상황을 도시하는 대상FP피크특징량(21)을 설명한 도면이다.
도 70 ~ 도 74는 평가부(11)에 관련한 각종 대상FP와 그 평가치(MD값)을 도시한 도면이다.
[기준FP의 선정]
상기 기준FP선정부(5)의 기능을, 도 5 ~ 도 9를 이용하여 상세히 설명한다.
도 5는 대상FP 및 기준FP의 리텐션타임을 도시하는 도면, 도 6은 기준FP의 리텐션타임출현패턴을 도시하는 도면, 도 7은 기준FP의 리텐션타임출현패턴을 도시하는 도면, 도 8은 대상FP와 기준FP의 리텐션타임출현거리의 일치수를 도시하는 도면, 도 9는 대상FP와 기준FP의 리텐션타임출현패턴의 일치도를 도시하는 도면이다.
도 5에서는, 대상FP(33) 및 기준FP(55) 각각의 리텐션타임을 도시한다. 도 6, 도 7에서는 대상FP(33) 및 기준FP(55) 각각의 리텐션타임으로부터 모든 리텐션타임간 거리를 산출하고, 그들 거리를 표형식으로 정리한 리텐션타임출현패턴을 도시한다. 도 8에서는, 이들 출현패턴으로부터 리텐션타임출현거리의 일치수를 산출하고, 이들 일치수를 표형식으로 정리한 리텐션타임출현거리의 일치수를 도시한다. 도 9에서는, 이 일치수를 기초로 하여 리텐션타임출현패턴의 일치도를 산출하고, 이들 일치도를 표형식으로 정리한 리텐션타임출현패턴의 일치도를 도시한다.
대상FP(33)의 피크귀속처리에 있어서, 대상FP(33)와 가능한한 FP패턴이 유사한 기준FP로 대상FP(33)의 각 피크를 귀속한다. 이 대상FP(33)에 유사한 기준FP를 복수의 기준FP로부터 선정하는 것이 정도가 높은 귀속을 행하는 것에 중요한 포인트이다.
그리고, 대상FP(33)의 FP패턴과의 유사성을 객관적으로 또한 간이적(簡易的)으로 평가하는 방법으로 하여, 리텐션타임출현패턴의 일치도에 의해 FP패턴의 유사성을 평가한다.
예를 들면, 대상FP(33) 및 기준FP(55)의 리텐션타임이 도 5와 같은 경우, 대상FP(33) 및 기준FP(55) 각각의 리텐션타임출현패턴은 도 6, 도 7과 같이 된다. 도 6, 도 7에서는 상단의 대상FP(33) 및 기준FP(55)에 대하여, 하단의 도표와 같이, 각 셀의 값이 리텐션타임간 거리로 구성된 표형식의 패턴으로서 작성되어 있다.
도 6에서, 대상FP(33)의 각 피크(35, 37, 39, 41, 43, 45, 47, 49, 51, 53)의 리텐션타임은 (10, 2), (10, 5), (10, 8), (11, 1), (11, 6), (12, 1), (12, 8), (13, 1), (13, 6), (14, 0)으로 되어 있다.
따라서, 피크(35) 및 피크(37) 사이의 리텐션타임간 거리는 (10, 5) - (10, 2) = (0, 3)이 된다. 마찬가지로, 피크(35) 및 피크(39) 사이는 (0, 6), 피크(37) 및 피크(39) 사이는 (0, 3) 등이 된다. 이하, 마찬가지이며, 도 6의 하단도표의 대상FP출현패턴(79)으로 된다.
도 7에서, 기준FP(55)의 각 피크(57, 59, 61, 63, 65, 67, 69, 71, 73, 75, 77)의 리텐션타임은 (10, 1), (10, 4), (10, 7), (11, 1), (11, 7), (12, 3), (12, 7), (13, 1), (13, 6), (14, 1), (14, 4)로 되어 있다.
따라서, 마찬가지로 리텐션타임간 거리는 도 7의 하단도표의 기준FP출현패턴(81)으로 된다.
이 도 6, 도 7과 같이 패턴화한 각 피크를 전체개별대비방식으로 비교하여 일치수를 구한다. 예를 들면, 도 6 하단도표의 대상FP출현패턴의 각 셀의 값과 도 7 하단도표의 기준FP출현패턴의 각 셀의 값을 비교하여, 도 8과 같이 일치수(83)를 얻는다.
즉, 대상FP(33)과 기준FP(55)의 리텐션타임출현패턴의 모든 리텐션타임간 거리를 행단위로 순번으로 전체개별대비방식으로 비교하여, 설정한 범위 내에서 거리가 일치한 수를 산출했다.
예를 들면, 도 6, 도 7의 대상 및 기준FP 리텐션타임출현패턴(79, 81)의 1행을 비교하면, 일치수는 7개이다. 이 7개의 일치수가 도 8의 대상 및 기준FP리텐션타임출현패턴의 1행째에 기재된다. 도 6, 도 7 중의 그 외의 행에 대해서도 마찬가지이며, 대상FP리텐션타임출현패턴의 1행 ~ 9행까지와, 기준FP리텐션타임출현패턴의 1행 ~ 10행까지를 전체개별대비방식으로 비교하여, 각각 일치수가 얻어진다.
도 8에 결과를 도시했다. 이 도 8에 있어서, 동그라미로 둘러싸인 좌단의 7의 수치는 대상 및 기준FP 리텐션타임출현패턴의 1행째를 비교한 결과이고, 그 옆의 7의 수치는 대상FP리텐션타임출현패턴의 1행째와 기준FP리텐션타임출현패턴의 2행째를 비교한 결과이다. 설정의 범위에 한정은 없지만, 바람직하게는 0.05분 ~ 0.2분의 범위이다. 실시예 1은 0.1분으로 했다.
리텐션타임출현패턴의 일치도를 RP로 하면, 대상FP(33)의 f행째의 리텐션타임출현패턴과 기준FP(55)의 g행째의 리텐션타임출현패턴의 일치도(RPfg)는 Tanimoto 계수를 이용하여,
RPfg = {1 - ( m / (a + b - m ) ) } × ( a - m + 1 )
로 산출한다.
또한, 식중 a는 대상FP(33)의 피크수(대상FP피크수), b는 기준FP(55)의 피크수(기준FP피크수), m은 출현거리의 일치수이다(도 8 참조). 도 8의 일치수(83)를 기초로 상기 식에 의해 각 리텐션타임출현패턴의 일치도(RP)를 산출했다(도 9의 일치도(85) 참조).
이들 RP의 최소치(RP_min)를 대상FP(33)와 기준FP(55)의 리텐션타임출현패턴의 일치도로 한다. 도 9에서는 (0.50)이 대상FP(33)의 기준FP에 대한 일치도가 된다.
이 일치도를 모든 기준FP에 대하여 산출하고, 가장 작은 일치도의 기준FP가 선정되고, 이 기준FP에 대하여 대상FP의 피크귀속을 행하게 한다.
기준FP선정부(5)는 대상FP(33) 및 기준FP(55)를 피크높이비로 패턴화할 수도 있다.
피크높이비로 패턴화한 각 피크를 전체개별대비방식으로 비교하여, 설정한 범위 내에서 높이비가 일치한 수를 산출한다. 이 산출에 의해 도 8과 마찬가지로 일치수를 얻을 수가 있다.
또한, 피크높이비로 패턴화하는 경우는, 1행 중에 같은 값이 복수 존재하는 케이스가 있으며, 이들을 복수회 카운트하지 않도록 해야한다.
일치도는 Tanimoto계수를 '높이비의 일치수/(대상FP피크수 + 기준FP피크수 - 높이비의 일치수)'로 하고, (1 - Tanimoto계수)가 영에 가까운 것으로 상기 일치도를 구할 수가 있다.
또한, (1 - Tanimoto계수)에 (대상FP피크수 - 높이비의 일치수 + 1)의 가중치부여를 하여, '(1 - Tanimoto계수) × (대상FP피크수 - 출현거리 또는 높이비의 일치수 + 1)'로 하여, 가중치에 의해 대상FP(33)의 피크(35, 37, ...)가 보다 많이 일치하는 기준FP를 선택할 수가 있다.
[피크패턴에 따른 특징량화]
상기 피크패턴작성부(7) 및 피크귀속부(9)의 기능을 도 10 ~ 도 69를 이용하여 상세히 설명한다.
도 10과 같이 귀속대상피크(45)를 기준FP(55)의 어느 하나의 피크에 귀속할 때, 어느 피크에 귀속시켜야 하는가라는 것이 된다. 가령, 이 피크귀속을 피크, 리텐션타임 또는 UV스펙트럼의 어느 하나의 정보만으로 행했다고 하면, 이들 3개의 정보의 어느 하나가 상기 약제간 오차와 분석 오차에 기인하는 오차를 포함하고 있기 때문에, 단독정보에 따른 피크귀속으로는 충분한 정도를 얻을 수 없다.
또한, 도 13, 도 14와 같이 귀속대상피크(45)와 기준FP(55)의 각 피크간에 리텐션타임의 어긋남의 허용폭을 설정하여, 그 허용폭 내에 존재하는 기준FP(55)의 피크와 UV스펙트럼정보의 2개의 정보에 따른 피크귀속에서는 모든 정보를 총합(總合)하여 귀속처를 판정하고 있기 때문에, 상기 단독정보에 따른 피크귀속에 비교하여 정도(精度)는 향상된다.
그러나, 3개의 정보를 사용한 피크귀속이었다고 해도, UV스펙트럼의 특성으로서, 유사성분의 UV스펙트럼은 거의 동일하게 되기 때문에, 귀속후보피크에 복수의 유사성분이 포함되어 있는 경우는 결국 피크정보만으로의 귀속이 되어버려서, 충분한 정도는 얻어질 수 없다. 그 때문에, 보다 정도가 높은 피크귀속을 행하기 위해서는, 이들 3개의 정보에 플러스하는 정보가 필요하다.
그래서, 도 11, 도 12와 같은 주변 피크의 정보를 포함시킨 피크패턴을 작성하여, 이 피크패턴의 비교에 의해 피크를 귀속하는 것으로 했다.
주변의 피크를 포함시킨 피크패턴으로 하면, 지금까지의 3개의 정보에 주변정보가 플러스되어, 4개의 정보에 따른 피크귀속이 가능하게 되어, 보다 높은 귀속정도가 얻어질 수 있다.
그 결과로, 한 번의 귀속처리로 대량의 피크를 고정도로 또한 효율적으로 일제(一齊)귀속할 수 있다.
또한, 피크귀속에서 사용하는 데이타를 주변정보를 포함한 4개의 정보로 하는 것으로, 기존의 피크귀속시에 설정하는 제약조건(피크정의 등)도 필요 없게 되었다.
도 11에서는, 귀속대상피크(45)에 대하여, 시간축방향의 쌍방에 존재하는 피크(43, 47)를 포함한 피크패턴(87)을 작성했다.
도 12에서는 귀속대상피크(45)에 대하여, 시간축방향의 쌍방에 존재하는 피크(41, 43, 47, 49)를 포함한 피크패턴(97)을 작성했다.
도 13, 도 14에서는, 귀속대상피크(45)와 기준FP(55)의 각 피크간에 리텐션타임의 어긋남의 허용폭을 설정하여, 그 허용폭 내에 존재하는 기준FP(55)의 피크를 귀속대상피크(45)에 대응하는 후보피크(이하, 귀속후보피크)로 했다.
도 15에서는 귀속대상피크(45)의 피크패턴(87)과 비교하는 피크패턴으로서 귀속후보피크(65)에 대하여 시간축방향 전후의 쌍방에 존재하는 피크(63, 67)를 포함시킨 피크패턴(89)을 작성했다.
도 16 ~ 도 18에서는 귀속대상피크(45)의 피크패턴(87)과 비교하는 피크패턴으로서 별도의 귀속후보피크(67, 69, 71)에 대하여 각각 시간축방향 전후의 쌍방에 존재하는 피크를 포함시킨 피크패턴(91, 93, 95)을 작성했다.
이 피크패턴의 비교를 보다 높은 정도로 행하게 하기 위해서는, 도 19 ~ 도 22와 같이, 대상FP 및 기준FP의 양방에서 주변의 피크수를 증가시킨 피크패턴을 작성하는 것이 긴요(肝要)하다.
예를 들면, 주변의 피크 4개를 포함시킨 합계 5개의 피크에 따른 피크패턴을 비교하면, 보다 높은 귀속정도가 얻어질 수 있다.
도 19에서는, 귀속대상피크(45)의 피크패턴(97)과 비교하는 피크패턴으로서 귀속후보피크(65)에 대하여 시간축방향 전후의 쌍방에 존재하는 피크(61, 63, 67, 69)를 포함시킨 피크패턴(99)을 작성했다.
도 20 ~ 도 22에서는 귀속대상피크(45)의 피크패턴(97)과 비교하는 피크패턴으로서 별도의 귀속후보피크(67, 69, 71)에 대하여 각각 시간축방향 전후의 쌍방에 존재하는 피크를 포함시킨 피크패턴(101, 103, 105)을 작성했다.
또한, 이 피크패턴에 따른 귀속을 보다 높은 정도(精度)로 행하게 하기 위해서는, 대상FP와 기준FP의 피크수가 틀리는(즉, 어느 쪽이 일방에 존재하지 않는 피크가 있는) 경우의 대응이 필요하다. 그 때문에, 도 23 ~ 도 25와 같이 귀속대상피크 및 귀속후보피크의 양방에서 피크패턴구성피크를 망라적으로 변화시킨 피크패턴을 작성하는 것이 긴요하다.
구체적으로는 대상FP의 귀속대상피크의 주변피크 중에서 미리 피크패턴구성피크의 후보가 되는 피크(이하, 피크패턴구성후보피크)를 설정하고, 이 피크패턴구성후보피크를 순번으로 피크패턴구성피크로 하여 피크패턴을 작성한다. 기준FP의 귀속후보피크에 대해서도 마찬가지로 피크패턴구성후보피크를 설정하고, 이 피크패턴구성후보피크를 순번으로 피크패턴구성피크로 하여 피크패턴을 작성한다.
예를 들면, 도 23과 같이 귀속대상피크(45)의 피크패턴구성후보피크로서 시간축방향 주변 4개(41, 43, 47, 49), 귀속후보피크(65)의 피크패턴구성후보피크로서 시간축방향 주변 4개(61, 63, 67, 69)로 하고, 피크패턴구성피크를 임의의 2개로 각각 설정한다. 이 경우, 도 24, 도 25와 같이 귀속대상피크(45) 및 귀속후보피크(65) 각각에서 피크패턴이 4C2 (=6) 패턴 작성된다.
더욱이, 피크패턴구성후보피크를 10개로 하고, 피크패턴구성피크를 임의의 2개로 설정하면, 귀속대상피크 및 귀속후보피크 각각에서 10C2 (=45) 패턴의 피크패턴이 작성된다. 피크패턴구성피크를 임의의 4개로 설정하면, 귀속대상피크 및 귀속후보피크 각각에서 10C4 (=210) 패턴의 피크패턴이 작성된다.
상기 피크귀속부(9)의 기능을, 도 26 ~ 67을 이용하여 상세히 설명한다.
피크귀속부(9)에서는 피크패턴작성부(7)에서 작성한 귀속대상피크와 귀속후보피크의 모든 피크패턴 사이에, 대응하는 피크 및 리텐션타임의 차를 기초로 피크패턴의 일치도(이하, P_Sim)를 산출한다. 피크귀속부(9)는 P_Sim의 최소치(이하, P_Sim_min)를 귀속대상피크와 귀속후보피크의 피크패턴의 일치도로 한다.
예를 들면, 도 26 ~ 도 61과 같이 귀속대상피크(45) 및 귀속후보피크(65)의 각각에서 피크패턴구성후보피크를 시간축방향 주변 4개로 하고, 피크패턴구성피크를 임의의 2개로 설정한다. 이 설정에서는, 귀속대상피크 및 귀속후보피크 각각에서 4C2 (=6) 패턴의 피크패턴이 작성된다. 따라서, 귀속대상피크(45) 및 귀속후보피크(65)의 P_Sim은 6 패턴 × 6 패턴 (= 36) 조로 계산되고, 이들 P_Sim의 최소치인 P_Sim_min을 귀속대상피크(45)와 귀속후보피크(65)의 일치도로 한다.
덧붙여, 귀속대상피크(45) 및 귀속후보피크(65)의 각각에서 피크패턴구성후보피크를 시간축방향 전후 주변 10개로 하고, 피크패턴구성피크를 임의의 2개로 설정하면, 귀속대상피크 및 귀속후보피크 각각에서 4C2 (=45) 패턴의 피크패턴이 작성된다. 따라서, 귀속대상피크(45) 및 귀속후보피크(65)의 P_Sim은 45 패턴 × 45 패턴 (= 2025) 조로 계산되고, 이들 P_Sim의 최소치인 P_Sim_min을 귀속대상피크(45)와 귀속후보피크(65)의 일치도로 한다. 또한, 피크패턴구성피크를 임의의 4개로 설정하면, 귀속대상피크 및 귀속후보피크 각각에서 피크패턴이 10C4 (=210) 패턴의 피크패턴이 작성된다. 따라서, 귀속대상피크(45) 및 귀속후보피크(65)의 P_Sim은 210 패턴 × 210 패턴 (= 44100)으로 계산되고, 이들 P_Sim의 최소치인 P_Sim_min을 귀속대상피크(45)와 귀속후보피크(65)의 일치도로 한다.
이 P_Sim은 귀속대상피크(45)의 모든 귀속후보피크에 대하여 마찬가지로 산출한다.
도 62, 도 63에서, 3개의 피크로 구성된 피크패턴을 비교하기 위한 피크패턴의 일치도의 산출방법을 설명한다. 이 경우, 귀속대상피크(45)의 피크패턴(87)과 귀속후보피크(67)의 피크패턴(91)을 예로 한다.
귀속대상피크(45)의 피크패턴(87)에서, 귀속대상피크(45)의 피크데이타 및 리텐션타임을 p1 및 r1, 피크패턴구성피크(43)의 피크데이타 및 리텐션타임을 dn1 및 cn1, 피크패턴구성피크(47)의 피크데이타 및 리텐션타임을 dn2 및 cn2로 한다.
귀속후보피크(67)의 피크패턴(91)에서, 귀속후보피크(67)의 피크데이타 및 리텐션타임을 p2 및 r2, 피크패턴구성피크(65)의 피크데이타 및 리텐션타임을 fn1 및 en1, 피크패턴구성피크(69)의 피크데이타 및 리텐션타임을 fn2 및 en2로 한다.
피크패턴의 일치도를 P_Sim으로 하면, 귀속대상피크(45)와 귀속후보피크(67)의 3개의 피크로 구성하는 피크패턴의 일치도(P_Sim (45 - 67))는
P_Sim (45 - 67) = (│p1 - p2│+1)×(│r1 - (r2 + d)│+ 1)
+ (│dn1 - fn1│+ 1)×(│cn1 - r1) - (en1 - r2│ + 1)
+ (│dn2 - fn2│+ 1)×(│cn2 - r1) - ( en2 - r2│ + 1)로 산출한다.
또한, 식 중의 d는 리텐션타임의 어긋남을 보정하는 값이다.
도 64에서, 5개의 피크로 구성된 피크패턴을 비교하기 위한 피크패턴의 일치도의 계산방법을 설명한다. 이 경우, 귀속대상피크(45)의 피크패턴(97)과 귀속후보피크(67)의 피크패턴(101)을 예로 한다.
귀속대상피크(45)의 피크패턴(97)에서, 귀속대상피크(45)의 피크데이타 및 리텐션타임을 p1 및 r1, 피크패턴구성피크(41, 43, 47, 49)의 피크데이타 및 리텐션타임을 각각 dn1 및 cn1, dn2 및 dn2, dn3 및 cn3, dn4 및 cn4로 한다.
귀속후보피크(67)의 피크패턴(101)에서, 귀속후보피크(67)의 피크데이타 및 리텐션타임을 p2 및 r2, 피크패턴구성피크(63, 65, 69, 71)의 피크데이타 및 리텐션타임을 각각 fn1 및 en1, fn2 및 en2, fn2 및 en3, fn3 및 en3, fn4 및 en4로 한다.
피크대상피크(45)와 귀속후보피크(67)의 5개의 피크로 구성하는 피크패턴의 일치도(P_Sim (45 - 67))는
P_Sim (45 - 67) = (│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)
로 산출한다.
또한, 식 중의 d는 리텐션타임의 어긋남을 보정하는 값이다.
피크귀속부(9)에서는, 도 65, 도 66과 같이 귀속대상피크와 귀속후보피크에서 UV스펙트럼의 일치도를 산출한다.
도 65는 귀속대상피크(45) 및 귀속후보피크(67)의 UV스펙트럼(107과 111)도이고, 도 66과 같이 이들 2개의 UV스펙트럼의 일치도(UV_Sim (45 - 67))는
UV_Sim (45 - 67) = RMSD (107 vs 111)
로 산출한다.
RMSD는 평균이승(平均二乘)편차이며, 대응하는 2점의 거리(dis)를 각각 이승하고, 그 상가(相加)평균의 평균근으로서 정의된다. 즉,
RMSD = √ { dis2 / n }
으로 산출한다.
n은 dis의 수이다.
여기서, UV스펙트럼의 파형은 극대파장 및 극소파장을 포함하고 있고, 이 극대파장 및 극소파장 혹은 어느 하나를 비교하는 것으로 일치도를 산출하는 것도 가능하다. 그러나, 흡수특성이 없는 화합물 혹은 흡수특성이 유사한 화합물 등에서는 극대파장, 극소파장은 동일하지만, 전체의 파장이 상당히 상이한 경우도 있어서, 극대파장과 극소파장의 비교로는 파형의 일치도를 산출할 수 없을 우려가 있다.
이것에 대하여, UV스펙트럼의 파형에 의해 RMSD룰 이용한 경우에는, 파형전체의 비교가 되기 때문에, UV스펙트럼의 파형의 일치도의 산출을 보다 정확하게 행하여, 흡수특성이 없는 화합물 혹은 흡수특성이 유사한 화합물에서도 정확하게 식별할 수 있게 된다.
이 UV스펙트럼의 일치도는 귀속대상피크(45)의 전체의 귀속후보피크에 대하여 마찬가지로 산출한다.
더욱이, 피크귀속부(9)에서는, 도 67과 같이 상기 2개의 일치도를 통합한 귀속후보피크의 일치도를 산출한다.
도 67과 같이 귀속후보피크의 일치도(SCORE (45 - 67))는 피크패턴과 UV스펙트럼의 각각의 일치도를 곱하여 산출한다. 피크패턴(45, 67)의 일치도를 표시하는 스코어는 P_Sim_min(45-67)이고, 대응하는 UV스펙트럼(107, 111)의 일치도를 표시하는 스코어가 UV_Sim(45-67)인 것으로 한다. 이 때, 귀속후보피크의 일치도 SOCORE(45-67)는
SCORE(45-67)
= P_Sim_min(45-67)×UV_Sim(45-67)
로 산출한다.
이 귀속후보피크의 일치도를 귀속대상피크(45)의 모든 귀속후보피크에 대하여 마찬가지로 계산한다.
그리고, 전(全) 귀속후보피크 사이에서 이 SCORE를 비교하여, SCORE가 최소가 되는 귀속후보피크를 귀속대상피크(45)의 귀속피크로서 판정한다.
피크귀속부(9)에서는 귀속대상피크가 귀속해야할 피크를 2개의 관점을 총합(總合)하여 판정하기 때문에, 정확한 피크귀속을 실현할 수 있다.
또한, 피크귀속부(9)에서는 대상FP의 기준FP로의 귀속결과를 기초로, 도 68과 같이, 대상FP(17)의 각 피크를 기준군FP(19)로 귀속한다.
대상FP(17)의 각 피크는 상기 귀속처리에 의해 기준군FP를 구성하는 기준FP에 귀속된다. 이 귀속결과를 기초로 최종적으로 기준군FP(19)의 피크에 귀속한다.
또한, 기준군FP(19)는 정상품(正常品)으로 정해진 복수의 기준FP 전체를 상기와 같이 귀속처리하여 작성한 것이고, 그 각 피크는 귀속된 피크의 평균치(흑점) ± 표준편차(종선)로 표시하고 있다.
도 69는 대상FP(17)를 기준군FP(19)로 귀속한 결과이고, 이 결과가 대상FP(17)의 귀속처리의 최종결과이다.
[MD값]
이 결과로부터, 상기와 같이 MT법에 의해 MD값(MD값 : 0.25, 2.99 등)을 구할 수가 있다. (도 70 ~ 도 74 참조)
[다성분 약제의 평가방법]
도 75는 본발명 실시예 1의 다성분 약제의 평가방법을 도시하는 공정도이다.
도 75와 같이, 패턴의 평가방법으로서의 다성분 약제의 평가방법은 패턴취득공정으로서의 FP작성공정(113)과, 기준패턴선정공정으로서의 기준FP선정공정(115)과, 피크패턴작성공정(117)과, 피크귀속공정(119)과, 평가공정(121)을 구비하고 있다. 이들 FP작성공정(113), 기준FP선정공정(115), 피크패턴작성공정(117), 피크귀속공정(119), 평가공정(121)은 본 실시예에 있어서 상기 다성분 약제의 평가장치(1)를 이용하여 행해지고, FP작성공정(113)은 FP작성부(3)의 기능에 의해 행해지고, 마찬가지로 기준FP선정공정(115), 피크패턴작성공정(117), 피크귀속공정(119), 평가공정(121)을 기준FP선정부(5), 피크패턴특정부(7), 피크귀속부(9), 평가부(11)의 기능에 의해 행하게 할 수 있다.
상기 FP작성공정(113)은 FP작성방법으로서 구비되어, 크로마토를 리텐션타임, 검출파장, 및 피크를 데이타로 하는 3D크로마토로 하여, 이 3D크로마토(15)로부터 특정 파장에서 검출된 피크와 그 리텐션타임 및 그 피크의 UV스펙트럼으로 FP(17)를 작성한다.
[다성분 약제의 평가프로그램]
도 76 ~ 도 91은 다성분 약제의 평가프로그램에 관한 플로우차트, 도 92는 3D크로마토의 데이타 예를 도시하는 도표, 도 93은 피크정보 데이타 예를 도시하는 도표, 도 94는 FP데이타 예를 도시하는 도표, 도 95는 스탭 S3에서 작성한 판정결과파일 예를 도시하는 도표, 도 96은 대상FP와 기준FP에서 대응하는 피크를 특정하는 과정에서 작성하는 2개의 중간파일(귀속후보피크스코어표, 귀속후보피크번호표) 예를 도시하는 도표, 도 97은 대상FP와 기준FP에 대응하는 피크를 특정한 결과에서의 조합(照合)결과파일 예를 도시하는 도표, 도 98은 기준군FP의 데이타 예를 도시하는 도표, 도 99는 대상FP귀속피크의 데이타인 대상FP의 피크데이타특징량파일 예를 도시하는 도표이다.
도 76은 평가대상약제를 평가하기 위한 처리 전체의 스탭을 도시하는 플로우차트이며, 시스템 기동에 의해 스타트하고, FP작성부(3)의 FP작성기능과, 기준FP선정부(5)의 기준FP선정기능과, 피크패턴작성부(7)의 피크패턴작성기능과, 피크귀속부(9)의 피크귀속기능과, 평가부(11)의 평가기능을 컴퓨터에 실현한다.
FP작성기능은 스탭 S1에서 실현되고, 기준FP선정기능은 스탭 S2에서 실현되고, 피크패턴작성기능은 스탭 S3에서 실현되고, 피크귀속기능은 스탭 S3 ~ S5에서 실현되고, 평가기능은 스탭 S6, S7에서 실현된다.
스탭 S1은 3D크로마토 및 특정의 검출파장에서의 피크정보를 입력데이타로 하여 'FP작성처리'가 실행된다.
3D크로마토는 평가대상약제를 HPLC로 분석하는 것에 의해 얻어지는 데이타이고, 도 92의 3D크로마토의 테이타 예(123)에서 도시하듯이, 리텐션타임, 검출파장, 피크(시그널강도)의 삼차원정보로 구성된 데이타이다. 피크정보는 동 HPLC분석에 의해 얻어지는 특정파장에 있어서의 크로마토 데이타를 HPLC데이타분석루트(예를 들면, ChemStation 등)로 처리하여 얻어지는 데이타이고, 도 93의 피크정보데이타 예(125)와 같이, 피크로서 검출된 전(全) 피크의 극대치 및 면적치와 그 시점의 리텐션타임 등으로 구성된 데이타이다.
스탭 S1에서는 컴퓨터의 FP작성부(3)(도 1)가 기능하여, 3D크로마토 및 피크정보로부터 상기 대상FP(17)(도 2)를 작성하여, 그 데이타를 파일로서 출력한다. 이 대상FP(17)는, 도 94의 FP데이타 예(127)와 같이, 리텐션타임과 피크높이와 피크높이마다의 UV스펙트럼으로 구성된 데이타이다.
스탭 S2는 스탭 S1에서 출력된 대상FP 및 전 기준FP를 입력으로 하여, '대상FP귀속처리1'이 실행된다.
스탭 S2에서는 컴퓨터의 기준FP선정부(5)가 기능하여, 전 기준FP에 대하여 대상FP(17)와의 리텐션타임출현패턴의 일치도를 산출하고, 대상FP(17)의 귀속에 적합한 기준FP를 선정한다.
기준FP는 정상품으로 정해진 약제의 3D크로마토와 피크정보로부터 상기 스탭 S1과 마찬가지의 처리에 의해 작성된 FP이다. 또한, 정상품은 안전성, 유효성이 확인된 약제로 정의되고, 제품 로트가 상이한 복수의 약제가 해당한다. 기준FP도, 도 94의 FP데이타 예(127)와 마찬가지로 구성된 데이타이다.
스탭 S3는 대상FP(117)와 스탭 S2에서 선정한 기준FP를 입력으로 하여, '대상FP귀속처리2'가 실행된다.
스탭 S3에서는, 컴퓨터의 피크패턴작성부(7)(도 1) 및 피크귀속부(9)(도 1)가 기능한다. 이 기능에 의해 대상FP(17)와 스탭 S2에서 선정한 기준FP의 전 피크에서, 도 23 ~ 도 61과 같이 망라적으로 피크패턴을 작성하고, 이어서 그들 피크패턴의 일치도(도 63 또는 도 64의 P_Sim)를 산출한다. 또한, 대상FP와 기준FP의 피크 사이에서 UV스펙트럼의 일치도(도 66의 UV-Sim)를 산출한다. 더욱이, 이들 2개의 일치도로부터 귀속후보피크의 일치도(도 67의 SCORE)를 산출하고, 그 산출결과를 파일(판정결과파일)로 출력한다.
스탭 S4는 스탭 S3에서 출력한 판정결과파일을 입력으로 하여, '대상FP귀속처리3'가 실행된다.
스탭 S4에서는 컴퓨터의 피크귀속부(7)가 기능하여, 대상FP(17)와 기준FP 사이에서, 귀속후보피크의 일치도(SCORE)를 기초로 대상FP의 각 피크에 대응하는 기준FP의 피크를 특정하고, 그 결과를 파일(조합결과파일)로 출력한다.
스탭 S5는 스탭 S4에서 출력한 조합(照合)결과파일과 기준군FP를 입력으로 하여, '대상FP귀속처리4'가 실행된다.
기준군FP는 전 기준FP로부터 상기 스탭 S2로부터 스탭 S4와 마찬가지의 처리에 의해 작성된 전 기준FP 사이의 피크대응데이타이다.
스탭 S5에서는 컴퓨터의 피크귀속부(7)가 기능하여, 대상FP(17)의 조합결과파일을 기초로, 도 68, 도 69와 같이, 대상FP(17)의 각 피크를 기준군FP의 피크에 귀속하여, 그 결과를 파일(피크데이타특정량파일)로 출력한다.
스탭 S6는 스탭 S5에서 출력한 피크데이타특정량파일과 기준군FP를 입력으로 하여, 'FP평가처리'가 실행된다.
스탭 S6에서는 컴퓨터의 평가부(11)가 기능하여, 스탭 S5에서 출력한 피크데이타특정량데이타와 기준FP의 동등성을 MT법에 의해 평가하고, 그 평가결과를 MD값으로서 출력한다(도 70 ~ 도 74).
스탭 S7은 스탭 S6에서 출력한 MD값을 입력으로 하여, '합부(合否)판정'이 실행된다.
스탭 S7에서는 컴퓨터의 평가부(11)가 기능하여, 스탭 S6에서 출력한 MD값과 미리 설정한 역치(MD값의 상한치)를 비교하여, 합부를 판정한다. (도 2의 그래프 23).
[S1 : FP작성처리 (단일 파장만 이용)]
도 77은 도 76 스탭 S1 'FP 작성처리'의 단일파장의 피크정보를 이용한 경우의 플로우차트이다.
도 77은 단일파장, 예를 들면 203nm로 평가대상의 FP를 작성하는 스탭의 상세이다. 이 처리에서는 3D크로마토와 검출파장이 203nm에서의 피크정보로부터 203nm로 검출된 피크에서의 리텐션타임과 피크 및 그들 피크의 UV스펙트럼으로 구성하는 FP를 작성한다.
스탭 S101에서는, '피크정보를 판독'의 처리가 실행된다.
이 처리에서는 FP의 작성에 필요한 2개의 데이타 중의 첫번째로 피크정보가 판독되고, 스탭 S102로 이행한다.
스탭 S102에서는 '피크의 리텐션타임(R1)'과 대응하는 피크데이타(P1)를 순번으로 취득'의 처리가 실행된다. 이 처리에서는 피크정보로부터, 피크의 리텐션타임(R1) 및 피크데이타(P1)를 1 피크씩 순번으로 취득하고, 스탭 S103으로 이행한다.
스탭 S103에서는 '3D크로마토를 판독'의 처리가 실행된다. 이 처리에서는 FP의 작성에 필요한 2개의 데이타 중의 2번째로 3D크로마토가 판독되고, 스탭 S104로 이행한다.
스탭 S104에서는 '피크의 리텐션타임(R2)과 대응하는 UV스펙트럼(U1)을 순번으로 취득'의 처리가 실행된다. 이 처리에서는 3D크로마토로부터 리텐션타임(R2) 및 UV스펙트럼(U1)을 HPLC분석시의 샘플링레이트/2마다 취득하고, 스탭 S105로 이행한다.
스탭 S105에서는 '│R1 - R2│ < = 역치 ?'의 판단처리가 실행된다. 이 처리에서는 스탭 S102 및 S104에서 판독된 R1과 R2가 역치의 범위에서 대응한 것인가의 여부가 판단된다. 대응하고 있는(YES) 경우는 2개의 리텐션타임은 같고, 리텐션타임이 R1의 피크의 UV스펙트럼은 U1이라고 판단하고, 스탭 S106으로 이행한다. 대응하지 않는(NO) 경우는 2개의 리텐션타임은 동일하지 않고, 리텐션타임이 R1의 피크의 UV스펙트럼은 U1은 아니라고 판단하고, 3D크로마토의 다음의 데이타와의 비교를 위해, 스탭 S104로 이행한다. 또한, 이 판단처리에서의 역치는 3D크로마토에 있어서의 '샘플링레이트'로 한다.
스탭 S106에서는 'U1을 최대치 1로 규격화'의 처리가 실행된다. 이 처리에서는, S105에서 R1의 UV스펙트럼으로 판단한 U1을 최대치 1로 규격화하고, 스탭 S107로 이행한다.
스탭 S107에서는 'R1과 P1 및 규격화한 U1을 출력(대상FP)'의 처리가 실행된다. 이 처리에서는, 피크정보로부터 취득한 R1과 P1 및 S106에서 규격화한 U1을 대상FP로 출력하고, 스탭 S108로 이행한다.
스탭 S108에서는 '전 피크의 처리종료 ?'의 판단처리가 실행된다. 이 처리에서는 피크정보 중의 모든 피크에 대하여 처리가 행해졌는지 여부가 판단되고, 전 피크에 대하여 처리가 종료하지 않았다면(NO), 미처리의 피크를 처리하기 위해, 스탭 S102로 이행한다. S102로부터 S108까지의 처리는 전 피크의 처리가 종료하기 까지 반복되며, 전 피크의 처리가 종료하면(YES), FP작성처리를 종료한다.
[S1 : FP작성처리(복수파장이용)]
도 78, 도 79는 도 76 스탭 S1 'FP작성처리'에서, 상기 단일파장의 피크정보 대신에, 복수파장의 피크정보를 이용한 경우의 플로우차트이다. 예를 들면, 203nm를 포함시켜서, 검출파장축방향에 복수(n 개)의 파장을 선택하여, FP를 작성하는 경우이다.
이 FP작성처리는 도 77과 같은 단일파장에서는 3D크로마토에서 검출되어 있는 전 피크를 망라할 수 없는 경우에 복수파장의 피크정보를 이용하여, 3D크로마토의 전 피크를 망라한 FP를 작성하기 위한 것이다.
또한, 도 78, 도 79는 상기 단일파장만을 이용한 FP작성처리에서 파장마다의 FP를 n 개 작성 후, 그들 FP로부터 복수파장에 따른 FP를 작성하는 스탭의 상세이다.
스탭 S110에서는 '파장마다 FP작성'의 처리가 실행된다. 이 처리에서는, 파장마다 상기 단일 파장만을 이용한 FP작성처리가 행해져서, n 개의 FP를 작성하고, 스탭 S111로 이행한다.
스탭 S111에서는 'FP를 피크수(내림차순)으로 리스트화'의 처리가 실행된다. 이 처리에서는 n 개의 FP를 피크수가 많은 순으로 리스트화하고, 스탭 S112로 이행한다.
스탭 S112에서는 n 개의 FP를 순번으로 처리하기 위한 카운터의 초기화로서 n에 1이 대입되고(n ← 1), 스탭 S113으로 이행한다.
스탭 S113에서는 '리스트 n 번째의 FP를 판독'의 처리가 실행된다. 이 처리에서는 리스트 n 번째의 FP를 판독하고, 스탭 S114로 이행한다.
스탭 S114에서는 '전 리텐션타임(X)을 취득'의 처리가 실행된다. 이 처리에서는 S113에서 판독한 FP의 리텐션타임정보를 모두 취득하고, 스탭 S115로 이행한다.
스탭 S115에서는 'n의 경신(n ← n+1)'의 처리가 실행된다. 이 처리에서는 처리를 다음의 FP로의 이행하기 위해, n의 경신으로서 n에 n+1을 대입하고, 스탭 S116으로 이행한다.
스탭 S116에서는, '리스트 n 번째의 FP의 판독'의 처리가 실행된다. 이 처리에서는 리스트 n 번째의 FP를 판독하고, 스탭 S117로 이행한다.
스탭 S117에서는 '전 리텐션타임(Y)을 취득'의 처리가 실행된다. 이 처리에서는 S116에서 판독한 FP의 리텐션타임정보를 모두 취득하고, 스탭 S118로 이행한다.
스탭 S118에서는 'X와 Y를 중복 없이 통합(Z)'의 처리가 실행된다. 이 처리에서는 S114에서 취득한 리텐션타임정보(X)와 S117에서 취득한 리텐션타임정보(Y)를 중복 없이 통합한 후, Z에 보존하고, 스탭 S119로 이행한다.
스탭 S119에서는 'X의 경신(X ← Z)'의 처리가 실행된다. 이 처리에서는 X의 경신으로서 X에 S118에서 보존한 Z를 대입하고, 스탭 S120으로 이행한다.
스탭 S120에서는 '전 FP처리종료 ?'의 판단처리가 실행된다. 이 처리에서는 S110에서 작성한 n 개의 FP 전체가 처리되었는지의 여부가 판단되고, 처리완료(YES)의 경우는 스탭 S121로 이행한다. 미처리의 FP가 있는(NO)의 경우는 미처리의 FP에 대하여 S115 ~ S120의 처리를 실행하기 위해, S115로 이행한다. 전 FP의 처리가 종료하기까지 S115 ~ S120의 처리를 반복한다.
스탭 S121에서는 n 개의 FP를 재차 순번으로 처리하기 위한 카운터의 초기화로서 n에 1이 대입(n ← 1)되고, 스탭 S122로 이행한다.
스탭 S122에서는 '리스트 n 번째의 FP를 판독'의 처리가 실행된다. 이 처리에서는 리스트 n 번째의 FP를 판독하고, 스탭 S123으로 이행한다.
스탭 S123에서는 '각 피크의 리텐션타임(R1), 피크데이타(P1) 및 UV스펙트럼(U1)을 순번으로 취득'의 처리가 실행된다. 이 처리에서는, S122에서 판독한 FP로부터 리텐션타임(R1), 피크데이타(P1) 및 UV스펙트럼(U1)을 1 피크씩 순번으로 취득하고, 스탭 S124로 이행한다.
스탭 S124에서는, 'X로부터 리텐션타임(R2)를 순번으로 취득'의 처리가 실행된다. 이 처리에서는 전 FP의 리텐션타임이 중복 없이 보존되어 있는 X로부터 1 리텐션타임(R2)씩 순번으로 취득하고, 스탭 S125로 이행한다.
스탭 S125에서는, 'R1 = R2 ?'의 판단처리가 실행된다. 이 처리에서는 S123에서 취득한 R1과 S124에서 취득한 R2가 같은가의 여부가 판단되고, 같은(YES) 경우는 스탭 S127로 이행한다. 같지 않은(NO) 경우는 스탭 S126으로 이행한다.
스탭 S126에서는 'X의 전 리텐션타임 비교종료 ?'의 판단처리가 실행된다. 이 처리에서는, S123에서 취득한 R1에 대하여, X의 전 리텐션타임과의 비교가 종료하고 있는지의 여부가 판단된다. 종료하고 있는(YES) 경우는, 리텐션타임이 R1의 피크는 처리완료라고 판단하고, 다음의 피크로 이행하기 위해, 스탭 S123으로 이행한다. 종료하지 않은(NO)의 경우는 X의 다음의 리텐션타임으로 이행하기 위해, 스탭 S124로 이행한다.
스탭 S127에서는 'R1에 (n - 1) × 분석시간(T)을 가산 (R1 ← R1 + (n - 1)×T )'의 처리가 실행된다. 이 처리에서는 가장 피크수가 많은 리스트 1번째의 FP에 존재하는 피크의 리텐션타임은 그대로이고, 리스트 1번째의 FP에는 존재하지 않고, 리스트 2번째의 FP에 존재하는 피크의 리텐션타임은 R1에 분석시간(T)이 가산되며, 더욱이, 리스트 1 ~ n-1 번째의 FP에는 존재하지 않고, 리스트 n 번째의 FP에 존재하는 피크의 리텐션타임은 R1에 (n - 1)×T 가 가산되고, 스탭 S128로 이행한다.
스탭 S128에서는 'R1, P1 및 U1을 출력 (대상FP)'의 처리가 실행된다. 이 처리에서는 S127에서 처리한 R1과 S123에서 취득한 P1 및 U1을 대상FP로 출력하고, 스탭 S129로 이행한다.
스탭 S129에서는 'X로부터 R2를 삭제'의 처리가 실행된다. 이 처리에서는 리텐션타임이 R1 (= R2)에 있어서의 처리가 S127, S128로 종료했기 때문에, X로부터 처리종료 리텐션타임(R2)을 삭제하고, S130으로 이행한다.
스탭 S130에서는 '전 피크 처리종료?'의 판단처리가 실행된다. 이 처리에서는 리스트 n 번째의 FP의 전 피크에 대하여 처리가 종료하였는지의 여부가 판단되고, 처리종료(YES)의 경우는 리스트 n 번째의 FP에서의 FP작성처리를 종료하고, 스탭 S131로 이행한다. 미처리의 피크가 있는(NO) 경우는 미처리의 피크를 처리하기 위해, 스탭 S123으로 이행한다. 전 피크의 처리가 종료하기까지 S123 ~ S130의 처리를 반복한다.
스탭 S131에서는 'n의 경신 (n ← n + 1)'의 처리가 실행된다. 이 처리에서는 처리를 다음의 FP로 이행하기 위해, n의 경신으로서 n에 n+1을 대입하고, 스탭 S132로 이행한다.
스탭 S132에서는 '전 FP처리종료?'의 판단처리가 실행된다. 이 처리에서는, S110에서 작성한 n개의 FP 모두가 처리되었는지의 여부가 판단되고, 처리완료(YES)의 경우는 FP작성처리를 종료한다. 미처리의 FP가 있는(NO) 경우는 미처리의 FP에 대하여 S122 ~ S132의 처리를 실행하기 위해, S122로 이행한다. 전 FP의 처리가 종료하기까지 S122 ~ S132의 처리를 반복한다.
[S2 : 대상FP귀속처리1]
도 80은 도 76 스탭 S2의 '대상FP귀속처리1'의 상세를 도시하는 폴로우차트이다. 이 처리는 귀속의 전(前)처리이고, 정상품으로 된 복수의 기준FP로부터 대상FP(17)의 귀속에 적합한 기준FP를 선정한다.
스탭 S201에서는 '대상FP를 판독'의 처리가 실행된다. 이 처리에서는, 귀속대상의 FP를 판독하고, 스탭 S202로 이행한다.
스탭 S202에서는 '전 리텐션타임(R1)을 취득'의 처리가 실행된다. 이 처리에서는 S201에서 판독한 대상FP의 리텐션타임정보를 모두 취득하고, 스탭 S203으로 이행한다.
스탭 S203에서는 '전 기준FP의 파일명을 리스트화'의 처리가 실행된다. 이 처리에서는 후에 전 기준FP를 순번으로 처리하기 위해 미리 전 기준FP의 파일명을 리스트화하고, 스탭 S204로 이행한다.
스탭 S204에서는 전 기준FP를 순번으로 처리하기 위한 카운터의 초기치로서 n에 1을 대입 (n ← 1)로 하고, 스탭 S205로 이행한다.
스탭 S205에서는, '리스트 n 번째의 기준FP(기준FPN)를 판독'의 처리가 실행된다. 이 처리에서는 S203에서 리스트화한 전 기준FP의 파일명 리스트의 n 번째 FP를 판독하고, 스탭 S206으로 이행한다.
스탭 S206에서는, '전 리텐션타임(R2)을 취득'의 처리가 실행된다. 이 처리에서는 S205에서 판독한 기준FP의 리텐션타임정보를 모두 취득하고, 스탭 S207로 이행한다.
스탭 S207에서는, 'R1과 R2의 리텐션타임출현패턴의 일치도를 산출(RPn_min)'의 처리가 실행된다. 이 처리에서는 S202에서 취득한 대상FP의 리텐션타임과 S206에서 취득한 기준FP의 리텐션타임으로부터 RPn_min을 산출하고, 스탭 S208로 이행한다. 또한, RPn_min의 상세한 계산플로우는 도 85의 섭루틴1에 의해 별로로 설명한다.
스탭 S208에서는, 'RPn_min의 보존(RPall_min)'의 처리가 실행된다. 이 처리에서는 S207에서 산출된 RPn_min을 RPall_min으로 보존하고, 스탭 S209로 이행한다.
스탭 S209에서는 'n의 경신 (n ← n+1)'의 처리가 실행된다. 이 처리에서는 처리를 다음의 FP로 이행하기 위해 n의 경신으로서 n에 n+1을 대입하고, 스탭 S210으로 이행한다.
스탭 S210에서는 '전 기준FP 처리종료?'의 판단처리가 실행된다. 이 처리에서는, 기준FP 모두가 처리되었는지의 여부가 판단되고, 처리완료(YES)의 경우는 스탭 S211로 이행한다. 미처리의 기준FP가 있는(NO) 경우는 미처리의 FP에 대하여 S205 ~ S210의 처리를 실행하기 위해, S205로 이행한다. 전 기준FP의 처리가 종료하기까지 S205 ~ S210의 처리를 반복한다.
스탭 S211에서는, 'RPALL_min으로부터 일치도가 최소로 되는 기준FP를 선정'의 처리가 실행된다. 이 처리에서는 전 기준FP에 대하여 산출한 RP1_min으로부터 RPn_min을 비교하여, 대상FP와의 리텐션타임출현패턴의 일치도가 최소로 되는 기준FP를 선정하고, 대상FP귀속처리1을 종료한다.
[S3 : 대상FP귀속처리2]
도 81은 도 76 스탭 S3의 '대상FP귀속처리2'의 상세를 도시하는 폴로우차트이다. 이 처리는 귀속의 본처리이고, 대상FP(17)와 스탭 S2에서 선정한 기준FP와의 사이에서, 상기와 같이 피크패턴 및 UV스펙트럼의 일치도로부터 각 귀속후보피크의 일치도(SCORE)를 산출한다.
스탭 S301에서는 '대상FP를 판독'의 처리가 실행된다. 이 처리에서는, 귀속대상의 FP를 판독하고, 스탭 S302로 이행한다.
스탭 S302에서는 '귀속대상피크의 리텐션타임(R1)과 피크데이타(P1) 및 UV스펙트럼(U1)을 순번으로 취득'의 처리가 실행된다. 이 처리에서는 S301에서 판독한 대상FP의 각 피크를 순번으로 귀속대상피크로 하여, R1과 P1 및 U1을 취득하고, 스탭 S303으로 이행한다.
스탭 S303에서는 '기준FP를 판독'의 처리가 실행된다. 이 처리에서는 도 80의 '대상FP귀속처리1'에서 선정된 기준FP를 판독하고, 스탭 S304로 이행한다.
스탭 S304에서는 '기준FP의 피크의 리텐션타임(R2)과 피크데이타(P2) 및 UV스펙트럼(U2)을 순번으로 취득'의 처리가 실행된다. 이 처리에서는 S303에서 판독한 기준FP로부터 R2와 P2 및 U2를 1피크씩 취득하고, 스탭 S305으로 이행한다.
스탭 S305에서는 '│R1 - (R2 + d)│ < 역치?'의 판단처리가 실행된다. 이 처리에서는 S302 및 S304에서 판독된 R1과 R2가 역치의 범위 내에 대응한 것인지의 여부가 판단된다. 대응하고 있는(YES) 경우는 리텐션타임이 R2의 피크는 리텐션타임이 R1의 피크의 귀속후보피크인 것으로 판단하고, 귀속후보피크의 일치도(SCORE)를 산출하기 위해 스탭 S306으로 이행한다. 대응하지 않는(NO) 경우는 리텐션타임이 R2의 피크와 리텐션타임이 R1의 피크에서는 리텐션타임이 너무 상이하기 때문에 귀속후보피크가 될 수 없다고 판단하고, 스탭 S309로 이행한다. 또한, 이 판단처리에서의 d는 대상FP와 기준FP의 피크의 리텐션타임을 보정하는 값으로, 초기치는 0으로 하고, 처리를 진행하는 중에 수시(隨時) 귀속한 피크 사이의 리텐션타임의 차를 구하여, 이값으로 d를 경신한다. 또한, 역치는 귀속후보피크로 해야하는지의 여부를 판단하기 위한 리텐션타임의 허용폭이다.
스탭 S306에서는 'UV스펙트럼의 일치도를 산출(UV_Sim)'의 처리가 실행된다. 이 처리에서는, S302에서 취득한 귀속대상피크의 U1과 S304에서 취득한 귀속후보피크의 U2로부터 UV_Sim을 산출하고, 스탭 S307로 이행한다. 또한, UV_Sim의 상세한 계산플로우는 도 86의 섭루틴2에서 별도로 기재한다.
스탭 S307에서는 '피크패턴의 일치도를 산출(P_Sim_min)'의 처리가 실행된다. 이 처리에서는 S302에서 취득한 귀속대상피크의 R1 및 P1과 S304에서 취득한 귀속대상피크의 R2 및 P2로부터 이들 피크에 대하여 망라적으로 피크패턴을 작성한다. 이들 피크패턴의 P_Sim_min을 산출하고, 스탭 S308로 이행한다. 또한, P_Sim_min의 상세한 계산플로우는 도 87의 섭루틴3에서 별도로 기재한다.
스탭 S308에서는 '귀속후보피크의 일치도를 산출(SOCORE)'의 처리가 실행된다. 이 처리에서는 S306에서 산출한 UV_Sim과 S307에서 산출한 P_Sim_min으로부터, 귀속대상피크와 귀속후보피크의 SCORE를
SCORE = UV_Sim × P_Sim_min
로 산출하고, 스탭 S310으로 이행한다.
스탭 S309에서는 'SCORE에 888888를 대입(SCORE ← 888888)'의 처리가 실행된다. 이 처리에서는 귀속대상피크의 귀속후보피크에 해당하지 않는 피크의 SCORE를 888888로 하고, 스탭 S310으로 이행한다.
스탭 S310에서는, 'SCORE의 보존(SCORE_all)'의 처리가 실행된다. 이 처리에서는 S308 혹은 S309에서 얻어진 SCOPRE을 SCORE_all로 보존하고, 스탭 S311로 이행한다.
스탭 S311에서는 '기준 전 피크의 처리종료?'의 판단처리가 실행된다. 이 처리에서는, 기준FP의 전 피크가 처리되었는지의 여부가 판단되고, 처리완료(YES)의 경우는 스탭 S312로 이행한다. 미처리의 피크가 있는(NO) 경우는 미처리의 피크에 대하여 S304 ~ S311의 처리를 실행하기 위해, S304로 이행한다. 전 피크의 처리가 종료하기까지 S304 ~ S311의 처리를 반복한다.
스탭 S312에서는, '판정결과파일에 SCORE_all을 출력하고, SCORE_all을 초기화(0으로 함)'의 처리가 실행된다. 이 처리에서는 판정결과파일에 SCORE_all을 출력한 후, SCORE_all을 초기화(영으로 함)하고, 스탭 S311로 이행한다.
스탭 S313에서는 '대상 전 피크의 처리종료?'의 판단처리가 실행된다. 이 처리에서는, 대상FP의 전 피크가 처리되었는지의 여부가 판단되고, 처리완료(YES)의 경우는 대상FP처리2를 종료한다. 미처리의 피크가 있는(NO) 경우는 미처리의 피크에 대하여 S302 ~ S313의 처리를 실행하기 위해, S302로 이행한다. 전 피크의 처리가 종료하기까지 S302 ~ S313의 처리를 반복한다.
도 95에 출력한 판정결과파일 예(129)를 도시한다.
[S4 : 대상FP귀속처리3]
도 82는 도 76 스탭 S4의 '대상FP귀속처리3'의 상세를 도시하는 플로우차트이다. 이 처리는 귀속의 후처리이고, 상기와 같이 산출한 귀속후보피크의 일치도(SCORE)로부터 대상FP의 각 피크에 대응하는 기준FP의 피크를 특정한다.
스탭 S401에서는, '판정결과파일을 판독'의 처리가 실행된다. 이 처리에서는 도 81의 '대상FP귀속처리2'에서 작성한 판정결과파일을 판독하고, 스탭 S402로 이행한다.
스탭 S402에서는, ''SCORE < 역치'의 조건을 만족한 데이타로 귀속후보피크스코어표를 작성'의 처리가 실행된다. 이 처리에서는, 판정결과파일의 SCORE를 기초로 도 96(윗 도면)의 귀속후보스코어표(131)를 작성하고, 스탭 S403으로 이행한다. 이 귀속후보피크스코어표는 기준FP의 피크마다에, 대상FP 전 피크에 대하여 산출된 SCORE로부터 역치보다 작은 SCORE만을 올림차순으로 나열한 표이다. 즉, 이 SCORE는 값이 작을수록 귀속해야할 피크의 가능성이 높다. 또한, 역치는 귀속후보로 해야할지의 여부를 판단하기 위한 SCORE의 상한치이다.
스탭 S403에서는, '귀속후보피크번호표를 작성'의 처리가 실행된다. 이 처리에서는, 귀속후보피크스코어표를 기초로 도 96(아래 도면)의 귀속후보번호표(133)를 작성하고, 스탭 S404로 이행한다. 이 귀속후보번호표는 귀속후보피크스코어표의 각 스코어를 그 스코어에 대응하는 대상FP의 피크번호에 치환한 표이다. 이 때문에 이 표는 기준FP의 피크마다에 대응시켜야할 대상FP의 피크번호가 순번으로 나열된 표로 되어 있다.
스탭 S404에서는, '귀속해야할 대상FP의 피크번호를 취득'의 처리가 실행된다. 이 처리에서는, S403에서 작성한 귀속후보피크번호표로부터 기준FP의 피크마다에 가장 상위에 위치하는 대상FP의 피크번호를 취득하고, 스탭 S405로 이행한다.
스탭 S405에서는, '취득한 피크번호가 내림차순(중복 없이)으로 나열되어 있는가?'의 판단처리가 실행된다. 이 처리에서는, S404에서 취득한 대상FP의 피크번호가 중복 없이 내림차순으로 나열되어 있는지의 여부가 판단된다. 나열되어 있는(YES) 경우는 기준FP의 각 피크에 대응하는 대상FP의 피크가 확정되었다고 판단하고, 스탭 S408로 이행한다. 나열되어 있지 않은(NO) 경우는 문제가 있었던 기준FP의 피크에 귀속해야할 대상FP의 피크를 수정하기 위해, 스탭 S406으로 이행한다.
스탭 S406에서는, '문제가 있었던 피크 사이에서 SCORE를 비교하고, 귀속후보피크번호표를 경신'의 처리가 실행된다. 이 처리에서는, 문제가 있었던 대상FP의 피크번호에 대응하는 SCORE를 귀속후보스코어표에서 비교하고, SCORE가 큰 쪽의 피크번호를 2번째에 위치하는 피크번호에 치환한 후보피크번호표로 경신하고, 스탭 S407로 이행한다.
스탭 S407에서는, '귀속후보피크스코어표를 경신'의 처리가 실행된다. 이 처리에서는, S406에서의 귀속후보피크번호표의 경신내용에 따라서, 귀속후보피크스코어표를 경신하고, 스탭 S408로 이행한다. 대상FP의 피크번호에 문제(중복 있음, 내림차순으로 나열되어 있지 않음)가 없어질 때까지 S404로부터 S407의 처리를 반복한다.
스탭 S408에서는, '귀속결과를 보존(TEMP)'의 처리가 실행된다. 이 처리에서는, 기준FP의 전 피크의 피크번호, 리텐션타임, 및 피크와 이들 피크에 대응하는 피크로서 특정한 대상FP의 피크데이타를 TEMP에 보존하고, 스탭 S409로 이행한다.
스탭 S409에서는, 'TEMP에 대상FP의 전 피크가 들어 있는가?'의 판단처리가 실행된다. 이 처리에서는, S408에서 보존한 TEMP 중에 대상FP의 전 피크의 피크데이타가 들어 있는지의 여부가 판단된다. 전부 들어 있는 경우(YES)는, 대상FP의 모든 피크에서 처리가 종료한 것으로 판단하고, S412로 이행한다. 들어 있지 않은 피크가 있는(NO) 경우는 들어 있지 않은 피크의 피크데이타를 TEMP에 추가하기 위해, 스탭 S410으로 이행한다. 스탭 S410에서는 'TEMP에 들어 있지 않은 대상FP의 피크의 리텐션타임을 보정'의 처리가 실행된다. 이 처리에서는, TEMP에 들어가 있지 않은 대상FP의 피크(보정이 필요한 대상FP의 피크)의 리텐션타임은
보정치 = k1 + (k2 - k1) * (t0 - t1) / (t2 - t1)
k1 : 보정이 필요한 대상FP의 피크근방에서 귀속된 2개의 기준FP측의 피크 중의 리텐션타임이 작은 피크의 리텐션타임
k2 : 보정이 필요한 대상FP의 피크근방에서 귀속된 2개의 기준FP측의 피크 중의 리텐션타임이 큰 피크의 리텐션타임
t0 : 보정이 필요한 대상FP의 피크의 리텐션타임
t1 : 보정이 필요한 대상FP의 피크 근방에서 귀속된 2개의 대상FP측의 피크 중의 리텐션타임이 작은 피크의 리텐션타임
t2 : 보정이 필요한 대상FP의 피크 근방에서 귀속된 2개의 대상FP측의 피크 중의 리텐션타임이 큰 피크의 리텐션타임
로 하여 기준FP에 있어서의 리텐션타임으로 보정하고, 스탭 S411로 이행한다.
스탭 S411에서는, '보정한 리텐션타임과 그 피크의 피크데이타를 TEMP에 추가하고, TEMP를 경신'의 처리가 실행된다. 이 처리에서는, S410에서 보정한 TEMP에 들어 있지 않은 대상FP의 피크의 리텐션타임을 TEMP 중의 기준FP의 리텐션타임과 비교하여, TEMP 중의 타당한 위치에 TEMP에 들어 있지 않은 대상FP의 피크가 보정한 리텐션타임 및 피크데이타를 추가하여, TEMP를 경신하고, S409로 이행한다. 대상FP의 전 피크가 추가되기까지, S409로부터 S411의 처리를 반복한다.
스탭 S412에서는 'TEMP를 조합결과파일에 출력'의 처리가 실행된다. 이 처리에서는, 기준FP의 전 피크와 대상FP의 전 피크에 대응관계를 특정한 TEMP를 조합결과파일로서 출력하고, 대상FP귀속처리3을 종료한다.
도 97에 상기와 같이 출력한 조합결과파일 예(135)를 도시한다.
[S4 : 대상FP귀속처리4]
도 83, 도 84는 도 76 스탭 S5의 '대상FP귀속처리4'의 상세를 도시하는 플로우차트이다. 이 처리는 귀속의 최종처리이고, 도 76 스탭 S4에서 작성한 조합결과파일을 기초로 대상FP의 각 피크를 기준군FP의 피크에 귀속한다. 또한, 기준군FP는 상기와 같이 전 기준FP 사이에서 피크의 대응관계를 특정한 FP이고, 도 98의 기준군FP데이타 예(137)와 같이, 기준군FP피크번호와 기준군리텐션타임과 피크높이로 구성된 데이타이다. 도 2의 기준군FP(19)에서 표시한 바와 같이, 각 피크는 평균값(흑점) ± 표준편차(종선)로 표시할 수 있다.
스탭 S501에서는, '조합결과파일을 판독'의 처리가 실행된다. 이 처리에서는 도 82의 S412에서 출력한 조합결과파일을 판독하고, 스탭 S502로 이행한다.
스탭 S502에서는, '기준군FP를 판독'의 처리가 실행된다. 이 처리에서는 대상FP의 각 피크의 최종적인 귀속상대인 기준군FP를 판독하고, 스탭 S503으로 이행한다.
스탭 S503에서는, '대상FP와 기준군FP를 통합하여 보존(TEMP)'의 처리가 실행된다. 이 처리에서는, 조합결과파일과 기준군FP에서 공통으로 존재하는 기준FP의 피크데이타를 기초로 2개의 파일을 통합하여, 그 결과를 TEMP로서 보존하고, 스탭 S504로 이행한다.
스탭 S504에서는, '기준FP에 대응하는 피크가 없는 대상FP의 피크의 리텐션타임을 보정'의 처리가 실행된다. 이 처리에서는, 조합결과파일에서 기준FP에 대응하는 피크가 없는 대상FP의 전 피크의 리텐션타임을 S503에서 보존한 TEMP의 리텐션타임으로 보정하고, 스탭 S505 이행한다. 또한, 리텐션타임의 보정은 상기 스탭 S4의 '대상FP귀속처리3'의 스탭 S410과 마찬가지의 방법으로 보정한다.
스탭 S505에서는, '보정한 리텬션타임(R1, R3)과 대응하는 피크데이타(P1)를 순번으로 취득'의 처리가 실행된다. 이 처리에서는, S504에서 보정한 리텐션타임을 R1 및 R3으로 하여, 대응하는 피크의 피크데이타를 P1으로서 순번으로 취득하고, 스탭 S507로 이행한다
스탭 S506에서는'TEMP로부터 대상FP의 귀속후보피크의 리텐션타임(R2)과 대응하는 피크데이타(P2)를 순번으로 취득'의 처리가 실행된다. 이 처리에서는, S503에서 보존한 TEMP로부터 대상FP가 귀속되어 있지 않은 리텐션타임을 R2로 하고, 대응하는 피크데이타를 P2로 하여 순번으로 취득하고, 스탭 S507로 이행한다.
스탭 S507에서는 '│R1 - R2│< 역치1 ?'의 판단처리가 실행된다. 이 처리에서는 S505 및 S506에서 취득한 R1과 R2의 차가 역치1보다 작은지의 여부가 판단된다. 작은(YES) 경우는 대상FP의 리텐션타임이 R1의 피크와 기준FP의 리텐션타임이 R2의 피크가 대응하는 가능성이 있는 것으로 판단하고, 스탭 S508로 이행한다. R1과 R2의 차가 역치1 이상(NO)의 경우는 대응할 가능성이 없는 것으로 판단하고, 스탭 S512로 이행한다.
스탭 S508에서는 'R1, R2에 대응하는 UV스펙트럼의 취득(U1, U2)'의 처리가 실행된다. 이 처리에서는 S507에서 대응할 가능성이 있다고 판단된 리텐션타임이 R1과 R2의 피크에 대응하는 UV스펙트럼을 각각의 FP로부터 취득하고, 스탭 S509로 이행한다.
스탭 S509에서는 'UV스펙트럼의 일치도를 산출(UV_Sim)'의 처리가 실행된다. 이 처리에서는, S508에서 취득한 UV스펙트럼 U1 및 UV2로부터 스탭 S3의 '대상FP귀속처리2'의 스탭 S306과 마찬가지의 방법으로 UV_Sim을 산출하고, 스탭 S510으로 이행한다. 또한, UV_Sim의 상세한 계산플로우는 도 86의 섭루틴2에서 별도로 설명한다.
스탭 S510에서는 'UV_Sim < 역치2 ?'의 판단처리가 실행된다. 이 처리에서는 S509에서 산출한 UV_Sim이 역치2 보다 작은지의 여부가 판단된다. 작은(YES) 경우는 UV스펙트럼이 U1의 피크와 U2의 피크가 대응하고 있는 것으로 판단하고, 스탭 S511로 이행한다. UV_Sim이 역치2 이상(NO)인 경우는 대응하지 않는 것으로 판단하고, 스탭 S507로 이행한다.
스탭 S511에서는 'R3 ← R2, 역치2 ← UV_Sim'의 처리가 실행된다. 이 처리에서는, S510에서 대응한다고 판단한 리텐션타임이 R3(즉, R1)을 대응상대의 리텐션타임인 R2로 경신한 후, 역치2를 UV_Sim의 값으로 경신하고, S507로 이행한다.
스탭 S512에서는 '모든 귀속후보피크의 리텐션타임을 비교종료?'의 판단처리가 실행된다. 이 처리에서는, R1과 모든 귀속후보피크의 리텐션타임의 비교가 종료했는지의 여부가 판단되고, 종료한(YES)의 경우는 스탭 S513으로 이행한다. 종료하지 않은(NO) 경우는 스탭 S507로 이행한다.
스탭 S513에서는 'R1, R3과 P1 및 역치2를 보존(TEMP2)'의 처리가 실행된다. 이 처리에서는, S510에서 대응한다고 판단한 리텐션타임(R1)과 대응상대의 리텐션타임(R2)으로 경신한 R3와 대응하는 피크(P1) 및 현시점의 역치2를 보존(TEMP2)하고, S507로 이행한다.
스탭 S514에서는 '모든 비대응피크의 리텐션타임을 비교종료?'의 판단처리가 실행된다. 이 처리에서는, 모든 비대응피크의 리텐션타임에서 귀속후보피크의 리텐션타임과의 비교가 종료했는지의 여부가 판단된다. 종료한(YES) 경우는 모든 비대응피크의 귀속처리가 종료한 것으로 판단하고, S516로 이행한다. 종료하지 않은(NO) 경우는 미처리의 비대응피크가 남아 있는 것으로 판단하고, 스탭 S515로 이행한다.
스탭 S515에서는 '역치2 ← 초기치'의 처리가 실행된다. 이 처리에서는, S511에서 UV_Sim에 경신되어 있는 역치2를 초기치로 복귀시키고, S505로 이행한다.
스탭 S516에서는 'TEMP2에 R3의 값이 같은 피크가 존재하는가?'의 판단처리가 실행된다. 이 처리에서는, TEMP 중의 동일한 피크에 복수의 비대응피크가 귀속되어 있는지의 여부가 판단된다. 같은 피크에 귀속된 비대응피크가 존재하는(YES) 경우는 스탭 S517로 이행한다. 존재하지 않는(NO) 경우는 스탭 S518로 이행한다.
스탭 S517에서는 'R3의 값이 같은 피크의 역치2를 비교하고, 값이 큰 피크의 R3를 원래의 값(R1)으로 되돌림'의 처리가 실행된다. 이 처리에서는 TEMP2 중의 R3의 값이 동일한 피크의 역치2를 비교하고, 값이 큰 피크의 R3의 값을 원래의 값(즉, R1)으로 복귀시키고, 스탭 S518로 이행한다.
스탭 S518에서는 'TEMP에 TEMP2의 피크를 추가(TEMP의 리텐션타임과 R3가 일치한 피크만)'의 처리가 실행된다. 이 처리에서는, TEMP의 리텐션타임과 R3가 일치한 피크만 TEMP에 R3에 대응하는 피크를 추가하고, 스탭 S519로 이행한다. R3가 TEMP의 리텐션타임과 일치하지 않는 피크는 기준군FP에 귀속상대가 되는 피크가 존재하지 않기 때문에 추가하지 않는다.
스탭 S519에서는 'TEMP 중의 대응FP의 피크를 출력(피크특징량파일)'의 처리가 실행된다. 이 처리에서는 기준군FP(137)에 귀속된 대상FP의 피크데이타를 피크데이타특징량파일로서 출력하고, 대상FP귀속처리4를 종료한다.
도 99에 상기와 같이 출력하는 피크데이타특징량의 파일 예(139)를 도시한다.
[섭루틴 1]
도 85는 도 80의 '기준FP선정처리'에 있어서의 '섭루틴1'의 상세를 도시하는 플로우차트이다. 이 처리는 FP간(예를 들면, 대상FP와 기준FP)의 리텐션타임출현패턴의 일치도를 계산한다.
스탭 S1001에서는 'x ← R1, y ← R2'의 처리가 실행된다. 이 처리에서는 도 80의 S202와 S206에서 취득한 R1 및 R2를 각각 x와 y에 대입하고, 스탭 S1002로 이행한다.
스탭 S1002에서는 'x, y의 데이타수를 취득(a,b)'의 처리가 실행된다. 이 처리에서는 x, y의 데이타수를 각각 a, b로서 취득하고, 스탭 S1003으로 이행한다.
스탭 S1003에서는 x의 리텐션타임을 순번으로 호출하기 위한 카운터의 초기치로서 i에 1을 대입 (i ← 1)하고, 스탭 S1004로 이행한다.
스탭 S1004에서는 'xi번째의 리텐션타임으로부터 전 거리를 취득(f)'의 처리가 실행된다. 이 처리에서는 xi번째의 리텐션타임과 그 이후의 전 리텐션타임간의 거리를 f로서 취득하고, 스탭 S1005로 이행한다.
스탭 S1005에서는 y의 리텐션타임을 순번으로 호출하기 위한 카운터의 초기치로서 j에 1을 대입(j ← 1)하고, 스탭 S1006으로 이행한다.
스탭 S1006에서는 'yj번째의 리텐션타임으로부터 전 거리를 취득(g)'의 처리가 실행된다. 이 처리에서는 yj번째의 리텐션타임과 그 이후의 전 리텐션타임간의 거리를 g로서 취득하고, 스탭 S1007로 이행한다.
스탭 S1007에서는 '│f의 각 리텐션타임 사이 거리 - g의 각 리텐션타임 사이 거리│< 역치'의 조건을 만족한 데이타수를 취득(m)]의 처리가 실행된다. 이 처리에서는 S1004 및 S1006에서 취득한 리텐션타임간 거리 f와 g를 전체개별대비방식으로 비교하여, '│f의 각 리텐션타임간 거리 - g의 각 리텐션타임간 거리│< 역치'의 조건을 만족한 데이타수를 m으로 취득하고, 스탭 S1008로 이행한다.
스탭 S1008에서는, 'f와 g의 리텐션타임출현패턴의 일치도를 산출(RPfg)'의 처리가 실행된다. 이 처리에서는, S1002에서 취득한 a, b와 S1007에서 취득한 m으로부터 RPfg를
RPfg = ( 1 - ( m / ( a + b - m ) ) ) × ( a - m + 1 )
로 산출하고, 스탭 S1009로 이행한다.
스탭 S1009에서는, 'RPfg를 보존(RP_all)'의 처리가 실행된다. 이 처리에서는, S1008에서 산출한 일치도를 RP_all로 보존하고, 스탭 S1010으로 이행한다.
스탭 S1010에서는 'j의 경신(j ← j + 1)'의 처리가 실행된다.이 처리에서는, y의 처리를 다음의 리텐션타임으로 이행하기 위한 j의 경신으로서 j에 j+1을 대입하고, 스탭 S1011로 이행한다.
스탭 S1011에서는 'y의 전 리텐션타임에서 처리종료?'의 판단처리가 실행된다. 이 처리에서는, y의 전 리텐션타임의 처리가 종료했는지의 여부가 판단된다. 종료한(YES) 경우는 y의 전 리텐션타임의 처리가 종료한 것으로 판단하고, S1012로 이행한다. 종료하지 않은(NO) 경우는 y 중에 미처리의 리텐션타임이 남아 있는 것으로 판단하고, 스탭 S1006으로 이행한다. 즉, S1006 ~ S1011의 처리는 y의 전 리텐션타임이 처리되기까지 반복한다.
스탭 S1012에서는 'i의 경신(i ← i + 1)'의 처리가 실행된다.이 처리에서는, x의 처리를 다음의 리텐션타임으로 이행하기 위한 i의 경신으로서 i에 i+1을 대입하고, 스탭 S1013으로 이행한다.
스탭 S1013에서는 'x의 전 리텐션타임에서 처리종료?'의 판단처리가 실행된다. 이 처리에서는, x의 모든 리텐션타임의 처리가 종료했는지의 여부가 판단된다. 종료한(YES) 경우는 x의 전 리텐션타임의 처리가 종료했다고 판단하고, 스탭 S1014로 이행한다. 종료하지 않은(NO) 경우는 x 중에 미처리의 리텐션타임이 남아 있는 것으로 판단하고, 스탭 S1004으로 이행한다. 즉, S1004 ~ S1013의 처리는 x의 모든 리텐션타임이 처리되기까지 반복한다.
스탭 S1014에서는 'RP_all로부터 최소치를 취득(RP_min)'의 처리가 실행된다. 이 처리에서는 대상FP와 기준FP의 리텐션타임출현패턴의 전체 조합의 RP가 보존된 RP_all 중의 최소치를 RP_min로서 취득하고, 그 RP_min을 도 80의 S207로 넘기고, 리텐션타임출현패턴의 일치도계산처리를 종료한다.
[섭루틴2]
도 86은 도 81의 '대상FP귀속처리2'에 있어서의 '섭루틴2'의 상세를 도시하는 플로우차트이다. 이 처리는 UV스펙트럼의 일치도를 계산한다.
스탭 S2001에서는 'x ← U1, y ← U2, z ← 0'의 처리가 실행된다. 이 처리에서는 도 81의 S302와 S304에서 취득한 UV스펙트럼 U1 및 U2를 각각 x와 y에 대입하고, 더욱이 UV스펙트럼간의 거리의 제곱합(z)의 초기치로서 0을 대입하고, 스탭 S2002로 이행한다.
스탭 S2002에서는 'x 데이타수를 취득(a)'의 처리가 실행된다. 이 처리에서는 x의 데이타수를 각각 a로서 취득하고, 스탭 S2003로 이행한다.
스탭 S2003에서는 'i ← 1'의 처리가 실행된다. 이 처리에서는, UV스펙트럼 U1 및 U2를 구성하는 각 검출파장에 있어서의 흡광도를 x 및 y로부터 순번으로 호출하기 위한 초기치로서 i에 1을 대입하고, 스탭 S2004로 이행한다.
스탭 S2004에서는 'xi번째의 데이타를 취득(b)'의 처리가 실행된다. 이 처리에서는 UV스펙트럼 U1을 대입한 x의 i번째의 흡광도 데이타를 b로서 취득하고, 스탭 S2005로 이행한다.
스탭 S2005에서는 'yi번째의 데이타를 취득(c)'의 처리가 실행된다. 이 처리에서는 UV스펙트럼 U2를 대입한 y의 i번째의 흡광도 데이타를 c로서 취득하고, 스탭 S2006로 이행한다.
스탭 S2006에서는 'UV스펙트럼간 거리(d)와 UV스펙트럼간 거리의 제곱합(z)을 산출'의 처리가 실행된다. 이 처리에서는, UV스펙트럼간 거리(d)와 UV스펙트럼 간 거리의 제곱합(z)을
d = b - c
z = z + d2
으로 산출하고, 스탭 S2007로 이행한다.
스탭 S2007에서는 'i의 경신(i ← i + 1)'의 처리가 실행된다.이 처리에서는, i의 경신으로서 i에 i+1을 대입하고, 스탭 S2008로 이행한다.
스탭 S2008에서는 'x의 전 데이타에서 처리종료?'의 판단처리가 실행된다. 이 처리에서는, x와 y의 모든 데이타의 처리가 종료했는지의 여부가 판단된다. 종료한(YES) 경우는 x와 y의 전 데이타의 처리가 종료했다고 판단하고, S2009로 이행한다. 종료하지 않은(NO) 경우는 x와 y에 미처리의 데이타가 남아 있는 것으로 판단하고, 스탭 S2004으로 이행한다. 즉, S2004 ~ S2008의 처리는 x와 y의 모든 흡광도 데이타가 처리되기까지 반복한다.
스탭 S2009에서는 'x와 y의 UV스펙트럼의 일치도를 계산(U_Sim)'의 처리가 실행된다. 이 처리에서는 UV_Sim을 상기 UV스펙트럼간 거리의 제곱합(z)과 x의 데이타수 a로부터
UV_Sim = √ ( z / a )
로 산출하고, 이 UV_Sim을 도 81의 스탭 S306으로 넘기고, UV스펙트럼의 일치도계산처리를 종료한다.
[섭루틴3]
도 87은 도 81의 '대상FP귀속처리2'에 있어서의 '섭루틴3'의 상세를 도시하는 플로우차트이다. 이 처리는 피크패턴의 일치도를 계산한다.
스탭 S3001에서는 '피크패턴구성후보수(m)와 피크패턴구성피크수(n)를 설정'의 처리가 실행된다. 이 처리에서는 피크패턴을 망라적으로 작성하기 위한 설정으로서, 피크패턴구성후보수(m)와 피크패턴구성피크수(n)를 각각 설정하고, 스탭 S3002로 이행한다.
스탭 S3002에서는 'x ← 대상FP명(名), r1 ← R1, p1 ← P1, y ← 기준FP명, r2 ← R2, p2 ← P2)'의 처리가 실행된다.이 처리에서는, 처리에 필요한 대상FP 및 기준FP의 파일명 및 도 81의 S302와 S304에서 취득한 리텐션타임 및 피크데이타를 각각 x, r1, p1과 y, r2, p2에 대입하고, 스탭 S3003으로 이행한다.
스탭 S3003에서는 'x의 전 리텐션타임을 취득(a)'의 처리가 실행된다. 이 처리에서는 S3002에서 x에 대입한 이름의 파일(대상FP)을 판독하여, 그 파일의 전 리텐션타임을 a로서 취득하고, 스탭 S3004로 이행한다.
스탭 S3004에서는 'y의 전 리텐션타임을 취득(b)'의 처리가 실행된다. 이 처리에서는 S3002에서 y에 대입한 이름의 파일(기준FP)을 판독하여, 그 파일의 전 리텐션타임을 b로서 취득하고, 스탭 S3005로 이행한다.
스탭 S3005에서는 'a로부터 r1의 피크패턴구성후보피크 m개의 리텐션타임(cm) 및 피크데이타(dm)를 취득'의 처리가 실행된다. 이 처리에서는 a로부터 귀속대상피크의 리텐션타임인 r1의 피크패턴구성후보피크 m개의 리텐션타임을 cm, 피크데이타를 dm으로서 각각 취득하고, 스탭 S3006으로 이행한다. 또한, 피크패턴구성후보피크 m개는 r1과 리텐션타임이 가까운 m개이다.
스탭 S3006에서는 'b로부터 r2의 피크패턴구성후보피크 m개의 리텐션타임(em) 및 피크데이타(fm)를 취득'의 처리가 실행된다. 이 처리에서는 b로부터 귀속후보피크의 리텐션타임인 r2의 피크패턴구성후보피크 m개의 리텐션타임을 em, 피크데이타를 fm으로서 각각 취득하고, 스탭 S3007로 이행한다. 또한, 피크패턴구성후보피크 m개는 r2와 리텐션타임이 가까운 m개이다.
스탭 S3007에서는 'cm, dm을 리텐션타임순(올림차순)으로 배열함'의 처리가 실행된다. 이 처리에서는, S3005에서 취득한 cm과 dm을 리텐션타임이 올림차순이 되도록 교체배열하고, 스탭 S3008로 이행한다.
스탭 S3008에서는 'em, fm을 리텐션타임순(올림차순)으로 배열함'의 처리가 실행된다. 이 처리에서는, S3006에서 취득한 em과 fm을 리텐션타임이 올림차순이 되도록 교체배열하고, 스탭 S3009로 이행한다.
스탭 S3009에서는 'cm, dm으로부터 피크패턴구성피크 n개의 리텐션타임(cn) 및 피크데이타(dn)를 순번으로 취득'의 처리가 실행된다. 이 처리에서는 피크패턴구성후보피크 m개의 cm 및 dm으로부터 피크패턴구성피크 n개의 리텐션타임을 cn, 피크데이타를 dn으로 순번으로 취득하고, 스탭 S3010으로 이행한다.
스탭 S3010에서는 'em, fm으로부터 피크패턴구성피크 n개의 리텐션타임(en) 및 피크데이타(fn)를 순번으로 취득'의 처리가 실행된다. 이 처리에서는 피크패턴구성후보피크 m개의 em 및 fm으로부터 피크패턴구성피크 n개의 리텐션타임을 en, 피크데이타를 fn으로 순번으로 취득하고, 스탭 S3011로 이행한다.
스탭 S3011에서는 '피크패턴의 일치도를 계산(P_Sim)'의 처리가 실행된다. 이 처리에서는, 지금까지 취득한 귀속대상피크의 r1과 p1 및 그 피크패턴구성피크 n개의 cn과 dn, 및 귀속후보피크의 r2와 p2 및 그 피크패턴구성피크 n개의 en과 fn으로부터 피크패턴의 일치도(P_Sim)를 n = 4의 경우를 예로 들면, 도 64 중과 같이
P_Sim = (│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)
로 산출하고, 스탭 S3012로 이행한다.
스탭 S3012에서는 'P_Sim을 보존(P_Sim_all)'의 처리가 실행된다. 이 처리에서는 S3011에서 산출된 P_Sim을 순차 P-Sim_all에 보존하고, 스탭 S3013으로 이행한다.
스탭 S3013에서는 'em 중의 m개로부터 취출하는 전체 조합 종료?'의 판단처리가 실행된다. 이 처리에서는 귀속후보피크의 피크패턴구성후보피크 m개로부터 피크패턴구성피크 n개를 취출하는 전체 조합에서 처리가 종료했는지의 여부가 판단된다. 종료한(YES) 경우는 귀속후보피크에서 망라적인 피크패턴의 작성과 그 패턴에 있어서의 일치도의 계산이 종료했다고 판단하고, 스탭 S3014로 이행한다. 종료하지 않은(NO) 경우는 m개로부터 n개를 취출하는 조합이 종료하지 않은 것으로 판단하고, 스탭 S3010으로 이행한다. 즉, S3010 ~ S3013의 처리는 m개로부터 n개를 취출하는 전체 조합에서 처리가 종료하기까지 반복된다.
스탭 S3014에서는 'cm 중의 m개로부터 취출하는 전체 조합 종료?'의 판단처리가 실행된다. 이 처리에서는 귀속대상피크의 피크패턴구성후보피크 m개로부터 피크패턴구성피크 n개를 취출하는 전체 조합에서 처리가 종료했는지의 여부가 판단된다. 종료한(YES) 경우는 귀속대상피크에서 망라적인 피크패턴의 작성과 그 패턴에 있어서의 일치도의 계산이 종료했다고 판단하고, 스탭 S3015로 이행한다. 종료하지 않은(NO) 경우는 m개로부터 n개를 취출하는 조합이 종료하지 않은 것으로 판단하고, 스탭 S3009로 이행한다. 즉, S3009 ~ S3014의 처리는 m개로부터 n개를 취출하는 전체 조합에서 처리가 종료하기까지 반복한다.
스탭 S3015에서는 'P_Sim_all로부터 최소치를 취득(P_Sim_min)'의 처리가 실행된다. 이 처리에서는 S3012에서 보존한 P_Sim_all의 최소치를 P_Sim_min으로 취득하고, 이 P_Sim_min을 도 81의 스탭 S307로 넘기고, 피크패턴의 일치도계산처리를 종료한다.
[기준군FP의 작성]
대상FP특징량데이타를 기준FP특징량데이타와 비교하기 위한 기준FP특징량파일은 도 88 ~ 도 91과 같이 작성된다.
도 88은 기준FP특징량통합파일을 작성하기 위한 플로우차트이고, 기준FP작성부의 FP작성기능과, 기준FP피크귀속부의 기준FP피크귀속기능과, 기준FP귀속결과통합부의 기준FP귀속결과통합기능과, 기준FP피크특징량작성부의 기준FP피크특징량기능을 컴퓨터에 실현시킨다.
기준FP작성기능은 스탭 S10001에서 실현된다. 기준FP피크귀속기능은 스탭 S10002, S10003, S10004에서 실현된다. 기준FP귀속결과통합기능은 S10005에서 실현된다. 기준FP피크특징량작성기능은 S10006에서 실현된다.
S10001 ~ S10004는, 도 76의 대상FP특징량통합파일의 작성에 관한 S1 ~ S4에 대응한다.
스탭 S10001은 3D크로마토 및 특정의 검출파장에 있어서의 피크정보를 입력데이타로 하여 'FP작성처리'가 실행된다.
3D크로마토 및 피크데이타와 함께 평가기준이 되는 복수의 평가기준 약제(기준한방약)마다 구비되어 있다.
스탭 S10001에서는, 컴퓨터의 기준FP작성부가 기능하여, 3D크로마토 및 피크정보로부터 기준FP가 대상FP(17)(도 2)와 마찬가지로 작성되고, 기준FP의 데이타가 파일로서 출력된다.
스탭 S10002는 스탭 S10001에서 출력한 전 기준FP를 입력으로 하여 '기준FP귀속처리1'이 실행된다.
스탭 S10002에서는 컴퓨터의 기준FP피크귀속부가 기능하여, 전 기준FP를 대상으로 하여, 선택한 조합 또한 순번으로 귀속스코어를 산출하기 위해 전 기준FP로부터 조합을 선택하여 스탭 S10003으로 이행한다.
스탭 S10003은 선택한 기준FP의 조합을 입력으로 하여, '기준FP귀속처리2'가 실행된다.
스탭 S10003에서는, 스탭 S2에서 선정한 기준FP의 조합의 전 피크에서, 도 23 ~ 도 61과 같이 망라적으로 피크패텬을 작성하고, 이어서 이들 피크패턴의 일치도(도 63 또는 도 64의 UV_Sim)를 산출한다. 또한, 선정한 기준FP의 조합의 피크 사이에서 UV스펙트럼의 일치도(도 66의 P_Sim)를 산출한다. 또한 이들 2개의 일치도로부터 귀속후보피크의 일치도(도 67의 SCORE)를 산출한다. 그 산출결과는 판정결과파일(도 95의 판정결과파일 예 129 참조)로 출력된다.
스탭 S10004는 스탭 S10003에서 출력한 판정결과파일을 입력으로 하여, '기준FP귀속처리3'가 실행된다.
스탭 S10004에서는 선정한 기준FP의 조합 사이에서 귀속후보피크의 일치도(SCORE)을 기초로 선정한 기준FP의 조합 사이에 대응하는 피크를 특정한다. 그 결과는 기준FP마다 귀속데이타로서 출력된다.
스탭 S10005는 스탭 S10004에서 출력한 전 기준FP귀속데이타를 입력으로 하여, '기준FP귀속결과통합처리'가 실행된다.
스탭 S10005에서는 컴퓨터의 기준FP귀속결과통합부가 기능하여, 기준FP피크귀속부에서 특정한 각 기준FP의 피크대응관계를 참조하여, 전 기준FP귀속데이타를 통합하여 기준FP대응표를 작성하고, 스탭 S10006으로 이행한다.
스탭 S10006에서는 컴퓨터의 기준FP피크특징량작성부가 기능하여, 기준FP귀속결과통합부에서 작성된 기준FP대응표를 기초로 전 기준FP특징량(기준군FP)을 작성한다. 이 기준FP피크특징량작성부에서의 처리는 기준FP대응표의 각 피크(예)에서 통계량(최대치, 최소치, 중앙값, 평균치 등)을 산출하고, 그 정보를 기초로 피크(예)를 선정한다. 선정한 피크(예)는 기준군FP(도 98의 기준군FP 예 137 참조)로서 출력된다.
[S10005 : 기준FP대응표의 작성]
도 89, 90은 스탭 S10005의 '기준FP귀속결과통합처리(기준FP대응표의 작성)'의 상세를 도시하는 플로우차트이다.
스탭 S10101에서는 '귀속순번 1번의 귀속데이타를 통합데이타로서 판독'의 처리가 실행된다. 이 처리에서는 스탭 S10004에서 첫번째로 귀속처리하여 피크의 대응관계를 특정한 기준FP귀속데이타를 통합데이타로서 판독하고, 스탭 S10102로 이행한다.
스탭 S10102에서는 '귀속데이타 2번째 이후를 순번으로 판독'의 처리가 실행된다. 이 처리에서는 먼저 스탭 S10004에서 2번째로 귀속처리하여 피크의 대응관계를 특정한 기준FP귀속데이타를 통합데이타로서 판독하고, 스탭 S10103으로 이행한다.
스탭 S10103에서는 '통합데이타와 귀속데이타를 공통으로 하는 피크데이타로 통합'의 처리가 실행된다. 이 처리에서는, 통합데이타와 귀속데이타에서 공통으로 존재하는 기준FP의 피크데이타를 기초로 2개의 파일을 통합하여, 그 결과로 통합데이타를 경신하고, 스탭 S10104로 이행한다.
스탭 S10104에서는 '귀속데이타 중의 모든 피크를 통합데이타에 추가?'의 판단처리가 실행된다. 이 처리에서는 귀속데이타의 모든 피크가 통합데이타에 추가되었는지의 여부가 판단된다. 추가된(YES) 경우는 스탭 S10105로 이행한다. 추가되지 않은 피크(흠락피크)가 있는(NO) 경우는 이 흠락(欠落)피크를 통합데이타에 추가처리하기 위해, 스탭 S10107로 이행한다. 또한, 흠락피크의 통합데이타로의 추가처리(S10107 - S10120)는 상기 S5(대상FP귀속처리4)에 있어서의 스탭 S504 - S517와 마찬가지의 처리가 행해진다.
스탭 S10121에서는 '통합데이타에 TEMP2의 데이타를 추가(모든 리텐션타임과 피크)'의 처리가 실행된다. 이 처리에서는 TEMP2의 모든 리텐션타임(R3)과 피크(P1)를 통합데이타의 해당 개소에 추가하고, 스탭 S10122로 이행한다.
스탭 S10122에서는 '역치2 ← 초기치, TEMP2 내의 데이타를 모두 삭제'의 처리가 실행된다. 이 처리에서는, UV_Sim에 경신되어 있는 역치2를 초기치로 복귀시키고, 전 흠락피크의 리텐션타임이나 피크 등의 데이타가 들어간 TEMP2로부터 모든 데이타를 삭제하고, 스탭 S10104로 복귀한다.
스탭 S10104로부터 이행하는 스탭 S10105에서 '모든 귀속데이타의 처리종료 ?'의 판단처리가 실행된다. 이 처리에서는 전 기준데이타의 처리가 종료했는지의 여부가 판단된다. 처리가 종료(YES)의 경우는 전 귀속데이타의 통합결과인 기준FP대응표를 출력하기 위해, 스탭 S10106으로 이행한다. 모든 처리가 종료하지 않은(NO)의 경우는 S10102로 복귀하여, 남은 귀속데이타를 순차처리한다.
스탭 S10106에서는 '통합데이타를 출력(기준FP대응표)'의 처리가 실행된다. 이 처리에서는 전 귀속데이타를 통합한 결과를 기준FP대응표로 출력하고, 기준FP대응표의 작성처리를 종료한다.
[S10006 : 피크특징량화처리]
도 91은 도 88 스탭 S10006의 '피크특징화처리(기준군FP의 작성)'의 상세를 도시하는 플로우차트이다.
스탭 S10201에서는 '기준FP대응표를 판독'의 처리가 실행된다. 이 처리에서는 스탭 S10005에서 작성한 기준FP대응표를 판독하고, 스탭 S10202로 이행한다.
스탭 S10202에서는 '각 피크(예)에서 통계량을 산출'의 처리가 실행된다. 이 처리에서는 기준FP대응표의 각 피크(예)에서 통계량(최대치, 최소치, 중앙값, 평균치, 분산, 표준편차, 존재수, 존재율)을 산출하고, 스탭 S10203으로 이행한다.
스탭 S10203에서는 '산출한 통계량을 참고로 피크(예)를 선정'의 처리가 실행된다. 이 처리에서는 스탭 S10202에서 산출한 통계량을 참고로 피크를 선정하고, 스탭 S10204로 이행한다.
스탭 S10204에서는 '선정한 피크(예)를 출력(기준군FP)'의 처리가 실행된다. 이 처리에서는 통계량에 의해 피크(예)의 선정결과를 기준군FP로 출력하고, 기준군FP의 작성처리를 종료한다.
도 98에 상기와 같이 출력하는 기준FP대응표 예(137)을 도시한다.
[실시예 1의 효과]
본 발명의 실시예 1에서는 평가대상이 되는 다성분 약제의 3D크로마토(15)로부터 특정의 파장, 예를 들면 203nm에서 검출된 피크와 그 리텐션타임, 및 그 피크의 UV스펙트럼으로 구성되는 대상FP(17)를 작성하는 FP작성공정(113)과, 대상FP(17)의 피크귀속에 적합한 기준FP를 복수의 기준FP로부터 선정하는 기준FP선정공정(115)과, 대상FP로 선정된 기준FP의 피크 각각에 대하여, 예를 들면 시간축방향 전후의 적어도 일방에 존재하는 2개의 피크를 포함한 3개의 피크로 구성되는 피크패턴을 작성하는 피크패턴작성공정(117)과, 상기 피크패턴 및 피크의 UV스펙트럼을 비교하여 대응하는 피크를 특정하는 피크귀속공정(119)과, 귀속한 피크와 복수의 기준FP의 피크를 예를 들면 MT법에 의해 비교평가하는 평가공정(121)을 구비했다.
평가대상이 되는 다성분 약제의 3D크로마토(15)를 이들 5개 공정(113, 115, 117, 119, 121)으로 처리하는 것으로, 평가대상 약제의 품질평가의 정도 및 효율을 보다 향상시킬 수 있다.
FP작성공정(113)에서 작성된 대상FP(17)는 3D크로마토(15)와 마찬가지로 삼차원의 정보(피크, 리텐션타임 및 UV스펙트럼)로 구성되어 있다. 그 때문에, 그 약제 특유의 정보를 그대로 계승한 데이타이다. 그럼에도 불구하고, 데이타용량은 약 1/70으로 압축되어 있기 때문에, 3D크로마토(15)에 비교하여, 처리해야할 정보량을 대폭으로 감소시킬 수 있어서 처리속도를 빠르게 할 수 있다.
FP작성공정(113)은 검출파장이 상이한 복수의 FP를 합성한 FP를 작성한다. 이것에 의해 1개의 파장으로 모든 성분을 검출할 수 없는 성분이 조합된 다성분 약제라고 해도, 복수의 검출파장의 FP를 합성하는 것으로 모든 성분을 포함한 품질평가가 가능하게 된다.
FP작성공정(113)은 3D크로마토에서 검출된 전 피크를 포함한 FP를 작성한다. 이 때문에, 다성분 약제인 한방약의 품질평가에 적합하다.
기준FP선정공정(115)에서는 대상FP에 귀속에 적합한 기준FP를 FP간의 리텐션타임출현패턴을 비교하여, 패턴의 일치도가 좋은 기준FP를 선정한다. 이것에 의해, 피크귀속처리공정(119)에서 패턴이 유사한 FP 사이에서 귀속처리가 가능하기 때문에, 정도가 높은 귀속이 가능하게 된다.
피크패턴작성공정(117)에서는 귀속대상피크 및 귀속후보피크 각각에서 복수의 주변피크를 사용하여 망라적으로 피크패턴을 작성한다. 이것에 의해, 대상FP와 기준FP에서 FP 전체의 패턴이 다소 상이하여도, 피크귀속처리공정(119)에서 정도가 높은 귀속이 가능하게 된다.
피크귀속공정(119)에서는 피크패턴작성공정(117)에서 작성된 피크패턴의 일치도에 추가하여, 귀속대상피크와 귀속후보피크의 UV스펙트럼의 일치도도 가미하여 귀속해야할 피크를 특정하고 있다. 그 때문에, 정도가 높은 귀속이 가능하게 된다.
피크귀속공정(119)에서는 대상FP의 전 피크를 기준FP의 피크로 일제 귀속한다. 그 때문에, 효율이 좋은 귀속처리가 가능하게 된다.
평가공정(121)에서는 다차원 데이타인 다성분으로 구성된 FP를 MT법으로 MD값으로서 1차원으로 집약하고, 복수의 평가대상로트를 간단히 비교평가한다. 이 때문에 복수의 성분으로 구성되는 다성분계 약제의 평가에 적합하다.
본 발명의 실시예의 다성분 약제의 평가프로그램은 각 기능을 컴퓨터에 실현시켜서, 평가의 정도 및 효율을 보다 향상시킬 수 있다.
본 발명의 실시예의 다성분 약제의 평가장치는 각 부(3, 5, 7, 9, 11)를 작용시켜서, 평가의 정도 및 효율을 보다 향상시킬 수 있다.
[피크패턴의 일치도계산(P_Sim)의 변형례]
도 63, 도 64, 도 87에서의 피크패턴의 일치도계산(P_Sim)은 FP를 피크높이로 작성한 상기 실시예의 경우에 대하여 적용하여, 비교대상의 피크높이의 차에 기초하여 계산했다.
한편, 본 발명의 평가방법, 평가프로그램, 및 평가장치에 있어서의 피크에 대해서는 상기와 같이 시그널강도(높이)의 극대치를 의미하는 경우와, 시그널강도의 면적값(피크면적)을 높이로 표현한 것을 의미하는 경우의 어느 것도 포함할 수 있다.
즉, FP를 피크면적으로 작성할 때도, 면적값을 높이로 표현하여 FP를 작성하기 위해, FP로서는 상기 실시예의 피크높이로 작성하는 경우와 마찬가지의 표현이 된다. 이 때문에, FP를 피크높이로 작성한 경우와 마찬가지로 상기 실시예의 처리에 의해 마찬가지로 평가할 수 있다.
단, FP를 피크면적으로 작성했을 때는 비교대상의 피크값의 차가 크게 되기 때문에, 비율에 기초한 계산으로 하여, 취급을 용이하게 하는 것이 적합하다.
이하, 비율에 기초하여 계산하는 피크패텬의 일치도(P_Sim)를 n = 2, n = 4의 경우를 예로 하여 표시한다.
n = 2의 경우
P_Sim = (p1 / p2＃1) × (│r1 - (r2 + d)│ + 1)
+ (dn1 / fn1＃1) × (│cn1 - r1) - (en1 - r2│ + 1)
+ (dn2 / fn2＃1) × (│cn2 - r1) - (en2 - r2│ + 1)
n = 4의 경우
P_Sim = (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 은 비교대상의 2개의 값의 비(큰 값 / 작은 값)인 것을 표시하고 있다.
또한, FP를 피크높이로 작성한 경우에도, 비율에 기초하여 피크패턴의 일치도(P_Sim)를 계산할 수 있고, FP를 피크면적으로 작성한 경우에도, 상기 높이의 차와 마찬가지로 피크면적값의 차에 기초하여 피크패턴의 일치도(P_Sim)를 얻을 수 있다.
[섭루틴2의 변형례]
도 100은 도 86 대신에 적용하는 섭루틴2의 변형례에 관한 것으로, 도 81의 '대상FP귀속처리2'에 있어서의 '섭루틴2'의 변형례의 상세를 도시하는 플로우차트이다. 이 변형례에 관한 처리에 의해 UV스펙트럼의 일치도를 계산한다.
본 섭루틴2의 변형례에서는, 도 86의 섭루틴2에서 RMSD에 UV패턴의 이동평균의 경사정보(DNS)를 가미하는 처리가 가능하도록 했다. DNS는 후술의 식으로 표현되고, UV패턴에 있어서의 이동평균치의 이동경사를 2개의 패턴으로 비교한 때의 경사부호( +/- )의 불일치수로서 정의된다. 즉, DNS는 UV패턴의 극대, 극소치의 위치의 일치상태를 평가하는 값이다.
이 DNS의 정보를 상기 RMSD에 가미하는 것으로 UV스펙트럼의 파형의 일치도를 보다 정확하게 산출할 수가 있다.
도 100의 변형례에 관한 섭루틴2에 있어서, 스탭 S2001 ~ S2008는 도 86의 섭루틴2와 거의 동일하다. 단, 스탭 S2001에서는 구간1 ← w1, 구간2 ← w2의 초기설정이 추가되어 행해지고, 후술하는 이동평균, 이동경사의 계산의 구간이 이용된다.
본 변형례의 섭루틴2에서는 DNS가미를 위해 스탭 S2010 ~ S2013을 추가하고, 스탭 S2009A에서 DNS를 가미한 일치도의 계산을 가능하게 했다.
스탭 S2010에서는 'DNS를 가미함?'의 판단처리가 실행되고, DNS를 가미한다고 판단된 때는(YES), 스탭 S2011로 이행하고, DNS를 가미하지 않는다고 판단된 때는(NO), 스탭 S2009A로 이행한다. DNS를 가미하는지 여부의 기인(起因)은 예를 들면 초기의 설정에 따른다. 예를 들면, FP를 피크면적으로 작성했을 때는 DNS가미, FP를 피크높이로 작성한 때는 DNS비가미로 설정한다.
단, FP를 피크높이로 작성한 상기 실시예의 경우에도, DNS를 가미하는 처리에서 UV패턴일치도를 계산할 수 있고, FP를 피크면적으로 작성한 경우에도 DNS를 가미하지 않고 상기 실시예의 처리에서 UV패턴일치도를 계산할 수 있다.
스탭 S2011에서는 '구간1(w1)에서 x와 y의 이동평균을 계산'의 처리가 실행되어, 구간1(w1)에 있어서의 이동평균이 구해진다. 구간1(w1)은 UV데이타의 파장에 관한 구간이고, 스탭 S2001의 초기설정에 있어서 w1 = 3이라면 구간1(3)로 되고, 3개의 파장에 있어서의 UV강도의 평균이 구해진다. 구체적으로는 도 101의 도표에서 후술한다.
스탭 S2012에서는 '구간2(w2)에서 x와 y의 이동경사를 계산'의 처리가 실행되고, 구간2(w2)에서의 이동경사가 구해진다. 구간2(w2)는 스탭 S2011에서 구해진 이동평균에 관한 구간이고, 스탭 S2001의 초기설정에 있어서 w2 = 3이라면 구간2(3)로 되고, 스탭 S2011에서 계산한 이동평균에 기초하여, 3개의 이동평균에 걸친 경사의 (±)가 구해진다. 구체적으로는 도 101의 도표에서 후술한다.
스탭 S2013에서는 'x와 y의 이동경사의 부호의 불일치를 계산(DNS)'의 처리가 실행되고, 스탭 S2012에서 계산된 이동경사로부터 경사의 (±)의 일치수가 계산된다. 이동경사 +는 도 66에서 상승하는 것을 표시하고, 이동경사 -는 하강하는 것을 표시한다.
스탭 S2013으로부터 스탭 S2009A로 이행하면, 이 스탭 S2009A의 처리에 있어서 DNS를 가미한 일치도의 계산이 행해진다.
스탭 S2009A에서는 'x와 y의 UV스펙트럼의 일치도를 계산(U_Sim)'의 처리가 실행되고, DNS를 가미한 일치도의 계산에서는 UV_Sim을 상기 UV스펙트럼간 거리의 제곱합(z)과 x의 데이타수 a와 DNS로부터
UV_Sim = √ ( z / a ) × 1. 1DNS
로 산출하여, 이 UV_Sim을 도 81의 스탭 S306으로 넘기고, UV스펙트럼의 일치도계산처리를 종료한다.
또한, 스탭 S2010으로부터 스탭 S2009A로 이행한 경우의 처리는 도 86의 스탭 S2009와 동일하다.
도 101은 이동평균 및 이동경사의 계산 예를 도시하는 도표이다.
도 101의 상단은 UV데이타 예, 중단은 이동평균의 계산 예, 하단은 이동경사의 계산 예를 도시한다. UV데이타 예는 구체적인 수치로 대입하고, UV강도를 a1 ~ a7으로 표기하고 있다. 예를 들면, 220nm의 UV강도가 a1, 221nm의 UV강도가 a2 등으로 된다. 이동평균의 계산 예 및 이동경사의 계산 예도 구체적인 수치로 대입하여, UV강도 a1 ~ a7을 사용하고 있다.
이동평균은 구간1(w1 = 3)을 예로 하여, 스탭 S2012(도 100)에서 구간(a1, a2, a3), 구간(a2, a3, a4), … 마다 계산한 값으로 하여 m1, m2, …가 계산된다. 이동경사도 구간2(3)를 예로 하여, 스탭 S2013(도 100)에서 구간(m1, m2, m3), 구간(m2, m3, m4), … 마다 계산한 값으로 하여 s1, …가 계산된다. 예를 들면, 이동평균의 차 m3 - m1이 이동경사가 되고, 그의 (±)를 취출한다.
이렇게 하여, FP를 피크면적으로 작성한 때는 기준군FP로의 귀속처리 및 기준FP귀속결과통합처리에서 DNS를 가미한 처리로 UV패턴일치도를 계산할 수 있다. 이 계산에 의해, 도 66에 도시하는 대응하는 2점의 거리(dis)가 피크높이로 작성한 FP에 비교하여 크게 되어도, 취급을 용이하게 하고, UV패턴일치도를 정확히 계산할 수 있다.
[기타]
본 발명의 실시예는 다성분 약제로서 한방약의 평가에 대해서 적용했지만, 그 외의 다성분 물질의 평가에도 적용할 수 있다.
상기 실시예의 FP는 3D크로마토 상에서의 전 피크를 대상으로 했지만, 미세한 데이타, 예를 들면 3D크로마토 상에서 피크면적이 5% 미만의 피크를 제외하고 FP를 작성할 수도 있다.
상기 실시예의 FP는 피크높이에 기초하여 작성하여, 도 70 ~ 도74의 평가를 얻었지만, 피크면적에 기초하여 FP를 작성한 경우에 대해서도, 피크높이에 기초하여 작성한 상기 실시예와 마찬가지의 수순에 의해 MT법에 의해 MD값을 구하여, 도 70 ~ 도 74와 마찬가지로 평가를 얻을 수가 있다.
1 다성분 약제의 평가장치
3 FP작성부(FP작성장치)
13 한방약(평가대상)
15 3D크로마토
17 대상FP(대상패턴)
19 기준군FP
21 기준군FP에 귀속한 대상FP
25 대상FP에 포함되는 피크의 UV스펙트럼
27 약제A의 FP
29 약제B의 FP
31 약제C의 FP
33 대상FP(리텐션타임 10.0 - 14.5분)
35, 37, 39, 41, 43, 45, 47, 49, 51, 53 대상FP(리텐션타임 10.0 - 14.5분) 중의 각 피크
55 기준FP(리텐션타임 10.0 - 14.5 분)
107 귀속대상피크의 UV스펙트럼
111 귀속후보피크의 UV스펙트럼
113 FP작성공정
123 3D크로마토 데이타 예
125 피크정보 데이타 예
127 FP 데이타 예
3 FP작성부(FP작성장치)
13 한방약(평가대상)
15 3D크로마토
17 대상FP(대상패턴)
19 기준군FP
21 기준군FP에 귀속한 대상FP
25 대상FP에 포함되는 피크의 UV스펙트럼
27 약제A의 FP
29 약제B의 FP
31 약제C의 FP
33 대상FP(리텐션타임 10.0 - 14.5분)
35, 37, 39, 41, 43, 45, 47, 49, 51, 53 대상FP(리텐션타임 10.0 - 14.5분) 중의 각 피크
55 기준FP(리텐션타임 10.0 - 14.5 분)
107 귀속대상피크의 UV스펙트럼
111 귀속후보피크의 UV스펙트럼
113 FP작성공정
123 3D크로마토 데이타 예
125 피크정보 데이타 예
127 FP 데이타 예
Claims (16)
- 평가대상의 크로마토그램 데이타로부터 검출된 피크와 그 리텐션타임으로 구성하는 핑거프린트 데이타(FP)를 작성하는 FP작성공정을 구비한 FP작성방법에 있어서,
상기 FP작성공정은 상기 크로마토그램 데이타를 리텐션타임, 검출파장, 및 피크를 데이타로 하는 삼차원 크로마토그램 데이타로 하고, 이 삼차원 크로마토그램 데이타로부터 특정파장에서 검출된 피크의 리텐션타임과 이 리텐션타임마다 나타나는 피크의 크기를 이차원으로 표시하고, 상기 이차원으로 표시된 피크마다 상기 삼차원 크로마토그램 데이타의 검출파장에 기초하는 리텐션타임마다의 이차원의 UV스펙트럼을 가진 핑거프린트 데이타를 작성하는 것을 특징으로 하는 FP작성방법. - 청구항 1에 있어서,
상기 평가대상은 다성분 약제인 것을 특징으로 하는 FP작성방법. - 청구항 2에 있어서,
상기 다성분 약제는 생약, 생약의 조합, 그들의 추출물, 한방약의 어느 하나인 것을 특징으로 하는 FP작성방법. - 청구항 1 내지 청구항 3중 어느 하나의 항에 있어서,
상기 FP작성공정은 상기 특정파장이 상이한 복수의 핑거프린트 데이타로부터 이들을 합성한 핑거프린트 데이타를 작성하는 것을 특징으로 하는 FP작성방법. - 청구항 1 내지 청구항 3중 어느 하나의 항에 있어서,
상기 FP작성공정은 상기 삼차원 크로마토그램 데이타의 전 피크를 추출하여 상기 핑거프린트 데이타를 작성하는 것을 특징으로 하는 FP작성방법. - 평가대상의 크로마토그램 데이타로부터 검출된 피크와 그 리텐션타임으로 구성하는 핑거프린트 데이타(FP)를 작성하는 FP작성기능을 컴퓨터에 실현시키는 평가대상의 FP작성프로그램을 기록한 컴퓨터로 읽을 수 있는 매체에 있어서,
상기 FP작성기능은 상기 크로마토그램 데이타를 리텐션타임, 검출파장, 및 피크를 데이타로 하는 삼차원 크로마토그램 데이타로 하고, 이 삼차원 크로마토그램 데이타로부터 특정파장에서 검출된 피크의 리텐션타임과 이 리텐션타임마다 나타나는 피크의 크기를 이차원으로 표시하고, 상기 이차원으로 표시된 피크마다 상기 삼차원 크로마토그램 데이타의 검출파장에 기초하는 리텐션타임마다의 이차원의 UV스펙트럼을 가진 핑거프린트 데이타를 작성하는 것을 특징으로 하는 FP작성프로그램을 기록한 컴퓨터로 읽을 수 있는 매체. - 청구항 6에 있어서,
상기 평가대상은 다성분 약제인 것을 특징으로 하는 FP작성프로그램을 기록한 컴퓨터로 읽을 수 있는 매체. - 청구항 7에 있어서,
상기 다성분 약제는 생약, 생약의 조합, 그들의 추출물, 한방약의 어느 하나인 것을 특징으로 하는 FP작성프로그램을 기록한 컴퓨터로 읽을 수 있는 매체. - 청구항 6 내지 청구항 8중 어느 하나의 항에 있어서,
상기 FP작성기능은 상기 특정파장이 상이한 복수의 핑거프린트 데이타로부터 이들을 합성한 핑거프린트 데이타를 작성하는 것을 특징으로 하는 FP작성프로그램을 기록한 컴퓨터로 읽을 수 있는 매체. - 청구항 6 내지 청구항 8중 어느 하나의 항에 있어서,
상기 FP작성기능은 상기 삼차원 크로마토그램 데이타의 전 피크를 추출하여 상기 핑거프린트 데이타를 작성하는 것을 특징으로 하는 FP작성프로그램을 기록한 컴퓨터로 읽을 수 있는 매체. - 평가대상의 크로마토그램 데이타로부터 검출된 피크와 그 리텐션타임으로 구성하는 핑거프린트 데이타(FP)를 작성하는 FP작성부를 구비한 평가대상의 FP작성장치에 있어서,
상기 FP작성부는 상기 크로마토그램 데이타를 리텐션타임, 검출파장, 및 피크를 데이타로 하는 삼차원 크로마토그램 데이타로 하고, 이 삼차원 크로마토그램 데이타로부터 특정파장에서 검출된 피크의 리텐션타임과 이 리텐션타임마다 나타나는 피크의 크기를 이차원으로 표시하고, 상기 이차원으로 표시된 피크마다 상기 삼차원 크로마토그램 데이타의 검출파장에 기초하는 리텐션타임마다의 이차원의 UV스펙트럼을 가진 핑거프린트 데이타를 작성하는 것을 특징으로 하는 FP작성장치. - 청구항 11에 있어서,
상기 평가대상은 다성분 약제인 것을 특징으로 하는 FP작성장치. - 청구항 12에 있어서,
상기 다성분 약제는 생약, 생약의 조합, 그들의 추출물, 한방약의 어느 하나인 것을 특징으로 하는 FP작성장치. - 청구항 11 내지 청구항 13중 어느 하나의 항에 있어서,
상기 FP작성부는 상기 특정파장이 상이한 복수의 핑거프린트 데이타로부터 이들을 합성한 핑거프린트 데이타를 작성하는 것을 특징으로 하는 FP작성장치. - 청구항 11 내지 청구항 13중 어느 하나의 항에 있어서,
상기 FP작성부는 상기 삼차원 크로마토그램 데이타의 전 피크를 추출하여 상기 핑거프린트 데이타를 작성하는 것을 특징으로 하는 FP작성장치. - 청구항 1 내지 청구항 3중 어느 하나의 항에 기재된 FP작성방법에 따라서 작성된 핑거프린터 데이타를 기록한 컴퓨터로 읽을 수 있는 매체에 있어서,
상기 삼차원 크로마토그램 데이타로부터 특정파장에서 검출된 피크와 그 리텐션타임 및 그 피크의 UV스펙트럼으로 작성되고, 상기 리텐션타임과 이 리텐션타임마다 나타나는 피크의 크기를 이차원으로 표시하고, 상기 이차원으로 표시된 피크마다 상기 삼차원 크로마토그램 데이타의 검출파장에 기초하는 리텐션타임마다의 이차원의 UV스펙트럼을 가지는 것을 특징으로 하는 핑거프린트 데이타를 기록한 컴퓨터로 읽을 수 있는 매체.
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2011123845 | 2011-06-01 | ||
| JPJP-P-2011-123845 | 2011-06-01 | ||
| PCT/JP2012/003607 WO2012164949A1 (ja) | 2011-06-01 | 2012-05-31 | Ｆｐ作成方法、作成プログラム、作成装置、及びｆｐ |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| KR20130037682A KR20130037682A (ko) | 2013-04-16 |
| KR101412655B1 true KR101412655B1 (ko) | 2014-06-26 |
Family
ID=47258819
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| KR1020127032269A KR101412655B1 (ko) | 2011-06-01 | 2012-05-31 | Fp작성방법, 작성프로그램, 작성장치, 및 fp |
Country Status (8)
| Country | Link |
|---|---|
| US (1) | US20140123736A1 (ko) |
| EP (1) | EP2717047A4 (ko) |
| JP (1) | JP5895847B2 (ko) |
| KR (1) | KR101412655B1 (ko) |
| CN (1) | CN102971625B (ko) |
| HK (1) | HK1181115A1 (ko) |
| TW (1) | TWI609181B (ko) |
| WO (1) | WO2012164949A1 (ko) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002214215A (ja) * | 2001-01-12 | 2002-07-31 | Council Scient Ind Res | クロマトグラフフィンガープリントならびに単一の医薬および処方物の標準化のための新規な方法 |
Family Cites Families (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2658344B2 (ja) * | 1989-01-26 | 1997-09-30 | 株式会社島津製作所 | クロマトグラフ用データ処理装置 |
| JP3094921B2 (ja) * | 1996-09-26 | 2000-10-03 | 株式会社島津製作所 | クロマトグラフ用データ処理装置 |
| UA78194C2 (en) * | 2000-12-08 | 2007-03-15 | Council Scient Ind Res | Method for revealing and identification of active components in extracts by chromatographic fingerprinting of "fingerprints" and processor of program data processing chromatograms |
| JP2005091068A (ja) * | 2003-09-16 | 2005-04-07 | Sumitomo Chemical Co Ltd | 着色試料の分析方法および分析装置 |
| KR101150231B1 (ko) * | 2004-01-28 | 2012-05-25 | 카운실 오브 사이언티픽 엔드 인더스트리얼 리서치 | 움직이는 크로마토그래픽 핑거프린팅을 이용한 식품 및의약품의 화학적 및 치료 가치의 표준화를 위한 방법 |
| CN1965233B (zh) * | 2004-06-07 | 2013-03-13 | 株式会社津村 | 多成分药剂的评价方法 |
| WO2006010117A2 (en) * | 2004-07-10 | 2006-01-26 | The Research Foundation Of State University Of New York | Production of flavonoids by recombinant microorganisms |
| US7447597B2 (en) * | 2005-05-06 | 2008-11-04 | Exxonmobil Research And Engineering Company | Data processing/visualization method for two (multi) dimensional separation gas chromatography xmass spectrometry (GCxMS) technique with a two (multiply) dimensional separation concept as an example |
| 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 |
-
2012
- 2012-05-31 CN CN201280001664.9A patent/CN102971625B/zh active Active
- 2012-05-31 KR KR1020127032269A patent/KR101412655B1/ko active IP Right Grant
- 2012-05-31 EP EP12792398.5A patent/EP2717047A4/en active Pending
- 2012-05-31 WO PCT/JP2012/003607 patent/WO2012164949A1/ja active Application Filing
- 2012-05-31 JP JP2012549178A patent/JP5895847B2/ja active Active
- 2012-05-31 US US13/805,294 patent/US20140123736A1/en not_active Abandoned
- 2012-06-01 TW TW101119811A patent/TWI609181B/zh active
-
2013
- 2013-07-18 HK HK13108442.2A patent/HK1181115A1/zh unknown
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002214215A (ja) * | 2001-01-12 | 2002-07-31 | Council Scient Ind Res | クロマトグラフフィンガープリントならびに単一の医薬および処方物の標準化のための新規な方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102971625B (zh) | 2015-11-25 |
| US20140123736A1 (en) | 2014-05-08 |
| EP2717047A4 (en) | 2014-12-03 |
| WO2012164949A1 (ja) | 2012-12-06 |
| HK1181115A1 (zh) | 2013-11-01 |
| KR20130037682A (ko) | 2013-04-16 |
| JPWO2012164949A1 (ja) | 2015-02-23 |
| TW201303290A (zh) | 2013-01-16 |
| CN102971625A (zh) | 2013-03-13 |
| JP5895847B2 (ja) | 2016-03-30 |
| TWI609181B (zh) | 2017-12-21 |
| EP2717047A1 (en) | 2014-04-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| KR101492898B1 (ko) | 피크귀속방법, 귀속프로그램, 및 귀속장치 | |
| KR101492895B1 (ko) | 패턴의 평가방법, 다성분 물질의 평가방법, 평가프로그램, 및 평가장치 | |
| KR101442117B1 (ko) | 패턴 또는 fp의 특징량작성방법, 작성프로그램, 및 작성장치 | |
| KR101436532B1 (ko) | 패턴의 평가방법, 다성분 물질의 평가방법, 평가프로그램, 및 평가장치 | |
| KR101412655B1 (ko) | Fp작성방법, 작성프로그램, 작성장치, 및 fp |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A201 | Request for examination | ||
| E902 | Notification of reason for refusal | ||
| E701 | Decision to grant or registration of patent right | ||
| GRNT | Written decision to grant | ||
| FPAY | Annual fee payment |
Payment date: 20170609 Year of fee payment: 4 |