JP7798677B2 - 制御装置、分離分析装置、分離分析方法、及び分離分析プログラム - Google Patents

Description

本開示は、制御装置、分離分析装置、分離分析方法、及び分離分析プログラムに関する。

例えば、特許文献１には、クロマトグラフィよりのクロマトグラムを受け、Ａ／Ｄ変換器でデジタルデータに変換して取込み、所定の分析処理を実行するクロマトグラフ用データ処理装置が記載されている。このクロマトグラフ用データ処理装置は、取込まれるデータをプロットする記録手段と、データよりクロマトグラムのピークを検出するピーク検出手段と、成分同定のためのピークの許容時間幅を予め記憶する許容時間幅記憶手段と、検出されたピーク値の相当時間が、許容時間幅内であるか否かを判別する手段と、判別の結果許容時間内であると、その成分名を記録手段の記録クロマトグラムのピーク値近傍に記録させる手段と、を備える。

また、特許文献２には、ガスクロマトグラフィ分析装置を用いてあらかじめ定めた対象物質の被測定物中への含有を判定する方法が記載されている。この方法は、目標とする検出下限を確保できるカラムを使用し、分析対象物質の標準物質を測定することによって得た結果から、分析対象物質に関して試料負荷量に応じたピーク面積と保持時間との相関式を算出し、被測定物を分析することにより得たクロマトグラム中にある分析対象物質のピークが検出されうる保持時間近傍の検出ピークの面積を得るとともに、前記相関式を活用して検出ピークの面積から分析対象物質のピークが出現する保持時間を推定し、その推定した保持時間と被測定物のクロマトグラム中の検出ピークの保持時間とを比較することにより被測定物中の分析対象物質の含有判定を行う。

血中の各種ヘモグロビン（以下、「Ｈｂ」ともいう。）についても上記のような分離分析方法を用いた成分の同定、測定が行われている。特に、ＨｂＤ、ＨｂＥ、ＨｂＳ、ＨｂＣ等の主要なヘモグロビン変異体は同定のニーズが高く、ピークの保持時間から同定が可能とされる。

特開平１－０７９６５６号公報 特開２０１２－１６３４７６号公報

しかしながら、昨今のＨｂの分離分析装置は、測定時間の短縮化によりＨｂ種の各ピーク保持時間の間は狭くなっている。このため、ピークの保持時間が近いＨｂ種のピークを区別することは困難である。例えば、陽イオン交換クロマトグラフィでＨｂＤを分析すると、ＨｂＤのピーク保持時間は、ＨｂＥのピーク保持時間と近い時間を示す。ＨｂＤのピーク保持時間の許容時間幅内にＨｂＥのピークが検出される場合もあり、上記特許文献１のように許容時間幅を設定してもＨｂＤのピークとＨｂＥのピークとの識別は困難である。また、事前に分析者によって試料濃度を調整することは可能だが、作業が煩雑であることから、血液の状態のまま測定できることが望ましい。

本開示は、上記の点に鑑みてなされたものであり、クロマトグラムのピークの保持時間が近い複数のピークについて、試料濃度を調整することなく、ピークの識別を行うことができる制御装置、分離分析装置、分離分析方法、及び分離分析プログラムを提供することを目的とする。

上記目的を達成するために、本開示の一態様に係る制御装置は、クロマトグラフィにより試料を分離、分析する分離分析装置の制御装置であって、予め準備された複数の試料に対してクロマトグラフィを実施して得られた、基準成分の基準ピークの保持時間、第１成分の第１ピークの保持時間と第２成分の第２ピークの保持時間との閾値、及び、試料中の成分の量を表す指標値の相関関係が記憶された記憶部と、前記基準成分と、前記第１成分又は前記第２成分である分析対象成分とを含む試料に対してクロマトグラフィを実施して前記基準ピークの保持時間、前記指標値、及び前記分析対象成分のピーク保持時間を取得する取得部と、前記記憶部に記憶された前記相関関係と、前記取得部が取得した前記基準ピークの保持時間及び前記指標値とに基づいて前記閾値を求め、前記閾値と前記分析対象成分のピーク保持時間とを比較して前記分析対象成分のピークが前記第１ピークであるか前記第２ピークであるかを識別する識別部と、を備える。

更に、上記目的を達成するために、本開示の一態様に係る分離分析装置は、分析カラムに導入する試料を調製するための試料調製ユニットと、前記分析カラムを有し、前記分析カラムの充填剤に対する前記試料中の成分の吸着及び脱着を制御するための分析ユニットと、前記分析カラムからの脱着液に含まれる成分を光学的に検出するための測光ユニットと、前記試料調製ユニット、前記分析ユニット、及び前記測光ユニットが各々接続された、上記制御装置と、を備える。

更に、上記目的を達成するために、本開示の一態様に係る分離分析方法は、予め準備された複数の試料に対してクロマトグラフィを実施して基準成分の基準ピークの保持時間、第１成分の第１ピークの保持時間と第２成分の第２ピークの保持時間との閾値、及び、試料中の成分の量を表す指標値の相関関係を取得し、前記基準成分と、前記第１成分又は前記第２成分である分析対象成分とを含む試料に対してクロマトグラフィを実施して前記基準ピークの保持時間、前記指標値、及び前記分析対象成分のピーク保持時間を取得し、前記相関関係と、前記取得した前記基準ピークの保持時間及び前記指標値とに基づいて前記閾値を求め、前記閾値と前記分析対象成分のピーク保持時間とを比較して前記分析対象成分のピークが前記第１ピークであるか前記第２ピークであるかを識別する。

更に、上記目的を達成するために、本開示の一態様に係る分離分析プログラムは、予め準備された複数の試料に対してクロマトグラフィを実施して基準成分の基準ピークの保持時間、第１成分の第１ピークの保持時間と第２成分の第２ピークの保持時間との閾値、及び、試料中の成分の量を表す指標値の相関関係を取得し、前記基準成分と、前記第１成分又は前記第２成分である分析対象成分とを含む試料に対してクロマトグラフィを実施して前記基準ピークの保持時間、前記指標値、及び前記分析対象成分のピーク保持時間を取得し、前記相関関係と、前記取得した前記基準ピークの保持時間及び前記指標値とに基づいて前記閾値を求め、前記閾値と前記分析対象成分のピーク保持時間とを比較して前記分析対象成分のピークが前記第１ピークであるか前記第２ピークであるかを識別することを、コンピュータに実行させる。

以上説明したように、本開示によれば、クロマトグラムのピークの保持時間が近い複数のピークについて、試料濃度を調整することなく、ピークの識別を行うことができる、という効果が得られる。

高速液体クロマトグラフィを利用したＨＰＬＣ装置の概略構成を示す図である。 ＨＰＬＣ装置の制御系の構成例を示すブロック図である。 クロマトグラムの一例を示す図である。 ＨｂＡ０、ＨｂＥ、及びＨｂＤを含むクロマトグラムの一例を示す図である。 クロマトグラムを用いて作成したピーク保持時間の相関の一例を示す散布図である。 実施形態に係る制御装置の機能的な構成の一例を示すブロック図である。 実施形態に係る分離分析プログラムによる相関関係導出処理の流れの一例を示すフローチャートである。 （Ａ）～（Ｄ）は、全ピーク面積毎のピーク保持時間の相関の一例を示す散布図である。 （Ａ）～（Ｄ）は、図８（Ａ）～図８（Ｄ）の各散布図から得られる各近似式の一例を示す図である。 （Ａ）は、傾きと全ピーク面積（指標値）との近似式の一例を示す図である。（Ｂ）は、切片と全ピーク面積（指標値）との近似式の一例を示す図である。 実施形態に係る分離分析プログラムによるピーク同定処理の流れの一例を示すフローチャートである。

