(57)【特許請求の範囲】
【請求項１】
一方の軸を試料の種別とし他方の軸を成分の種別とした２次元マップの各領域に、対応する試料の成分量を濃淡又はカラーで表示するクロマトグラフィー質量分析の分析結果表示方法において、
表示処理装置に、複数の試料に関する情報、前記複数の試料に対するクロマトグラフィー質量分析の結果から求められた前記２次元マップの各領域に表示すべき濃淡又はカラーに関する情報、及び前記クロマトグラフィー質量分析によって得られた前記２次元マップの各領域に対応する分析データ及び分析結果に関する情報を取り込む工程と、
前記２次元マップの各領域に表示すべき濃淡又はカラーに関する情報を用いて、表示装置に、２次元マップを表示する工程と、
前記２次元マップ上で拡大表示すべき試料の範囲及び成分の範囲の指定を受け付ける工程と、
前記２次元マップ上で、前記指定された範囲に含まれる領域を拡大表示し、残りの領域を縮小表示する工程と、
前記クロマトグラフィー質量分析によって得られた前記２次元マップの各領域に対応する分析データ及び分析結果に関する情報を用いて、前記拡大表示された各領域に同定情報の有無及び成分の対応付けの際の誤差に関する情報を重ねて表示する工程と
を有することを特徴とするクロマトグラフィー質量分析の分析結果表示方法。
【請求項２】
請求項１記載のクロマトグラフィー質量分析の分析結果表示方法において、前記２次元マップ上で拡大表示すべき試料の範囲又は成分の範囲の指定は、拡大表示すべき試料又は成分の開始位置の指定と拡大表示すべき試料又は成分の数の指定によって行われることを特徴とするクロマトグラフィー質量分析の分析結果表示方法。
【請求項３】
請求項１記載のクロマトグラフィー質量分析の分析結果表示方法において、前記成分の対応付けの際の誤差に関する情報には、質量電荷比の誤差に関する情報、及び保持時間の誤差に関する情報が含まれることを特徴とするクロマトグラフィー質量分析の分析結果表示方法。
【請求項４】
請求項１記載のクロマトグラフィー質量分析の分析結果表示方法において、前記成分の対応付けの際の誤差の基準となるデータが含まれる領域にマークを表示することを特徴とするクロマトグラフィー質量分析の分析結果表示方法。
【請求項５】
請求項１記載のクロマトグラフィー質量分析の分析結果表示方法において、前記拡大表示されている領域内の試料が選択されたとき、前記クロマトグラフィー質量分析によって得られた前記２次元マップの各領域に対応する分析データ及び分析結果に関する情報を用いて、当該試料の拡大表示されている領域の成分に関する詳細情報を前記２次元マップの外に前記成分の種別を表す軸方向に並べて表示することを特徴とするクロマトグラフィー質量分析の分析結果表示方法。
【請求項６】
請求項１記載のクロマトグラフィー質量分析の分析結果表示方法において、前記２次元マップ上の特定の領域が選択されたとき、当該領域に対応する測定データを保持時間と質量電荷比に対するイオン強度のマップとして表示し、選択された成分をそのマップ上にマーク表示することを特徴とするクロマトグラフィー質量分析の分析結果表示方法。
【請求項７】
請求項１記載のクロマトグラフィー質量分析の分析結果表示方法において、保持時間によるソートの指示が入力されたとき、前記２次元マップにおける成分の表示順序を保持時間の順に並べ替えて表示することを特徴とするクロマトグラフィー質量分析の分析結果表示方法。
【請求項８】
請求項１記載のクロマトグラフィー質量分析の分析結果表示方法において、質量電荷比によるソートの指示が入力されたとき、前記２次元マップにおける成分の表示順序を質量電荷比の順に並べ替えて表示することを特徴とするクロマトグラフィー質量分析の分析結果表示方法。
【請求項９】
請求項１記載のクロマトグラフィー質量分析の分析結果表示方法において、試料のクラスタリングの指示が入力されたとき、前記２次元マップにおける試料の表示順序を最も類似した試料が隣りに表示されるように入れ替えて表示することを特徴とするクロマトグラフィー質量分析の分析結果表示方法。
【請求項１０】
請求項１記載のクロマトグラフィー質量分析の分析結果表示方法において、成分のクラスタリングの指示が入力されたとき、前記２次元マップにおける成分の表示順序を最も類似した成分が隣りに表示されるように入れ替えて表示することを特徴とするクロマトグラフィー質量分析の分析結果表示方法。
【請求項１１】
複数の試料に関する情報、前記複数の試料に対するクロマトグラフィー質量分析の結果から求められた一方の軸を試料の種別とし他方の軸を成分の種別とした２次元マップの各領域に表示すべき濃淡又はカラーに関する情報、及びクロマトグラフィー質量分析によって得られた前記２次元マップの各領域に対応する分析データ及び分析結果に関する情報を取り込んで表示処理を行う表示処理装置と、
表示装置と、
入力装置とを有し、
前記表示処理装置は、前記２次元マップの各領域に表示すべき濃淡又はカラーに関する情報を用いて、前記表示装置に２次元マップを表示し、前記入力装置から入力された前記２次元マップ上で拡大表示すべき試料の範囲及び成分の範囲の指定に従って、前記２次元マップ上で、前記指定された範囲に含まれる領域を拡大表示すると共に残りの領域を縮小表示し、前記クロマトグラフィー質量分析によって得られた前記２次元マップの各領域に対応する分析データ及び分析結果に関する情報を用いて、前記拡大表示された各領域に同定情報の有無及び成分の対応付けの際の誤差に関する情報を重ねて表示することを特徴とするクロマトグラフィー質量分析の分析結果表示装置。
【請求項１２】
請求項１１記載のクロマトグラフィー質量分析の分析結果表示装置において、前記入力装置は、前記２次元マップ上で拡大表示すべき試料の範囲又は成分の範囲の開始位置を入力するためのスクロールバー、及び前記拡大表示すべき試料又は成分の数を入力するためのスクロールバーを前記表示装置の表示画面上に備えることを特徴とするクロマトグラフィー質量分析の分析結果表示装置。
【請求項１３】
請求項１１記載のクロマトグラフィー質量分析の分析結果表示装置において、前記成分の対応付けの際の誤差に関する情報には、質量電荷比の誤差に関する情報、及び保持時間の誤差に関する情報が含まれることを特徴とするクロマトグラフィー質量分析の分析結果表示装置。
【請求項１４】
請求項１１記載のクロマトグラフィー質量分析の分析結果表示装置において、前記成分の対応付けの際の誤差の基準となるデータが含まれる領域にマークを表示することを特徴とするクロマトグラフィー質量分析の分析結果表示装置。
【請求項１５】
請求項１１記載のクロマトグラフィー質量分析の分析結果表示装置において、前記拡大表示されている領域内の試料が選択されたとき、前記クロマトグラフィー質量分析によって得られた前記２次元マップの各領域に対応する分析データ及び分析結果に関する情報を用いて、当該試料の拡大表示されている領域の成分に関する詳細情報を前記２次元マップの外に前記成分の種別を表す軸方向に並べて表示することを特徴とするクロマトグラフィー質量分析の分析結果表示装置。
【請求項１６】
請求項１１記載のクロマトグラフィー質量分析の分析結果表示装置において、前記２次元マップ上の特定の領域が選択されたとき、当該領域に対応する測定データを保持時間と質量電荷比に対するイオン強度のマップとして表示し、選択された成分をそのマップ上にマーク表示することを特徴とするクロマトグラフィー質量分析の分析結果表示装置。
【請求項１７】
請求項１１記載のクロマトグラフィー質量分析の分析結果表示装置において、前記入力装置は保持時間によるソートを指示する入力部を有し、当該指示が入力されたとき、前記２次元マップにおける成分の表示順序を保持時間の順に並べ替えて表示することを特徴とするクロマトグラフィー質量分析の分析結果表示装置。
【請求項１８】
請求項１１記載のクロマトグラフィー質量分析の分析結果表示装置において、前記入力装置は質量電荷比によるソートを指示する入力部を有し、当該指示が入力されたとき、前記２次元マップにおける成分の表示順序を質量電荷比の順に並べ替えて表示することを特徴とするクロマトグラフィー質量分析の分析結果表示装置。
【請求項１９】
請求項１１記載のクロマトグラフィー質量分析の分析結果表示装置において、前記入力装置は試料のクラスタリングの指示を入力する入力部を有し、当該指示が入力されたとき、前記２次元マップにおける試料の表示順序を最も類似した試料が隣りに表示されるように入れ替えて表示することを特徴とするクロマトグラフィー質量分析の分析結果表示装置。
【請求項２０】
請求項１１記載のクロマトグラフィー質量分析の分析結果表示装置において、前記入力装置は成分のクラスタリングの指示を入力する入力部を有し、当該指示が入力されたとき、前記２次元マップにおける成分の表示順序を最も類似した成分が隣りに表示されるように入れ替えて表示することを特徴とするクロマトグラフィー質量分析の分析結果表示装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、クロマトグラフィー質量分析の分析結果を表示する方法及び表示装置に関する。
【背景技術】
【０００２】
液体クロマトグラフィー質量分析など、分離装置と質量分析装置を連結し、複数成分から成る試料を測定した場合、試料成分がイオンとして順次観測される。その結果、質量電荷比に対するイオン強度から成る質量スペクトルを時間（所謂、保持時間）ごとに集積した、３次元データを得ることができる。さらに、複数の試料の測定結果から、同一成分ごとにイオン強度の値を整理することにより、例えば試料を列に、成分を行に割り振ったイオン強度の表を作成することができる。この表は試料間での成分の量的変化を示しており、各種の標準化、クラスタリング、検定などの統計的な解析手法などを適用することによって、特徴的な成分を抽出する等の研究がなされている。
【０００３】
ここでは、このような複数の試料と複数の成分に対応する量を図示したもの、すなわち、横軸と縦軸に試料と成分を配置し、対応する量を色又は濃淡で示したものをヒートマップと呼ぶ。例えば、健常者と特定の疾患の患者に由来する試料の場合、この表を解析することにより、疾患の有無に関係すると予想される成分を発見できる可能性がある。また、投薬後の代謝産物の量的な変動を時間を追って解析することにより、代謝経路の研究など、応用分野が広い。
【０００４】
このような、試料と成分に対する存在量を示す図の一つの表現形態を、特許文献１の図8に見ることができる。この例では、当該データに対してクラスタリングを実施し、類似する変動パターンを示すものを隣接させて表示している。特許文献２は、個々のデータに関連するデータを、この表の中に含めることによって、有意な情報を抽出する方法を提示している。
【０００５】
なお、複数の試料をクロマトグラフィー質量分析した結果から、同一成分と予想されるイオンを選別する場合、少なくとも保持時間と質量電荷比の値に基づいて判断する。もし、観測した保持時間が測定ごとにある程度の変動を示す場合、補正を講じる必要がある。そのような補正の例が非特許文献１に示されている。
【０００６】
【特許文献１】特表2005-527904号
【特許文献２】特開2004-133903号公報
【非特許文献１】Xiang Zhang. et al. (2005) Data pre-processing in liquid chromatography-mass spectrometry based proteomics. Bioinformatics, 21, 4054-4059.
【発明の開示】
【発明が解決しようとする課題】
【０００７】
保持時間を補正したとしても、試料に含まれる成分が非常に多種類であり、保持時間や質量電荷比において適切な分離が得られない場合、同一成分と予想されるものが別であったり、別とされたものが同一であるような過誤に対するリスクが常に存在する。もし、着目する成分を抽出したとしても、同一成分であることの判定に誤りがあった場合には、その結果は意味のないものになってしまう。これまでの、ヒートマップの表示方法においては、このようなリスクに対し、簡便な評価手段を提供してはいない。
【０００８】
がんと正常組織、投薬や手術前後の血液や尿、遺伝子破壊株と野性株、表現型が異なる個体由来の組織など、多成分が混在する試料の成分を比較するケースは非常に多い。例えば、健常者と患者から血液や尿などを採取し、その成分を比較することによって、特定の疾患の初期段階においてのみ観測される成分を発見し、診断に利用しようとする試みや、特定の成分と予後の病態の因果関係を立証し、治療方針の立案に役立てようとするなど、様々なアプローチが成されている。
【０００９】
このようなアプローチにおいて、着目する成分を統計的に抽出したとしても、それまでの解析結果に過誤がないことを確認する必要がある。その検証には、同一成分であることの評価や、同定結果、各種の統計解析結果の確認などを含んでおり、多くの時間と労力を必要としている。
【００１０】
本発明は、試料と成分に対する量的な変化の解析結果を尊重しつつ、個々のデータを効率よく検証する手段の提供を目的とする。
【課題を解決するための手段】
【００１１】
本発明では、クロマトグラフィー質量分析の結果から求めた、試料と成分に対応する成分量を濃淡又はカラーで示したマップ表示において、観測したイオンの質量電荷比、保持時間、イオン量、同一成分の判定に関わる質量誤差、保持時間誤差、成分の同定に関わる同定結果やその有無、同定結果識別文字、同定の信頼性指標、同定結果に対するキーワード検索の結果などを、記号、数値、文字を用いて重ねて表示する。
【００１２】
また、試料と成分に対応する成分量を表示する個々の領域を文字の大きさによって規定することにより、観測したイオンの質量電荷比、保持時間、イオン量、同一成分の判定に関わる質量誤差、保持時間誤差、成分の同定に関わる同定結果やその有無、同定の信頼性指標、同定結果に対するキーワード検索の結果などを、試料と成分で特定される個々の領域の上下又は左右の、マップ外の領域に表記する。
【００１３】
このような表示において、一部の領域の拡大表示と、それ以外の領域の縮小表示を併用し、マウスによるスライダーの操作などの数値指定に連動し、リアルタイムに表示を更新する。
【００１４】
さらに、試料と成分を選択することにより、該当する測定データを保持時間と質量電荷比に対するイオン強度のマップとして表示し、選択された成分をそのマップ上にマーク表示し、その近傍を順次拡大表示する。
【００１５】
本発明によるクロマトグラフィー質量分析の分析結果表示方法は、一方の軸を試料の種別とし他方の軸を成分の種別とした２次元マップの各領域に、対応する試料の成分量を濃淡又はカラーで表示するであり、表示処理装置に、複数の試料に関する情報、複数の試料に対するクロマトグラフィー質量分析の結果から求められた２次元マップの各領域に表示すべき濃淡又はカラーに関する情報、及びクロマトグラフィー質量分析によって得られた２次元マップの各領域に対応する分析データ及び分析結果に関する情報を取り込む工程と、２次元マップの各領域に表示すべき濃淡又はカラーに関する情報を用いて、表示装置に、２次元マップを表示する工程と、２次元マップ上で拡大表示すべき試料の範囲及び成分の範囲の指定を受け付ける工程と、２次元マップ上で、指定された範囲に含まれる領域を拡大表示し、残りの領域を縮小表示する工程と、クロマトグラフィー質量分析によって得られた２次元マップの各領域に対応する分析データ及び分析結果に関する情報を用いて、拡大表示された各領域に同定情報の有無及び成分の対応付けの際の誤差に関する情報を重ねて表示する工程とを有する。
【００１６】
また、本発明によるクロマトグラフィー質量分析の分析結果表示装置は、複数の試料に関する情報、複数の試料に対するクロマトグラフィー質量分析の結果から求められた一方の軸を試料の種別とし他方の軸を成分の種別とした２次元マップの各領域に表示すべき濃淡又はカラーに関する情報、及びクロマトグラフィー質量分析によって得られた２次元マップの各領域に対応する分析データ及び分析結果に関する情報を取り込んで表示処理を行う表示処理装置と、表示装置と、入力装置とを有し、表示処理装置は、２次元マップの各領域に表示すべき濃淡又はカラーに関する情報を用いて、表示装置に２次元マップを表示し、入力装置から入力された２次元マップ上で拡大表示すべき試料の範囲及び成分の範囲の指定に従って、２次元マップ上で、指定された範囲に含まれる領域を拡大表示すると共に残りの領域を縮小表示し、クロマトグラフィー質量分析によって得られた２次元マップの各領域に対応する分析データ及び分析結果に関する情報を用いて、拡大表示された各領域に同定情報の有無及び成分の対応付けの際の誤差に関する情報を重ねて表示する。
【発明の効果】
【００１７】
本発明では生体に由来する試料なども想定している。このような試料においては、数千以上の成分を観測する可能性がある。このような膨大なデータを解析した場合、着目成分の数も数百に上る場合も考えられる。それらの検証作業において、次の効果を得ることができる。
【００１８】
（１）測定結果の評価
複数の試料において非常に多くの成分の量的な変動などを解析し、着目すべきものを絞り込む過程において、測定データの確認作業は非常に煩雑かつ困難な作業となる。例えば、従来から一般に採用されているマスクロマトグラムや質量スペクトルによる確認作業では、そもそも保持時間、質量電荷比、及びイオン強度の３次元で構成されるデータを、保持時間とイオン強度、又は質量電荷比とイオン強度の２次元としてしか表現できず、関連する情報を俯瞰することすら容易ではない。
【００１９】
本発明により、着目する試料及び成分に対し、該当する測定データの領域を、個々に拡大表示する操作を簡便に実施することが可能となる。さらに、異なる試料間で、同時に比較することが可能となる。
【００２０】
このことは、特定の試料にのみ観測された夾雑物の存在や、他の成分の影響によりその成分のイオン化が抑制されている状態や、装置のトラブルなどに起因する測定データの局所的な感度低下など、特異なケースをも発見できる可能性が高い。すなわち、そのような状態で下された、それが特徴的な成分であるとの判断そのものを再考する機会を与えることができる。
【００２１】
（２）同一成分の対応付けの評価
表示されたヒートマップにおいて、異なる成分を同一と判断している場合、その解析結果の信頼性は大きく損なわれる。そのような場合、あるものは特徴的な成分として認識され、あるものはその特徴が見えなくなるなど、多くの弊害が生まれる。
【００２２】
本発明では、そのような判断の際に、多くの情報を提供することができる。すなわち、質量分析装置から得た成分の対応付けとして最も基本的な、保持時間の誤差、質量電荷比の誤差などの情報を、矢印などで視覚的に表現することができるため、誤差の傾向の把握を容易にしている。保持時間の誤差や質量電荷比の誤差の傾向が異なる成分については、先に示した測定結果の評価などを適用することが可能となる。その他の情報についても、文字列などで明示するため、詳細なデータを検討することができる。
【００２３】
（３）同定結果の評価
保持時間と質量電荷比によって対応付けられた成分に対し、同定結果の有無や、同定結果そのもの、及び同定の信頼性に関する指標を表示することにより、同定又は成分の対応付けの過誤の可能性を示すことができる。特に、保持時間と質量電荷比の値から同一と判断された成分において、異なる同定結果が対応する場合、同定結果そのものの過誤を検討する必要が生じる。
【００２４】
同様にして、異なる成分と判断されていても、同定結果が同じ場合には、本来は同じ成分としなければならない可能性がある。
【００２５】
本発明による同定結果識別文字によって、同一成分とされるものに対して異なる同定結果が与えられていないか等を効率よく判断することができる。
【００２６】
（４）ヒートマップの精査
本発明では、ヒートマップの全域を表示しつつ、その一部分を拡大表示することができる。すなわち、全体の傾向を確認しつつ、着目すべき部分を検証する作業を、限られた画面上のスペースにおいて実施することができる。さらに、拡大表示する位置や表示する領域の大きさを、マウスによるスライダなどの操作と連動させて変えることによって、ヒートマップ全体の中での部分拡大の操作を効率よく実施できる。
【００２７】
このように、本発明では、複数の試料と成分の比較解析などにおいて、特徴的な成分などを検討する際に、全体の情報を維持しつつ、細部の具体的な情報を提供することが可能となる。そのため、全体を俯瞰して得た特徴的な領域について、様々な過誤の可能性を否定し、正しい解析結果を得るための有効な手段を提供することができる。
【発明を実施するための最良の形態】
【００２８】
本発明は、ヒートマップの表現方法に関するものである。以下、図面を参照しつつ、本発明を詳細に説明する。
【００２９】
１．試料特性
ヒートマップの性格上、そこに表現される試料は、多くの成分を含んでいることが前提となる場合が多い。また、成分としては、タンパク質やペプチド、及び代謝物などの他、環境中に含まれる低分子化合物なども想定されるが、限定されるものではない。
【００３０】
２．測定装置
本発明によるクロマトグラフィー質量分析装置の構成例を図１に示す。本例の装置は、液体クロマトグラフと質量分析装置を備える。
【００３１】
液体クロマトグラフ１０１で分離された試料成分は、質量分析装置本体１０２に導入され、イオン源１０３、質量分離装置１０４を経て検出器１０５で検出される。検出信号は、制御装置１０６を経て、データ処理装置１０８において質量スペクトルとして順次収集される。なお、それぞれの機器及び装置は、信号線１０７で接続され、制御及び情報収集に利用される。また、データ処理装置１０８は、表示装置１０９とキーボード１１０、及びマウス１１１を有している。
【００３２】
ここでは質量分離装置１０４として、イオントラップ形式のものを示したが、時系列的に試料成分に由来する質量スペクトルを得ることができるものであれば、磁場や四重極などを採用した装置であってもよい。
【００３３】
この構成において、液体クロマトグラフによる試料成分の分離開始に同期して、質量スペクトルを連続して取得する。さらに、質量電荷比に対するイオン強度から成る質量スペクトルを、それを観測した保持時間ごとに集積する。
【００３４】
質量分析装置の中には、生成したイオンをガス分子と衝突させることによって活性化させるなどして、イオン分子内の結合を開裂させ、より小さいイオンを生成させた後で質量スペクトル（所謂ＭＳ／ＭＳスペクトル）を得ることができるものもある。このＭＳ／ＭＳスペクトルは、成分イオンが開裂して生成されたものであるため、成分イオンの構造などに関連する情報を得ることができる。以下、ＭＳ／ＭＳスペクトルに対比し、成分イオンを観測した質量スペクトルを明示する必要がある場合は、ＭＳスペクトルと呼ぶ。また、ＭＳ／ＭＳスペクトルを得るために開裂させたイオンは、前駆イオンと呼ばれる。
【００３５】
さらに、取得したＭＳスペクトルから前駆イオンを選択し、リアルタイムにそのＭＳ／ＭＳスペクトルを測定できるものがある。このような場合、混在するＭＳスペクトルとＭＳ／ＭＳスペクトルを区別し、ＭＳ／ＭＳスペクトルから成分の定性的な情報を得ることができる場合がある。また、その近傍の保持時間において測定したＭＳスペクトルに含まれる前駆イオンの情報から、量的知見を得ることができる。
【００３６】
本発明は、ＭＳ／ＭＳスペクトルを測定するような測定にも対応できる。
なお、必ずしも液体クロマトグラムを用いる必要はなく、ガスクロマトグラムなどの分離手段と質量分析装置を連結した構成であれば、本発明を適用できる。
【００３７】
また、データ処理装置１０８の内部に実装される処理の一例として、対応成分識別処理１１２、同定処理１１３、解析処理１１４、表示処理１１５を挙げる。特に、本発明に関わるのが表示処理であり、それ以外の処理は、表示されるデータの作成に関与する。
【００３８】
３．処理の概要
前述した装置構成において、データ処理装置１０８内で実行される各処理１１２〜１１５のフローを図２に示す。この図において、前述した装置構成によって得た測定データから、本発明を適用するためのヒートマップ等を表示するまでの処理内容の一例を説明する。
【００３９】
（１）測定データファイル（２０１）
各試料を測定した結果は、測定データファイル２０１として集積される。本発明は、複数の試料を成分ごとに比較する技術に関するものであるため、ここでは複数の測定データファイルを想定する。
【００４０】
（２）対応成分識別処理（２０２）
個々の測定データファイルには、保持時間と質量電荷比によって分離された複数の成分の情報が含まれている。ここでは、異なる測定データから、同じ成分を抽出し、後の比較を可能とする。なお、同一成分であることの判定は、基本的には、イオンを観測した保持時間と質量電荷比が、指定した誤差範囲に収まることを以って行うことができる。
【００４１】
また、判定するイオンが有する価数の一致や、同位体イオンの強度分布の類似性を評価することや、後述する同定処理の結果が一致することなど、成分が対応していると見なす判定要素は複数存在する。
【００４２】
（３）同定処理（２０３）
同定処理は、測定データにおいて個々の成分に相当するＭＳスペクトル中のイオンの質量電荷比の値や、そのイオンを前駆イオンとして得たＭＳ／ＭＳスペクトルなどから、それがどのようなものであったかを予測するものである。予測結果は、化学式や、アミノ酸の配列などとして取得する。
【００４３】
（４）統合データ（２０４）
上述した対応成分識別処理２０２や同定処理２０３の結果は、この例においては、統合データ２０４として一元的にまとめている。統合データの構造は、後に詳述する図３のデータ構造に近い。図３は、ヒートマップの各点を説明するための情報を示しているが、その原形がここで作成される。
【００４４】
すなわち、統合データ２０４には、図３における個々の試料（Ｓ）と成分（Ｃ）の情報が存在し、さらに特定の試料と成分に対応する具体的に観測されたイオンの情報等（Ｐ）と同様の基本構造を有する。
【００４５】
（５）解析処理（２０５）
統合データは、試料と成分で特定することができるイオン量の値が基本となる。ここでの解析処理は、そのようなデータから特徴を見出すための処理であり、一般的には、試料単位でのイオン量の標準化、成分単位でのイオン量の標準化、各イオン量の多変量解析などが該当する。特に標準化では、平均を１とする方法や、平均を０、分散を１とする方法などが採用される。また、多変量解析では、主成分分析や、成分及び試料に関するクラスタリング解析が採用されることが多い。また、試料が健常者と患者由来などのように、識別できる場合には、成分ごとにそれらの識別に従って差異を認めるかどうかを検定するなどの方法も選択される。
【００４６】
これらの処理によって、観測したイオン量からヒートマップを表示するためのデータが生成され、さらにクラスタリングなどの結果に応じて試料や成分を表示する順番が入れ替わるなどの影響を受ける。
【００４７】
（６）表示データ（２０６）
表示データ２０６には、本発明における表示内容に関わるデータを収載しているものとする。ただし、クラスタリングの結果については省略した。
【００４８】
表示データは具体的に図３及び図４に示した。また、その内容については、この後で詳述する。
【００４９】
なお、先の解析処理２０５において、標準化などを何もしていない場合、統合データ２０４と表示データ２０６は実質的に同一内容となる。
【００５０】
（７）表示処理（２０７）
表示処理２０７では、表示データ２０６を図６に示すヒートマップとして表示する。ここで、図６（Ａ）が通常のヒートマップの例、図６（Ｂ）が一部を拡大表示した例である。この処理が本発明に直接対応している。その詳細は後述する。
【００５１】
４．表示データ
本発明を適用する際に、想定されるデータ構造の一例を図３及び図４に示す。表示データは、ヒートマップに関係するデータを補足する基本データと、ヒートマップデータ、表示用データに分けている。
【００５２】
基本データは、図３に示すように、ヒートマップデータの補足情報配列３０１、成分情報配列３０２、試料情報配列３０３からなる。ヒートマップデータは、図４（Ａ）に示すように、ヒートマップデータの配列４０１からなる。また、表示データは、図４（Ｂ）に示すように、試料拡大表示用データ４０２、成分拡大表示用データ４０３、ヒートマップの画面座標Ｘ（試料方向）４０４、ヒートマップの画面座標Ｙ（成分方向）４０５、画面表示座標での文字高さ４０６からなる。
【００５３】
それぞれの記号の意味は次のとおりである。
（１）基本データ（図３）
Ｐj,k 添え字ｊ，ｋで限定されるヒートマップの補足情報。
測定データにおいて観測した値や、同定処理の結果等で構成される。
一例を次に挙げる。
・測定データにおいて観測したイオン量
・測定データにおいて観測した保持時間
・測定データにおいて観測した質量電荷比
・測定データにおいて観測したイオンの価数
・測定データにおいて観測したイオンの同定結果
Ｃk   添え字ｋで限定される、当該成分を代表する情報。
個々の試料に依存しない値をとり、一例は次のとおりである。
・成分を代表する保持時間
・成分を代表する質量電荷比
・成分を代表する価数
・成分の同定結果
Ｓj   添え字ｊで限定される試料を代表する情報。
個々の成分に依存しない値をとり、一例は次のとおりである。
・試料の名称
・試料に対応する測定データファイルの所在
・試料に対応する同定結果ファイルの所在
【００５４】
（２）ヒートマップデータ（図４（Ａ））
Ｈj,k 添え字ｊ，ｋで限定されるヒートマップの各値。
【００５５】
（３）表示用データ（図４（Ｂ））
ＳＰ  試料を表示する際の、拡大開始位置
ＳＺ  試料を表示する際の、拡大表示する試料数
ＣＰ  成分を表示する際の、拡大開始位置
ＣＺ  成分を表示する際の、拡大表示する成分数
Ｘ１  試料を表示する際の、表示位置左端（画面上の座標）  (固定値)
Ｘ２  試料を表示する際の、拡大表示位置左端（画面上の座標）
Ｘ３  試料を表示する際の、拡大表示位置右端（画面上の座標）
Ｘ４  試料を表示する際の、表示位置右端（画面上の座標）  (固定値)
Ｙ１  成分を表示する際の、表示位置上端（画面上の座標）  (固定値)
Ｙ２  成分を表示する際の、拡大表示位置上端（画面上の座標）
Ｙ３  成分を表示する際の、拡大表示位置下端（画面上の座標）
Ｙ４  成分を表示する際の、表示位置下端（画面上の座標）  (固定値)
ＴＨ  文字列の高さ(画面上の座標)  (固定値)
また、添え字の意味は次のとおりである。
ｊ 試料を特定する数
ｋ 成分を特定する数
Ｊ 試料数（ｊの最大値）
Ｋ 成分数（ｋの最大値）
【００５６】
ここで、基本データとヒートマップデータを表示する際に、画面上の位置を特定するための情報や、表示する範囲を規定する情報が表示用データとなる。この場合、基本データ及びヒートマップデータでは、試料を列に、成分を行に設定している。また、Ｐj,k、Ｃk、Ｓjで表される基本データは、それぞれが複数の内容を含んでいる。それらの詳細を示す場合、例えば、「Ｐj,kの保持時間」、「Ｃkの質量電荷比」、「Ｓjの試料名称」のように表記する。配列そのものを指す場合には、単にＰやＨと記述する。
【００５７】
ここで、基本データとヒートマップデータは、共通する添え字ｊとｋで限定される。すなわち、Ｈj,kとＰj,kは、共通する成分Ｃkと試料Ｓjに対応している。また、配列ＨとＰは、１対１に対応する。
【００５８】
従って、表示したヒートマップの一つの点の情報は、Ｈの配列上の一点に相当し、さらにその点に対応するデータがＰの配列上の１点に対応する。Ｐの配列上の１点は、観測したイオン強度や保持時間など、複数のデータから構成されている。
【００５９】
５．ヒートマップ表示
本発明は、ヒートマップの表示方法に関するものである。ここでは、ヒートマップを表示するための処理内容を詳述する前に、ここで表示されたヒートマップについて、図６を用いて説明する。
【００６０】
図６（Ａ）は、一般的なヒートマップの例を示す図である。この例では、横方向が試料、縦方向が成分を示している。図３の基本データ、図４のヒートマップデータとしては、試料数の最大値Ｊが８、成分数の最大値Ｋが１５の、８×１５点から成るＰ及びＨの配列が対応する。すなわち、ヒートマップの左上、右下の点に対応するデータは、それぞれＨ1,1とＨ8,15となる。
【００６１】
図６（Ｂ）は、ｋ＝７からｋ＝９の成分を拡大した例を示す図である。ここで、（Ｘ１，Ｙ１）のような表記は、ヒートマップを表示した際の画面上の座標を示す。横方向をＸ、縦方向をＹで示し、ヒートマップの上端がＹ１、下端がＹ４、拡大表示の上端がＹ２、拡大表示の下端をＹ３とした。同様に、左端をＸ１、拡大表示の左端をＸ２、拡大表示の右端をＸ３、ヒートマップの右端をＸ４としている。すなわち、左上の隅が（Ｘ１，Ｙ１）、右下の隅が（Ｘ４，Ｙ４）となる。また、この例では、全ての試料を拡大した例となっている。すなわち、Ｘ１とＸ２は同値、Ｘ３とＸ４も同値となる。
【００６２】
また、拡大表示した部分に表示される○印や矢印などの記号は、成分の対応付けの状況や、同定結果の有無などを示している。例えば○印６０１が付与された場所は、該当する試料の該当する成分に対して、それが何であるかを示す同定情報が存在していることを示している。成分が正しく対応付けられているとすれば、横方向に並んだ○印の成分において、その同定結果が一致する可能性が高い。もし、横方向で同定結果が異なる場合には、同定の過誤や成分の対応付けが正しく行われていない可能性を示唆している。
【００６３】
矢印６０３は、成分の対応付けの際の誤差を示している。ここでは、成分の対応付けについては特に限定するものではないが、少なくとも保持時間と質量電荷比がある誤差範囲に収まることを以って対応していることを判断するものとする。すなわち、基準とする保持時間と質量電荷比に対して、ヒートマップ上の各点は保持時間の誤差と質量電荷比の誤差の値を保持している。もし、該当するイオンが存在しない場合には、そのような情報や同定結果が存在しない。
【００６４】
図６（Ｂ）の例では、左右方向に質量電荷比の誤差、上下方向に保持時間の誤差を示す座標軸を設定し、ベクトルで示している。なお、矢印の根が原点であり、双方の誤差０を意味する。
【００６５】
ここでは、基準位置に相当する点があれば、●印６０２で示した。例えば、その成分において最も強く観測されたデータを基準として●印を付与するなどが想定される。
なお、基準とする位置は、当該成分において観測した保持時間と質量電荷比の平均値なども採用することができる。
【００６６】
矢印、及びマーク類の一例を図７に示す。ここでは二つの数値を図示する場合の説明と例（７０１〜７０８）、及び同定結果の有無と信頼度などを示す例（７０９〜７１１）を示す。
【００６７】
保持時間の誤差や質量電荷比の誤差などを示すためには、ヒートマップ上の１点を拡大した領域７０１の中に、それぞれを示す座標軸７０２，７０３を設定する。ここでは、それらの座標軸の交点を原点とし、保持時間は上から下、質量電荷比は左から右方向に増加する座標軸を想定した。この座標軸を用いて、該当する数値を表現した例を、７０４〜７０８に示した。７０４は、独立した二つの数値を座標軸に垂直な線分で示している。７０５では二つの数値から成るベクトルを、７０６ではベクトルの矢の位置にマークを、７０７では二つの数値を楕円の横方向の幅と縦方向の幅で示した。７０８はさらにもう二つの数を矢印のデザインが異なる二つのベクトルで示した。
【００６８】
同定結果の有無と信頼度などを示す例の丸印７０９は同定結果がある場合に表示できる。例えば、輪を示す線の幅と同定における信頼性のスコアを連動させ、信頼性の高いものを強調させる表示などが考えられる。また、同結果においてリン酸など特定の修飾が示唆されている場合、左上に“Ｐ”などのマークを付与すると、より効果的な表示７１０となる。右上の数値“７５”は、同定の信頼性スコアを示す。また、表示７１１のように、信頼性のスコアが一定以上の場合に星印を付与してもよい。
【００６９】
文字を表示するその他の例として、表示７１２のように、同定結果を識別するための１〜数文字を表示することも考えられる。この例で示した“Ｑ４Ａ”は、当該成分の同定結果が、多の同定結果と同一かどうかを識別するために利用される。すなわち、複数の試料について、同定結果を比較する際に、それらが同じかどうかを、個々の同定結果を参照することなく確認することができる。
【００７０】
このように、矢印やマーク類のデザインは多くの態様を想定できる。その元になるデータは、図３に示す基本データにおけるＰの中に含める。
【００７１】
ヒートマップの全体は、その傾向を把握する目的において有効であるが、個々の点を厳密に評価するためには殆ど情報を与えない。しかし、このように保持時間や質量電荷比の誤差、及び同定に関する情報を効果的に提示することによって、ヒートマップを作成した際の正当性などについて評価することができる。
【００７２】
６．表示処理
（１）ヒートマップ表示
前述したヒートマップ表示を行うための処理について、成分の拡大や表示の位置指定を中心に、図５を用いて説明する。なお、試料についての扱いも同様となるため、ここでは割愛した。
【００７３】
図５における、ＣＺ、ＣＰ、ＴＨ、Ｙ１、Ｙ２、Ｙ３、Ｙ４、及びＨ、Ｐ、Ｃについては、図４で説明した。ここで、Ｙ１、Ｙ４、及びＴＨは、それぞれ画面上の表示位置の上端、下端と、文字の高さであり、予め設定されているものとする。
【００７４】
拡大表示数指定操作５０１によって、成分の中で拡大表示する数を指定する。また、表示位置指定操作５０２は、成分のどの位置から拡大表示するかを指定する。それぞれの操作の結果は、ＣＺとＣＰに格納される。なお、図６の例では、拡大表示数ＣＺは３、拡大開始位置ＣＰは７となる。
【００７５】
さらに、表示座標計算処理５０３では、ＣＺ、ＣＰ、ＴＨの値と、Ｙ１、Ｙ４の値を用いて拡大表示の上端Ｙ２と下端Ｙ３を決定する。例えば拡大表示の高さを文字の２倍と固定した場合、拡大表示位置の座標には次の関係がある。
Ｙ３−Ｙ２＝２×ＴＨ×ＣＺ                …(1)
表示する成分数の最大値Ｋから、縮小して表示する成分数は次のとおりである。
縮小表示成分数＝Ｋ−ＣＺ                  …(2)
拡大表示の上に位置する縮小表示成分数は次のとおりである。
拡大表示の上の縮小表示成分数＝ＣＰ−１    …(3)
【００７６】
ヒートマップにおいて、縮小表示する領域の画面表示における合計高さは、式(1)を用いて求めることができる。
縮小表示する成分の高さ合計＝（Ｙ４−Ｙ１）−（Ｙ３−Ｙ２）
＝（Ｙ４−Ｙ１）−２×ＴＨ×ＣＺ    …(4)
さらに、縮小表示する成分あたりの高さは、次のようになる。
縮小表示する成分あたりの高さ＝縮小表示する成分の高さ合計／縮小表示成分数
＝（（Ｙ４−Ｙ１）−２×ＴＨ×ＣＺ）／（Ｋ−ＣＺ）  …(5)
従って、Ｙ２は式(2)と式(5)から、次の式で示すことができる。
Ｙ２＝Ｙ１＋縮小表示する成分あたりの高さ×拡大表示の上の縮小表示成分数
＝Ｙ１＋（（Ｙ４−Ｙ１）−２×ＴＨ×ＣＺ）／（Ｋ−ＣＺ）×（ＣＰ−１）
また、Ｙ３は式(1)を変形して次のように求めることができる。
Ｙ３＝Ｙ２＋２×ＴＨ×ＣＺ
このようにして、既知の数値から、求めるべきＹ２及びＹ３を求める。
【００７７】
ヒートマップ等表示処理５０４では、拡大表示の上の縮小表示であれば、Ｙ１とＹ２の間を、縮小して表示する成分数で案分し、等間隔に表示すればよい。同様に、拡大表示は、Ｙ２とＹ３の間を案分して表示し、拡大表示の下の縮小表示は、Ｙ３とＹ４について案分して表示する。
【００７８】
図６のヒートマップでは、ｋの値が１〜６まではＹ１とＹ２の間に等間隔に表示し、ｋが７〜９まではＹ２とＹ３の間に、ｋが１０〜１５まではＹ３からＹ４の間に表示する。
【００７９】
以上、成分の表示について説明した。試料についても、同様にして拡大する範囲や拡大開始位置、及び拡大する座標を求めることができる。
【００８０】
このようにして決定されたＸ１、Ｘ２、Ｘ３、Ｘ４、及びＹ１、Ｙ２、Ｙ３、Ｙ４の各座標に対して、同様に決定した各縮小及び拡大領域のヒートマップデータを案分し、ｊとｋによって特定された試料及び成分について、対応するＨを濃淡やカラーで表現する。
【００８１】
（２）矢印等の表示
さらに、拡大表示した部分に矢印などを表示する方法については次のとおりである。
ヒートマップ上の１点を拡大した領域７０１の中に、例えば、左から右が質量電荷比の誤差７０２、上から下が保持時間の誤差７０３の座標系を設定する。さらに領域７０１の範囲を、それぞれの座標系で−３０秒から＋３０秒、及び０．０５Ｄａから＋０．０５Ｄａなどと予め決めておくことにより、ベクトル等の描画が可能となる。また、誤差の値そのものは、次のようにして計算することができる。
保持時間の誤差＝Ｐj,kの保持時間−Ｃkの保持時間
質量電荷比の誤差＝Ｐj,kの質量電荷比−Ｃkの質量電荷比
【００８２】
従って、誤差の座標系において、矢印の根元が（０，０）、矢の先が（質量電荷比の誤差，保持時間の誤差）などとし、ベクトルを表示する。
なお、この処理は、ヒートマップ等表示処理５０４で実施する。
【００８３】
（３）同定結果識別文字の表示
同定結果の表示において、その内容を識別するための１〜数文字を表示する。その文字列の作成方法については、全ての同定結果を整理し、その種類に対して与えた序数を文字として表示する方法や、同定結果の文字列から生成する方法などが考えられる。先の方法は、同定結果が少ない場合に有効であるが、多い場合に数文字で表記しきれないなどの課題がある。一方後者は、異なる同定結果から同じ識別文字を生成する可能性がある。試料の数が数十〜数百の場合を考慮し、ここでは、後者の方法について図１３を用いて一例を示す。
【００８４】
図１３では、同定結果の文字列のコード配列Ａ（１３０１）を仮定する。なお、ここで扱う文字は７ビットで構成されるＡＳＣＩＩコードとする。文字コードからなる配列Ａの値の合計値（整数）をΣＡ（１３０２）に登録し、さらにその値を同定結果識別文字のコード配列Ｂ（整数）に展開する。
【００８５】
たとえば、同定結果が“ＩＤ７・・・０”のＸ文字からなる場合、“Ｉ”の文字コードがＡ1に相当する。同様にして“Ｄ”がＡ2、最後の“０”がＡXに対応する。さらに、これらの値の合計値であるΣＡは、Ａ1＋Ａ2＋Ａ3＋・・・＋ＡXの結果として得ることができる。このΣＡのビットパターンの右側から６ビットずつ区切り、それぞれをＢ3、Ｂ2、Ｂ1にコピーする。さらに、表示可能な文字に変換するため、例えば４８（“０”（ゼロ）の文字に対応するＡＳＣＩＩコード）を加算する。これによって、６ビットで表すことのできる数値、０〜６３が、“０”〜“ｏ”として表現される。このようにして、Ｂ1、Ｂ2、Ｂ3からなる同定結果識別文字Ｂ（１３０３）を生成する。ここでは、図７の例（７１２）に示した“Ｑ４Ａ”を、Ｂ1、Ｂ2、Ｂ3に対応する文字として例示した。
なお、この処理も、ヒートマップ等表示処理５０４で実施する。
【００８６】
前述した表示方法に基づき、実施例を挙げてさらに説明する。図８の例では、試料と成分について拡大表示を実施しているヒートマップと、試料や成分の情報を合わせて表示している。さらに、図９において関連する処理のためのボタン等を配置した画面の一例を示す。これらの画面は常に同時に表示されているものとする。
【００８７】
１．  画面の構成要素
図８及び図９の表示における主な構成要素について、個々に説明する。
（１）成分の拡大表示数指定用のスクロールバー（８０１）
ヒートマップ等の表示において、拡大表示して示す成分の数を指定する。スライダーをマウス等で上下に移動させることにより、図５におけるＣＺを更新する。さらに、表示座標計算処理５０３とヒートマップ等表示処理５０４をリアルタイムに連動させる。そうすることによって、全体を概観するヒートマップから、特定の部分の拡大表示をスムーズに連動させることができる。
【００８８】
このスライダーを上に移動した場合の例を図１０に示す。図８では４成分が拡大表示されていたが、図１０では２成分になっていることが分かる。拡大表示の開始位置を示すスライダー１００２は影響を受けないものの、デンドログラム１００６、ヒートマップ１００５、各成分の関連データのグラフ１００８、及び各成分の関連データ１００９は、拡大表示数の減少に対応して更新されている。
【００８９】
（２）成分の拡大表示位置指定用のスクロールバー（８０２）
拡大表示する位置を、成分の全体から指定する。スライダーの上下をマウス等で指定することによって、図５におけるＣＰを更新する。さらに、表示座標計算処理とヒートマップ表示処理をリアルタイムに連動させる。すなわち、上下の指定によって、リアルタイムに拡大領域を上下させることができる。
【００９０】
図１０の状態からこのスライダーを下に移動した例を図１１に示す。この例においては、拡大表示を開始する位置が下方向にずれている。この場合、拡大表示数１１０１は同じだが、デンドログラム１１０６、ヒートマップ１１０５、各成分の関連データのグラフ１１０８、及び各成分の関連データ１１０９は、拡大表示の開始位置の変更に対応して更新されている。
【００９１】
（３）試料の拡大表示数指定用のスクロールバー（８０３）
ヒートマップ等の表示において、拡大表示して示す試料の数を指定する。この指定も、ヒートマップ表示処理までリアルタイムに連動させる。
【００９２】
（４）試料の拡大表示位置指定用のスクロールバー（８０４）
拡大表示する位置を、試料の全体から指定する。これも、同様にヒートマップ表示処理までリアルタイムに連動させ、拡大表示の位置を自由にさせることを可能とする。
【００９３】
（５）ヒートマップ表示（８０５）
試料と成分で限定される個々の成分の存在量をカラー又は濃淡で表示する。この例では、試料と成分の双方について一部を拡大しているが、ヒートマップ（Ｈ）のカラー又は濃淡表示については省略している。また、拡大表示された内容は、対応するヒートマップの基本データ（Ｐ）を用いて描画している。
【００９４】
なお、拡大表示位置をマウスでドラッグ（特定のボタンを押下しつつマウスのカーソルを移動）することにより、成分と試料の双方の拡大表示位置を変更することも考えられる。これは、成分の拡大表示位置指定用のスクロールバー８０２と試料の拡大表示位置指定用のスクロールバー８０４を、マウスの一回の操作で同時に修正する操作に相当する。
【００９５】
（６）成分ツリー表示（８０６）
ヒートマップの個々の成分に関し、個々の類似性をデンドログラムで示している。なお、デンドログラムの表示と、ヒートマップの各成分は対応する。デンドグラムの詳細は省略するが、ウォード法などのクラスタリングによって、ヒートマップを成分ごとに見た場合、その内容が一致するもの、又は最も似ているものが隣どうしになるように成分の位置を定め、その距離をツリー表示に反映させたものである。こうすることにより、ヒートマップの全体を見た場合、類似する特性を示す成分が近傍に配置されるため、より全体の特徴を把握しやすくなる。
【００９６】
（７）試料ツリー表示（８０７）
ヒートマップの個々の試料に関し、個々の類似性をデンドログラムで示している。なお、デンドログラムの表示と、ヒートマップの各試料は対応する。
【００９７】
（８）各成分の関連データのグラフ（８０８）
ヒートマップの成分に対応する情報であり、選択されている試料ｊに対応する、Ｈj,kやＰj,k又はＣkの数値をプロットしたものである。ここではＨj,kの例を示しており、右側ほど値が大きい。なお、選択されている試料は、イタリック体で示しており、ここでは“Sample Ｄ”が対応する。
【００９８】
（９）各成分の関連データ（８０９）
ヒートマップの拡大表示した成分に対応する情報であり、選択されている試料ｊに対応する、Ｈj,kやＰj,k又はＣkの値を表記したもの。
【００９９】
この例では、「Ｃkの保持時間」、「Ｃkの質量電荷比」、「Ｐj,kの同定結果」を表示している。なお、選択されている試料は、イタリック体で示しており、ここでは“Sample Ｄ”が対応する。
【０１００】
このように、より具体的な情報を提供することによって、同一成分と考えていたものでも、同定結果が異なる場合を簡便に示すことができる。また、同定結果のみならず、ＰやＣとして保持している全ての内容をここに表示することができる。
【０１０１】
（１０）各試料の関連データのグラフ（８１０）
ヒートマップの試料に対応する情報であり、選択されている成分ｋに対応する、Ｈj,kやＰj,k又はＳｋの数値をプロットしたもの。ここではＨj,kの例を示しており、上側ほど値が大きい。なお、ここで選択されている成分とは、イタリック体で示した“Protein Ｂ”を意味している。
【０１０２】
（１１）各試料の関連データ（８１１）
ヒートマップの拡大表示した試料に対応する情報であり、選択されている成分ｋに対応する、Ｈj,kやＰj,k又はＳｋの値を表記したものである。この例では、「Ｓjの試料名称」を示している。
【０１０３】
以上が、図８における構成要素の説明となる。さらに、図９の画面が同時に表示されている状態を想定し、個々の構成要素について説明する。
【０１０４】
（１２）試料のクラスタリングボタン（９０１）
試料についてクラスタリング処理を実施するためのボタン。このボタンの指定により、試料についてのクラスタリングを行い、ヒートマップやデンドログラム等を描画する。クラスタリングによって、Ｈ、Ｐ、Ｓの各データは、対応する添え字ｊごとに入れ替えを実施し、整合性を維持する。
なお、試料のクラスタリングは、図２における解析処理に相当する。
【０１０５】
（１３）成分のクラスタリングボタン（９０２）
成分についてクラスタリング処理を実施するためのボタン。このボタンの指定により、成分についてのクラスタリングを行い、ヒートマップやデンドログラム等を描画する。クラスタリングによって、Ｈ、Ｐ、Ｃの各データは、対応する添え字ｋごとに入れ替えを実施し、整合性を維持する。
なお、成分のクラスタリングは、図２における解析処理に相当する。
【０１０６】
（１４）保持時間ソートボタン（９０３）
成分について、保持時間の順に整列するためのボタン。
「Ｃkの保持時間」が昇順となるように、Ｈ、Ｐ、Ｃの各データを対応する添え字ｋごとに入れ替える。なお、成分についてのデンドログラムは消去される。
なお、保持時間によるソートは、図２における解析処理に相当する。
【０１０７】
また、保持時間ごとに成分をソートすることによって、ある保持時間の領域に依存する成分の対応付けの特性などを確認し易くなる。
【０１０８】
（１５）質量電荷比ソートボタン（９０４）
成分について、質量電荷比の順に整列するためのボタン。
「Ｃkの質量電荷比」が昇順となるように、Ｈ、Ｐ、Ｃの各データを対応する添え字ｋごとに入れ替える。なお、成分についてのデンドログラムは消去される。
なお、質量電荷比によるソートは、図２における解析処理に相当する。
【０１０９】
また、このようなソート機能によって、質量電荷比に依存する成分の対応付けの特性などを確認し易くなる。
【０１１０】
（１６）キーワード入力領域（９０５）
キーワード検索する文字列を入力するための領域。次に説明するキーワード検索で利用する。
【０１１１】
（１７）キーワード検索ボタン（９０６）
設定されたキーワードが同定結果の中に含まれるか検索するためのボタン。
このボタンがマウス等で指定されると、キーワード入力領域に設定された文字列が、全ての同定情報に含まれているかを検索し、該当するものがある場合に、そのデータを選択状態とする。選択状態は、イタリック体で示したが、色を変える、マークを表示するなど、他の方法で明示してもよい。
【０１１２】
（１８）測定データ表示ボタン（９０７、９０８）
選択された一つ又は複数の成分について、選択された試料の測定データと合わせて表示するためのボタン。該当するボタンは、２Ｄ（Ｌ）（９０７）と、２Ｄ（Ｒ）（９０８）の二種類を示している。これは、図１２の表示例における、左側と右側の更新を意味する。
【０１１３】
なお、図１２の画面では、縦軸が保持時間、横軸が質量電荷比から成る２次元マップ１２０１において、測定データを概観することができる。ここで、観測されたイオンは濃淡又はカラー表示１２０４で明示する。また、図８の画面で選択された成分に対応する保持時間と質量電荷比は、○印１２０５で示している。すなわち、測定データ全域の中で、着目する成分が具体的にどこに位置しているのかを一望することができる。
【０１１４】
また、２次元マップ１２０１内の矩形１２０６は、拡大表示している範囲を示している。この矩形の中央に位置する保持時間に相当する質量スペクトル１２０２と、マスクロマトグラム又はトータルイオンクロマトグラム等（１２０３）が同時に表示されている。なお、それぞれ（１２０２、１２０３）の表示範囲は、２次元マップ１２０１の表示範囲に一致させる。
【０１１５】
図１２の下側には、拡大表示位置１２０６に対応する２次元マップ１２０７を表示している。拡大表示位置の中央に相当する質量スペクトル１２１０と、マスクロマトグラム１２０８により、具体的なイオンの増減や夾雑物の存在などを確認することができる。また、それぞれ（１２０８、１２１０）の表示範囲は、２次元マップ１２０７に一致させる。
【０１１６】
文字表示領域１２０９には、拡大表示位置に対応するタンパク質の同定結果などを表示する。これにより、同定結果と各イオンの観測状況などを合わせて確認することが可能となる。
【０１１７】
また、図１２は、左右対称に構成されている。すなわち、二つの試料について、それぞれ全体を俯瞰するマップ１２０１，１２１２や拡大したマップ１２０７，１２１３などがそれぞれ対応する。また、それらを繋ぐ線分１２１１は、左右に表示された成分の対応関係を示し、同一成分がそれぞれ対応する保持時間を、直線で結んだものとなる。
【０１１８】
このように、ボタンが指定された時点で選択されている試料と成分について、図１２の画面の左側、又は右側の情報、及び中央の対応関係を更新する。
【０１１９】
また、キーボードの矢印キーの押下などによって、拡大表示する位置を、全ての○印に対応し、順次移動させる。すなわち、双方の試料の選択された成分の測定状況を、簡便な操作で確認することができる。
【０１２０】
２．  操作の一例
図８、図９及び図１２の画面の例を用いて、基本的な操作と処理について説明する。
（１）試料、成分の並べ替え
試料のクラスタリングボタン９０１、成分のクラスタリング９０２、保持時間のソートボタン９０３、質量電荷比のソートボタン９０４をマウス等で指定することによって、ヒートマップ内の試料や成分を並べ替える。クラスタリングを指定した場合、該当するデンドログラムが表示される。
【０１２１】
ここで、例えば成分のクラスタリングを採用した場合、類似する傾向を示す成分を近傍に集めることができる。また、質量電荷比に基づいてソートされたヒートマップでは、質量電荷比の誤差の傾向を把握しやすくし、さらに保持時間が近い同一成分の可能性のあるものを確認することなどに寄与する。保持時間に基づいてソートした場合、成分の対応付けにおいて保持時間に依存する現象を確認し易い。
【０１２２】
（２）試料や成分の選択
表示されている成分や試料の情報をマウスで直接指定することにより、当該成分又は当該試料を選択状態として、強調表示する。また、ヒートマップ内をクリックすることによってその成分や試料を選択状態にする。同様にして、ドラッグした場合、当該範囲に相当する試料や成分が、複数選択される。
【０１２３】
また、キーワードを設定（９０５）し、検索ボタン９０６を指定することによって、指定された文字列を「Ｐj,kの同定結果」に含む成分が選択状態となる。
この例では、太字のイタリック体で表示している。
【０１２４】
選択された試料や成分は、図１２に示すような測定データの確認手段によって、検証することができる。さらに、選択された成分を、上下の矢印キーなどの操作によって順次拡大表示すると、より効果的な検証が可能となる。
【０１２５】
（３）試料や成分の拡大表示
拡大表示は、表示する試料や成分の数を指定するスクロールバー８０１，８０３と、表示位置を指定するスクロールバー８０２，８０４で指定する。これらのスクロールバーの操作によって、リアルタイムに画面を更新する。その際、ツリー表示６０６，６０７、グラフ８０８，８１０、文字情報８０９，８１１も、ヒートマップの更新に同期する。
なお試料や成分の拡大表示は、この例の場合２行の文字領域に固定している。
【０１２６】
このような拡大表示の連携は、前述した座標系の計算によって実現できる。すなわち、Ｘ１、Ｘ２、Ｘ３、Ｘ４、及びＹ１、Ｙ２、Ｙ３、Ｙ４で示す座標は、ツリー表示やグラフ、文字情報においても、それぞれの表示の位置を定めるために利用できる。
【０１２７】
（４）測定データの表示
測定データ表示ボタン９０７、９０８は、選択された試料と成分について、測定データを表示する。
【０１２８】
例えば、特定の成分が選択された状態において、“Sample Ｄ”の部分をクリックし当該サンプルを選択する。続いて“２Ｄ（Ｌ）”ボタンをクリックし、図１２の左側のマップを表示する。さらに、“Sample Ｆ”の部分をクリックし、“２Ｄ（Ｒ）”ボタンのクリックで図１２の右側のマップ表示及び対応する成分を直線で結ぶ。
【０１２９】
なお、図１２のマップとして表示されるデータは、選択された試料jにおいて、「Ｓjの試料に対応する測定データファイルの所在」にファイルとして登録されており、図２の２０１に相当する。
【０１３０】
また、選択された試料ｊ１とｊ２において対応する成分ｋは、「Ｐj1,kの保持時間」と「Ｐj1,kの質量電荷比」によって特定できる成分と、「Ｐj2,kの保持時間」と「Ｐj2,kの質量電荷比」によって特定できる成分の対応として規定できる。これらの値を、選択された全ての成分について抽出し、１２０５の○印を与える。さらに、「Ｐj1,kの保持時間」と「Ｐj2,kの保持時間」を直線で結んだものが、１２１１で示した保持時間の対応表示に相当する。対応する成分が多い場合には、「Ｐj1,kのイオン量」や「Ｐj2,kのイオン量」が大きいものに着目して表示するなど、一部の対応関係を間引いて表示することも考えられる。
【０１３１】
このようにして、図１２のような表示を、異なる試料について表示し、対応する成分について比較することにより、成分の対応付けの良し悪しを評価することが可能となる。すなわち、このような測定時の状況を直接観測する手段に戻ることによって、ケミカルノイズの影響や、他のイオンの影響、夾雑物の影響などによる過誤を、効率よく判断することができる。
【０１３２】
ここで、図１２を用いて、夾雑物の影響と判断できる一例を示す。図１２の左右のマップは比較すべき成分のものとする。ここで拡大表示されているマップ１２０７，１２１３を比較すると、濃淡で示す領域が異なることが分かる。また、質量スペクトル表示１２１０，１２１４も異なっている。この場合、右の試料において、保持時間と質量がわずかに異なる成分も溶出している。さらに、その成分の影響により、イオン強度の値が実際のものよりも大きく観測されているものと推測できる。
【０１３３】
もし、図１２における○印１２０５が複数ある場合、例えば矢印キーの押下などによって拡大する位置を換える。連続して矢印キーを押下することにより、多くの箇所を連続して確認することを可能とする。
【０１３４】
このように、同定された成分を名称などから選択し、ヒートマップ上で該当する部分を拡大表示しながら、個々の成分の対応関係などを確認することができる。さらに、同一クラスタに含まれる同定されていない成分なども、保持時間の誤差や質量電荷比の誤差などの情報を提供することができる。従って、成分が同定されていないものについても、同定された成分から、それが着目されるべきものであることを判断する。また、クラスタリングなどによって特徴的な成分を検出した場合、その成分の測定データを簡便な操作で確認することができる。
【０１３５】
本発明によると、ヒートマップの全体を俯瞰し、特徴的な細部を効率よく確認することが可能となるため、着目する成分の検証作業を支援することができる。
【図面の簡単な説明】
【０１３６】
【図１】本発明によるクロマトグラフィー質量分析装置の構成例を示す図。
【図２】データ処理装置内で実行される処理のフローの例を示す図。
【図３】表示に用いられるデータ構造の一例を示す図。
【図４】ヒートマップのデータ構造の一例を示す図。
【図５】拡大表示の処理に関するフローの例を示す図。
【図６】本発明を適用した表示の一例を示す図。
【図７】本発明を適用した表示の他の例を示す図。
【図８】本発明による表示画面の一例を示す図。
【図９】処理を指示するボタン等を配置した画面の一例を示す図。
【図１０】本発明による表示画面の一例を示す図。
【図１１】本発明による表示画面の一例を示す図。
【図１２】測定データを確認する画面の一例を示す図。
【図１３】同定結果識別文字を説明する図。
【符号の説明】
【０１３７】
１０１    液体クロマトグラフ
１０２    質量分析装置本体
１０３    イオン源
１０４    質量分離装置
１０５    検出器
１０６    制御装置
１０７    信号線
１０８    データ処理装置
１０９    表示装置
１１０    キーボード
１１１    マウス
１１２    対応成分識別処理
１１３    同定処理
１１４    解析処理
１１５    表示処理
２０１    測定データファイル
２０２    対応成分識別処理
２０３    同定処理
２０４    統合データ
２０５    解析処理
２０６    表示データ
２０７    表示処理
３０１    ヒートマップデータの補足情報配列
３０２    成分情報配列
３０３    試料情報配列
４０１    ヒートマップデータの配列
４０２    試料拡大表示用データ
４０３    成分拡大表示用データ
４０４    ヒートマップの画面座標Ｘ（試料方向）
４０５    ヒートマップの画面座標Ｙ（成分方向）
４０６    画面表示座標での文字高さ
５０１    拡大表示数指定操作
５０２    表示位置指定操作
５０３    表示座標計算処理
５０４    ヒートマップ等表示処理
６０１    同定情報マーク
６０２    成分基準マーク
６０３    成分評価ベクトル
７０１    座標系説明図
７０２    質量電荷比の誤差の座標軸
７０３    保持時間の誤差の座標軸
７０４    表示例１
７０５    表示例２
７０６    表示例３
７０７    表示例４
７０８    表示例５
７０９    同定結果表示例１
７１０    同定結果表示例２
７１１    同定結果表示例３
７１２    同定結果表示例４
８０１    成分の拡大表示数指定用のスクロールバー
８０２    成分の拡大表示位置指定用のスクロールバー
８０３    試料の拡大表示数指定用のスクロールバー
８０４    試料の拡大表示位置指定用のスクロールバー
８０５    ヒートマップ表示
８０６    成分のクラスタリング結果
８０７    試料のクラスタリング結果
８０８    各成分の関連データのグラフ
８０９    各成分の関連データ
８１０    各試料の関連データのグラフ
８１１    各試料の関連データ
９０１    試料のクラスタリングボタン
９０２    成分のクラスタリングボタン
９０３    成分の保持時間ソートボタン
９０４    成分の質量電荷比ソートボタン
９０５    キーワード入力領域
９０６    キーワー検索ドボタン
９０７    測定データ表示ボタン（左側）
９０８    測定データ表示ボタン（右側）
１００１    成分の拡大表示数指定用のスクロールバー
１００２    成分の拡大表示位置指定用のスクロールバー
１００５    ヒートマップ表示
１００６    成分のクラスタリング結果
１００８    各成分の関連データのグラフ
１００９    各成分の関連データ
１１０１    成分の拡大表示数指定用のスクロールバー
１１０２    成分の拡大表示位置指定用のスクロールバー
１１０５    ヒートマップ表示
１１０６    成分のクラスタリング結果
１１０８    各成分の関連データのグラフ
１１０９    各成分の関連データ
１２０１    マップ表示
１２０２    質量スペクトル表示
１２０３    クロマトグラム表示
１２０４    観測イオン表示
１２０５    成分の位置表示
１２０６    拡大範囲表示
１２０７    拡大マップ表示
１２０８    拡大クロマトグラム表示
１２０９    文字情報表示
１２１０    質量スペクトル表示
１２１１    保持時間の対応表示
１２１２    マップ表示（右）
１２１３    拡大マップ表示（右）
１２１４    質量スペクトル表示（右）
１３０１    同定結果の文字列コード配列
１３０２    同定結果の文字列の合計値
１３０３    同定結果代表コード
【図１】
0004860575.tif 000002
【図２】
0004860575.tif 000003
【図３】
0004860575.tif 000004
【図４】
0004860575.tif 000005
【図５】
0004860575.tif 000006
【図６】
0004860575.tif 000007
【図７】
0004860575.tif 000008
【図８】
0004860575.tif 000009
【図９】
0004860575.tif 000010
【図１０】
0004860575.tif 000011
【図１１】
0004860575.tif 000012
【図１２】
0004860575.tif 000013
【図１３】
0004860575.tif 000014