(57)【要約】  代謝物や代謝酵素の測定値データについて、専門家による分析検討を支援する支援装置、支援方法、及びコンピュータプログラムを提供する。  座標取得部３０１２は、代謝パスウェイマップ上における代謝物、代謝酵素の位置座標を取得する。関連性強度算出部３０１５は、測定値データ格納部３０４に格納された代謝物等の測定値データから、各々の代謝パスウェイの関連性強度を算出する。表示内容決定部３０１７は、選択入力受付部３０１６が受け付けた代謝パスウェイの選択入力等に基づき、代謝パスウェイマップ上において、座標取得部３０１２が取得した座標の位置に存在する代謝物等の測定値の増減が識別可能であるように表示内容を決定し、表示装置１０２に表示する。
2011055820.tif 000002
【特許請求の範囲】
【請求項１】
複数種類の代謝物及び／又は代謝酵素の各々について、異なる条件で得られた複数の測定値データを取得する測定値データ取得手段と、
生体内の代謝パスウェイマップを表すデータを取得するパスウェイマップデータ取得手段と、
前記測定値データ取得手段が取得した測定値データに基づいて、前記代謝物及び／又は代謝酵素と、代謝パスウェイ各々との関連性の強度を算出する関連性強度算出手段と、
表示手段と、
前記表示手段への表示内容を決定する表示内容決定手段とを備え、
前記表示内容決定手段は、
前記表示手段に、前記関連性の強度が識別可能であるように、二種以上の代謝パスウェイを表す内容を表示させることが可能であり、
さらに、
前記表示手段に前記二種以上の代謝パスウェイを表す内容が表示された場合に、ユーザからの、いずれかの代謝パスウェイの選択入力を受け付けるパスウェイ選択入力受付手段を備え、
前記表示内容決定手段は、
前記パスウェイ選択入力受付手段が受け付けた選択入力により選択された代謝パスウェイについて、前記表示手段に代謝パスウェイマップを表示させるとともに、当該代謝パスウェイマップ上に存在する代謝物及び／又は代謝酵素について、前記測定値データが示す測定値の増減が識別可能であるように、前記表示手段への表示内容を決定する
ことを特徴とする支援装置。
【請求項２】
前記関連性強度算出手段は、
前記測定値データにより表される測定値の増減を加味して、前記代謝物及び／又は代謝酵素と、代謝パスウェイ各々との関連性の強度を算出する
ことを特徴とする請求項１に記載の支援装置。
【請求項３】
関連性強度の算出に加味される増減の判定に用いる測定値の範囲指定入力を受け付ける範囲選択入力受付手段を備える
ことを特徴とする請求項２に記載の支援装置。
【請求項４】
前記関連性強度算出手段は、
測定値の増減の起こりにくさに基づいて代謝パスウェイごとの関連性強度を算出する
ことを特徴とする請求項２又は３に記載の支援装置。
【請求項５】
前記関連性強度算出手段は、
測定値の増減から算出される代謝パスウェイごとの自己エントロピー又はシャノン・エントロピーを関連性の強度として算出する
ことを特徴とする請求項２から４のいずれかに記載の支援装置。
【請求項６】
前記関連性強度算出手段は、
前記測定値データ取得手段が測定値データを取得した代謝物及び／又は代謝酵素の種類の数に基づいて、前記代謝物及び／又は代謝酵素と、代謝パスウェイ各々との関連性の強度を算出する
ことを特徴とする請求項１に記載の支援装置。
【請求項７】
前記関連性の強度を、フィッシャーの正確確率検定により算出する
ことを特徴とする請求項６に記載の支援装置。
【請求項８】
前記関連性強度算出手段は、
フィッシャーの正確確率検定により算出された関連性の強度を、ＦＤＲを用いて、さらに補正して関連性強度を算出する
ことを特徴とする請求項７に記載の支援装置。
【請求項９】
前記支援装置は、さらに、
代謝パスウェイマップ上に表示される代謝物及び／又は代謝酵素のマップ上の位置情報を取得する座標取得手段を備える
ことを特徴とする請求項１から８のいずれかに記載の支援装置。
【請求項１０】
前記表示内容決定手段は、
測定値が増加したか減少したかにより、前記パスウェイマップ上において前記代謝物及び／又は代謝酵素を表す部分の色を変化させる
ことを特徴とする請求項１から９のいずれかに記載の支援装置。
【請求項１１】
前記表示内容決定手段は、
測定値の増加量、減少量に基づき、前記パスウェイマップ上において前記代謝物及び／又は代謝酵素を表す部分の大きさを変化させる
ことを特徴とする請求項１から１０のいずれかに記載の支援装置。
【請求項１２】
前記支援装置は、さらに、
前記代謝物及び／又は代謝酵素を表す部分の色、大きさを変化させるための増加量、減少量を算出する測定値の範囲を指定する第２の範囲選択入力受付手段を備える
ことを特徴とする請求項１０又は１１に記載の支援装置。
【請求項１３】
前記表示内容決定手段は、
前記パスウェイマップ上において前記代謝物及び／又は代謝酵素を表す部分の近傍に、測定値の増減の様子を表す増減パターンを表示させる
ことを特徴とする請求項１から１２のいずれかに記載の支援装置。
【請求項１４】
前記増減パターンにおいて、測定値の増減の様子の判定に閾値を用いる
ことを特徴とする請求項１３に記載の支援装置。
【請求項１５】
前記閾値の入力を受け付ける閾値入力受付手段を備える
ことを特徴とする請求項１４に記載の支援装置。
【請求項１６】
前記閾値として、
前記増減パターンを表す線の始点の位置決定に用いる第１の閾値と、始点の位置が決定された後の増減パターンの生成に用いる第２の閾値とが設定される
ことを特徴とする請求項１４又は１５に記載の支援装置。
【請求項１７】
前記閾値として、測定値との比較に用いる閾値が設定される
ことを特徴とする請求項１４又は１５に記載の支援装置。
【請求項１８】
前記閾値として、前記代謝物及び／又は代謝酵素ごとに異なった閾値を用いる
ことを特徴とする請求項１４から１７のいずれかに記載の支援装置。
【請求項１９】
前記増減パターンの生成に用いる測定値の範囲指定入力を受け付けるパターン表示範囲選択入力受付手段を備える
ことを特徴とする請求項１３から１８のいずれかに記載の支援装置。
【請求項２０】
前記代謝パスウェイマップ上に表示された前記代謝物及び／又は代謝酵素の選択入力を受け付ける代謝物選択入力受付手段を備え、
前記表示内容決定手段は、
選択された代謝物及び／又は代謝酵素について、測定値又は測定値を表すグラフを表示するように表示内容を決定する
ことを特徴とする請求項１３から１９のいずれかに記載の支援装置。
【請求項２１】
表示手段が接続されたコンピュータを用いて、複数種類の代謝物及び／又は代謝酵素の各々について、異なる条件で得られた複数の測定値データから、複数種類の代謝パスウェイ各々への影響の分析、検討を支援する支援方法であって、
生体内の代謝パスウェイマップを表すデータを取得するパスウェイマップデータ取得ステップと、
前記測定値データに基づいて、前記代謝物及び／又は代謝酵素と、代謝パスウェイ各々との関連性の強度を算出する関連性強度算出ステップと、
前記表示手段への表示内容を決定する第１の表示内容決定ステップとを含み、
前記第１の表示内容決定ステップにおいて、
前記表示手段に、前記関連性の強度が識別可能であるように、二種以上の代謝パスウェイを表す内容を表示させることが可能であり、
さらに、
前記表示手段に前記二種以上の代謝パスウェイを表す内容が表示された場合に、ユーザからの、いずれかの代謝パスウェイの選択入力を受け付けるパスウェイ選択入力受付ステップと、
前記パスウェイ選択入力受付ステップにおいて受け付けた選択入力により選択された代謝パスウェイについて、前記表示手段に代謝パスウェイマップを表示させるとともに、当該代謝パスウェイマップ上に存在する代謝物及び／又は代謝酵素について、前記測定値データが示す測定値の増減が識別可能であるように、前記表示手段への表示内容を決定する第２の表示内容決定ステップを含む
ことを特徴とする支援方法。
【請求項２２】
前記関連性強度算出ステップにおいて、
前記測定値データにより表される測定値の増減を加味して、前記代謝物及び／又は代謝酵素と、代謝パスウェイ各々との関連性の強度を算出する
ことを特徴とする請求項２１に記載の支援方法。
【請求項２３】
表示手段が接続されたコンピュータ上で動作するコンピュータプログラムであって、
前記コンピュータを、
複数種類の代謝物及び／又は代謝酵素の各々について、異なる条件で得られた複数の測定値データを取得する測定値データ取得手段、
生体内の代謝パスウェイマップを表すデータを取得するパスウェイマップデータ取得手段、
前記測定値データ取得手段が取得した測定値データに基づいて、前記代謝物及び／又は代謝酵素と、代謝パスウェイ各々との関連性の強度を算出する関連性強度算出手段、及び
前記表示手段への表示内容を決定する表示内容決定手段として動作させ、
前記表示内容決定手段は、
前記表示手段に、前記関連性の強度が識別可能であるように、二種以上の代謝パスウェイを表す内容を表示させることが可能であり、
さらに、前記コンピュータプログラムは、前記コンピュータを、
前記表示手段に前記二種以上の代謝パスウェイを表す内容が表示された場合に、ユーザからの、いずれかの代謝パスウェイの選択入力を受け付けるパスウェイ選択入力受付手段として動作させ、
前記表示内容決定手段は、
前記パスウェイ選択入力受付手段が受け付けた選択入力により選択された代謝パスウェイについて、前記表示手段に代謝パスウェイマップを表示させるとともに、当該代謝パスウェイマップ上に存在する代謝物及び／又は代謝酵素について、前記測定値データが示す測定値の増減が識別可能であるように、前記表示手段への表示内容を決定する
ことを特徴とするコンピュータプログラム。
【請求項２４】
前記関連性強度算出手段は、
前記測定値データにより表される測定値の増減を加味して、前記代謝物及び／又は代謝酵素と、代謝パスウェイ各々との関連性の強度を算出する
ことを特徴とする請求項２３に記載のコンピュータプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、支援装置、支援方法、及びコンピュータプログラムに関し、特に、生体内の代謝経路における代謝物、代謝酵素等の増減状態把握による病態診断、検査、薬物投与の影響判断などを支援する支援装置、支援方法、及び当該支援装置、支援方法を実現するコンピュータプログラムに関する。
【背景技術】
【０００２】
生体内における代謝経路（代謝パスウェイ）に関する情報を管理するデータベースの一例として、京都大学化学研究所が運用するＫＥＧＧ（生命システム情報統合データベース）がある。ＫＥＧＧでは、代謝経路の分類方法の一例として、糖質代謝、脂質代謝、アミノ酸代謝等に分類しており、現時点で３００種以上の代謝経路が知られている。代謝パスウェイは、生体反応の新たな解明により、今後も増加が予想される。【０００３】
薬物などの対象化合物の投与が生体に及ぼす作用を定量的に、統合的に記述してプロファイリングすることにより、標的分子を推定し、薬効プロファイルを推定し、副作用や毒性を予測し、化合物を選定しようとする技術が特許文献１に記載されている。【０００４】
特許文献２には、候補生理活性物質に曝されていない細胞に由来する第１のサンプルのＦＴＭＳ（フーリエ変換質量スペクトル）ピークプロファイルと、当該候補生理活性物質に曝された細胞に由来する第２のサンプルのＦＴＭＳピークプロファイルとを比較する技術が開示されている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００３−３２５０７６号公報
【０００６】
【特許文献２】国際公開第２００２／００４９５７号
【発明の概要】
【発明が解決しようとする課題】
【０００７】
創薬のための試験では、生体への薬物投与の影響調査など、例えば測定対象の血中、尿中から数十から数百、ときには千種類を超える生体内化合物の濃度測定を行う。疾患診断において精密検査を行う場合でも同様である。このような場合、測定値のみからでは、代謝パスウェイへの影響を詳細に見極めることは極めて困難であるし、従来の血液検査のように測定値が基準値の範囲内にあるか否かを判定するだけでは、重篤な疾患が発見されない可能性もある。特許文献１には、対象化合物を作用させたときの観測分子の変動の結果を比較分析する点について記載があるが、最終的には、医師、研究者などの専門家による判断は不可欠である。特許文献２には、ＦＴＭＳピークプロファイルをコンピュータで比較する点や、ピークプロファイルをデータベース化するといった内容が記載されているが、コンピュータによる比較だけでは発見できない異常の発生も充分に考えられるところである。ところが、上記したように、例えば血中、尿中から得られる測定値データが多種である一方、代謝パスウェイも３００種を超える現状では、医師、研究者といった専門家でも、どの代謝パスウェイから検討を加えれば良いのか、途方にくれることになる。【０００８】
本発明は、上記のような点に鑑みてなされたものであって、代謝物、代謝酵素の測定値データについて、専門家による分析検討を支援する支援装置、支援方法、及びコンピュータプログラムを提供することを目的としている。
【課題を解決するための手段】
【０００９】
上記の問題点を解決するために、本発明に係る支援装置は、複数種類の代謝物及び／又は代謝酵素の各々について、異なる条件で得られた複数の測定値データを取得する測定値データ取得手段と、生体内の代謝パスウェイマップを表すデータを取得するパスウェイマップデータ取得手段と、前記測定値データ取得手段が取得した測定値データに基づいて、前記代謝物及び／又は代謝酵素と、代謝パスウェイ各々との関連性の強度を算出する関連性強度算出手段と、表示手段と、前記表示手段への表示内容を決定する表示内容決定手段とを備え、前記表示内容決定手段は、前記表示手段に、前記関連性の強度が識別可能であるように、二種以上の代謝パスウェイを表す内容を表示させることが可能であり、さらに、前記表示手段に前記二種以上の代謝パスウェイを表す内容が表示された場合に、ユーザからの、いずれかの代謝パスウェイの選択入力を受け付けるパスウェイ選択入力受付手段を備え、前記表示内容決定手段は、前記パスウェイ選択入力受付手段が受け付けた選択入力により選択された代謝パスウェイについて、前記表示手段に代謝パスウェイマップを表示させるとともに、当該代謝パスウェイマップ上に存在する代謝物及び／又は代謝酵素について、前記測定値データが示す測定値の増減が識別可能であるように、前記表示手段への表示内容を決定することを特徴としている。【００１０】
本発明によれば、測定値データに基づいて、代謝物及び／又は代謝酵素と、代謝パスウェイ各々との関連性の強度を算出するとともに、当該関連性の強度が識別可能であるように、二種以上の代謝パスウェイを表す内容を表示させることが可能であり、ユーザは、関連性の強度を識別しつつ、いずれかの代謝パスウェイを選択することができ、選択された代謝パスウェイについて、表示手段に代謝パスウェイマップを表示させるとともに、当該代謝パスウェイマップ上に存在する代謝物及び／又は代謝酵素について、前記測定値データが示す測定値の増減が識別可能であるように表示手段への表示を行う。従って、医師、研究者などの専門家は、測定値データに基づいて関連性の高いと判断される代謝パスウェイを効率的に選択しつつ、代謝物等の測定値の増減状態分析、検討することができ、従来技術と比較すれば、専門家による分析、検討の負担は大きく軽減される。なお、関連性の強度は、測定値データに基づいて、各々の代謝パスウェイについて算出される値であって、算出方法は種々考えられるところである。具体的な算出方法のいくつかの例は、後に詳細に説明する。測定値データと関連しない代謝パスウェイについては算出を割愛し、算出負担の軽減を図ることもできる。【００１１】
なお、本発明の支援装置における支援の対象として、例えば、病態診断支援、検査支援、創薬分野における生体内での薬物代謝の調査支援などが挙げられるが、これらに限定されることはない。【００１２】
測定値データ取得手段は、予めハードディスク、半導体メモリ等の記憶媒体に記憶された測定値データを取得しても良いし、ネットワーク等を介して外部から測定値データを取得しても良い。ユーザがキーボード等の入力装置を用いて入力した測定値データを取得しても良い。【００１３】
パスウェイマップデータ取得手段も同様であり、予めハードディスク、半導体メモリ等の記憶媒体に記憶されたパスウェイマップを表すデータを取得しても良いし、ネットワーク等を介して外部からデータを取得しても良い。「パスウェイマップを表すデータ」は、パスウェイマップのイメージデータやＨＴＭＬ形式のデータなどに限定されず、動的にパスウェイマップを構成することができるようなデータ形式であっても良い。【００１４】
「関連性の強度が識別可能であるように」とは、必ずしも常に関連性強度が識別可能である場合に限定されず、スクロール等を経て関連性強度が表示されるような形態も含む。「代謝パスウェイを表す内容」は、下記実施の形態では、パスウェイ一覧表示領域６０１２に、パスウェイ名称、関連性強度を一覧表示しているが、これに限定されず、アイコンを複数表示するような形態も含まれる。アイコンの形状が各々の代謝パスウェイを表すような実施形態も含まれる。さらに、アイコンをクリックすれば関連性強度が表示されるような実施形態も含まれる。クリックで関連性強度が表示され、ダブルクリックで、当該ダブルクリックで選択された代謝パスウェイマップが表示されるような実施形態も含まれる。もっとも、パスウェイ選択入力受付手段が受け付ける代謝パスウェイの選択入力が、ポインティングデバイスによる選択入力に限定されないことは勿論である。【００１５】
前記関連性強度算出手段は、前記測定値データにより表される測定値の増減を加味して、前記代謝物及び／又は代謝酵素と、代謝パスウェイ各々との関連性の強度を算出することが好ましい。この場合、関連性強度の算出に加味される増減の判定に用いる測定値の範囲指定入力を受け付ける範囲選択入力受付手段を備えることができる。指定された範囲により、測定値の増減は変化し得ることになり、算出される関連性強度も変化し得る。【００１６】
前記関連性強度算出手段は、測定値の増減の起こりにくさ（確率）に基づいて代謝パスウェイごとの関連性強度を算出することができる。例えば、測定値の増減から算出される代謝パスウェイごとの自己エントロピー又はシャノン・エントロピーを関連性の強度として算出することができる。【００１７】
前記関連性強度算出手段は、前記測定値データ取得手段が測定値データを取得した代謝物及び／又は代謝酵素の種類の数に基づいて、前記代謝物及び／又は代謝酵素と、代謝パスウェイ各々との関連性の強度を算出する構成とすることもできる。この場合、例えば、前記関連性の強度を、フィッシャーの正確確率検定により算出する構成とすることができるし、フィッシャーの正確確率検定により算出された関連性の強度を、ＦＤＲを用いて、さらに補正して関連性強度を算出する構成とすることもできる。【００１８】
前記支援装置は、さらに、代謝パスウェイマップ上に表示される代謝物及び／又は代謝酵素のマップ上の位置情報を取得する座標取得手段を備える構成とすることができる。【００１９】
前記表示内容決定手段は、測定値が増加したか減少したかにより、前記パスウェイマップ上において前記代謝物及び／又は代謝酵素を表す部分の色を変化させる構成とすることができる。【００２０】
前記表示内容決定手段は、測定値の増加量、減少量に基づき、前記パスウェイマップ上において前記代謝物及び／又は代謝酵素を表す部分の大きさを変化させる構成とすることもできる。【００２１】
ここで、代謝物及び／又は代謝酵素を表す部分の色、大きさを変化させるための増加量、減少量を見る測定値の範囲を指定する第２の範囲選択入力受付手段を備えることもできる。この第２の範囲選択入力受付手段は、前記した、関連性強度の算出に加味される増減の判定に用いる測定値の範囲指定入力を受け付ける範囲選択入力受付手段、と兼用することもできるが、両範囲を個別に指定することが可能なようにしておけば、より実用に効果的である。【００２２】
前記表示内容決定手段は、前記パスウェイマップ上において前記代謝物及び／又は代謝酵素を表す部分の近傍に、測定値の増減の様子を表す増減パターンを表示させる構成とすることができる。この増減パターンは、表示／非表示を切り替えることが可能な構成とすることができ、一例として、パスウェイマップと異なるレイヤに増減パターンを形成する構成とすることができる。【００２３】
前記増減パターンにおいて、測定値の増減の様子の判定に閾値を用いる構成とすることができる。閾値については、下記に詳細に説明するが、前記閾値の入力を受け付ける閾値入力受付手段を備えておくことが好ましく、例えば、前記閾値として、前記増減パターンを表す線の始点の位置決定に用いる第１の閾値と、始点の位置が決定された後の増減パターンの生成に用いる第２の閾値とが設定される構成や、前記閾値として、測定値との比較に用いる閾値が設定される構成が考えられる。前記閾値として、前記代謝物及び／又は代謝酵素ごとに異なった閾値を用いる構成としても良い。【００２４】
前記増減パターンの生成に用いる測定値の範囲指定入力を受け付けるパターン表示範囲選択入力受付手段を備える構成としても良い。この範囲選択は、前記した、関連性強度の算出に加味される増減の判定に用いる測定値の範囲指定入力を受け付ける範囲選択入力受付手段や、前記第２の範囲選択入力受付手段の双方、またはいずれかと兼用することもできるが、個別に設けても良い。【００２５】
前記代謝パスウェイマップ上に表示された前記代謝物及び／又は代謝酵素の選択入力を受け付ける代謝物選択入力受付手段を備え、前記表示内容決定手段は、選択された代謝物及び／又は代謝酵素について、測定値又は測定値を表すグラフを表示するように表示内容を決定する構成とすることができる。測定値を表すグラフとして、例えば棒グラフを表示させる場合、例えばＪＡＶＡ（登録商標）の機能を利用して棒グラフを生成することが可能であるが、これに限定されない。この機能は、増減パターンのみでは、具体的な増減の様子を正確に把握しにくい場合もあり得るために設けられるものであり、下記に説明するような棒グラフに限定されず、測定値を、より正確に表す折れ線グラフとすることも考えられるし、測定値データの具体的な数値を表示するような実施形態も可能であろう。【００２６】
本発明に係る支援方法は、表示手段が接続されたコンピュータを用いて、複数種類の代謝物及び／又は代謝酵素の各々について、異なる条件で得られた複数の測定値データから、複数種類の代謝パスウェイ各々への影響の分析、検討を支援する支援方法であって、生体内の代謝パスウェイマップを表すデータを取得するパスウェイマップデータ取得ステップと、前記測定値データに基づいて、前記代謝物及び／又は代謝酵素と、代謝パスウェイ各々との関連性の強度を算出する関連性強度算出ステップと、前記表示手段への表示内容を決定する第１の表示内容決定ステップとを含み、前記第１の表示内容決定ステップにおいて、前記表示手段に、前記関連性の強度が識別可能であるように、二種以上の代謝パスウェイを表す内容を表示させることが可能であり、さらに、前記表示手段に前記二種以上の代謝パスウェイを表す内容が表示された場合に、ユーザからの、いずれかの代謝パスウェイの選択入力を受け付けるパスウェイ選択入力受付ステップと、前記パスウェイ選択入力受付ステップにおいて受け付けた選択入力により選択された代謝パスウェイについて、前記表示手段に代謝パスウェイマップを表示させるとともに、当該代謝パスウェイマップ上に存在する代謝物及び／又は代謝酵素について、前記測定値データが示す測定値の増減が識別可能であるように、前記表示手段への表示内容を決定する第２の表示内容決定ステップを含むことを特徴としている。【００２７】
前記関連性強度算出ステップにおいて、前記測定値データにより表される測定値の増減を加味して、前記代謝物及び／又は代謝酵素と、代謝パスウェイ各々との関連性の強度を算出することが好ましい。【００２８】
本発明に係るコンピュータプログラムは、表示手段が接続されたコンピュータ上で動作するコンピュータプログラムであって、前記コンピュータを、複数種類の代謝物及び／又は代謝酵素の各々について、異なる条件で得られた複数の測定値データを取得する測定値データ取得手段、生体内の代謝パスウェイマップを表すデータを取得するパスウェイマップデータ取得手段、前記測定値データ取得手段が取得した測定値データに基づいて、前記代謝物及び／又は代謝酵素と、代謝パスウェイ各々との関連性の強度を算出する関連性強度算出手段、及び、前記表示手段への表示内容を決定する表示内容決定手段として動作させ、前記表示内容決定手段は、前記表示手段に、前記関連性の強度が識別可能であるように、二種以上の代謝パスウェイを表す内容を表示させることが可能であり、さらに、前記コンピュータプログラムは、前記コンピュータを、前記表示手段に前記二種以上の代謝パスウェイを表す内容が表示された場合に、ユーザからの、いずれかの代謝パスウェイの選択入力を受け付けるパスウェイ選択入力受付手段として動作させ、前記表示内容決定手段は、前記パスウェイ選択入力受付手段が受け付けた選択入力により選択された代謝パスウェイについて、前記表示手段に代謝パスウェイマップを表示させるとともに、当該代謝パスウェイマップ上に存在する代謝物及び／又は代謝酵素について、前記測定値データが示す測定値の増減が識別可能であるように、前記表示手段への表示内容を決定することを特徴としている。【００２９】
前記関連性強度算出手段は、前記測定値データにより表される測定値の増減を加味して、前記代謝物及び／又は代謝酵素と、代謝パスウェイ各々との関連性の強度を算出する構成とすることが好ましい。
【発明の効果】
【００３０】
本発明に係る支援装置等によると、代謝物、代謝酵素の測定値データについて、専門家による分析検討を支援することが可能になるという効果を奏する。
【図面の簡単な説明】
【００３１】
【図１】本発明の支援装置の一実施の形態における装置全体構成の概略について説明するための図である。
【図２】コンピュータ１０１等のハードウェア構成の一例について説明するためのブロック図である。
【図３】本発明の実施の形態における支援装置の機能的構成の一例について説明するための機能ブロック図である。
【図４】代謝パスウェイの一例としてのクエン酸回路について、代謝パスウェイマップの一例を示す図である。
【図５】代謝経路ＤＢ格納部３０３に格納される代謝パスウェイ情報の内容の一例を示す図である。
【図６】測定値データ格納部３０４に格納される測定値データの内容の一例を示す図である。
【図７】選択情報量（自己エントロピー）を利用した関連性強度算出について説明するための図である。
【図８】測定値データの変動確率算出の一例について説明するための図である。
【図９】シャノン・エントロピーを利用した関連性強度算出について説明するための図である。
【図１０】関連性強度を代謝物の数から算出する場合の代謝パスウェイ情報の一例を示す図である。
【図１１】ｐ値の算出について説明するための図である。
【図１２】座標取得部３０１２の処理について説明するための図である
【図１３】測定値の増減を示す円、増減パターンをパスウェイマップに重ねて表示した様子の一例を示す図である。
【図１４】前処理の内容について説明するためのフローチャートである。
【図１５】座標情報の内容の一例を示す図である。
【図１６】代謝装置起動時の処理内容について説明するためのフローチャートである。
【図１７】表示内容生成処理の内容の一例について説明するためのフローチャートである。
【図１８】表示装置１０２に表示される表示画面の一例を示す図である。
【図１９】パスウェイ一覧表示領域６０１２に表示され得る画面の全体像の一例を示す図である。
【図２０】測定値データの棒グラフ６０１１１を表示した様子の一例を示す図である。
【図２１】増減判定範囲指定の入力を受け付けるために表示画面に表示されるウィンドウの一例を示す図である。
【図２２】関連性強度の算出を行う範囲の開始を「条件１」、終了を「条件４」と設定し、増減パターンを表示する範囲の開始を「条件１」、終了を「条件３」、円の色を決定する増減判定の範囲の開始を「条件２」、終了を「条件３」と設定した場合の表示画面の一例を示す図である。
【図２３】増減表示・パターン生成処理の内容の一例を示すフローチャートである。
【図２４】開始時閾値、開始以後閾値を設定するために表示されるウィンドウの例を示す図である。
【図２５】閾値を用いた増減パターン生成について説明するための図である。
【図２６】測定値の閾値を設定した場合の増減パターンの例を示す図である。
【発明を実施するための形態】
【００３２】
以下、本発明の実施の形態について、図面を参照しながら説明する。
図１は、本発明の支援装置の一実施の形態における装置全体構成の概略について説明するため図である。
【００３３】
本実施の形態の支援装置は、表示装置１０２が接続されたコンピュータ１０１上でコンピュータプログラムが動作することにより実現される。コンピュータ１０１には、入力装置の例としてキーボード１０３やマウス１０４が接続されており、例えばＬＡＮケーブル１０５を介して、ＬＡＮ、インターネット等のネットワーク１０６と通信を行うことができる。
【００３４】
図２は、コンピュータ１０１等のハードウェア構成の一例について説明するためのブロック図である。コンピュータ１０１は、ＣＰＵ２０１を中心として構成され、ＣＰＵ２０１には、記憶素子であるＲＯＭ２０２やＲＡＭ２０３が接続されている。また、他のデータ記憶装置としてＨＤ（ハードディスク）ドライブ２０４や、ＣＤ／ＤＶＤドライブ２０５が接続されており、ハードディスクやＣＤ（コンパクト・ディスク）、ＤＶＤなどのディスク媒体からデータの読み出しを行うことができる。
【００３５】
さらにＣＰＵ２０１には、ＮＩＣ（ネットワーク・インタフェース・カード）２０６が接続され、コネクタに前記ＬＡＮケーブル１０５を挿入してネットワーク１０６と通信を行うことができる。ネットワーク１０６には、構内のＬＡＮ、インターネット等の広域ネットワークなどが含まれる。本実施の形態の支援装置を実現するコンピュータプログラムは、例えばＤＶＤに記憶されたものをハードディスク上にインストールして使用することができる。支援装置は、ハードディスクドライブ２０４から、必要なプログラムモジュールを適宜読み出してＲＡＭ２０３へと格納し、さらにＣＰＵ２０１が読み出して動作することにより実現される。
【００３６】
図３は、本実施の形態における支援装置の機能的構成の一例について説明するための機能ブロック図である。本実施の形態の支援装置は、支援装置制御部３０１、代謝経路マップ格納部３０２、代謝経路ＤＢ（データベース）格納部３０３、測定値データ格納部３０４を含んでいる。入力装置３０５は、例えば図１に示されるキーボード１０３やマウス１０４を含むが、これらに限定されず、他のポインティングデバイスやソフトウェアキーボード、音声入力装置などを備えて利用することもできる。本実施の形態では、代謝経路マップ格納部３０２、代謝経路ＤＢ（データベース）格納部３０３、測定値データ格納部３０４としては、ＨＤドライブ２０４が用いられる。
【００３７】
支援装置制御部３０１は、代謝経路マップ管理部３０１１、座標取得部３０１２、代謝経路データ管理部３０１３、測定値データ管理部３０１４、関連性強度算出部３０１５、選択入力受付部３０１６、表示内容決定部３０１７、表示制御部３０１８を含んでいる。本実施の形態では、これらの各部は、ＣＰＵ２０１等においてコンピュータプログラムが動作して、その機能が実現される。なお、支援装置制御部３０１は、図２では詳細な図示が省略された入出力インタフェースやビデオＲＡＭ等、通常のコンピュータに備えられえるハードウェア資源を含んでおり、それら各部とＣＰＵ２０１とが協働して、上記各部の機能が実現される。
【００３８】
代謝経路マップ管理部３０１１は、代謝経路マップ格納部３０２からのデータ読み出し、データ内容修正等を行う部分である。代謝経路データ管理部３０１３は、代謝経路ＤＢ格納部３０３からのデータ読み出し、データ内容追加修正等を行う。測定値データ管理部３０１４は、測定値データ格納部３０４からのデータ読み出し、データ入力、データ内容追加修正等を行う。
【００３９】
座標取得部３０１２、関連性強度算出部３０１５の処理内容については後に詳細に説明する。選択入力受付部３０１６は、マウス１０４、キーボード１０３等を含む入力装置３０５からの選択入力を受け付ける。選択の対象は、表示装置１０２に表示されたパスウェイマップ上の代謝物、代謝酵素であったり、パスウェイ一覧として表示された代謝パスウェイのいずれかであったりする。表示内容決定部３０１７は、表示装置１０２への表示内容を決定する。表示制御部３０１８は、実際に表示装置１０２の表示画面への表示制御を行う部分である。
【００４０】
代謝経路マップ格納部３０２には、代謝経路マップの情報が格納される。代謝経路マップとは、代謝経路（代謝パスウェイ）を図示したものである。図４に代謝パスウェイの一例としてのクエン酸回路について、代謝パスウェイマップの一例を示す。同図に示されるように、代謝経路マップとは、代謝物、代謝酵素、及びそれらの生体反応経路を図に表したものである。同図において丸印は代謝物を示し、その近傍に代謝物名称が表記されている。四角は代謝酵素を表し、その中に、代謝酵素の機能等を表す識別子が記入されている。同図４に示されるような代謝経路マップ上において、代謝物等の測定値の増減をビジュアルに表示することは、医師、研究者などの専門家による分析、検討に大きな負担軽減となる。
【００４１】
なお、代謝パスウェイマップは、例えば上記したＫＥＧＧからダウンロードしたデータを代謝経路マップ格納部３０２としてのＨＤドライブ２０４に格納して利用することができる。もっとも、ネットワーク１０６を介して必要に応じて外部から取得するようにしても良い。なお、代謝経路マップを表すデータは、必ずしも図４に示したようなイメージ形式でなくても良い。例えば、代謝パスウェイの名称に加えて代謝物、代謝酵素、代謝経路などの情報をデータベースとして格納しておき、必要に応じて動的に、図４に例示したようなマップイメージを作成するような実施形態も可能である。
【００４２】
図４の中で、点線で表示された部分４０１は、連続する（関連する）他の代謝パスウェイを示す。マウス１０４などのポインティングデバイスで当該部分４０１を選択してクリック等を行うことで、当該他の代謝パスウェイマップの表示に切り替えることができる。
【００４３】
代謝経路ＤＢ格納部３０３には、代謝パスウェイについての情報が格納される。図５は、代謝経路ＤＢ格納部３０３に格納される代謝パスウェイ情報の内容の一例を示す図である。本実施の形態では、代謝パスウェイの情報として代謝パスウェイ名称、測定値データに含まれる代謝物の数、代謝パスウェイに含まれる代謝物の数、及び関連性強度を示す値を格納している。
【００４４】
測定値データに含まれる代謝物の数は、測定値データ格納部３０４の内容からカウントできる。図６は、測定値データ格納部３０４に格納される測定値データの内容の一例を示す図である。同図に示されるように、測定値データとしては、代謝物名称、代謝物に付与された識別子（ＩＤ）、及び異なる条件における各代謝物の測定値が保持される。同図に示されるＩＤは上記したＫＥＧＧにおいて割り振られた化合物ＩＤであるが、ＩＤの付与の仕方はこれに限定されない。なお、異なる条件としては、例えば、ある化合物を投与してからの経過時間が異なるものや、投与濃度が異なるものなどが考えられるが、それらに限定されない。また、測定値データ格納部３０４に格納される測定値データは、必ずしもユーザによる実測値でなくてもよく、外部から取得した参考用データや文献等で公表されているようなデータを入力してもよい。適当な値を入力し、本実施の形態の支援装置をシミュレータ的に利用することもできる。
【００４５】
図５に示した代謝パスウェイ情報に含まれる「測定値データに含まれる代謝物の数」は、図６の測定値データにデータが格納されている代謝物の中で、当該代謝パスウェイに含まれる化合物の数である。「代謝パスウェイに含まれる代謝物の数」は、例えば図４に示したクエン酸回路であれば丸印の数、即ち１２となる。ここでは代謝物のみを対象としたが、代謝酵素の数としてもよいし、代謝物の数と代謝酵素の数の合計としてもよい。なお、図５で示したように、本実施の形態では、「関連性強度」として、測定値データから算出される各代謝パスウェイとの関連性の強度を示す数値を保持している。関連性強度は、図３に示した関連性強度算出部３０１５において算出される。
【００４６】
関連性強度の算出方法は種々考えられるが、代謝物及び／又は代謝酵素の測定値データから得られる測定値の増減（変動量）を加味して関連性強度を算出することが好ましい。測定値の変動量が大きいほど、当該代謝物及び／又は代謝酵素を含む代謝パスウェイと関連性が高いと判断できるからである（例外は存在し得るが、概ね成立する）。図７は、関連性強度の算出に自己エントロピーを用いた例について説明するための図である。同図の例では、代謝物数１９の代謝パスウェイの例を示している。
【００４７】
自己エントロピーＨ（同図の例では単位はｂｉｔｓ）は、ｉ番目（ｉは１から１９）の代謝物の各々についてｐｉを算出し、−ｌｏｇ２（ｐｉ）の総和として算出することができる。測定値データがない場合、ｉ番目の代謝物について、変動量は０となるが、この場合には、−ｌｏｇ２（ｐｉ）＝０と定義する。図８は、ｐｉの算出方法の一例について説明するための図である。
【００４８】
同図８に示される代謝物Ｘの測定値データが条件１から条件４まで、図示された如くであるとする。「変動量」は、条件１の測定値と条件２の測定値との差分、条件１と条件３、条件１と条件４とのそれぞれの変動量である。「変動の大きさ」は、変動量の絶対値である。同図の例では、変動の大きさの合計は５４．４となる。
【００４９】
次に、変動の大きさの単純割合を算出する。同図に示されるように、単純割合は、「変動の大きさ」を「変動の大きさの合計」で割った数値である。そして、その逆数を求め、当該「単純割合の逆数の合計（同図の例では３２．７８）」で「単純割合の逆数」を割った値を代謝物Ｘについてのｐ（Ｘ）の値とする。条件１と条件２、条件１と条件３、条件１と条件４とについて、それぞれ算出されたｐ（Ｘ）の値の計は１．００となる。
【００５０】
同図８の例からもわかるように、ｐ（Ｘ）の値は代謝物測定値の変動の起こりにくさ（確率）を示す指標となる（大きな変動ほど起こりにくいため、ｐ（Ｘ）の値は小さくなる）。このようにして得られたｐ（Ｘ）の値について、ｉ番目の代謝物についての値をｐｉとし、図７に示されるように−ｌｏｇ２（ｐｉ）の総和として算出することにより代謝パスウェイごとの自己エントロピーが求められる。この自己エントロピーを関連性強度とする。
【００５１】
もっとも、測定値の増減を加味する関連性強度の算出方法の他の一例として、条件１と条件２、条件２と条件３、条件３と条件４との差分（増減量の絶対値）を順次加算してトータルの変動量を算出し、この結果から自己エントロピー（もしくは後述するシャノン・エントロピー）を算出するような実施形態も可能である。
【００５２】
関連性強度の算出方法の他の一例としてシャノン・エントロピーを利用する方法がある。図９は、シャノン・エントロピーを利用した関連性強度算出について説明するための図である。同図に示されるように、代謝物測定値データの変動量（同図の例では変動量の絶対値）から、代謝パスウェイごとのシャノン・エントロピーを算出することで、関連性強度を求めることができる。同図の例は、３１種の代謝物を含む代謝パスウェイ（例えばbeta-Alanine metabolism）の場合の例であるが、他の代謝パスウェイでも同様の考え方で関連性強度を算出することができる。シャノン・エントロピーＨ（同図の例では単位はｂｉｔｓ）の値が大きいほど関連性は高いと判定できる。測定値データに含まれない代謝物については、変動量は０となるが、自己エントロピーと同様に、ｌｏｇ２（ｐｉ）＝０と定義する。なお、測定値の増減を加味した関連性強度の算出方法が、自己エントロピーやシャノン・エントロピーに限定されるわけではなく、他の関連性強度算出方法も存在し得る。
【００５３】
測定値の増減を加味した自己エントロピーやシャノン・エントロピー等を関連性強度とした場合、１つのパスウェイについて、複数の条件での測定値データがあれば、複数の関連性強度を算出することができる。図５に示される「関連性強度１」は、条件１と条件２との間の増減から算出した関連性強度、「関連性強度２」は、条件１と条件３との間の増減から算出した関連性強度、「関連性強度３」は、条件１と条件４との間の増減から算出した関連性強度の例である。これらの関連性強度は、測定値データから予め算出しておいても良いし、増減判定範囲の指定を受けたときに算出するようにしても良い。
【００５４】
本実施の形態のように、関連性強度という指標が無ければ、１つずつ代謝パスウェイを切換え表示させ、後述する円等の図形の色、大きさ、測定値の増減パターンなどを手掛かりにして変動の大きそうなものを目視で探索していかなければならない。関連性強度の大きなパスウェイを優先して探索することにより、変動の原因に近い箇所を容易に見つけることが可能となる。また、どの実験条件（あるいは実験条件の差）がパスウェイに大きな影響を与えたのかが、パスウェイ毎、あるいはパスウェイ全体について判断しやすくなるため、専門家によって、薬効や毒性がパスウェイに与える影響を比較することが容易となる。
【００５５】
なお、関連性強度としては、測定値の変動を加味する算出方法に限らず、代謝物の数から算出する方法も可能である。図１０は、関連性強度を代謝物の数から算出する場合の代謝パスウェイ情報の一例を示す図である。同図に示される「ｐ値」は、例えばフィッシャーの正確確率検定と称される方法を用いて、代謝パスウェイごとに算出される値である。
【００５６】
図１１は、ｐ値の算出について説明するための図である。例えば全ての代謝パスウェイ上にある代謝物の種類の合計をｎ、ある代謝パスウェイＡ上にある代謝物の種類の合計をｘとする。現在発見されている３００種以上の代謝パスウェイにおいて、代謝物の種類ｎの値は約６千近い値となる。もっとも、新たな代謝パスウェイ、新たな代謝物が増加しない限りｎの値は一定である。ｘの値は、予め代謝パスウェイごとに求めておくことができる。ｘの値は、図１０に、「パスウェイに含まれている代謝物数」として格納されている。
【００５７】
「パスウェイに含まれている代謝物数」ｘは、測定値データに含まれる代謝物の種類の数ａと、測定値データに含まれない代謝物の種類の数ｃとの合計となる。図１０に示されるようにａの値は「測定値データに含まれる代謝物数」として、代謝パスウェイの情報に含まれている。ｃの値は、ｘからａを減算して求めることができる。
【００５８】
図１１において、「測定値データに含まれる代謝物の種類」の合計（ｙ＝ａ＋ｂ）は、図６に示した測定値データから求めることができる。測定値データが取得されている代謝物の種類の数がｙである。従って、ｂの値は、ｙからａを減算して求めることができる。
【００５９】
前記したように全代謝物の種類ｎの値は一定であるから、「測定値データに含まれない代謝物の種類」の合計（ｃ＋ｄ）は、ｎからｙを減算して求めることができる。ｃの値は上記に算出されているので、ｄの値も、上記（ｃ＋ｄ）からｃの値を減算して求めることができる。
【００６０】
上記のように、ａ、ｂ、ｃ、ｄの値が算出されると、フィッシャーの正確確率検定を利用してｐ値が算出される。下記（数１）は、代謝パスウェイＡに対するｐ値Ｐａを求める式である。
【００６１】
（数１）
Pa
= ａ＋ｂCａ × ｃ＋ｄCｃ
／ ｎCａ＋ｃ 
= ((a＋b)！(c＋d)！(a＋c)！(b＋d)！) ／ (n！a！b！c！d！)  
ただし、n=a＋b＋c＋d  
【００６２】
ｐ値は代謝パスウェイごとに求められ、その数値が小さいほど、測定値データに含まれる代謝物との関連性が強いことを示す値となる。フィッシャーの正確確率検定を用いて算出されたｐ値は、その値が小さいほど、偶然には生じにくいことを示すものだからである。この手法は、単に、測定値データに含まれる代謝物が代謝パスウェイ上に存在する割合を示す数値ではなく（そのような方法であれば、ａ＝５、ｂ＝１０の場合と、ａ＝１０、ｂ＝２０の場合とでは関連性強度が同一と算出されることになる。）、関連性強度を、より好適に算出する方法では有り得るが、代謝物の種類の数から代謝パスウェイへの関連性強度を算出する方法は、上記（数１）の手法に限定されるわけではない。
【００６３】
なお、ｐ値の補正は、いわゆるＦＤＲ（Ｆａｌｓｅ  Ｄｉｓｃｏｖｅｒｙ  Ｒａｔｅ）を利用して行うことができる。「補正後ｐ値」は、例えばＦＤＲを０．０５に設定して補正して得ることができる。例えば、この補正後ｐ値が０．０５以上の代謝パスウェイを、測定値データと関連性の高い代謝パスウェイであると判定することができる。
【００６４】
次に、代謝パスウェイマップ上における測定値の増減の表示について説明する。前記したように、本実施の形態の支援装置では、測定値データから判別される代謝物ごとの測定値の増減を、代謝パスウェイマップ上にビジュアルに表示することを特徴としている。
【００６５】
このようなビジュアルな表示を可能とする処理を行うのが、座標取得部３０１２である。代謝経路マップ格納部３０２に格納されている代謝パスウェイマップには、測定値の増減を示す表示は付加されていないため、代謝パスウェイマップ上に測定値の増減をビジュアルに表示するためには、代謝パスウェイマップ上の代謝物の位置の座標を算出する必要がある。図１２は、座標取得部３０１２の処理について説明するための図である。
【００６６】
例えば代謝パスウェイマップの左上隅を基準点（０，０）とすると、この代謝パスウェイ（ＴＣＡサイクル）のパスウェイマップ上のコハク酸（Succinate）の表示位置の座標は（ｘ，ｙ）と検出できる。図６に示した例の場合、コハク酸の測定値は、条件１と条件４とで比較すると０．４２減少している。この場合、図１３に示すように、コハク酸の表示位置の座標（ｘ，ｙ）を中心として、測定値の減少を示す、例えば青色の円５０１を代謝パスウェイマップに重ねて表示する。
【００６７】
本実施の形態では、さらに、円５０１の近傍に測定値の増減を示す折れ線グラフ形状の図形（以下、単に「増減パターン」という。）５０２を表示する。増減パターン５０２は、パスウェイマップとは別のレイヤに作成して、表示手段１０２への表示、非表示を切り替えるようにすることができる。図６の例では、さらにスクシニル補酵素Ａについても、測定値データが得られており、条件１と条件４とで比較すると２．６０増加しているので、図１３に示すように、スクシニル補酵素Ａの表示位置の座標を中心として、測定値の増加を示す、例えば赤色の円５０３を代謝パスウェイマップに重ねて表示する。さらに、円５０３の近傍に測定値の増減パターン５０４を表示する。
【００６８】
以下、本実施の形態の支援装置について、その処理内容について説明する。もっとも、本実施の形態においては、支援装置を動作させるために、若干の前処理を行う。まず、前処理の内容について説明する。図１４は、前処理の内容について説明するためのフローチャートである。
【００６９】
前処理では、まず、全ての代謝パスウェイマップの情報を取得する（Ｓ１０１）。この情報は、例えば前記したＫＥＧＧから取得することができるが、ＫＥＧＧに限定されない。また、この情報はイメージデータである必要はなく、動的にマップイメージを構成できるような情報であっても良いが、本実施の形態では、画面表示を迅速にする等のために、予めイメージデータを構成、格納しておくようにしている。取得した情報は、例えば代謝パスウェイごとのＨＴＭＬ形式ファイルなどとして、代謝経路マップ格納部３０２に格納しておく。
【００７０】
次に、各代謝パスウェイに含まれる代謝物数を取得する（Ｓ１０２）。この代謝物数は、上記代謝パスウェイマップの情報から計数することも可能ではあるが、別に用意されている代謝パスウェイの管理テーブルから取得する方が、より迅速処理が可能となる。管理テーブルとして、ＫＥＧＧの場合には、ファイル名「ｃｐｄ＿ｍａｐ．ｔａｂ」が用意されており、当該ファイルから取得することができる。ＫＥＧＧ以外の場合でも、同様の管理テーブルがあれば、当該管理テーブルから取得することもできる。取得した値は、上記したように、図５（若しくは図１０）に例示した代謝パスウェイ情報（パスウェイに含まれる代謝物数）として、代謝経路ＤＢ格納部３０３に格納される。
【００７１】
次に、全ての代謝パスウェイに含まれる代謝物数ｎを計数する（Ｓ１０３）。このｎの値も、上記した管理テーブルに格納されていれば、迅速に数値を取得することができる。もっとも、代謝経路マップ格納部３０２に格納されている情報から計数することも可能である。取得、若しくは計数した値は、例えば代謝経路ＤＢ格納部３０３に格納しておくことができる。
【００７２】
次に、各代謝物について、代謝パスウェイマップ上における位置情報を取得する（Ｓ１０４）。具体的な処理については、上記した座標取得部３０１２の処理内容として説明した。図１２に例示したような、各代謝物のマップ上の位置を示す座標情報を、全ての代謝パスウェイの全代謝物について、前処理で予め算出しておく。この位置（座標）情報は、例えば代謝経路マップ管理部３０１１を介して、代謝経路マップ格納部３０２に格納しておくことができる。図１５は、座標情報の内容の一例を示す図である。図１５の例では代謝物を例示したが、代謝酵素の座標情報、若しくは代謝物と代謝酵素の両方の座標情報を保持することもできる。
【００７３】
以上のように前処理を終了した後、本実施の形態の代謝装置が起動された場合の具体的な処理内容、表示装置１０２への表示内容等について、以下に説明する。図１６は、代謝装置起動時の処理内容について説明するためのフローチャートである。
【００７４】
支援装置が起動されると、まず、測定値データ格納部３０４から測定値データを読み出す（Ｓ２０１）。起動前には、測定値データ格納部３０４に測定値データが格納されている必要がある。測定値データは、例えば入力装置３０５を用いてユーザが入力することができるが、ネットワーク１０６を介して外部から取得することもできる。例えばＬＣ−ＭＳ（液体クロマトグラフィ−質量分析装置）やＧＣ−ＭＳ（ガスクロマトグラフィ−質量分析装置）にて取得された分析データを取得することもできる。
【００７５】
ステップＳ２０２では、増減判定の範囲変更指示がなされたか否かを判定する。装置起動時は、このステップの判定結果はＹＥＳとなるようにしておき、代謝パスウェイごとの関連性強度の算出へと進む（Ｓ２０３）。本実施の形態では、上記した自己エントロピーやシャノン・エントロピーを用いた方法のように、代謝物の数だけでなく、測定値の増減を加味して関連性強度を算出する。関連性強度の算出のための増減として、初期的には、測定値データの最初と最後（図６の例では「条件１」と「条件４」）との間で増減を算出する。後述するように、この範囲は変更することができ、一度ステップＳ２０３にて関連性強度が算出された場合には、範囲変更の入力がなされた場合に改めて関連性強度を算出する。算出された関連性強度は、図５（若しくは図１０）に例示した代謝パスウェイ情報として、代謝経路ＤＢ格納部３０３に格納される。なお、本実施の形態では、初期設定を、測定値データの最初と最後としたが、予め他に初期的な範囲を指定しておくようにすることもできる。また、範囲変更の入力がなされた場合に関連性強度を算出するのではなく、予め全てのケースについて関連性強度を算出しておくような実施形態も可能である。
【００７６】
次に、表示装置１０２への表示内容生成処理を行う（Ｓ２０４）。図１７は、表示内容生成処理の内容の一例について説明するためのフローチャートである。図１８は、表示装置１０２に表示される表示画面の一例を示す図である。以下、表示画面の表示内容について説明する。
【００７７】
表示装置１０２の表示画面６０１には、パスウェイマップ表示領域６０１１、パスウェイ一覧表示領域６０１２、測定値データ表示領域６０１３が設けられる。これらの各領域は、図１８に示したように並べて表示するのではなく、重ねて表示し、最前面の領域を切り替えるようにすることもできる。各々の領域は、いずれもスクロールバーを操作して、画面に表示される範囲を移動することができる。
【００７８】
パスウェイマップ表示領域６０１１には、図１３に例示したようなパスウェイマップが表示される。測定値データの増減を示す円５０１、５０３、増減パターンを示す折れ線グラフ５０２、５０４も併せて表示される。なお、円５０１、５０３の色、大きさを決定するための増減は、測定値データにおける条件範囲、例えば条件１と条件２との間で比較するか、条件１と条件４との間で比較するかで異なることもある。また、増減パターンは、どれだけの数値変化があった場合に増加、減少と判断するかにより変化し得る。そのような条件設定、数値変化の閾値設定については後述する。
【００７９】
パスウェイ一覧表示領域６０１２には、パスウェイの一覧が表示される。図１９は、この領域に表示される画面の全体像の一例を示す図であるが、実際には、測定値データの内容に応じて多数の代謝パスウェイが一覧表示される。パスウェイ一覧表示領域６０１２には、測定値データと関連のある代謝パスウェイのみを一覧表示しても良いし、全代謝パスウェイを一覧表示するようにしても良い。もっともいずれの場合も、最初は、関連性強度により示される関連性の強い順に上から一覧表示する。なお、図１９の例では、関連性強度１（条件１と条件２との間の増減から算出した関連性強度）、関連性強度２（条件１と条件３との間の増減から算出した関連性強度）、関連性強度３（条件１と条件４との間の増減から算出した関連性強度）を並べて表示しているが、設定されている条件範囲で算出された関連性強度のみを表示するようにしてもよい。
【００８０】
初期画面では、パスウェイ一覧表示領域６０１２には、関連性強度が大きい順に上から代謝パスウェイ名称が一覧表示されるとともに、パスウェイマップ表示領域６０１１には、測定値データから算出された関連性強度がもっとも大きいＴＣＡサイクルのパスウェイマップが表示される。表示画面６０１上には、マウスカーソル１０４１が表示されていて、パスウェイ一覧表示領域６０１２上で任意の代謝パスウェイを選択入力することができる。マウス１０４（キーボード１０３でも良い）等の入力装置３０５を用いたパスウェイ選択入力は選択入力受付部３０１６（図３参照）にて受け付けられて、表示内容決定部３０１７に送られ、表示画面６０１の表示内容に反映される。具体的には、選択された代謝パスウェイのパスウェイマップがパスウェイマップ表示領域６０１１に表示される。本実施の形態では、この際にも、測定値データ等に基づいて、円や増減パターンが併せて表示される。
【００８１】
パスウェイ一覧表示領域６０１２の下部に設けられたスクロールバー６０１２１によりパスウェイ一覧表示領域６０１２に表示される部分を左右に移動させると、図１９に例示された関連性強度の数値が表示される。マウスカーソル１０４１で「関連性強度１」等と表示された領域をクリックすれば、パスウェイ一覧を表示する際の、関連性強度の降順、昇順の切り替えも可能である。もっとも、昇順に切り替えた場合に、全く関連性の無い代謝パスウェイが最上段に表示されると操作が煩雑となる可能性もあるから、上記したように、３００種以上の代謝パスウェイを全て一覧表示するよりは、多少でも測定値データと関連のある代謝パスウェイのみを一覧表示する方が好ましい。そのようにしても、図４に例示したように関連するパスウェイはマップ上に表示され、当該パスウェイに表示を変更することの可能であるから、使い勝手が悪くなることは少ないと考えられる。パスウェイ一覧表示領域６０１２の上下方向には、スクロールバー６０１２２を操作することにより、パスウェイ一覧表示領域６０１２に表示される部分を移動させることができる。
【００８２】
測定値データ表示領域６０１３には測定値データが表示される。スクロールバーの操作により、左右方向、上下方向に移動させることが可能であるのは、上記したパスウェイ一覧表示領域６０１２と同様である。なお、パスウェイマップ表示領域６０１１に表示されるパスウェイマップの大きさは、必ずしも、代謝経路マップ格納部３０２に格納されたサイズと等しくする必要はなく、拡大、縮小も可能である。その際には、円５０１、５０３や増減パターン５０２、５０４の位置も修正する必要があるが、パスウェイマップ上の座標情報が取得されているので、マップの表示サイズに合わせて、円５０１等の表示も容易に修正することができる。パスウェイマップ表示領域６０１１にパスウェイマップの全体を表示するようにすることも可能であり、その場合には、パスウェイマップ表示領域６０１１にはスクロールバーが不要となる。
【００８３】
なお、本実施の形態では、例えば円５０１、５０３の位置にマウスカーソル１０４１を移動させてクリックすると、図２０に例示されるように、測定値データの棒グラフ６０１１１が表示される。後述するように、本実施の形態の増減パターン５０２は、所定の閾値によりパターン化されているため、増減状態の、より詳細な把握には、棒グラフ６０１１１の表示は好適である。もっとも、棒グラフに限定されず、測定値を、より正確に表す折れ線グラフとすることも考えられるし、測定値データの具体的な数値を表示するような実施形態も可能であろう。選択された代謝物の測定値データが、測定値データ表示領域６０１３の先頭に表示されるようにすることもできる。
【００８４】
前記したように、関連性強度を算出するための増減を判定する範囲は変更することが可能である。図２１は、増減判定範囲指定の入力を受け付けるために表示画面に表示されるウィンドウの一例を示す図である。同図に例示されたようなプルダウンメニューにおいて増減を判定する範囲を選択することができる。範囲指定の方法がプルダウンメニューに限定されないことは勿論であり、開始が「条件４」で、終了が「条件１」といった逆方向の指定を可能とすることもできる。なお、本実施の形態では、同様の方法で、「増減パターンを表示する範囲」、「円の色等を決定するための増減判定範囲」を個別に指定することができる。もっとも、関連性強度を算出するための範囲を、そのまま増減パターンや円の色等に反映させるような実施形態も可能であるし、いずれかの範囲を兼用するような実施形態も可能である。
【００８５】
図２２は、関連性強度の算出を行う範囲の開始を「条件１」、終了を「条件４」と設定し、増減パターンを表示する範囲の開始を「条件１」、終了を「条件３」、円の色を決定する増減判定の範囲の開始を「条件２」、終了を「条件３」と設定した場合の表示画面の一例を示す図である。同図に示すように、増減パターン５０２、５０４の形態が図１８とは変化するとともに、円５０１が、測定値の増加を表す赤色に変わり、円５０３が、測定値の減少を表す青色に変化する。
【００８６】
このように、範囲を個別に設定することができるようにする利点について、以下、簡単に説明する。例えば、生体に何らかの刺激を与える（例えば投薬）ことを考える。ここで、条件１から条件４が以下のようであったとする。
条件１：コントロール（刺激無し）
条件２：刺激後８時間
条件３：刺激後２４時間
条件４：刺激後７２時間
【００８７】
関連性強度により、与えた刺激に対して変化が大きい（大きな影響を受けた）代謝パスウェイを知ることができる。このパスウェイ上もしくは近傍に作用点があると考えられ、作用点を探す手がかりとなる。
【００８８】
一方、増減パターン（折れ線）により、与えた刺激に対して亢進されたのか、それとも抑制されたのかを知ることができる。１つのパスウェイ上で亢進・抑制は共に起こり得るが、折れ線のパターンで、パスウェイ上の分子の挙動を容易に区別できる。また、どの分子が同じような制御を受けているかが判別できる。
【００８９】
さらに、円の色・大きさにより、特に注目する条件間での変化を理解することが容易となる。例えば条件１−２間では刺激直後の変化を、条件２−３、２−４、３−４間では、刺激後に持続する変化の大きさと方向に着目して眺めることができる。早い時間に変化して元に戻る分子は、後から変化して元に戻る、あるいは変化し続ける分子よりも作用点に近いと推測できるし、逆に後者から結果となる現象に関係している分子を知ることができる。原因あるいは結果のいずれに近い分子かを知ることで、薬効や毒性のメカニズムを考察することが容易となる。
【００９０】
投与量の違いを調べるような事例でも同様に考えることができる。例えば、以下のような投与量の異なる条件で実験を行ったとする。
条件１：コントロール（投与量０）
条件２：投与量１単位
条件３：投与量３単位
条件４：投与量１０単位
【００９１】
関連性強度、および増減パターン（折れ線）については、上記の事例と同様であるが、条件１−２間に着目すると、投与によって初めて生じる変化を知ることができる。条件２−３間、あるいは条件３−４間では、投与量の違いによる差異を検討することができる。以上に説明したように、関連性強度を算出する条件の範囲と、増減パターンの範囲、円の色（ないし大きさ）を決定する増減の範囲を個別に設定することができるようにしておくことは、専門家による分析、検討に大きく寄与すると言える。
【００９２】
以上に説明したような表示内容を生成する処理について、以下、図１７に戻って説明を続ける。表示内容生成処理では、まず、測定値データ表示領域６０１３の表示内容を生成する（Ｓ３０１）。この処理は、測定値データ格納部３０４に格納されている測定値データを読み出して、配置する処理である。
【００９３】
次に、増減判定範囲が変更されたか否かを判定する（Ｓ３０２）。装置の起動時には、図１６のステップＳ２０２と同様、判定結果がＹＥＳとなるようにしておく。増減判定範囲が変更された場合（Ｓ３０２：ＹＥＳ）、パスウェイ一覧表示領域６０１２の表示内容を生成する（Ｓ３０３）。上記したように、増減判定範囲が変更された場合、各々の代謝パスウェイの関連性強度が変化するためである。
【００９４】
次に選択されたパスウェイマップを、代謝経路マップ格納部３０２から読み出す（Ｓ３０４）。装置の起動時には、関連性強度のもっとも強い代謝パスウェイが選択されるように設定する。装置の起動後は、前記したようにマウスカーソル１０４１の移動、選択により、代謝パスウェイの選択入力が可能であり、装置の起動後には、このステップでは、選択されたパスウェイのマップを読み出す処理を行う。
【００９５】
次に座標情報を取得する（Ｓ３０５）。座標情報は前処理で取得され、格納済みであるので、ここでは、格納された座標情報を読み出す。そして、増減表示・パターン生成処理（Ｓ３０６）へと進む。図２３は、増減表示・パターン生成処理の内容の一例を示すフローチャートである。
【００９６】
増減表示・パターン生成処理は、読み出した座標情報で表される位置に円５０１等を配置したり、増減パターン５０２等を生成する処理である。増減表示・パターン生成処理では、まず、変数ｉ、ｊにそれぞれ１を設定する（Ｓ４０１）。変数ｉは、測定値データに含まれる代謝物の種類のカウンタである。変数ｊは、各々の代謝物について、測定値データに含まれる条件の数のカウンタである。
【００９７】
次のステップＳ４０２では、ｉ番目の代謝物について、ｊ番目の条件の測定値を読み出す。そして、ｊ＝１であるか否かを判定する（Ｓ４０３）。ここでｊ＝１か否かを判定するのは、増減判定開始の条件として「条件１」が設定されている場合の例である。増減判定開始の条件が「条件１」以外の場合で、順方向の場合には、開始条件として設定された条件となるまで、特に処理を行わず、カウンタｊの値を１ずつ加算すれば良い。上記したように、例えば、開始が「条件４」で、終了が「条件２」といった逆方向の増減判定の場合には、増減表示・パターン生成処理は図２３に例示された内容とは異なった処理となる（カウンタｊの値を順次減算するような処理となる）。
【００９８】
図２３に戻ってステップＳ４０３でＹＥＳの場合、ｊ番目の測定値を開始時閾値と比較する（Ｓ４０４）。以下、閾値について説明する。本実施の形態では、増減パターン５０２等の折れ線グラフ形態のパターンを生成する際に、測定値を閾値と比較する処理を行う。図２４（ａ）は、開始時閾値、開始以後閾値を設定するために表示装置１０２に表示されるウィンドウの一例を示す図である。同図のウィンドウには、開始時閾値を入力するボックス７０１と、開始以後閾値を入力するボックス７０２とが設けられている。同図の例では、開始時閾値として「０．５」、開始以後閾値として「０．２」が設定されている。
【００９９】
本実施の形態では、図１８等に例示されたように、増減パターン５０２等は、四角形の上辺と下辺との間での折れ線グラフ形状として生成される。増減判定開始の条件として「条件１」が設定された場合、「条件１」の測定値が開始時閾値と比較され、「条件１」の測定値が開始時閾値よりも大きい場合、前記四角形の上辺から折れ線グラフがスタートする。開始以後閾値については、後にも説明するが、一つ前の条件の測定値との差分を開始以後閾値と比較して、開始以後閾値よりも大きな測定値の変動があった場合に、折れ線グラフが逆側（上辺側から下辺側、または、下辺側から上辺側）に移動する。
【０１００】
なお、開始時閾値、開始以後閾値は代謝物ごとに異なる値を設定することもできる。図２４（ｂ）は、代謝物ごとに閾値を設定する場合のウィンドウの一例を示す図である。同図のように代謝物、代謝酵素ごとに閾値が設定された場合には、設定された閾値は、代謝物ごとに、例えばＨＤドライブ２０４に格納しておくことができる。
【０１０１】
図２５は、閾値を用いた増減パターン生成について説明するための図である。同図（ａ）の例は、図６に例示した測定値データの中で「alpha-keto-glutarate」の測定値データを用い、開始時閾値０．５、開始以後閾値０．２とした場合の増減パターンである。条件１の測定値は、０，４９であり、開始時閾値を下回るため、増減パターン四角形の下辺側からスタートとなる。条件２の測定値は０．０３と減少しているが、下辺側からスタートしているため、上辺側へのプロットのシフトはない。
【０１０２】
条件３の測定値は０．７５であり、開始以後閾値０．２を上回る増加を示しているため、増減パターンが上限側にシフトする。条件４の測定値は０．０９であり、条件３の測定値と比較すると、開始以後閾値０．２を上回る減少を示しているため、増減パターンが下辺側にシフトする。
【０１０３】
図２３に戻り、処理内容について説明を続ける。ステップＳ４０５において、ｊ＝１の場合の測定値が開始時閾値を上回る場合（Ｓ４０５：ＹＥＳ）には、増減パターンは上辺側からスタートする（Ｓ４０６）。Ｓ４０５の判定結果がＮＯの場合には、増減パターンは下辺側からスタートする（Ｓ４０７）。そしてｊの値に１を加算し（Ｓ４０８）、ｊ番目の測定値が無い場合（増減判定範囲終了に至った場合）までは（Ｓ４０９：ＮＯ）、ステップＳ４０２に戻る。
【０１０４】
ステップＳ４０３での判定結果がＮＯの場合、Ｓ４０２で取得した測定値と、一つ手前（ｊ−１番目）の測定値とを比較する（Ｓ４１０）。増減量が開始以後閾値よりも大きい場合（Ｓ４１１：ＹＥＳ）、原則として逆側（下辺側から上辺側、または、上辺側から下辺側）にシフトする（Ｓ４１２）。もっとも、図２３に図示は省略したが、一つ手前が下辺側で、さらに減少した場合、または、一つ手前が上辺側で、さらに増加した場合には、逆側にシフトしない。逆側にシフトする場合もしない場合も、増減パターンの横軸方向には条件一つ分進行させる。従って、測定値と閾値との関係によっては、増減パターンの枠である四角形（四角形に限定されないことは勿論である。）の範囲内に折れ線が出現しないような場合も有り得る。図２５の（ｂ）及び（ｃ）は、閾値設定により増減パターンの形態が変化する様子を示すものである。
【０１０５】
以上のように、ｉ番目の代謝物（または、代謝酵素）について増減パターンの生成を終了すると（Ｓ４０９：ＹＥＳ）、図の例では、１番目の測定値（条件１）と、ｊ−１番目の測定値（条件４）とを比較する（Ｓ４１３）。これは、増減判定範囲の最初と最後の比較であって、必ずしも１番目ではない場合もあることは前述の通りである。比較結果に応じて円５０１等の色、及び大きさを決定する（Ｓ４１４）。増減量が大きい場合に円の直径を大きく、増減量が小さい場合に円の直径を小さく、というようにすることもできるし、円に限定されないことは勿論である。円の配置は、取得した代謝物、または代謝酵素のマップ上の座標情報を参照して決定される。その後、増減パターンの大きさ及び位置を調整する（Ｓ４１５）。増減パターンの配置位置は、例えば、円５０１等の直径を考慮して決定することができる。増減パターンの大きさは、等しく表示されるように調整することが好ましいが、大きさの調整は必須ではない。
【０１０６】
以上の処理を終了すると、カウンタｉの値に１を加算する（Ｓ４１６）。ｉ番目の代謝物が存在しないか否かを判定しているが（Ｓ４１７）、これは、測定値データが必ずしも代謝パスウェイ上の全ての代謝物について存在するわけではないからである。ｉ番目の代謝物の測定値データが存在する場合（Ｓ４１７：ＮＯ）には、カウンタｊに１（若しくは、増減判定範囲の開始に相当する値）を設定して（Ｓ４１８）、ステップＳ４０２へと進む。ｉ版目の代謝物の測定値データが無い場合（Ｓ４１７：ＹＥＳ）、代謝パスウェイ上の全ての代謝物（及び／又は、代謝酵素）について、増減パターンの生成等の処理を終了したか否かを判定し（Ｓ４１９）、終了していれば（Ｓ４１９：ＹＥＳ）、増減表示・パターン生成処理を終了して図１７のフローチャートに戻る。終了していなければ（Ｓ４１９：ＮＯ）、ステップＳ４１６に進んでカウンタｉに１を加算する。
【０１０７】
図１７に戻り、増減表示・パターン生成処理（Ｓ３０６）の後は、上記に詳細に説明したような閾値の変更がなされたか否かを判定し（Ｓ３０７）、閾値が変更された場合には（Ｓ３０７：ＹＥＳ）、再度、増減表示・パターン生成処理を行う（Ｓ３０６）。上記したように、閾値が変更されると、増減パターンも変化する場合があるからである。
【０１０８】
閾値変更がない場合（Ｓ３０７：ＮＯ）、代謝物の選択入力を受け付けたか否かを判定する（Ｓ３０８）。この代謝物選択は、パスウェイマップ表示領域６０１１において、代謝物を表す丸、増減を表す円５０１等、代謝酵素を表す四角などの上にマウスカーソル１０４１を移動させ、マウスクリック等を行うことで入力される。代謝物選択入力を受け付けた場合には（Ｓ３０８；ＹＥＳ）、図２０に例示したような棒グラフを生成する（Ｓ３０９）。棒グラフの生成は、例えばＪＡＶＡ（登録商標）の機能を利用して実現することができるが、実現方法がＪＡＶＡ（登録商標）に限定されないことは勿論である。
【０１０９】
なお、パスウェイ選択入力を受け付けた場合（Ｓ３１０：ＹＥＳ）には、ステップＳ３０４へと戻る。このパスウェイ選択は、パスウェイ一覧表示領域６０１２において、代謝パスウェイの表示の上にマウスカーソル１０４１を移動させ、マウスクリック等を行うことで入力される。この場合にステップＳ３０４に戻るのは、パスウェイマップ表示領域６０１１に表示されている代謝パスウェイと異なる代謝パスウェイが選択された場合、パスウェイマップ表示領域６０１１の表示内容（増減パターン等を含む。）が変更されるためである。パスウェイ選択入力も無い場合には（Ｓ３１０：ＮＯ）、表示内容生成処理を終了して、図１６のフローチャートに戻る。
【０１１０】
図１６に戻り、表示内容生成処理（Ｓ２０４）の後には、生成された表示内容を実際に表示装置１０２に表示させる表示処理を実行し（Ｓ２０５）、何の選択入力も無い場合（代謝物、代謝パスウェイの選択入力等が無い場合）には（Ｓ２０６：ＮＯ）、支援装置の動作終了指示を受け付けた場合には（Ｓ２０７：ＹＥＳ）、動作を終了する。動作終了指示が無い場合（Ｓ２０７：ＮＯ）、ステップＳ２０６へと戻り、何らかの選択入力が検出された場合に（Ｓ２０６：ＹＥＳ）、ステップＳ２０２へと戻る。ステップＳ２０６の選択入力には、増減判定範囲の選択入力も含まれており、範囲が変更された場合にはステップＳ２０３の判定結果がＹＥＳとなる。
【０１１１】
以上に説明したように、本実施の形態に係る支援装置により、測定値データと代謝パスウェイとの関連性強度を参照しつつ、代謝パスウェイマップ上において代謝物、代謝酵素の増減の様子をビジュアルに表現することができるため、医師、研究者等の専門家による測定値データの分析、検討を支援することができる。測定値を単に基準値と比較したり、コンピュータに格納されたデータベースと比較するような手法では発見できないような疾患等の発見にも寄与することが可能となるものであり、実用的効果は極めて大きい。
【０１１２】
（変形例）
以上、本発明の実施の形態について詳細に説明したが、本発明の範囲が上記実施の形態に限定されることはなく、例えば、以下に例示するような変形例を考えることもできる。
【０１１３】
（１）上記実施の形態では、代謝物について図示したが、代謝酵素についても同様の図示が可能である。例えば代謝酵素を示す四角の近辺に増減パターンを表示するようにすることもできる。なお、図１８等に例示したように、代謝酵素の識別子がマップ上に表示される場合、増減を表す円は、当該識別子を囲む四角の、さらに周辺に表示するようにして、例えば円を表示した場合でも代謝酵素の識別子が判読可能なようにすることもできる。
【０１１４】
（２）上記実施の形態では、代謝物の増減を表現するのに円５０１等を表示したが、形状が円に限定されないことは勿論である。また、当該図形の色も、赤、青に限定されず、任意の色彩を用いることができる。
【０１１５】
（３）上記実施の形態では、「開始時閾値」と「開始以後閾値」とを設定し、測定値変動量が「開始以後閾値」よりも大きい場合に、増減パターンの上辺側、下辺側について逆側にプロットするようにした。もっとも、変動量の閾値ではなく、単に測定値の閾値を設定するような実施形態も可能である。測定値の閾値も、代謝物、代謝酵素ごとに設定するようにすることもできる。図２６は、測定値の閾値を設定した場合の増減パターンの例を示す図である。この場合には、各々の条件ごとに閾値を上回るか否かを判定し、閾値を上回る場合には上辺側、それ以外の場合には下辺側として、変動パターンを生成することができる。この場合にも、設定される閾値により増減パターンの形態が変化し得る。
【０１１６】
（４）上記実施の形態では、増減パターンとして、閾値を用いて上辺、下辺のいずれかにシフトするようなパターンとしたが、このような２段階のパターンだけでなく、３段階、４段階の増減パターンとすることも可能である。例えば、測定値が０．３以下は下辺、０．７以上は上辺とし、その間の測定値の場合には、上辺と下辺との中央にプロットするようにして増減パターンを生成すれば３段階の増減パターンが生成できる。
【０１１７】
（５）本実施の形態の支援装置は、さらに、代謝物、代謝酵素、代謝パスウェイ等の名称の入力を受け付ける名称入力受付手段を備えることができる。名称は、例えばキーボード、ソフトウェアキーボード等の文字列入力装置を用いて入力する他、音声入力、手書き入力等の任意の文字列入力手法を利用することができる。本願発明者らは、名称入力、検索の際に、Ｊａｒｏ−Ｗｉｎｋｌｅｒ距離を利用することが好適であることを確認している。
【０１１８】
（６）本実施の形態の支援装置は、専門化による分析、検討を支援するものであるが、補助的に、コンピュータによるパターン比較を用いることを否定するわけではない。例えば、測定値の増減パターンをデータベース化しておき、類似のパターンが現れた場合に、その旨を表示して、専門家による分析、検討の補助に利用するような実施形態も可能である。より具体的には、健常状態、病態ごとの測定値パターンをデータベース化しておくような形態が可能である。
【０１１９】
（７）上記実施の形態の処理を実現するプログラム言語は何ら限定されないし、必ずしも、上記した処理を行うコードを含んでいないプログラムでも、上記に説明した処理を実現することは可能である。例えば、プログラム言語に用意された機能を活用して処理を実現するような実施形態も考えられる。また、上記フローチャートは実施形態の一例について説明するためのものであって、この処理内容、処理順序等に限定されることはない。
【０１２０】
（８）上記実施の形態では、パスウェイ一覧表示領域に、パスウェイ名称及び関連性強度を一覧表示したが、表示形態は名称の一覧表示に限定されず、例えば、パスウェイ名称と関連性強度が識別可能に付加表示されたアイコンを並べるような表示形態も考えられる。
【産業上の利用可能性】
【０１２１】
本発明は、例えば、診断、検査、薬物投与の影響判断等について、専門家による分析検討を支援する支援装置等に適用することができる。
【符号の説明】
【０１２２】
１０１                  コンピュータ
１０２                  表示装置（ディスプレイ）
１０３                  キーボード
１０４                  マウス
１０４１                マウスカーソル
１０５                  ＬＡＮケーブル
１０６                  ネットワーク
２０１                  ＣＰＵ
２０２                  ＲＯＭ（リード・オンリー・メモリ）
２０３                  ＲＡＭ（ランダム・アクセス・メモリ）
２０４                  ＨＤ（ハードディスク）ドライブ
２０５                  ＣＤ／ＤＶＤドライブ
２０６                  ＮＩＣ（ネットワーク・インタフェース・カード）
３０１                  支援装置制御部
３０１１                代謝経路マップ管理部
３０１２                座標取得部
３０１３                代謝経路データ管理部
３０１４                測定値データ取得部
３０１５                関連性強度算出部
３０１６                選択入力受付部
３０１７                表示内容決定部
３０１８                表示制御部
３０２                  代謝経路マップ格納部
３０３                  代謝経路ＤＢ格納部
３０４                  測定値データ格納部
３０５                  入力装置
４０１                  連続するパスウェイ
５０１、５０３          測定値の増減を示す円
５０２、５０４          測定値増減パターン
６０１                  表示画面
６０１１                パスウェイマップ表示領域
６０１２                パスウェイ一覧表示領域
６０１３                測定値データ表示領域
【図１】
2011055820.tif 000003
【図２】
2011055820.tif 000004
【図３】
2011055820.tif 000005
【図４】
2011055820.tif 000006
【図５】
2011055820.tif 000007
【図６】
2011055820.tif 000008
【図７】
2011055820.tif 000009
【図８】
2011055820.tif 000010
【図９】
2011055820.tif 000011
【図１０】
2011055820.tif 000012
【図１１】
2011055820.tif 000013
【図１２】
2011055820.tif 000014
【図１３】
2011055820.tif 000015
【図１４】
2011055820.tif 000016
【図１５】
2011055820.tif 000017
【図１６】
2011055820.tif 000018
【図１７】
2011055820.tif 000019
【図１８】
2011055820.tif 000020
【図１９】
2011055820.tif 000021
【図２０】
2011055820.tif 000022
【図２１】
2011055820.tif 000023
【図２２】
2011055820.tif 000024
【図２３】
2011055820.tif 000025
【図２４】
2011055820.tif 000026
【図２５】
2011055820.tif 000027
【図２６】
2011055820.tif 000028
【国際調査報告】
2011055820.tif 000029
2011055820.tif 000030
2011055820.tif 000031
2011055820.tif 000032