以下、図面を参照して、本開示の技術を実施するための形態の一例について詳細に説明する。なお、動作、作用、機能が同じ働きを担う構成要素及び処理には、全図面を通して同じ符号を付与し、重複する説明を適宜省略する場合がある。各図面は、本開示の技術を十分に理解できる程度に、概略的に示してあるに過ぎない。よって、本開示の技術は、図示例のみに限定されるものではない。また、本実施形態では、本開示と直接的に関連しない構成や周知な構成については、説明を省略する場合がある。

図１は、高速液体クロマトグラフィ（Ｈｉｇｈ Ｐｅｒｆｏｒｍａｎｃｅ Ｌｉｑｕｉｄ Ｃｈｒｏｍａｔｏｇｒａｐｈｙ：ＨＰＬＣ）を利用したＨＰＬＣ装置Ｘの概略構成を示す図である。このＨＰＬＣ装置Ｘは、分離分析装置の一例である。

ＨＰＬＣ装置Ｘは、採血管１１をセットして、全血中のグリコヘモグロビン（ＨｂＡ１ｃ）の濃度を自動で測定するように構成されたものである。このＨＰＬＣ装置Ｘは、複数の溶離液ボトル１２Ａ、１２Ｂ、１２Ｃ、１２Ｄ、１２Ｅ（図１では５個）等を含む装置本体２を備えている。

各溶離液ボトル１２Ａ～１２Ｅは、後述する分析カラム６０に供給すべき溶離液Ａ～Ｅを各々保持したものである。各溶離液は、用途に応じて例えば組成、成分比、ｐＨ、浸透圧等が異なる。

装置本体２は、試料調製ユニット５、分析ユニット６、及び測光ユニット７を有している。

採血管１１は、例えば、ラック（図示省略）に収納され、後述する試料調製ユニット５におけるノズル５１により採取可能な位置に移動するように構成されている。

試料調製ユニット５は、採血管１１から採取した血液から、分析カラム６０に導入する試料を調製するためのものである。この試料調製ユニット５は、ノズル５１、及び希釈槽５３を有している。

ノズル５１は、採血管１１の血液試料１３をはじめとする各種の液体を採取するためのものであり、液体の吸引・吐出が可能であるとともに、上下方向及び水平方向に移動可能とされている。このノズル５１の動作は、後述する制御装置１００によって制御される。

分析ユニット６は、分析カラム６０の充填剤に対する生体成分の吸着・脱着をコントロールし、各種の生体成分を測光ユニット７に供するためのものである。分析ユニット６における設定温度は、例えば４０℃程度とされる。分析カラム６０は、試料中のヘモグロビンを選択的に吸着させるための充填剤を保持させたものである。充填剤としては、例えばメタクリル酸－メタクリル酸エステル共重合体が使用される。

分析ユニット６は、分析カラム６０の他に、マニホールド６１、送液ポンプ６２、及びインジェクションバルブ６３を有している。

マニホールド６１は、複数の溶離液ボトル１２Ａ～１２Ｅのうちの特定の溶離液ボトルから、分析カラム６０に選択的に溶離液を供給させるためのものである。このマニホールド６１は、配管８０Ａ～８０Ｅを介して溶離液ボトル１２Ａ、１２Ｂ、１２Ｃ、１２Ｄ、１２Ｅに各々接続され、配管８４を介してインジェクションバルブ６３に接続されている。

送液ポンプ６２は、溶離液をインジェクションバルブ６３に移動させるための動力を付与するためのものであり、配管８４の途中に設けられている。

インジェクションバルブ６３は、一定量の導入用試料を採取するとともに、その導入用試料を分析カラム６０に導入可能とするものであり、複数の導入ポート及び排出ポート（図示省略）を備えている。このインジェクションバルブ６３には、インジェクションループ６４が接続されている。このインジェクションループ６４は、一定量（例えば数μＬ）の液体を保持可能なものであり、インジェクションバルブ６３を適宜切り替えることにより、インジェクションループ６４が希釈槽５３と連通して希釈槽５３からインジェクションループ６４に導入用試料が供給される状態、インジェクションループ６４がプレフィルターＰＦ及び配管８５を介して分析カラム６０と連通してインジェクションループ６４から導入用試料が分析カラム６０に導入される状態を選択することができる。このようなインジェクションバルブ６３としては、例えば六方バルブを使用することができる。なお、プレフィルターＰＦは、試料や溶離液を濾過するためのフィルターである。

測光ユニット７は、分析カラム６０からの脱着液に含まれるヘモグロビンを光学的に検出するためのものであり、配管８７を介して、分析カラム６０からの脱着液を排出するための廃液槽８８に接続されている。

図２は、ＨＰＬＣ装置Ｘの制御系の構成例を示すブロック図である。

図２に示すように、ＨＰＬＣ装置Ｘは、制御装置１００を備えている。制御装置１００は、ＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）１００Ａ、ＲＯＭ（Ｒｅａｄ Ｏｎｌｙ Ｍｅｍｏｒｙ）１００Ｂ、ＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）１００Ｃ、及び入出力インターフェース（Ｉ／Ｏ）１００Ｅがバス１００Ｆを介して各々接続された構成となっている。また、ＨＰＬＣ装置Ｘは、オペレータからの入力を受け付ける操作部（図示省略）を備えている。

Ｉ／Ｏ１００Ｅには、試料調製ユニット５、分析ユニット６、及び測光ユニット７が接続されている。

分析カラム６０から排出される各種ヘモグロビンを含む溶離液は、配管８６を介して測光ユニット７に供給される。この溶離液は、配管８７を介して廃液槽８８に導かれる。

測光ユニット７においては、溶離液に対して連続的に光が照射され、その受光結果（吸光度）が制御装置１００に出力される。そして、制御装置１００においてクロマトグラムが演算される。

図３は、クロマトグラムの一例を示す図である。図３において、横軸は測定開始からの経過時間を示し、縦軸は吸光度を示す。

図３に示すように、クロマトグラムは、測定を開始してからの経過時間と受光結果としての吸光度との関係を表すグラフとして表される。このクロマトグラムのピークがどの位置に現れるかによって、どのヘモグロビン（図３の例では、ＨｂＡ０、ＨｂＤ、ＨｂＳ、ＨｂＣ）が検出されたかを知ることができると共に、ピーク部分（山なりの部分）における吸光度の積算値、すなわちピーク部分の面積の大きさによってヘモグロビンの量（濃度）を知ることができる。ＨｂＡ０は、各種ヘモグロビンの中で最もピークが大きく、識別が容易であるため、基準となる成分である。なお、測定開始時点から各ヘモグロビンのピークが出現する時点までの経過時間をピーク保持時間（リテンションタイム）という。

なお、本実施形態に係る液体クロマトグラフィ測定処理を実行するための分離分析プログラムは、本実施形態では一例としてＲＯＭ１００Ｂに予め記憶される。ＣＰＵ１００Ａは、ＲＯＭ１００Ｂに記憶された分離分析プログラムをＲＡＭ１００Ｃに書き込んで実行する。また、ＣＤ－ＲＯＭ（Ｃｏｍｐａｃｔ Ｄｉｓｃ Ｒｅａｄ Ｏｎｌｙ Ｍｅｍｏｒｙ）等の記憶媒体に分離分析プログラムを記憶し、これをＣＤ－ＲＯＭドライブ等で読み込むことにより実行するようにしてもよい。

分離分析プログラムは、例えば、制御装置１００に予めインストールされていてもよい。分離分析プログラムは、不揮発性の非一時的記憶媒体に記憶して、又はネットワーク回線を介して配布して、制御装置１００に適宜インストールしたり、アップグレードしたりすることで実現してもよい。なお、不揮発性の非一時的記憶媒体の例としては、ＣＤ-ＲＯＭ、光磁気ディスク、ＨＤＤ（Ｈａｒｄ Ｄｉｓｋ Ｄｒｉｖｅ）、ＤＶＤ-ＲＯＭ（Ｄｉｇｉｔａｌ Ｖｅｒｓａｔｉｌｅ Ｄｉｓｃ Ｒｅａｄ Ｏｎｌｙ Ｍｅｍｏｒｙ）、フラッシュメモリ、メモリカード等が想定される。

ところで、陽イオン交換クロマトグラフィでＨｂＤを分析すると、ＨｂＡ０のピーク保持時間、ＨｂＥのピーク保持時間、及び、ＨｂＤのピーク保持時間は、一例として、図４に示すようになる。

図４は、ＨｂＡ０、ＨｂＥ、及びＨｂＤを含むクロマトグラムの一例を示す図である。図４において、横軸は時間を示し、縦軸は吸光度を示す。

図４に示すように、ＨｂＡ０のピークは容易に識別可能であるが、ＨｂＥのピーク保持時間と、ＨｂＤのピーク保持時間とは近いことが分かる。このため、ピーク保持時間だけでＨｂＤのピークとＨｂＥのピークとを識別するのは難しい。

例えば、ＨｂＡ０とＨｂＥを含む試料、及び、ＨｂＡ０とＨｂＤを含む試料を準備し、各試料を様々な希釈濃度に調整した上で各調整試料を測定し、測定して得られた各クロマトグラムに現れるＨｂＡ０ピークのピーク保持時間とＨｂＤピークのピーク保持時間、及び、ＨｂＡ０ピークのピーク保持時間とＨｂＥピークのピーク保持時間をクロマトグラムごとにプロットした散布図を作成すると、図５に示すようになる。図５から、ＨｂＡ０のピーク保持時間とＨｂＥのピーク保持時間、及び、ＨｂＡ０のピーク保持時間とＨｂＤのピーク保持時間が相関することが読み取れる。

図５において、横軸はＨｂＡ０のピーク保持時間であるＴｒＡ０を示し、縦軸はＨｂＥ又はＨｂＤのピーク保持時間であるＴｒＸを示す。図５は、クロマトグラムの各ピークの面積を表す値を全て合計した合計値（全ピークの合計面積を表す値。以下、「Ｔａｒｅａ」ともいう。）が所定範囲（例えば、０～１３０，０００）を示すクロマトグラムを用いて作成された散布図である。また、図５において、ＨｂＡ０のピーク保持時間とＨｂＥのピーク保持時間との相関を第１相関といい、ＨｂＡ０のピーク保持時間とＨｂＤのピーク保持時間との相関を第２相関という。

図５に示すように、ＨｂＥとＨｂＤとでは一定の傾向はあるものの、矢印で示すように重複する部分があり、ＨｂＡ０の保持時間との相関のみでは識別することは困難である。換言すると、ＨｂＡ０ピークのピーク保持時間に基づいてＨｂＥピークのピーク保持時間とＨｂＤピークのピーク保持時間との閾値を定めても、ＨｂＥピークとＨｂＤピークを識別することは困難である。原因としては、試料中に含まれるヘモグロビンの量によってピーク保持時間がずれるため、ピーク保持時間が近いＨｂＥとＨｂＤとでは識別が困難になると考えられる。

これに対して、本発明者は、試料中のヘモグロビンの量とＨｂＡ０ピークのピーク保持時間に基づいてＨｂＥピークのピーク保持時間とＨｂＤピークのピーク保持時間との閾値を定めれば、ピーク保持時間が近いＨｂＥのピークとＨｂＤのピークを識別できることを見出した。以下、一例として、ＨｂＡ０を基準成分とし、ＨｂＥを分析対象成分の第１成分とし、ＨｂＤを分析対象成分の第２成分として説明する。この場合、ＨｂＡ０のピークを基準ピークといい、ＨｂＥのピークを第１ピークといい、ＨｂＤのピークを第２ピークという。なお、分析対象とする試料は、ヘモグロビンを含む試料に限定されるものではなく、ピーク保持時間が比較的近い成分を含む試料であれば、適用することが可能とされる。

本実施形態に係る制御装置１００のＣＰＵ１００Ａは、ＲＯＭ１００Ｂに記憶されている分離分析プログラムをＲＡＭ１００Ｃに書き込んで実行することにより、図６に示す各部として機能する。

図６は、本実施形態に係る制御装置１００の機能的な構成の一例を示すブロック図である。

図６に示すように、本実施形態に係る制御装置１００のＣＰＵ１００Ａは、導出部１０１、取得部１０２、及び識別部１０３として機能する。

導出部１０１は、予め準備された複数の試料に対してクロマトグラフィを実施して得られた、ＨｂＡ０の基準ピークの保持時間である基準ピーク保持時間、ＨｂＥの第１ピークの保持時間である第１ピーク保持時間とＨｂＤの第２ピークの保持時間である第２ピーク保持時間との閾値、及び、試料中のヘモグロビンの量を表す指標値により定まる相関関係を導出する。相関関係は、相関式として表してもよいし、相関テーブルとして表してもよい。また、指標値には、クロマトグラムにおける、全ピークの面積（＝Ｔａｒｅａ）、ＨｂＥの第１ピークの面積、ＨｂＤの第２ピークの面積、ＨｂＡ０の基準ピークの面積、ＨｂＥの第１ピークの高さ、ＨｂＤの第２ピークの高さ、及び、ＨｂＡ０の基準ピークの高さのいずれかが用いられる。これらの指標値は、クロマトグラムから得られる。

相関式又は相関テーブルは、例えば、基準ピーク保持時間と第１ピーク保持時間との相関を表す第１相関、基準ピーク保持時間と第２ピーク保持時間との相関を表す第２相関、及び、指標値により定められる。つまり、導出部１０１は、指標値毎にクロマトグラムを分類し、各クロマトグラムの第１ピーク保持時間と第２ピーク保持時間をプロットして作成した各指標値毎の散布図に表れる第１相関及び第２相関から、基準ピーク保持時間、指標値、及び分析対象成分のピークの識別に用いる閾値を変数とし、分析対象成分のピークの識別に用いる閾値を導出するための相関式又は相関テーブルを導出する。導出部１０１により導出された相関式又は相関テーブルはＲＯＭ１００Ｂに記憶される。なお、制御装置１００が基準ピーク保持時間、指標値、及び分析対象成分のピークの識別に用いる閾値を変数とする相関式又は相関テーブルを記憶したＲＯＭ１００Ｂを備えていれば、導出部１０１を備える必要は無い。例えば、相関式又は相関テーブルを予め記憶させたＲＯＭを用いて制御装置１００を製造する場合、制御装置１００の構成から導出部１０１を省くことができる。ここで、相関式の一例を下記の式（１）（以下、「相関式（１）」という。）に示す。なお、この相関式（１）を生成する具体的な処理については後述する。

ｙ＝（ｃ×Ｔａｒｅａ＋ｄ）×ＴｒＡ０＋（ｅ×Ｔａｒｅａ＋ｆ） ・・・（１）

但し、ｙは閾値、Ｔａｒｅａは全ピーク面積（指標値）、ＴｒＡ０は基準ピークであるＨｂＡ０のピーク保持時間（基準ピーク保持時間）を示す。ｃ、ｄ、ｅ、ｆは、全ピーク面積（指標値）毎の第１相関及び第２相関から得られる定数である。

取得部１０２は、分析対象とする試料に対してクロマトグラフィを実施してクロマトグラムを取得する。具体的には、試料調製ユニット５、分析ユニット６、及び測光ユニット７を制御してクロマトグラフィを実施する。そして、測光ユニット７が出力する光学測定値である吸光度を取得し、取得した吸光度と測定開始からの経過時間を用いてクロマトグラムを作成する。取得部１０２は、クロマトグラムから、基準ピーク保持時間、全ピーク面積（指標値）、及び分析対象成分のピーク保持時間を取得する。ここでいう試料は、基準成分の一例であるＨｂＡ０と、分析対象成分の一例であるＨｂＥ又はＨｂＤを含むものとする。但し、ＨｂＥ及びＨｂＤのどちらを含むかは未知である。

識別部１０３は、ＲＯＭ１００Ｂに予め記憶された相関関係から、取得部１０２で取得されたクロマトグラムから得られる分析対象成分のピークがＨｂＥの第１ピークであるかＨｂＤの第２ピークであるかを識別する。

具体的に、識別部１０３は、一例として、上記相関式（１）を用いて、取得部１０２で取得されたクロマトグラムから得られる基準ピーク保持時間及び全ピーク面積（指標値）に対応する閾値を求め、求めた閾値と、クロマトグラムから得られる分析対象成分のピーク保持時間とを比較することにより、分析対象成分のピークが第１ピークであるか第２ピークであるかを識別する。ＨｂＥとＨｂＤの場合、分析対象成分のピーク保持時間が閾値より小さい場合、ＨｂＥと同定され、分析対象成分のピーク保持時間が閾値以上である場合、ＨｂＤと同定される。

次に、図７～図１１を参照して、本実施形態に係る制御装置１００の作用を説明する。

図７は、本実施形態に係る分離分析プログラムによる相関関係導出処理の流れの一例を示すフローチャートである。なお、図７の例では、ＨｂＥとＨｂＤの２つのピークを識別するための上記相関式（１）を生成する処理について具体的に説明する。

分離分析プログラムによる相関関係導出処理は、制御装置１００のＣＰＵ１００Ａが、ＲＯＭ１００Ｂに記憶されている分離分析プログラムをＲＡＭ１００Ｃに書き込むことで実行される。

図７のステップＳ１０１では、ＣＰＵ１００Ａが、上記相関式（１）の生成に用いる複数のクロマトグラムを取得する。

ここで、クロマトグラムの作成に用いる溶離液及び検体の一例について具体的に説明する。溶離液としては、以下の溶離液Ａ、Ｂ、Ｃを準備した。ヘモグロビンの溶出力は、溶離液Ａ＜溶離液Ｃ＜溶離液Ｂ、である。

溶離液Ａ：原料として、濃度１．４４ｗｔ％のリン酸二水素ナトリウム二水和物、濃度０．１６ｗｔ％のリン酸水素二ナトリウムを含み、ｐＨ５．０８に調整。
溶離液Ｂ：原料として、濃度０．０２ｗｔ％のリン酸二水素ナトリウム二水和物、濃度０．５０ｗｔ％のリン酸水素二ナトリウムを含み、ｐＨ８．０に調整。
溶離液Ｃ：原料として、濃度０．１２ｗｔ％のリン酸二水素ナトリウム二水和物、濃度０．３３ｗｔ％のリン酸水素二ナトリウムを含み、ｐＨ６．８２に調整。

また、検体としては、ＨｂＡ０及びＨｂＤを含みＨｂＥを含まない血液検体を２２検体準備し、ＨｂＡ０及びＨｂＥを含みＨｂＤを含まない血液検体を１９検体準備した。

次に、上述のＨＰＬＣ装置Ｘを用いて、上記で準備した血液検体を以下の手順で複数回測定する。

（Ｓ１）装置本体２に溶離液Ａを流して分析カラム６０を平衡化する。
（Ｓ２）溶血した所定量の血液検体を分析カラム６０に導入する。
（Ｓ３）溶離液Ａを所定時間（例えば、１３秒）流す。
（Ｓ４）溶離液Ａと溶離液Ｃを所定割合（例えば、１：９）で混合した液を所定時間（例えば、５秒）流してＨｂＡ０を溶出する。
（Ｓ５）溶離液Ｃを所定時間（例えば、１７秒）流してＨｂＡ２、ＨｂＥ、ＨｂＤを溶出する。また、溶離液Ｂを所定時間（例えば、２秒）流して分析カラム６０に残ったヘモグロビンを全て溶出する。
（Ｓ６）溶離液Ａを所定時間（例えば、５秒）流す。
（Ｓ７）検出波長が例えば４２０ｎｍである測光ユニット７（光学検出器）により得られた吸光度からクロマトグラムを作成する。

（Ｓ８）上記クロマトグラムに現れたピークの保持時間と、ピークの大きさから、ＨｂＤに由来する第２ピーク、ＨｂＥに由来する第１ピーク、ＨｂＡ０に由来する基準ピークを同定する。具体的には、ＨｂＤに由来する第２ピークは約２５秒に出現するピークである。ＨｂＥに由来する第１ピークは約２３秒に出現するピークである。ＨｂＡ０に由来する基準ピークは約１９秒に出現するピークであり、ピークの高さ又は面積が最も大きいピークである。
（Ｓ９）ＨｂＤに由来する第２ピーク保持時間（ＴｒＨｂＤ）、ＨｂＡ０に由来する基準ピーク保持時間（ＴｒＡ０）、クロマトグラムの全ピーク面積（Ｔａｒｅａ）をＨｂＤ及びＨｂＡ０を含む検体毎に得る。
（Ｓ１０）同様に、ＨｂＥに由来する第１ピーク保持時間（ＴｒＨｂＥ）、ＨｂＡ０に由来する基準ピーク保持時間（ＴｒＡ０）、クロマトグラムの全ピーク面積（Ｔａｒｅａ）をＨｂＥ及びＨｂＡ０を含む検体毎に得る。

なお、分析カラム６０に導入されたヘモグロビン量の指標としてクロマトグラムの全ピーク面積（Ｔａｒｅａ）を用いたが、これに限定されるものではなく、分析カラム６０に導入されたヘモグロビン量の指標となる値であればよい。例えば、総ヘモグロビン量の相当量を占めるＨｂＡ０の基準ピークの面積、総ヘモグロビン量の一定量を占める変異ヘモグロビンであるＨｂＤの第２ピークの面積又はＨｂＥの第１ピークの面積を用いることができる。また、ピークの面積の代わりにピークの高さを用いてもよい。また、全ピーク面積（Ｔａｒｅａ）と相関性のある特定のピークの面積から全ピーク面積（Ｔａｒｅａ）を推測することで間接的に得るようにしてもよい。

ステップＳ１０２では、ＣＰＵ１００Ａが、ステップＳ１０１で取得した複数のクロマトグラムを指標値（例えば、全ピーク面積Ｔａｒｅａ）の大きさごとに分類する。具体的には、一例として、図８（Ａ）～図８（Ｄ）に示すように、全ピーク面積Ｔａｒｅａ毎にクロマトグラムを分類し各クロマトグラムの第１ピーク保持時間と第２ピーク保持時間を全ピーク面積Ｔａｒｅａ毎にプロットすることで、全ピーク面積Ｔａｒｅａ毎に散布図を作成する。

図８（Ａ）～図８（Ｄ）は、全ピーク面積Ｔａｒｅａ毎のピーク保持時間の相関の一例を示す散布図である。図８（Ａ）～図８（Ｄ）において、横軸はＴｒＡ０を示し、縦軸はＴｒＸ（ＴｒＨｂＤ及びＴｒＨｂＥ）を示す。

図８（Ａ）は、全ピーク面積Ｔａｒｅａが第１範囲（例えば、９０，００１～１３０，０００）を示すクロマトグラムの第１相関及び第２相関を分類した散布図（試料１）である。図８（Ｂ）は、全ピーク面積Ｔａｒｅａが第２範囲（例えば、６０，００１～９０，０００）を示すクロマトグラムの第１相関及び第２相関を分類した散布図（試料２）である。図８（Ｃ）は、全ピーク面積Ｔａｒｅａが第３範囲（例えば、３０，００１～６０，０００）を示すクロマトグラムの第１相関及び第２相関を分類した散布図（試料３）である。図８（Ｄ）は、全ピーク面積Ｔａｒｅａが第４範囲（例えば、０～３０，０００）を示すクロマトグラムの第１相関及び第２相関を分類した散布図（試料４）である。なお、全ピーク面積Ｔａｒｅａの範囲の例示は、クロマトグラムから得られた面積を表す値である。また、図８（Ａ）の散布図に用いたクロマトグラムのＴａｒｅａの平均値は１１７，５１７を示し、図８（Ｂ）の散布図に用いたクロマトグラムのＴａｒｅａの平均値は７７，８２６を示し、図８（Ｃ）の散布図に用いたクロマトグラムのＴａｒｅａの平均値は３８，９３６を示し、図８（Ｄ）の散布図に用いたクロマトグラムのＴａｒｅａの平均値は１１，０２７を示す。

図８（Ａ）～図８（Ｄ）に示すように、いずれの散布図においてもＨｂＥとＨｂＤとを明確に識別可能であることが分かる。つまり、全ピーク面積Ｔａｒｅａと基準ピーク保持時間ＴｒＡ０に応じて、第１ピーク保持時間ＴｒＨｂＥと第２ピーク保持時間ＴｒＨｂＤとの閾値を設定すれば、ＨｂＥに由来する第１ピークと、ＨｂＤに由来する第２ピークとを識別できると考えられる。

ステップＳ１０３では、ＣＰＵ１００Ａが、基準ピーク保持時間ＴｒＡ０を変数とする近似式を導出する。具体的には、一例として、上述の図８（Ａ）～図８（Ｄ）の散布図毎に第１相関と第２相関とを識別可能な近似式を導出する。この近似式の一例を下記の式（２）に示す。この近似式は、第１相関式の一例である。

ｙ＝ａ×ＴｒＡ０＋ｂ ・・・（２）

但し、ｙは第１ピーク保持時間ＴｒＨｂＥと第２ピーク保持時間ＴｒＨｂＤとの閾値、ＴｒＡ０は基準ピーク保持時間を示す。ａは傾きを示す定数、ｂは切片を示す定数である。

傾きａ、切片ｂは、一例として、図９（Ａ）～図９（Ｄ）に示すように、閾値ｙと基準ピーク保持時間ＴｒＡ０との一次関数の式である式（２）を適宜調整することにより設定される。

図９（Ａ）～図９（Ｄ）は、図８（Ａ）～図８（Ｄ）の各散布図から得られる各近似式の一例を示す図である。

図９(Ａ)に示す近似式Ｄ１は、図８（Ａ）の散布図（Ｔａｒｅａの平均値＝１１７，５１７）から得られ、例えば、以下の通りである。

ｙ＝２．４×ＴｒＡ０－２２

図９(Ｂ)に示す近似式Ｄ２は、図８（Ｂ）の散布図（Ｔａｒｅａの平均値＝７７，８２６）から得られ、例えば、以下の通りである。

ｙ＝２．３×ＴｒＡ０－２０．１

図９(Ｃ)に示す近似式Ｄ３は、図８（Ｃ）の散布図（Ｔａｒｅａの平均値＝３８，９３６）から得られ、例えば、以下の通りである。

ｙ＝２．２×ＴｒＡ０－１８．５

図９(Ｄ)に示す近似式Ｄ４は、図８（Ｄ）の散布図（Ｔａｒｅａの平均値＝１１，０２７）から得られ、例えば、以下の通りである。

ｙ＝２．１２×ＴｒＡ０－１７．３５

上述の近似式Ｄ１～Ｄ４について、Ｔａｒｅａの平均値、傾きａ、及び切片ｂの対応関係をまとめると以下の表１のようになる。

（表１）

なお、第１ピーク保持時間ＴｒＨｂＥと第２ピーク保持時間ＴｒＨｂＤとを識別可能な関数であれば一次関数でなくてもよい。例えば、二次関数、三次関数の係数を適宜調整することで、第１ピーク保持時間ＴｒＨｂＥと第２ピーク保持時間ＴｒＨｂＤとを識別可能な関数を導出してもよい。

また、本例では、第１ピーク保持時間ＴｒＨｂＥと第２ピーク保持時間ＴｒＨｂＤとを識別可能な一次関数を得るために、傾きａ、切片ｂを適宜調整することで設定した。別の方法として、例えば、ＨｂＥの第１ピーク保持時間ＴｒＨｂＥとＨｂＡ０の基準ピーク保持時間ＴｒＡ０との第１近似式を、下記の式で求める。

ＴｒＨｂＥ＝ｇ×ＴｒＡ０＋ｈ

また、ＨｂＤの第２ピーク保持時間ＴｒＨｂＤとＨｂＡ０の基準ピーク保持時間ＴｒＡ０との第２近似式を、下記の式で求める。

ＴｒＨｂＤ＝ｊ×ＴｒＡ０＋ｋ

そして、第１近似式の傾きｇと第２近似式の傾きｊとの平均値を、傾き「（ｇ＋ｊ）／２」とし、第１近似式の切片ｈと第２近似式の傾きｋとの平均値を、切片「（ｈ＋ｋ）／２」とする一次関数を得るようにしてもよい。この方法でも、同様に、第１ピーク保持時間ＴｒＨｂＥと第２ピーク保持時間ＴｒＨｂＤとを識別可能な一次関数を導出することができる。

次に、ステップＳ１０４では、ＣＰＵ１００Ａが、全ピーク面積Ｔａｒｅａ（指標値）、基準ピーク保持時間ＴｒＡ０、及び第１ピークと第２ピークの識別に用いる保持時間の閾値を変数とする相関式（１）を導出する。具体的には、一例として、図１０（Ａ）及び図１０（Ｂ）に示すように、傾きａと全ピーク面積Ｔａｒｅａ（指標値）との近似式、切片ｂと全ピーク面積Ｔａｒｅａ（指標値）との近似式をそれぞれ作成する。これらの近似式は、第２相関式の一例である。

図１０（Ａ）は、傾きａと全ピーク面積Ｔａｒｅａ（指標値）との近似式の一例を示す図である。図１０（Ｂ）は、切片ｂと全ピーク面積Ｔａｒｅａ（指標値）との近似式の一例を示す図である。

図１０（Ａ）に示す近似式は、上述の図９（Ａ）～図９（Ｄ）の散布図の各々から得られる（Ｔａｒｅａの平均値、傾きａ）の組をプロットして得られたものである。つまり、点Ｐ１は（１１７５１７、２．４）、点Ｐ２は（７７８２６、２．３）、点Ｐ３は（３８９３６、２．２）、点Ｐ４は（１１０２７、２．１２）を示す。図１０（Ａ）に示す近似式は、以下の式（３）により表される。

ａ＝ｃ×Ｔａｒｅａ＋ｄ ・・・（３）

但し、傾きｃは定数であり、例えば、ｃ＝２．６×１０^－６である。切片ｄは定数であり、例えば、ｄ＝２．１である。

同様に、図１０（Ｂ）に示す近似式は、上述の図９（Ａ）～図９（Ｄ）の散布図の各々から得られる（Ｔａｒｅａの平均値、切片ｂ）の組をプロットして得られたものである。つまり、点Ｑ１は（１１７５１７、－２２）、点Ｑ２は（７７８２６、－２０．１）、点Ｑ３は（３８９３６、－１８．５）、点Ｑ４は（１１０２７、－１７．３５）を示す。図１０（Ｂ）に示す近似式は、以下の式（４）により表される。

ｂ＝ｅ×Ｔａｒｅａ＋ｆ ・・・（４）

但し、傾きｅは定数であり、例えば、ｅ＝－４．３×１０^－５である。切片ｆは定数であり、例えば、ｆ＝－１６．８である。

上述の式（２）に、式（３）及び式（４）を代入すると、上述の相関式（１）が得られる。つまり、

ｙ＝（ｃ×Ｔａｒｅａ＋ｄ）×ＴｒＡ０＋（ｅ×Ｔａｒｅａ＋ｆ）
＝（２．６×１０^－６×Ｔａｒｅａ＋２．１）×ＴｒＡ０＋（－４．３×１０^－５×Ｔａｒｅａ－１６．８） ・・・（１）

が導出される。ｙは第１ピーク保持時間ＴｒＨｂＥと第２ピーク保持時間ＴｒＨｂＤとの閾値、Ｔａｒｅａは全ピーク面積の面積値、ＴｒＡ０は基準ピーク保持時間を示す。ｃ、ｄ、ｅ、ｆは定数である。

なお、上記では、分析カラム６０へのヘモグロビンの導入量の違い以外の原因（例えば、環境温度の変化など）による閾値の変動の指標となる基準ピークとして血液検体に含まれるＨｂＡ０成分に由来する基準ピークを用いて、基準ピーク保持時間ＴｒＡ０を閾値ｙの相関式の関数に用いた。

クロマトグラフィにより検出可能であり、かつ、検体に含まれる成分に由来するピークであれば、ＨｂＡ０以外のピークを利用してもよい。

また、本例では、全ピーク面積Ｔａｒｅａの大きさごとにクロマトグラムを分けて、基準ピーク保持時間ＴｒＡ０を横軸にし、第１ピーク保持時間ＴｒＨｂＥ及び第２ピーク保持時間ＴｒＨｂＤであるＴｒＸを縦軸にした散布図を作成した。そして、閾値ｙ、基準ピーク保持時間ＴｒＡ０、及び全ピーク面積Ｔａｒｅａ（指標値）を関数とした相関式を求めた。

別の方法として、基準ピーク保持時間ＴｒＡ０の大きさごとにクロマトグラムを分けて、全ピーク面積Ｔａｒｅａを横軸にし、第１ピーク保持時間ＴｒＨｂＥ及び第２ピーク保持時間ＴｒＨｂＤであるＴｒＸを縦軸にした散布図を作成してもよい。この場合、閾値ｙ、基準ピーク保持時間ＴｒＡ０、及び全ピーク面積Ｔａｒｅａ（指標値）を関数とした相関式を求めることができる。

具体的に、相関式は、全ピーク面積Ｔａｒｅａ（指標値）と第１ピーク保持時間ＴｒＨｂＥとの相関を表す第１相関、全ピーク面積Ｔａｒｅａ（指標値）と第２ピーク保持時間ＴｒＨｂＤとの相関を表す第２相関、及び、基準ピーク保持時間ＴｒＡ０により定められる。つまり、基準ピーク保持時間ＴｒＡ０の大きさごとに第１相関及び第２相関を分類した結果から、第１ピーク保持時間ＴｒＨｂＥと第２ピーク保持時間ＴｒＨｂＤとの閾値ｙと全ピーク面積Ｔａｒｅａ（指標値）との近似式（例えば、ｙ＝ｌ×Ｔａｒｅａ＋ｍ）を基準ピーク保持時間ＴｒＡ０の大きさごとに求める。基準ピーク保持時間ＴｒＡ０とｌの相関式（例えば、ｌ＝ｎ×ＴｒＡ０＋ｏ）、基準ピーク保持時間ＴｒＡ０とｍの相関式（例えば、ｍ＝ｐ×ＴｒＡ０＋ｑ）、及び、閾値ｙと全ピーク面積Ｔａｒｅａ（指標値）との近似式（例えば、ｙ＝ｌ×Ｔａｒｅａ＋ｍ）から、基準ピーク保持時間ＴｒＡ０、全ピーク面積Ｔａｒｅａ（指標値）、及び分析対象成分のピークの識別に用いる閾値ｙを変数とし、分析対象成分のピークの識別に用いる閾値ｙを導出するための相関式（ｙ＝（ｎ×ＴｒＡ０＋ｏ）×Ｔａｒｅａ＋（ｐ×ＴｒＡ０＋ｑ））を導出する。なお、上述の図１０（Ａ）及び図１０（Ｂ）に示す近似式を得るためには、例えば、基準ピーク保持時間ＴｒＡ０の平均値を採用すればよい。

また、全ピーク面積Ｔａｒｅａ（指標値）と、基準ピーク保持時間ＴｒＡ０と、第１ピーク保持時間ＴｒＨｂＥ及び第２ピーク保持時間ＴｒＨｂＤであるＴｒＸとの３軸からなる三次元プロットを作成し、閾値ｙ、基準ピーク保持時間ＴｒＡ０、及び全ピーク面積Ｔａｒｅａ（指標値）を関数とした相関式を求めてもよい。

次に、ステップＳ１０５では、ＣＰＵ１００Ａが、ステップＳ１０４で導出した相関式（１）をＲＯＭ１００Ｂに記憶し、本分離分析プログラムによる一連の相関関係導出処理を終了する。

次に、図１１を参照して、本実施形態に係る制御装置１００によるピーク同定処理を説明する。なお、当然ながら、ピーク同定処理で用いる指標値は、上述の相関関係導出処理で用いた指標値と同一である。つまり、相関関係導出処理で全ピーク面積Ｔａｒｅａが用いられていれば、ピーク同定処理でも全ピーク面積Ｔａｒｅａを用いる。

図１１は、本実施形態に係る分離分析プログラムによるピーク同定処理の流れの一例を示すフローチャートである。なお、ピーク同定処理では、ＨｂＥ又はＨｂＤが溶出される位置に未知のピークＸｐが溶出された場合に、未知のピークＸｐがＨｂＥの第１ピークであるかＨｂＤの第２ピークであるかを識別する。図１１の例では、ＨｂＥの第１ピークとＨｂＤの第２ピークを識別するための上記相関式（１）を用いて、ピークを同定する処理について具体的に説明する。

分離分析プログラムによるピーク同定処理は、制御装置１００のＣＰＵ１００Ａが、ＲＯＭ１００Ｂに記憶されている分離分析プログラムをＲＡＭ１００Ｃに書き込むことで実行される。

図１１のステップＳ１１１では、ＣＰＵ１００Ａが、分析対象とする試料に対してクロマトグラフィを実施して得られたクロマトグラムを取得する。ここでいう試料は、基準成分の一例であるＨｂＡ０と、分析対象成分の一例であるＨｂＥ又はＨｂＤを含むものとする。但し、ＨｂＥ及びＨｂＤのどちらを含むかは未知である。

ステップＳ１１２では、ＣＰＵ１００Ａが、ステップＳ１１１で取得したクロマトグラムから、基準ピーク保持時間ＴｒＡ０、分析対象成分（ＨｂＥ又はＨｂＤ）のピーク保持時間（つまり、未知のピークＸｐの保持時間）、及び、指標値として全ピーク面積Ｔａｒｅａを取得する。

ステップＳ１１３では、ＣＰＵ１００Ａが、ステップＳ１１２で取得した基準ピーク保持時間ＴｒＡ０及び全ピーク面積Ｔａｒｅａ（指標値）を、上述の相関式（１）に代入し、閾値ｙを導出する。

ステップＳ１１４では、ＣＰＵ１００Ａが、ステップＳ１１２で取得した分析対象成分（ＨｂＥ又はＨｂＤ）のピーク保持時間が閾値ｙよりも小さいか否かを判定する。分析対象成分（ＨｂＥ又はＨｂＤ）のピーク保持時間が閾値ｙよりも小さいと判定した場合（肯定判定の場合）、ステップＳ１１５に移行し、分析対象成分（ＨｂＥ又はＨｂＤ）のピーク保持時間が閾値ｙ以上であると判定した場合（否定判定の場合）、ステップＳ１１６に移行する。

ステップＳ１１５では、ＣＰＵ１００Ａが、分析対象成分のピークをＨｂＥの第１ピークと同定し、本分離分析プログラムによる一連のピーク同定処理を終了する。

一方、ステップＳ１１６では、ＣＰＵ１００Ａが、分析対象成分のピークをＨｂＤの第２ピークと同定し、本分離分析プログラムによる一連のピーク同定処理を終了する。

上記では指標値として全ピーク面積Ｔａｒｅａを用いたが、他の指標値を用いても同様にピークを識別することが可能である。

例えば、指標値としてＨｂＥの第１ピーク又はＨｂＤの第２ピークの面積を用いる場合、上述の相関関係導出処理において、基準ピーク保持時間ＴｒＡ０、及び、ＨｂＥの第１ピーク又はＨｂＤの第２ピークの面積に基づき定数ｃ、ｄ、ｅ、ｆを求め、基準ピーク保持時間ＴｒＡ０、閾値、及び、ＨｂＥの第１ピーク又はＨｂＤの第２ピークの面積を相関式の変数とすればよい。

この場合、ピーク同定処理では、クロマトグラムから基準ピーク保持時間ＴｒＡ０、分析対象成分（ＨｂＥ又はＨｂＤ）のピーク保持時間（つまり、未知のピークＸｐの保持時間）、及び、指標値として分析対象成分（ＨｂＥ又はＨｂＤ）のピークの面積（つまり、未知のピークＸｐの面積）を取得する。取得した基準ピーク保持時間ＴｒＡ０、及び、分析対象成分（ＨｂＥ又はＨｂＤ）のピークの面積（指標値）を相関式に代入し、閾値を導出する。そして、分析対象成分（ＨｂＥ又はＨｂＤ）のピーク保持時間と閾値とを比較する。

また、指標値としてＨｂＡ０の基準ピークの面積を用いる場合、上述の相関関係導出処理において、基準ピーク保持時間ＴｒＡ０、及び、ＨｂＡ０の基準ピークの面積に基づき定数ｃ、ｄ、ｅ、ｆを求め、基準ピーク保持時間ＴｒＡ０、閾値、及び、ＨｂＡ０の基準ピークの面積を相関式の変数とすればよい。

この場合、ピーク同定処理では、クロマトグラムから基準ピーク保持時間ＴｒＡ０、分析対象成分（ＨｂＥ又はＨｂＤ）のピーク保持時間（つまり、未知のピークＸｐの保持時間）、及び、指標値としてＨｂＡ０の基準ピークの面積を取得する。取得した基準ピーク保持時間ＴｒＡ０、及び、ＨｂＡ０の基準ピークの面積（指標値）を相関式に代入し、閾値を導出する。そして、分析対象成分（ＨｂＥ又はＨｂＤ）のピーク保持時間と閾値とを比較する。

また、指標値としてＨｂＥの第１ピーク又はＨｂＤの第２ピークの高さを用いる場合、上述の相関関係導出処理において、基準ピーク保持時間ＴｒＡ０、及び、ＨｂＥの第１ピーク又はＨｂＤの第２ピークの高さに基づき定数ｃ、ｄ、ｅ、ｆを求め、基準ピーク保持時間ＴｒＡ０、閾値、及び、ＨｂＥの第１ピーク又はＨｂＤの第２ピークの高さを相関式の変数とすればよい。

この場合、ピーク同定処理では、クロマトグラムから基準ピーク保持時間ＴｒＡ０、分析対象成分（ＨｂＥ又はＨｂＤ）のピーク保持時間（つまり、未知のピークＸｐの保持時間）、及び、指標値として分析対象成分（ＨｂＥ又はＨｂＤ）のピークの高さ（つまり、未知のピークＸｐの高さ）を取得する。取得した基準ピーク保持時間ＴｒＡ０、及び、分析対象成分（ＨｂＥ又はＨｂＤ）のピークの高さ（指標値）を相関式に代入し、閾値を導出する。そして、分析対象成分（ＨｂＥ又はＨｂＤ）のピーク保持時間と閾値とを比較する。

また、指標値としてＨｂＡ０の基準ピークの高さを用いる場合、上述の相関関係導出処理において、基準ピーク保持時間ＴｒＡ０、及び、ＨｂＡ０の基準ピークの高さに基づき定数ｃ、ｄ、ｅ、ｆを求め、基準ピーク保持時間ＴｒＡ０、閾値、及び、ＨｂＡ０の基準ピークの高さを相関式の変数とすればよい。

この場合、ピーク同定処理では、クロマトグラムから基準ピーク保持時間ＴｒＡ０、分析対象成分（ＨｂＥ又はＨｂＤ）のピーク保持時間（つまり、未知のピークＸｐの保持時間）、及び、指標値としてＨｂＡ０の基準ピークの高さを取得する。取得した基準ピーク保持時間ＴｒＡ０、及び、ＨｂＡ０の基準ピークの高さ（指標値）を相関式に代入し、閾値を導出する。そして、分析対象成分（ＨｂＥ又はＨｂＤ）のピーク保持時間と閾値とを比較する。

このように本実施形態によれば、試料中のヘモグロビンの量とＨｂＡ０ピークのピーク保持時間に基づいてＨｂＥピークのピーク保持時間とＨｂＤピークのピーク保持時間との閾値を定めれば、ピーク保持時間が近いＨｂＥのピークとＨｂＤのピークを識別できることに着目し、ＨｂＡ０の基準ピーク保持時間、ＨｂＥの第１ピーク保持時間、ＨｂＤの第２ピーク保持時間、及び、試料中のヘモグロビンの量を表す指標値により定まる相関関係が用いられる。これにより、ピーク保持時間が近く識別が困難なＨｂＥのピークとＨｂＤのピークを、試料濃度を調整することなく、識別することができる。

なお、本実施形態は、上述した実施の形態には限定されず、種々に変更可能である。例えば、血液中のヘモグロビン濃度を測定するためのＨＰＬＣ装置に限らず、血液以外の検体を用いる場合、ヘモグロビン濃度以外の成分を測定する場合、あるいはＨＰＬＣ装置以外の液体クロマトグラフィ装置についても本実施形態を適用することができる。

以上、実施形態に係る制御装置を例示して説明した。実施形態は、制御装置が備える各部の機能をコンピュータに実行させるためのプログラムの形態としてもよい。実施形態は、このプログラムを記憶したコンピュータが読み取り可能な非一時的記憶媒体の形態としてもよい。

その他、上記実施形態で説明した制御装置の構成は、一例であり、主旨を逸脱しない範囲内において状況に応じて変更してもよい。

また、上記実施形態で説明したプログラムの処理の流れも、一例であり、主旨を逸脱しない範囲内において不要なステップを削除したり、新たなステップを追加したり、処理順序を入れ替えたりしてもよい。

また、上記実施形態では、プログラムを実行することにより、実施形態に係る処理がコンピュータを利用してソフトウェア構成により実現される場合について説明したが、これに限らない。実施形態は、例えば、ハードウェア構成や、ハードウェア構成とソフトウェア構成との組み合わせによって実現してもよい。

Ｘ ＨＰＬＣ装置
５ 試料調製ユニット
６ 分析ユニット
７ 測光ユニット
１１ 採血管
１３ 血液試料
６０ 分析カラム
１００ 制御装置
１０１ 導出部
１０２ 取得部
１０３ 識別部

Claims (12)

1. クロマトグラフィにより試料を分離、分析する分離分析装置の制御装置であって、
   予め準備された複数の試料に対してクロマトグラフィを実施して得られた、基準成分の基準ピークの保持時間、第１成分の第１ピークの保持時間と第２成分の第２ピークの保持時間との閾値、及び、試料中の成分の量を表す指標値の相関関係が記憶された記憶部と、
   前記基準成分と、前記第１成分又は前記第２成分である分析対象成分とを含む試料に対してクロマトグラフィを実施して前記基準ピークの保持時間、前記指標値、及び前記分析対象成分のピーク保持時間を取得する取得部と、
   前記記憶部に記憶された前記相関関係と、前記取得部が取得した前記基準ピークの保持時間及び前記指標値とに基づいて前記閾値を求め、前記閾値と前記分析対象成分のピーク保持時間とを比較して前記分析対象成分のピークが前記第１ピークであるか前記第２ピークであるかを識別する識別部と、
   を備えた制御装置。
2. 前記相関関係は、前記基準ピークの保持時間と前記第１ピークの保持時間との相関を表す第１相関、前記基準ピークの保持時間と前記第２ピークの保持時間との相関を表す第２相関、前記閾値、及び、前記指標値により定められる
   請求項１に記載の制御装置。
3. 前記相関関係は、前記指標値毎に、前記第１相関及び前記第２相関から、前記基準ピークの保持時間、前記閾値、及び前記指標値を変数とし、前記分析対象成分のピークの識別に用いる閾値を導出するための相関式又は相関テーブルとして表される
   請求項２に記載の制御装置。
4. 前記相関関係は、前記指標値と前記第１ピークの保持時間との相関を表す第１相関、前記指標値と前記第２ピークの保持時間との相関を表す第２相関、前記閾値、及び、前記基準ピークの保持時間により定められる
   請求項１に記載の制御装置。
5. 前記相関関係は、前記基準ピークの保持時間毎に、前記第１相関及び前記第２相関から、前記基準ピークの保持時間、前記閾値、及び前記指標値を変数とし、前記分析対象成分のピークの識別に用いる閾値を導出するための相関式又は相関テーブルとして表される
   請求項４に記載の制御装置。
6. 前記指標値は、前記クロマトグラフィを実施して得られる、全ピークの面積、前記第１ピークの面積、前記第２ピークの面積、前記基準ピークの面積、前記第１ピークの高さ、前記第２ピークの高さ、及び、前記基準ピークの高さのいずれかである
   請求項１に記載の制御装置。
7. 前記試料は、ヘモグロビンを含み、
   前記基準成分はＨｂＡ０であり、前記第１成分はＨｂＥであり、前記第２成分はＨｂＤである
   請求項１に記載の制御装置。
8. 分析カラムに導入する試料を調製するための試料調製ユニットと、
   前記分析カラムを有し、前記分析カラムの充填剤に対する前記試料中の成分の吸着及び脱着を制御するための分析ユニットと、
   前記分析カラムからの脱着液に含まれる成分を光学的に検出するための測光ユニットと、
   前記試料調製ユニット、前記分析ユニット、及び前記測光ユニットが各々接続された、請求項１～請求項７の何れか１項に記載の制御装置と、
   を備えた分離分析装置。
9. 予め準備された複数の試料に対してクロマトグラフィを実施して基準成分の基準ピークの保持時間、第１成分の第１ピークの保持時間と第２成分の第２ピークの保持時間との閾値、及び、試料中の成分の量を表す指標値の相関関係を取得し、
   前記基準成分と、前記第１成分又は前記第２成分である分析対象成分とを含む試料に対してクロマトグラフィを実施して前記基準ピークの保持時間、前記指標値、及び前記分析対象成分のピーク保持時間を取得し、
   前記相関関係と、前記取得した前記基準ピークの保持時間及び前記指標値とに基づいて前記閾値を求め、前記閾値と前記分析対象成分のピーク保持時間とを比較して前記分析対象成分のピークが前記第１ピークであるか前記第２ピークであるかを識別する、
   分離分析方法。
10. 前記予め準備された複数の試料に対してクロマトグラフィを実施して得られる複数のクロマトグラムを前記指標値の大きさごとに分類し、
    前記第１ピークの保持時間と前記第２ピークの保持時間との閾値、及び、前記基準ピークの保持時間の第１相関式を前記指標値の大きさごとに求めることで複数の前記第１相関式を取得し、
    前記複数の第１相関式から、前記指標値、及び、前記指標値に対応する第１相関式の係数の第２相関式を求め、
    前記第１相関式及び前記第２相関式に基づいて、前記相関関係を表す相関式を取得する、
    請求項９に記載の分離分析方法。
11. 前記予め準備された複数の試料に対してクロマトグラフィを実施して得られる複数のクロマトグラムを前記基準ピークの保持時間の大きさごとに分類し、
    前記第１ピークの保持時間と前記第２ピークの保持時間との閾値、及び、前記指標値の第１相関式を前記基準ピークの保持時間の大きさごとに求めることで複数の前記第１相関式を取得し、
    前記複数の第１相関式から、前記基準ピークの保持時間、及び、前記基準ピークの保持時間に対応する第１相関式の係数の第２相関式を求め、
    前記第１相関式及び前記第２相関式に基づいて、前記相関関係を表す相関式を取得する、
    請求項９に記載の分離分析方法。
12. 予め準備された複数の試料に対してクロマトグラフィを実施して基準成分の基準ピークの保持時間、第１成分の第１ピークの保持時間と第２成分の第２ピークの保持時間との閾値、及び、試料中の成分の量を表す指標値の相関関係を取得し、
    前記基準成分と、前記第１成分又は前記第２成分である分析対象成分とを含む試料に対してクロマトグラフィを実施して前記基準ピークの保持時間、前記指標値、及び前記分析対象成分のピーク保持時間を取得し、
    前記相関関係と、前記取得した前記基準ピークの保持時間及び前記指標値とに基づいて前記閾値を求め、前記閾値と前記分析対象成分のピーク保持時間とを比較して前記分析対象成分のピークが前記第１ピークであるか前記第２ピークであるかを識別することを、
    コンピュータに実行させるための分離分析プログラム。