(57)【要約】【課題】情報間の参照構造に基づく関係性を解析することができる、情報間関係性解析装置、情報間関係性解析方法、プログラムおよび記録媒体を提供することを課題とする。【解決手段】本発明にかかる情報間関係性解析装置は、概念に対応する情報であって他の概念に対応する情報を参照するための参照情報を少なくとも含む概念情報を複数含む概念情報群を取得し、概念情報群の中から、特定概念情報を設定し、概念情報群について、各概念情報を状態とみて各参照情報を通して概念情報を遷移するモデルで捉え、当該モデルをマルコフ過程と考えたときの特定概念情報の存在確率を算出し、概念情報に基づいて、予め設定された評価概念情報を当該評価概念情報に含まれる参照情報と共に無効にしたときの特定概念情報の無効後存在確率を算出し、算出された存在確率および無効後存在確率に基づいて、評価概念情報が特定概念情報に及ぼす影響度情報を算出する。【選択図】        図１
2005275794.tif 000002
【特許請求の範囲】
【請求項１】
概念に対応する情報であって他の上記概念に対応する情報を参照するための参照情報を少なくとも含む概念情報を複数含む概念情報群を取得する概念情報群取得手段と、
上記概念情報群取得手段により取得された上記概念情報群の中から、特定の上記概念情報である特定概念情報を設定する特定概念情報設定手段と、
上記概念情報群取得手段により取得された上記概念情報群について、各々の上記概念情報を単位に当該概念情報を相互に参照する上記参照情報があるとき、各々の上記概念情報を状態とみて各々の上記参照情報を通して上記概念情報を遷移するモデルで捉え、当該モデルをマルコフ過程と考えたときの上記特定概念情報の存在確率を算出する存在確率算出手段と、
上記概念情報群取得手段により取得された上記概念情報群に含まれる上記概念情報に基づいて、評価対象として予め設定された上記概念情報である評価概念情報を当該評価概念情報に含まれる上記参照情報と共に無効にしたときの上記特定概念情報の上記存在確率である無効後存在確率を算出する無効後存在確率算出手段と、
上記存在確率算出手段により算出された上記存在確率および上記無効後存在確率算出手段により算出された上記無効後存在確率に基づいて、上記評価概念情報が上記特定概念情報設定手段により設定された上記特定概念情報に及ぼす影響度に関する情報である影響度情報を算出する影響度算出手段と、
を備えたことを特徴とする情報間関係性解析装置。
【請求項２】
上記存在確率算出手段および上記無効後存在確率算出手段は、数式１に基づいて、それぞれ、上記存在確率および上記無効後存在確率を算出すること、
【数１】
2005275794.tif 000022
（ここで、「Ｒ（ｐ）」は概念情報識別情報「ｐ」に対応する概念情報の存在確率であり、「Ｒ（ｑ）」は概念情報識別情報「ｐ」に対応する概念情報への参照情報をもつ概念情報識別情報「ｑ」に対応する概念情報の存在確率であり、「Ｃ（ｑ）」は概念情報識別情報「ｑ」に対応する概念情報から他の概念情報への参照情報の数である。また、「Ｇ」は概念情報をノード、参照情報をエッジとしたグラフであり、「ｎ」はグラフ「Ｇ」のノードである概念情報の総数である。「ε」は予め設定される定数である。）
を特徴とする請求項１に記載の情報間関係性解析装置。
【請求項３】
上記影響度算出手段は、数式２に示すように、上記存在確率と上記無効後存在確率との差に基づいて、上記影響度情報を算出すること、
【数２】
2005275794.tif 000023
（ここで、「ｄＲ（ｐ）(s)」は評価概念情報識別情報「ｓ」に対応する評価概念情報が注目概念情報識別情報「ｐ」対応する注目概念情報に及ぼす影響度情報である。「Ｒ（ｐ）(s)」は、概念情報識別情報「ｓ」に対応する概念情報を参照情報と共に無効にしたときの、概念情報識別情報「ｐ」に対応する概念情報の存在確率である無効後存在確率であり、「Ｒ（ｐ）」は、概念情報識別情報「ｐ」に対応する概念情報の存在確率である。）
を特徴とする請求項１または２に記載の情報間関係性解析装置。
【請求項４】
上記特定概念情報設定手段により設定された上記特定概念情報が複数の場合、上記影響度算出手段により算出された上記影響度情報に基づいて、上記特定概念情報の上記影響度情報からの第１主成分得点を算出する主成分得点算出手段、
をさらに備えたことを特徴とする請求項１〜３のいずれか一つに記載の情報間関係性解析装置。
【請求項５】
上記影響度算出手段により算出された上記影響度情報および／または上記主成分得点算出手段により算出された上記第１主成分得点に基づいて、上記特定概念情報および／または上記評価概念情報を序列する序列手段、
をさらに備えたことを特徴とする請求項１〜４のいずれか一つに記載の情報間関係性解析装置。
【請求項６】
概念に対応する情報であって他の上記概念に対応する情報を参照するための参照情報を少なくとも含む概念情報を複数含む概念情報群を取得する概念情報群取得ステップと、
上記概念情報群取得ステップにより取得された上記概念情報群の中から、特定の上記概念情報である特定概念情報を設定する特定概念情報設定ステップと、
上記概念情報群取得ステップにより取得された上記概念情報群について、各々の上記概念情報を単位に当該概念情報を相互に参照する上記参照情報があるとき、各々の上記概念情報を状態とみて各々の上記参照情報を通して上記概念情報を遷移するモデルで捉え、当該モデルをマルコフ過程と考えたときの上記特定概念情報の存在確率を算出する存在確率算出ステップと、
上記概念情報群取得ステップにより取得された上記概念情報群に含まれる上記概念情報に基づいて、評価対象として予め設定された上記概念情報である評価概念情報を当該評価概念情報に含まれる上記参照情報と共に無効にしたときの上記特定概念情報の上記存在確率である無効後存在確率を算出する無効後存在確率算出ステップと、
上記存在確率算出ステップにより算出された上記存在確率および上記無効後存在確率算出ステップにより算出された上記無効後存在確率に基づいて、上記評価概念情報が上記特定概念情報設定ステップにより設定された上記特定概念情報に及ぼす影響度に関する情報である影響度情報を算出する影響度算出ステップと、
を含むことを特徴とする情報間関係性解析方法。
【請求項７】
上記存在確率算出ステップおよび上記無効後存在確率算出ステップは、数式１に基づいて、それぞれ、上記存在確率および上記無効後存在確率を算出すること、
【数３】
2005275794.tif 000024
（ここで、「Ｒ（ｐ）」は概念情報識別情報「ｐ」に対応する概念情報の存在確率であり、「Ｒ（ｑ）」は概念情報識別情報「ｐ」に対応する概念情報への参照情報をもつ概念情報識別情報「ｑ」に対応する概念情報の存在確率であり、「Ｃ（ｑ）」は概念情報識別情報「ｑ」に対応する概念情報から他の概念情報への参照情報の数である。また、「Ｇ」は概念情報をノード、参照情報をエッジとしたグラフであり、「ｎ」はグラフ「Ｇ」のノードである概念情報の総数である。「ε」は予め設定される定数である。）
を特徴とする請求項６に記載の情報間関係性解析方法。
【請求項８】
上記影響度算出ステップは、数式２に示すように、上記存在確率と上記無効後存在確率との差に基づいて、上記影響度情報を算出すること、
【数４】
2005275794.tif 000025
（ここで、「ｄＲ（ｐ）(s)」は評価概念情報識別情報「ｓ」に対応する評価概念情報が注目概念情報識別情報「ｐ」対応する注目概念情報に及ぼす影響度情報である。「Ｒ（ｐ）(s)」は、概念情報識別情報「ｓ」に対応する概念情報を参照情報と共に無効にしたときの、概念情報識別情報「ｐ」に対応する概念情報の存在確率である無効後存在確率であり、「Ｒ（ｐ）」は、概念情報識別情報「ｐ」に対応する概念情報の存在確率である。）
を特徴とする請求項６または７に記載の情報間関係性解析方法。
【請求項９】
上記特定概念情報設定ステップにより設定された上記特定概念情報が複数の場合、上記影響度算出ステップにより算出された上記影響度情報に基づいて、上記特定概念情報の上記影響度情報からの第１主成分得点を算出する主成分得点算出ステップ、
をさらに含むことを特徴とする請求項６〜８のいずれか一つに記載の情報間関係性解析方法。
【請求項１０】
上記影響度算出ステップにより算出された上記影響度情報および／または上記主成分得点算出ステップにより算出された上記第１主成分得点に基づいて、上記特定概念情報および／または上記評価概念情報を序列する序列ステップ、
をさらに含むことを特徴とする請求項６〜９のいずれか一つに記載の情報間関係性解析方法。
【請求項１１】
概念に対応する情報であって他の上記概念に対応する情報を参照するための参照情報を少なくとも含む概念情報を複数含む概念情報群を取得する概念情報群取得ステップと、
上記概念情報群取得ステップにより取得された上記概念情報群の中から、特定の上記概念情報である特定概念情報を設定する特定概念情報設定ステップと、
上記概念情報群取得ステップにより取得された上記概念情報群について、各々の上記概念情報を単位に当該概念情報を相互に参照する上記参照情報があるとき、各々の上記概念情報を状態とみて各々の上記参照情報を通して上記概念情報を遷移するモデルで捉え、当該モデルをマルコフ過程と考えたときの上記特定概念情報の存在確率を算出する存在確率算出ステップと、
上記概念情報群取得ステップにより取得された上記概念情報群に含まれる上記概念情報に基づいて、評価対象として予め設定された上記概念情報である評価概念情報を当該評価概念情報に含まれる上記参照情報と共に無効にしたときの上記特定概念情報の上記存在確率である無効後存在確率を算出する無効後存在確率算出ステップと、
上記存在確率算出ステップにより算出された上記存在確率および上記無効後存在確率算出ステップにより算出された上記無効後存在確率に基づいて、上記評価概念情報が上記特定概念情報設定ステップにより設定された上記特定概念情報に及ぼす影響度に関する情報である影響度情報を算出する影響度算出ステップと、
を含む情報間関係性解析方法をコンピュータに実行させることを特徴とするプログラム。
【請求項１２】
上記存在確率算出ステップおよび上記無効後存在確率算出ステップは、数式１に基づいて、それぞれ、上記存在確率および上記無効後存在確率を算出すること、
【数５】
2005275794.tif 000026
（ここで、「Ｒ（ｐ）」は概念情報識別情報「ｐ」に対応する概念情報の存在確率であり、「Ｒ（ｑ）」は概念情報識別情報「ｐ」に対応する概念情報への参照情報をもつ概念情報識別情報「ｑ」に対応する概念情報の存在確率であり、「Ｃ（ｑ）」は概念情報識別情報「ｑ」に対応する概念情報から他の概念情報への参照情報の数である。また、「Ｇ」は概念情報をノード、参照情報をエッジとしたグラフであり、「ｎ」はグラフ「Ｇ」のノードである概念情報の総数である。「ε」は予め設定される定数である。）
を特徴とする請求項１１に記載のプログラム。
【請求項１３】
上記影響度算出ステップは、数式２に示すように、上記存在確率と上記無効後存在確率との差に基づいて、上記影響度情報を算出すること、
【数６】
2005275794.tif 000027
（ここで、「ｄＲ（ｐ）(s)」は評価概念情報識別情報「ｓ」に対応する評価概念情報が注目概念情報識別情報「ｐ」対応する注目概念情報に及ぼす影響度情報である。「Ｒ（ｐ）(s)」は、概念情報識別情報「ｓ」に対応する概念情報を参照情報と共に無効にしたときの、概念情報識別情報「ｐ」に対応する概念情報の存在確率である無効後存在確率であり、「Ｒ（ｐ）」は、概念情報識別情報「ｐ」に対応する概念情報の存在確率である。）
を特徴とする請求項１１または１２に記載のプログラム。
【請求項１４】
上記特定概念情報設定ステップにより設定された上記特定概念情報が複数の場合、上記影響度算出ステップにより算出された上記影響度情報に基づいて、上記特定概念情報の上記影響度情報からの第１主成分得点を算出する主成分得点算出ステップ、
をさらに含むことを特徴とする請求項１１〜１３のいずれか一つに記載のプログラム。
【請求項１５】
上記影響度算出ステップにより算出された上記影響度情報および／または上記主成分得点算出ステップにより算出された上記第１主成分得点に基づいて、上記特定概念情報および／または上記評価概念情報を序列する序列ステップ、
をさらに含むことを特徴とする請求項１１〜１４のいずれか一つに記載のプログラム。
【請求項１６】
上記請求項１１〜１５のいずれか一つに記載されたプログラムを記録したことを特徴とするコンピュータ読み取り可能な記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、情報間関係性解析装置、情報間関係性解析方法、プログラムおよび記録媒体に関し、特に、相互に参照される情報（例えば、リンクを張るＷｅｂページ）群においてターゲットとする情報を無効にしたときの他への影響を算出し、互いの関係性を定量化することができる、情報間関係性解析装置、情報間関係性解析方法、プログラムおよび記録媒体に関するものである。
【背景技術】
【０００２】
昨今の電子化の流れに伴って、情報量は増加し、これまで個別に扱われていた情報を結び付けて全体を系として捉えることが計算機性能の向上により可能な状況となってきている。
【０００３】
ここで、情報の効果的活用に向けて、それらの情報間の依存関係をはじめとした未知の関係性の発見が求められている。
【０００４】
例えば、大規模な情報ネットワークの解析には、因果関係に基づく現象の原因解明および問題に対する対策の効果測定への期待があり、例えば、疾患に関連する遺伝子の探索もその情報ネットワークの解析の典型例の１つと捉えられる。
【０００５】
また、事象全体を系として捉える解析の重要性が認識されつつあり（例えば、非特許文献１参照。）、ライフサイエンス分野（例えば、非特許文献２、非特許文献３、非特許文献４参照。）においても、多くの変数とパラメータを関係付け、生物学的現象の計算機処理を図るシステム生物学への期待が高まってきている。
【０００６】
また、創薬においても、個別の化合物の特定から生理化学的ネットワークの同定への重要性が認識され始めてきている（例えば、非特許文献５参照。）。
【０００７】
ところで、これまで、疾患関連遺伝子は、遺伝子の位置を疾患との相関に基づいて領域を狭めていくことにより同定する方法であるポジショナルクローニングから得られている。このポジショナルクローニングと呼ばれる方法は、家系が付与された試料の収集をはじめ、多くの時間と労力を要する生物学的実験に基づいている。
【０００８】
また、生物学的実験において、特定の遺伝子をノックアウト、即ち機能欠損させ、症状を観察して遺伝子推定を図る、ノックアウトスタディと呼ばれるアプローチもよく知られている（例えば、非特許文献６参照。）。具体的には、例えばマウスなどの実験動物を用いて、注目する遺伝子を破壊（ノックアウト）することにより機能欠損させ、症状を観察して遺伝子機能の推定を図る。そして、実験動物との間の遺伝子配列の相同性を通してヒトの遺伝子機能が推定される。
【０００９】
つまり、従来の遺伝子機能解析の手段としては、生物学的実験が主体で断片的なものであった（例えば、非特許文献７参照。）。例えば、疾患関連遺伝子の抽出は、生物学的実験により得られるのが一般的で、蓄積されてきている。また、ポジショナルクローニングにおいては、特に、疾患の原因が一箇所ではない多因子性疾患では莫大な数の試料が必要とされ（例えば、非特許文献８参照。）、疾患の原因となる遺伝子探索への適用性に課題がある。また、遺伝子機能は生成される酵素による制御からもたらされるコンテクストに基づく（例えば、非特許文献９参照。）ことから、特に高血圧や糖尿病、喘息などの多因子遺伝子病の遺伝的素因を理解するためには、複数ある疾患関連遺伝子の組み合わせとその相互作用が総合的に考慮されなければならない（例えば、非特許文献１０参照。）とされている。
【００１０】
この状況の下、疾患の原因となる遺伝子探索は、薬剤開発におけるターゲット候補選定など、産業上の理由からも、中心的な課題に位置付けられている。
【００１１】
そこで、一般的に、情報間の依存関係などの関係性を解析する従来の手法としては、ＣＲＭ（Ｃｕｓｔｏｍｏｒ  Ｒｅｌａｔｉｏｎｓｈｉｐ  Ｍａｎａｇｅｍｅｎｔ）分野において、購買履歴データからどの商品が一緒に買われるかを分析するバスケット分析がよく知られている（例えば、非特許文献１１参照。）。
【００１２】
また、情報処理による疾患関連遺伝子の抽出の方法としては、既知の疾患関連遺伝子を基に、塩基配列の相同性による推定が挙げられる。
【００１３】
また、例えば、従来技術である非特許文献１２は、医学生物学的知見が論文の形でまとめられているとみて、論文データベースの「ＰｕｂＭｅｄ（例えば、非特許文献１３参照。）」を用いて、遺伝子、化合物、疾患の３つの間の関係を論文中の単語共起に基づいて求める技術について記載しており、疾患関連遺伝子の探索に活用し得ると考えられる。
【００１４】
また、論文中における遺伝子名の共起頻度により、網羅的に遺伝子間の関係を可視化する試み（例えば、非特許文献１４参照。）はよく知られ、広く活用されている。
【００１５】
【非特許文献１】アルバート＝ラズロ・バラバシ，  “新ネットワーク思考−世界のしくみを読み解く”，  ＮＨＫ出版，  ２００２
【非特許文献２】“２１世紀の医療・福祉を支える科学技術特集／５．ゲノム情報とその利用”，  信学誌，  Ｖｏｌ．８４，  Ｎｏ．５，  ｐｐ．３４１−３６７，  ２００１
【非特許文献３】森下  真一，  久光  徹，  高木  利久，  “特集  ゲノム情報科学”，  情報処理，  Ｖｏｌ．４３，  Ｎｏ．１，  ｐｐ．１−４１，  ２００２
【非特許文献４】阿久津  達也，  麻生川  稔，  小長谷  明彦，  “特集「分子生物情報学の新展開」”，  人工知能学会誌，  Ｖｏｌ．１５，  Ｎｏ．１，  ｐｐ．２−５５，  ２０００
【非特許文献５】Ｓ．Ｈｕａｎｇ，  “Ｒａｔｉｏｎａｌ  ｄｒｕｇ  ｄｉｓｃｏｖｅｒｙ：  ｗｈａｔ  ｃａｎ  ｗｅ  ｌｅａｒｎ  ｆｒｏｍ  ｒｅｇｕｌａｔｏｒｙ  ｎｅｔｗｏｒｋｓ？”，  Ｄｒｕｇ  Ｄｉｓｃｏｖ．  Ｔｏｄａｙ，  ｖｏｌ．１５，  Ｎｏ．７，  ｐｐＳ１６３−９，  ２００２
【非特許文献６】Ｂ．Ａｌｂｅｒｔｓ  ｅｔ  ａｌ．，  “Ｍｏｌｅｃｕｌａｒ  Ｂｉｏｌｏｇｙ  ｏｆ  Ｔｈｅ  Ｃｅｌｌ”，  Ｇａｒｌａｎｄ  Ｐｕｂｌｉｓｈｉｎｇ  Ｎｅｗ  Ｙｏｒｋ，  １９９４
【非特許文献７】中村  祐輔，  “概論−個の遺伝子研究からゲノムワイドの遺伝子研究へ”，  実験医学，  Ｖｏｌ．１８，  Ｎｏ．１２，  ｐｐ．１８−２０，  ２０００
【非特許文献８】安部  道子，  名倉  潤，  三木  哲郎，  “高血圧感受性遺伝子を探して”，  医学のあゆみ，  Ｖｏｌ．１９７，  Ｎｏ．１３，  ｐｐ．１１１７−１１２１，  ２００１
【非特許文献９】Ｓ．Ｈｕａｎｇ，  “Ｔｈｅ  ｐｒａｃｔｉｃａｌ  ｐｒｏｂｌｅｍｓ  ｏｆ  ｐｏｓｔ−ｇｅｎｏｍｅ  ｂｉｏｌｏｇｙ”，  Ｎａｔ．  Ｂｉｏｔｅｃｈ．，  ｖｏｌ．１８，  ｐｐ４７１−４７２，  ２０００
【非特許文献１０】鈴木  洋一，  “疾患関連遺伝子の総合的解析と相互作用”，  メディカル・サイエンス・ダイジェスト，  Ｖｏｌ．２９，  Ｎｏ．７，  ｐｐ．２１−２４，  ２００３
【非特許文献１１】山西  健司，  “データ・テキストマイニング”，  計算工学，  ６，  ｐｐ．３８６−３９３，  ２００１
【非特許文献１２】株式会社ワールドフュージョンが開発したテキストマイニングシステム「構ｋｏｚｏ造  ｖｅｒ１」のホームページ「ｈｔｔｐ：／／ｗｗｗ．ｗ−ｆｕｓｉｏｎ．ｃｏ．ｊｐ／ｗｆ／ｋｏｚｏ．ｈｔｍｌ」
【非特許文献１３】論文データベースＰｕｂＭｅｄのホームページ「ｈｔｔｐ：／／ｗｗｗ．ｎｃｂｉ．ｎｌｍ．ｎｉｈ．ｇｏｖ／ｅｎｔｒｅｚ／ｑｕｅｒｙ．ｆｃｇｉ？ｄｂ＝ＰｕｂＭｅｄ」
【非特許文献１４】ＴＫ．Ｊｅｎｓｓｅｎ，  Ａ．Ｌａｅｇｒｅｉｄ，  Ｊ．Ｋｏｍｏｒｏｗｓｋｉ，  Ｅ．Ｈｏｖｉｇ，  “Ａ  ｌｉｔｅｒａｔｕｒｅ  ｎｅｔｗｏｒｋ  ｏｆ  ｈｕｍａｎ  ｇｅｎｅｓ  ｆｏｒ  ｈｉｇｈ−ｔｈｒｏｕｇｈｐｕｔ  ａｎａｌｙｓｉｓ  ｏｆ  ｇｅｎｅ  ｅｘｐｒｅｓｓｉｏｎ”，  Ｎａｔ．Ｇｅｎｅｔ．，  ｖｏｌ．２８，  Ｎｏ．５，  ｐｐ２１−２８，  ２００１
【発明の開示】
【発明が解決しようとする課題】
【００１６】
しかしながら、非特許文献１１に記載の手法は、事象の同時出現という視点によるものであり、背景にある構造が反映されていない点で潜在的な関係性の発見に限界がある、という問題点があった。
【００１７】
また、非特許文献１２に記載の従来技術では、遺伝子を個別にみていることから、抽出が局所的な範囲に限られており、全体を系とみて遺伝子間の相互作用が考慮された疾患関連遺伝子の抽出が果たされていない、という問題点があった。
【００１８】
また、非特許文献１４に記載の従来技術は、疾患情報は含まれないことから、疾患関連遺伝子の抽出に直接つながるものではない、という問題点があった。
【００１９】
本発明は上記問題点に鑑みてなされたもので、相互に参照される情報（例えば、リンクを張るＷｅｂページ）群においてターゲットとする情報を無効にしたときの他への影響を算出し、互いの関係性を定量化することができる、情報間関係性解析装置、情報間関係性解析方法、プログラムおよび記録媒体を提供することを目的としている。
【課題を解決するための手段】
【００２０】
本発明では、情報間の参照構造に基づく関係性の解析を提案する。具体的には、相互に参照される情報群において、ターゲットとする情報を無効にしたときの他への影響を算出し、互いの関係性を定量化する。
【００２１】
また、本発明は、生物学的実験においてよく知られるノックアウトスタディの基本的な考え方を情報処理手法として導入したものとみることができ、観察結果を得る手段には状態遷移モデル（例えば、「Ｌ．Ｐａｇｅ，  Ｓ．Ｂｒｉｎ，  Ｒ．Ｍｏｔｗａｎｉ  ａｎｄ  Ｔ．Ｗｉｎｏｇｒａｄ，  “Ｔｈｅ  ＰａｇｅＲａｎｋ  Ｃｉｔａｔｉｏｎ  Ｒａｎｋｉｎｇ：  Ｂｒｉｎｇｉｎｇ  Ｏｒｄｅｒ  ｔｏ  ｔｈｅ  Ｗｅｂ”，  Ｓｔａｎｆｏｒｄ  Ｄｉｇｉｔａｌ  Ｌｉｂｒａｒｙ  Ｔｅｃｈｎｏｌｏｇｉｅｓ，  Ｗｏｒｋｉｎｇ  Ｐａｐｅｒ  ＳＩＤＬ−ＷＰ−１９９９−０１２０，  １９９８」参照。）を採用する。
【００２２】
なお、大規模情報ネットワークにおける解析において、リンクによる参照構造を用いたＷｅｂページの重要度評価方法（例えば、「Ｌ．Ｐａｇｅ，  Ｓ．Ｂｒｉｎ，  Ｒ．Ｍｏｔｗａｎｉ  ａｎｄ  Ｔ．Ｗｉｎｏｇｒａｄ，  “Ｔｈｅ  ＰａｇｅＲａｎｋ  Ｃｉｔａｔｉｏｎ  Ｒａｎｋｉｎｇ：  Ｂｒｉｎｇｉｎｇ  Ｏｒｄｅｒ  ｔｏ  ｔｈｅ  Ｗｅｂ”，  Ｓｔａｎｆｏｒｄ  Ｄｉｇｉｔａｌ  Ｌｉｂｒａｒｙ  Ｔｅｃｈｎｏｌｏｇｉｅｓ，  Ｗｏｒｋｉｎｇ  Ｐａｐｅｒ  ＳＩＤＬ−ＷＰ−１９９９−０１２０，  １９９８」参照。）は、膨大なＷｅｂページの中から、有用なものを上位にランキングするための有効な手段としてよく知られ、インターネットにおけるサーチエンジンに適用されている（例えば、「山名  早人，  近藤  秀和，  “解説  サーチエンジンＧｏｏｇｌｅ”，  情報処理，  Ｖｏｌ．４２，  Ｎｏ．８，  ｐｐ．７７５−７８０，  ２００１」参照。）。
【００２３】
上述の方法は、「多くの重要なＷｅｂページからの被リンク数が多いＷｅｂページほど、重要なＷｅｂページである」という経験則に基づいており、全てのページについて注目の度合いが数値化される。ここで、「注目の度合い」とは、他ページからのリンクによる「支持票」の集まるページは注目に値するということであり、必ずしも“重要”や“良好”などの基準が導入されているわけではない。そして、ページ間の相互参照関係より、グラフ理論の応用から定常状態遷移の確率分布で定式化し、確率の量で算出される。
【００２４】
すなわち、本発明は、相互参照構造をもつ情報ネットワークの一部に摂動を与え、系における変動を状態遷移モデルからの確率値を通して定量化するものである。
【００２５】
そして、本発明によれば、疾患と遺伝子群が医学生物学的結びつきによりリンク構造をもつネットワークにおいて、遺伝子をノックアウト、即ち当該遺伝子をリンク構造から除外し、注目する疾患に関する変動から遺伝子の疾患に対する寄与を序列化することにより、疾患関連遺伝子が自動抽出される。つまり、本発明は、上述のノックアウトスタディに基づくアプローチを導入したものといえる。
【００２６】
また、本発明は、遺伝子の「票」の疾患への「支持」のつながり度合いから、その遺伝子の疾患との関連の深さを把握しようとするところに着眼点がある。なお、サーチエンジンでは、全ページのランキングがねらいであり、着眼点を異にしている。
【００２７】
本発明は、相互に参照される情報群において、ターゲットとする情報を無効にしたときの他への影響を算出し、互いの関係性の定量化を図るものである。まず、本発明の基本原理を以下に説明する。
【００２８】
ページを単位にそれらを相互に参照するリンクがあるとき、ページを状態とみてリンクを通してユーザがページを遷移するモデルで捉えることができ（例えば、「Ｌ．Ｐａｇｅ，  Ｓ．Ｂｒｉｎ，  Ｒ．Ｍｏｔｗａｎｉ  ａｎｄ  Ｔ．Ｗｉｎｏｇｒａｄ，  “Ｔｈｅ  ＰａｇｅＲａｎｋ  Ｃｉｔａｔｉｏｎ  Ｒａｎｋｉｎｇ：  Ｂｒｉｎｇｉｎｇ  Ｏｒｄｅｒ  ｔｏ  ｔｈｅ  Ｗｅｂ”，  Ｓｔａｎｆｏｒｄ  Ｄｉｇｉｔａｌ  Ｌｉｂｒａｒｙ  Ｔｅｃｈｎｏｌｏｇｉｅｓ，  Ｗｏｒｋｉｎｇ  Ｐａｐｅｒ  ＳＩＤＬ−ＷＰ−１９９９−０１２０，  １９９８」、「山名  早人，  近藤  秀和，  “解説  サーチエンジンＧｏｏｇｌｅ”，  情報処理，  Ｖｏｌ．４２，  Ｎｏ．８，  ｐｐ．７７５−７８０，  ２００１」参照。）、マルコフ過程と考えたときのページｐのユーザの存在確率Ｒ（ｐ）は、次式（１）を繰り返して適用することにより求められる。すなわち、ページ群からなるネットワークにて、次式（１）によりＲ（ｐ）を算出する。具体的には、状態間の推移確率に関する行列の固有値問題を解くことにより求められる（例えば、「Ｌ．Ｐａｇｅ，  Ｓ．Ｂｒｉｎ，  Ｒ．Ｍｏｔｗａｎｉ  ａｎｄ  Ｔ．Ｗｉｎｏｇｒａｄ，  “Ｔｈｅ  ＰａｇｅＲａｎｋ  Ｃｉｔａｔｉｏｎ  Ｒａｎｋｉｎｇ：  Ｂｒｉｎｇｉｎｇ  Ｏｒｄｅｒ  ｔｏ  ｔｈｅ  Ｗｅｂ”，  Ｓｔａｎｆｏｒｄ  Ｄｉｇｉｔａｌ  Ｌｉｂｒａｒｙ  Ｔｅｃｈｎｏｌｏｇｉｅｓ，  Ｗｏｒｋｉｎｇ  Ｐａｐｅｒ  ＳＩＤＬ−ＷＰ−１９９９−０１２０，  １９９８」、「ｈｔｔｐ：／／ｗｗｗ．ｋｕｓａｓｔｒｏ．ｋｙｏｔｏ−ｕ．ａｃ．ｊｐ／ｂａｂａ／ｗａｉｓ／ｐａｇｅｒａｎｋ．ｈｔｍｌ」参照。）。
【数１】
2005275794.tif 000003
【００２９】
ここで、Ｒ（ｑ）は、ページｐへのリンクをもつページｑのユーザの存在確率である。ｎは、対象とするグラフＧ、即ち、各ページをノード、ページ間のリンクをエッジとしたグラフのノード総数（ページ数）、Ｃ（ｑ）は、ページｑから出ていく外向きのリンク数である。また、存在確率Ｒ（ｐ）は、次式（２）を満たす。
【数２】
2005275794.tif 000004
【００３０】
すなわち、存在確率Ｒ（ｐ）は、ランダムにリンクをたどるユーザがページｐを訪れる確率に対応し、リンクを張っている他のページの値を基に決められる。ここでは（および後述する本実施例では）、ユーザは、１−εの確率で現在のページからのリンクをたどり、εの確率で全く無関係なページへジャンプするモデルを考えることとし、ε＝０．１５とする（例えば、「Ｌ．Ｐａｇｅ，  Ｓ．Ｂｒｉｎ，  Ｒ．Ｍｏｔｗａｎｉ  ａｎｄ  Ｔ．Ｗｉｎｏｇｒａｄ，  “Ｔｈｅ  ＰａｇｅＲａｎｋ  Ｃｉｔａｔｉｏｎ  Ｒａｎｋｉｎｇ：  Ｂｒｉｎｇｉｎｇ  Ｏｒｄｅｒ  ｔｏ  ｔｈｅ  Ｗｅｂ”，  Ｓｔａｎｆｏｒｄ  Ｄｉｇｉｔａｌ  Ｌｉｂｒａｒｙ  Ｔｅｃｈｎｏｌｏｇｉｅｓ，  Ｗｏｒｋｉｎｇ  Ｐａｐｅｒ  ＳＩＤＬ−ＷＰ−１９９９−０１２０，  １９９８」、「ｈｔｔｐ：／／ｗｗｗ．ｋｕｓａｓｔｒｏ．ｋｙｏｔｏ−ｕ．ａｃ．ｊｐ／ｂａｂａ／ｗａｉｓ／ｐａｇｅｒａｎｋ．ｈｔｍｌ」参照。）。なお、このモデルでは、全ての状態が相互到達可能な既約であるマルコフ連鎖を考えることにより解の一意性が確保される（例えば、「Ｌ．Ｐａｇｅ，  Ｓ．Ｂｒｉｎ，  Ｒ．Ｍｏｔｗａｎｉ  ａｎｄ  Ｔ．Ｗｉｎｏｇｒａｄ，  “Ｔｈｅ  ＰａｇｅＲａｎｋ  Ｃｉｔａｔｉｏｎ  Ｒａｎｋｉｎｇ：  Ｂｒｉｎｇｉｎｇ  Ｏｒｄｅｒ  ｔｏ  ｔｈｅ  Ｗｅｂ”，  Ｓｔａｎｆｏｒｄ  Ｄｉｇｉｔａｌ  Ｌｉｂｒａｒｙ  Ｔｅｃｈｎｏｌｏｇｉｅｓ，  Ｗｏｒｋｉｎｇ  Ｐａｐｅｒ  ＳＩＤＬ−ＷＰ−１９９９−０１２０，  １９９８」参照。）。そして、どこからもリンクをもたないページのユーザの存在確率はε／ｎとなる。
【００３１】
そして、特定のページｓをノックアウト、即ち、ページを参照関係とともに無効にしたときに得られるページｐのユーザの存在確率をＲ（ｐ）(s)とするとき、ページｓのページｐに対する関係は、次式（３）のｄＲ（ｐ）(s)により定量化される。そして、特定のページｓ分だけ、算出を繰り返す。
【数３】
2005275794.tif 000005
【００３２】
ここで、図５は、７つのページとそれらの相互参照関係の例（例えば、「ｈｔｔｐ：／／ｗｗｗ．ｋｕｓａｓｔｒｏ．ｋｙｏｔｏ−ｕ．ａｃ．ｊｐ／ｂａｂａ／ｗａｉｓ／ｐａｇｅｒａｎｋ．ｈｔｍｌ」参照。）であり、矢印がリンクによる参照を表している。そして、図５におけるＲと、ページ＃２、ページ＃３およびページ＃７のノックアウトによるｄＲ(#2)、ｄＲ(#3)およびｄＲ(#7)を図６に示す。そして、図６には、各ページから外向きに出るリンク数および他のページから入ってくる被リンク数のノックアウト前後の差をあわせて示す。
【００３３】
図６に示すように、ページ＃５の関係をみる場合、ページ＃２とは０．０５９、ページ＃３とはノックアウトで外向きリンクが１つ減ることから０．０７３、ページ＃７とは被リンクが１つなくなることから−０．０３３が得られている。すなわち、ｄＲの絶対値を通してページ間の関係の大きさを理解することができる。
【００３４】
したがって、上述した目的を達成するために、請求項１に記載の情報間関係性解析装置は、概念に対応する情報であって他の上記概念に対応する情報を参照するための参照情報を少なくとも含む概念情報を複数含む概念情報群を取得する概念情報群取得手段と、上記概念情報群取得手段により取得された上記概念情報群の中から、特定の上記概念情報である特定概念情報を設定する特定概念情報設定手段と、上記概念情報群取得手段により取得された上記概念情報群について、各々の上記概念情報を単位に当該概念情報を相互に参照する上記参照情報があるとき、各々の上記概念情報を状態とみて各々の上記参照情報を通して上記概念情報を遷移するモデルで捉え、当該モデルをマルコフ過程と考えたときの上記特定概念情報の存在確率を算出する存在確率算出手段と、上記概念情報群取得手段により取得された上記概念情報群に含まれる上記概念情報に基づいて、評価対象として予め設定された上記概念情報である評価概念情報を当該評価概念情報に含まれる上記参照情報と共に無効にしたときの上記特定概念情報の上記存在確率である無効後存在確率を算出する無効後存在確率算出手段と、上記存在確率算出手段により算出された上記存在確率および上記無効後存在確率算出手段により算出された上記無効後存在確率に基づいて、上記評価概念情報が上記特定概念情報設定手段により設定された上記特定概念情報に及ぼす影響度に関する情報である影響度情報を算出する影響度算出手段とを備えたことを特徴とする。
【００３５】
この装置によれば、概念（例えば、疾患や遺伝子など）に対応する情報であって他の概念に対応する情報を参照するための参照情報（例えば、Ｗｅｂページにおけるリンクなど）を少なくとも含む概念情報（例えば、Ｗｅｂページなど）を複数含む概念情報群を取得し、取得された概念情報群の中から、特定の概念情報である特定概念情報を設定し、取得された概念情報群について、各々の概念情報を単位に当該概念情報を相互に参照する参照情報があるとき、各々の概念情報を状態とみて各々の参照情報を通して概念情報を遷移するモデルで捉え、当該モデルをマルコフ過程と考えたときの特定概念情報の存在確率を算出し、取得された概念情報群に含まれる概念情報に基づいて、評価対象として予め設定された概念情報である評価概念情報を当該評価概念情報に含まれる参照情報と共に無効にしたときの特定概念情報の存在確率である無効後存在確率を算出し、算出された存在確率および算出された無効後存在確率に基づいて、評価概念情報が特定概念情報に及ぼす影響度に関する情報である影響度情報を算出するので、概念情報間の参照情報に基づく関係性の解析ができる。すなわち、この装置によれば、相互に参照される概念情報群において、ターゲットとする概念情報（評価概念情報）を無効にしたときの他への影響を算出するので、互いの関係性を定量化することができる。また、この装置によれば、参照情報に基づいて互いの概念情報の関係性を定量化しているので、背景にある構造を反映しており、情報間の潜在的な関係性を発見することができる。
【００３６】
また、この装置によれば、例えば疾患関連遺伝子の抽出において、遺伝子を個別にみず、全体を系とみているので、遺伝子間の相互作用が考慮された疾患関連遺伝子を局所的な範囲に限定せずに抽出することができる。具体的には、この装置によれば、例えば、疾患と遺伝子群が医学生物学的結び付きによりリンク構造をもつネットワークを状態遷移モデルで捉え、或る状態（例えば、遺伝子）を無効にしたときの、疾患に相当する状態へ与える変動から、遺伝子の疾患に対する寄与を定量化することができる。また、この装置によれば、例えば、潜在的な疾患関連遺伝子の発見の他、疾患の発症に強く関与する主効果遺伝子を抽出することができる。
【００３７】
また、この装置によれば、例えば、トラフィックに応じて地点間を結ぶ工事の際の最適な迂回路の設計や、人から人への情報の伝達、企業間の資金の流れなどのコネクションに応じたコミュニティ生成に関する予測が可能となる。
【００３８】
また、請求項２に記載の情報間関係性解析装置は、請求項１に記載の情報間関係性解析装置において、上記存在確率算出手段および上記無効後存在確率算出手段は、数式１に基づいて、それぞれ、上記存在確率および上記無効後存在確率を算出することを特徴とする。
【数４】
2005275794.tif 000006
（ここで、「Ｒ（ｐ）」は概念情報識別情報「ｐ」に対応する概念情報の存在確率であり、「Ｒ（ｑ）」は概念情報識別情報「ｐ」に対応する概念情報への参照情報をもつ概念情報識別情報「ｑ」に対応する概念情報の存在確率であり、「Ｃ（ｑ）」は概念情報識別情報「ｑ」に対応する概念情報から他の概念情報への参照情報の数である。また、「Ｇ」は概念情報をノード、参照情報をエッジとしたグラフであり、「ｎ」はグラフ「Ｇ」のノードである概念情報の総数である。「ε」は予め設定される定数である。）
【００３９】
これは存在確率算出手段および無効後存在確立算出手段の一例を一層具体的に示すものである。この装置によれば、存在確率算出手段および無効後存在確率算出手段は、数式１に基づいて、それぞれ、存在確率および無効後存在確率を算出するので、既存の状態遷移モデルを用いて存在確率および無効後存在確率を算出することができる。
【数５】
2005275794.tif 000007
（ここで、「Ｒ（ｐ）」は概念情報識別情報「ｐ」に対応する概念情報の存在確率であり、「Ｒ（ｑ）」は概念情報識別情報「ｐ」に対応する概念情報への参照情報をもつ概念情報識別情報「ｑ」に対応する概念情報の存在確率であり、「Ｃ（ｑ）」は概念情報識別情報「ｑ」に対応する概念情報から他の概念情報への参照情報の数である。また、「Ｇ」は概念情報をノード、参照情報をエッジとしたグラフであり、「ｎ」はグラフ「Ｇ」のノードである概念情報の総数である。「ε」は予め設定される定数である。）
【００４０】
また、請求項３に記載の情報間関係性解析装置は、請求項１または２に記載の情報間関係性解析装置において、上記影響度算出手段は、数式２に示すように、上記存在確率と上記無効後存在確率との差に基づいて、上記影響度情報を算出することを特徴とする。
【数６】
2005275794.tif 000008
（ここで、「ｄＲ（ｐ）(s)」は評価概念情報識別情報「ｓ」に対応する評価概念情報が注目概念情報識別情報「ｐ」対応する注目概念情報に及ぼす影響度情報である。「Ｒ（ｐ）(s)」は、概念情報識別情報「ｓ」に対応する概念情報を参照情報と共に無効にしたときの、概念情報識別情報「ｐ」に対応する概念情報の存在確率である無効後存在確率であり、「Ｒ（ｐ）」は、概念情報識別情報「ｐ」に対応する概念情報の存在確率である。）
【００４１】
これは影響度算出手段の一例を一層具体的に示すものである。この装置によれば、影響度算出手段は、数式２に示すように、存在確率と無効後存在確率との差に基づいて、影響度情報を算出するので、簡単な四則演算により影響度情報を算出することができる。
【数７】
2005275794.tif 000009
（ここで、「ｄＲ（ｐ）(s)」は評価概念情報識別情報「ｓ」に対応する評価概念情報が注目概念情報識別情報「ｐ」対応する注目概念情報に及ぼす影響度情報である。「Ｒ（ｐ）(s)」は、概念情報識別情報「ｓ」に対応する概念情報を参照情報と共に無効にしたときの、概念情報識別情報「ｐ」に対応する概念情報の存在確率である無効後存在確率であり、「Ｒ（ｐ）」は、概念情報識別情報「ｐ」に対応する概念情報の存在確率である。）
【００４２】
また、請求項４に記載の情報間関係性解析装置は、請求項１〜３のいずれか一つに記載の情報間関係性解析装置において、上記特定概念情報設定手段により設定された上記特定概念情報が複数の場合、上記影響度算出手段により算出された上記影響度情報に基づいて、上記特定概念情報の上記影響度情報からの第１主成分得点を算出する主成分得点算出手段をさらに備えたことを特徴とする。
【００４３】
この装置によれば、設定された特定概念情報が複数の場合、算出された影響度情報に基づいて、特定概念情報の影響度情報からの第１主成分得点を算出するので、例えば、複数の疾患を共に対象とした関連遺伝子を自動的に抽出することができ、新しい薬剤や治療法の開発におけるコスト削減を実現することができる。
【００４４】
また、請求項５に記載の情報間関係性解析装置は、請求項１〜４のいずれか一つに記載の情報間関係性解析装置において、上記影響度算出手段により算出された上記影響度情報および／または上記主成分得点算出手段により算出された上記第１主成分得点に基づいて、上記特定概念情報および／または上記評価概念情報を序列する序列手段をさらに備えたことを特徴とする。
【００４５】
この装置によれば、算出された影響度情報および／または算出された第１主成分得点に基づいて、特定概念情報および／または評価概念情報を、例えば昇順または降順に序列するので、例えば、影響度情報および／または第１主成分得点の絶対値の大きい遺伝子は疾患との深い関連性があるなどの知見を容易に得ることができる。
【００４６】
また、本発明は情報間関係性解析方法に関するものであり、請求項６に記載の情報間関係性解析方法は、概念に対応する情報であって他の上記概念に対応する情報を参照するための参照情報を少なくとも含む概念情報を複数含む概念情報群を取得する概念情報群取得ステップと、上記概念情報群取得ステップにより取得された上記概念情報群の中から、特定の上記概念情報である特定概念情報を設定する特定概念情報設定ステップと、上記概念情報群取得ステップにより取得された上記概念情報群について、各々の上記概念情報を単位に当該概念情報を相互に参照する上記参照情報があるとき、各々の上記概念情報を状態とみて各々の上記参照情報を通して上記概念情報を遷移するモデルで捉え、当該モデルをマルコフ過程と考えたときの上記特定概念情報の存在確率を算出する存在確率算出ステップと、上記概念情報群取得ステップにより取得された上記概念情報群に含まれる上記概念情報に基づいて、評価対象として予め設定された上記概念情報である評価概念情報を当該評価概念情報に含まれる上記参照情報と共に無効にしたときの上記特定概念情報の上記存在確率である無効後存在確率を算出する無効後存在確率算出ステップと、上記存在確率算出ステップにより算出された上記存在確率および上記無効後存在確率算出ステップにより算出された上記無効後存在確率に基づいて、上記評価概念情報が上記特定概念情報設定ステップにより設定された上記特定概念情報に及ぼす影響度に関する情報である影響度情報を算出する影響度算出ステップとを含むことを特徴とする。
【００４７】
この方法によれば、概念（例えば、疾患や遺伝子など）に対応する情報であって他の概念に対応する情報を参照するための参照情報（例えば、Ｗｅｂページにおけるリンクなど）を少なくとも含む概念情報（例えば、Ｗｅｂページなど）を複数含む概念情報群を取得し、取得された概念情報群の中から、特定の概念情報である特定概念情報を設定し、取得された概念情報群について、各々の概念情報を単位に当該概念情報を相互に参照する参照情報があるとき、各々の概念情報を状態とみて各々の参照情報を通して概念情報を遷移するモデルで捉え、当該モデルをマルコフ過程と考えたときの特定概念情報の存在確率を算出し、取得された概念情報群に含まれる概念情報に基づいて、評価対象として予め設定された概念情報である評価概念情報を当該評価概念情報に含まれる参照情報と共に無効にしたときの特定概念情報の存在確率である無効後存在確率を算出し、算出された存在確率および算出された無効後存在確率に基づいて、評価概念情報が特定概念情報に及ぼす影響度に関する情報である影響度情報を算出するので、概念情報間の参照情報に基づく関係性の解析ができる。すなわち、この方法によれば、相互に参照される概念情報群において、ターゲットとする概念情報（評価概念情報）を無効にしたときの他への影響を算出するので、互いの関係性を定量化することができる。また、この方法によれば、参照情報に基づいて互いの概念情報の関係性を定量化しているので、背景にある構造を反映しており、情報間の潜在的な関係性を発見することができる。
【００４８】
また、この方法によれば、例えば疾患関連遺伝子の抽出において、遺伝子を個別にみず、全体を系とみているので、遺伝子間の相互作用が考慮された疾患関連遺伝子を局所的な範囲に限定せずに抽出することができる。具体的には、この方法によれば、例えば、疾患と遺伝子群が医学生物学的結び付きによりリンク構造をもつネットワークを状態遷移モデルで捉え、或る状態（例えば、遺伝子）を無効にしたときの、疾患に相当する状態へ与える変動から、遺伝子の疾患に対する寄与を定量化することができる。また、この方法によれば、例えば、潜在的な疾患関連遺伝子の発見の他、疾患の発症に強く関与する主効果遺伝子を抽出することができる。
【００４９】
また、この方法によれば、例えば、トラフィックに応じて地点間を結ぶ工事の際の最適な迂回路の設計や、人から人への情報の伝達、企業間の資金の流れなどのコネクションに応じたコミュニティ生成に関する予測が可能となる。
【００５０】
また、請求項７に記載の情報間関係性解析方法は、請求項６に記載の情報間関係性解析方法において、上記存在確率算出ステップおよび上記無効後存在確率算出ステップは、数式１に基づいて、それぞれ、上記存在確率および上記無効後存在確率を算出することを特徴とする。
【数８】
2005275794.tif 000010
（ここで、「Ｒ（ｐ）」は概念情報識別情報「ｐ」に対応する概念情報の存在確率であり、「Ｒ（ｑ）」は概念情報識別情報「ｐ」に対応する概念情報への参照情報をもつ概念情報識別情報「ｑ」に対応する概念情報の存在確率であり、「Ｃ（ｑ）」は概念情報識別情報「ｑ」に対応する概念情報から他の概念情報への参照情報の数である。また、「Ｇ」は概念情報をノード、参照情報をエッジとしたグラフであり、「ｎ」はグラフ「Ｇ」のノードである概念情報の総数である。「ε」は予め設定される定数である。）
【００５１】
これは存在確率算出ステップおよび無効後存在確立算出ステップの一例を一層具体的に示すものである。この方法によれば、存在確率算出ステップおよび無効後存在確率算出ステップは、数式１に基づいて、それぞれ、存在確率および無効後存在確率を算出するので、既存の状態遷移モデルを用いて存在確率および無効後存在確率を算出することができる。
【数９】
2005275794.tif 000011
（ここで、「Ｒ（ｐ）」は概念情報識別情報「ｐ」に対応する概念情報の存在確率であり、「Ｒ（ｑ）」は概念情報識別情報「ｐ」に対応する概念情報への参照情報をもつ概念情報識別情報「ｑ」に対応する概念情報の存在確率であり、「Ｃ（ｑ）」は概念情報識別情報「ｑ」に対応する概念情報から他の概念情報への参照情報の数である。また、「Ｇ」は概念情報をノード、参照情報をエッジとしたグラフであり、「ｎ」はグラフ「Ｇ」のノードである概念情報の総数である。「ε」は予め設定される定数である。）
【００５２】
また、請求項８に記載の情報間関係性解析方法は、請求項６または７に記載の情報間関係性解析方法において、上記影響度算出ステップは、数式２に示すように、上記存在確率と上記無効後存在確率との差に基づいて、上記影響度情報を算出することを特徴とする。
【数１０】
2005275794.tif 000012
（ここで、「ｄＲ（ｐ）(s)」は評価概念情報識別情報「ｓ」に対応する評価概念情報が注目概念情報識別情報「ｐ」対応する注目概念情報に及ぼす影響度情報である。「Ｒ（ｐ）(s)」は、概念情報識別情報「ｓ」に対応する概念情報を参照情報と共に無効にしたときの、概念情報識別情報「ｐ」に対応する概念情報の存在確率である無効後存在確率であり、「Ｒ（ｐ）」は、概念情報識別情報「ｐ」に対応する概念情報の存在確率である。）
【００５３】
これは影響度算出ステップの一例を一層具体的に示すものである。この方法によれば、影響度算出ステップは、数式２に示すように、存在確率と無効後存在確率との差に基づいて、影響度情報を算出するので、簡単な四則演算により影響度情報を算出することができる。
【数１１】
2005275794.tif 000013
（ここで、「ｄＲ（ｐ）(s)」は評価概念情報識別情報「ｓ」に対応する評価概念情報が注目概念情報識別情報「ｐ」対応する注目概念情報に及ぼす影響度情報である。「Ｒ（ｐ）(s)」は、概念情報識別情報「ｓ」に対応する概念情報を参照情報と共に無効にしたときの、概念情報識別情報「ｐ」に対応する概念情報の存在確率である無効後存在確率であり、「Ｒ（ｐ）」は、概念情報識別情報「ｐ」に対応する概念情報の存在確率である。）
【００５４】
また、請求項９に記載の情報間関係性解析方法は、請求項６〜８のいずれか一つに記載の情報間関係性解析方法において、上記特定概念情報設定ステップにより設定された上記特定概念情報が複数の場合、上記影響度算出ステップにより算出された上記影響度情報に基づいて、上記特定概念情報の上記影響度情報からの第１主成分得点を算出する主成分得点算出ステップをさらに含むことを特徴とする。
【００５５】
この方法によれば、設定された特定概念情報が複数の場合、算出された影響度情報に基づいて、特定概念情報の影響度情報からの第１主成分得点を算出するので、例えば、複数の疾患を共に対象とした関連遺伝子を自動的に抽出することができ、新しい薬剤や治療法の開発におけるコスト削減を実現することができる。
【００５６】
また、請求項１０に記載の情報間関係性解析方法は、請求項６〜９のいずれか一つに記載の情報間関係性解析方法において、上記影響度算出ステップにより算出された上記影響度情報および／または上記主成分得点算出ステップにより算出された上記第１主成分得点に基づいて、上記特定概念情報および／または上記評価概念情報を序列する序列ステップをさらに含むことを特徴とする。
【００５７】
この方法によれば、算出された影響度情報および／または算出された第１主成分得点に基づいて、特定概念情報および／または評価概念情報を、例えば昇順または降順に序列するので、例えば、影響度情報および／または第１主成分得点の絶対値の大きい遺伝子は疾患との深い関連性があるなどの知見を容易に得ることができる。
【００５８】
また、本発明はプログラムに関するものであり、請求項１１に記載の情報間関係性解析方法をコンピュータに実行させることを特徴とするプログラムは、概念に対応する情報であって他の上記概念に対応する情報を参照するための参照情報を少なくとも含む概念情報を複数含む概念情報群を取得する概念情報群取得ステップと、上記概念情報群取得ステップにより取得された上記概念情報群の中から、特定の上記概念情報である特定概念情報を設定する特定概念情報設定ステップと、上記概念情報群取得ステップにより取得された上記概念情報群について、各々の上記概念情報を単位に当該概念情報を相互に参照する上記参照情報があるとき、各々の上記概念情報を状態とみて各々の上記参照情報を通して上記概念情報を遷移するモデルで捉え、当該モデルをマルコフ過程と考えたときの上記特定概念情報の存在確率を算出する存在確率算出ステップと、上記概念情報群取得ステップにより取得された上記概念情報群に含まれる上記概念情報に基づいて、評価対象として予め設定された上記概念情報である評価概念情報を当該評価概念情報に含まれる上記参照情報と共に無効にしたときの上記特定概念情報の上記存在確率である無効後存在確率を算出する無効後存在確率算出ステップと、上記存在確率算出ステップにより算出された上記存在確率および上記無効後存在確率算出ステップにより算出された上記無効後存在確率に基づいて、上記評価概念情報が上記特定概念情報設定ステップにより設定された上記特定概念情報に及ぼす影響度に関する情報である影響度情報を算出する影響度算出ステップとを含むことを特徴とする。
【００５９】
このプログラムによれば、概念（例えば、疾患や遺伝子など）に対応する情報であって他の概念に対応する情報を参照するための参照情報（例えば、Ｗｅｂページにおけるリンクなど）を少なくとも含む概念情報（例えば、Ｗｅｂページなど）を複数含む概念情報群を取得し、取得された概念情報群の中から、特定の概念情報である特定概念情報を設定し、取得された概念情報群について、各々の概念情報を単位に当該概念情報を相互に参照する参照情報があるとき、各々の概念情報を状態とみて各々の参照情報を通して概念情報を遷移するモデルで捉え、当該モデルをマルコフ過程と考えたときの特定概念情報の存在確率を算出し、取得された概念情報群に含まれる概念情報に基づいて、評価対象として予め設定された概念情報である評価概念情報を当該評価概念情報に含まれる参照情報と共に無効にしたときの特定概念情報の存在確率である無効後存在確率を算出し、算出された存在確率および算出された無効後存在確率に基づいて、評価概念情報が特定概念情報に及ぼす影響度に関する情報である影響度情報を算出するので、概念情報間の参照情報に基づく関係性の解析ができる。すなわち、このプログラムによれば、相互に参照される概念情報群において、ターゲットとする概念情報（評価概念情報）を無効にしたときの他への影響を算出するので、互いの関係性を定量化することができる。また、このプログラムによれば、参照情報に基づいて互いの概念情報の関係性を定量化しているので、背景にある構造を反映しており、情報間の潜在的な関係性を発見することができる。
【００６０】
また、このプログラムによれば、例えば疾患関連遺伝子の抽出において、遺伝子を個別にみず、全体を系とみているので、遺伝子間の相互作用が考慮された疾患関連遺伝子を局所的な範囲に限定せずに抽出することができる。具体的には、このプログラムによれば、例えば、疾患と遺伝子群が医学生物学的結び付きによりリンク構造をもつネットワークを状態遷移モデルで捉え、或る状態（例えば、遺伝子）を無効にしたときの、疾患に相当する状態へ与える変動から、遺伝子の疾患に対する寄与を定量化することができる。また、このプログラムによれば、例えば、潜在的な疾患関連遺伝子の発見の他、疾患の発症に強く関与する主効果遺伝子を抽出することができる。
【００６１】
また、このプログラムによれば、例えば、トラフィックに応じて地点間を結ぶ工事の際の最適な迂回路の設計や、人から人への情報の伝達、企業間の資金の流れなどのコネクションに応じたコミュニティ生成に関する予測が可能となる。
【００６２】
また、請求項１２に記載のプログラムは、請求項１１に記載のプログラムにおいて、上記存在確率算出ステップおよび上記無効後存在確率算出ステップは、数式１に基づいて、それぞれ、上記存在確率および上記無効後存在確率を算出することを特徴とする。
【数１２】
2005275794.tif 000014
（ここで、「Ｒ（ｐ）」は概念情報識別情報「ｐ」に対応する概念情報の存在確率であり、「Ｒ（ｑ）」は概念情報識別情報「ｐ」に対応する概念情報への参照情報をもつ概念情報識別情報「ｑ」に対応する概念情報の存在確率であり、「Ｃ（ｑ）」は概念情報識別情報「ｑ」に対応する概念情報から他の概念情報への参照情報の数である。また、「Ｇ」は概念情報をノード、参照情報をエッジとしたグラフであり、「ｎ」はグラフ「Ｇ」のノードである概念情報の総数である。「ε」は予め設定される定数である。）
【００６３】
これは存在確率算出ステップおよび無効後存在確立算出ステップの一例を一層具体的に示すものである。このプログラムによれば、存在確率算出ステップおよび無効後存在確率算出ステップは、数式１に基づいて、それぞれ、存在確率および無効後存在確率を算出するので、既存の状態遷移モデルを用いて存在確率および無効後存在確率を算出することができる。
【数１３】
2005275794.tif 000015
（ここで、「Ｒ（ｐ）」は概念情報識別情報「ｐ」に対応する概念情報の存在確率であり、「Ｒ（ｑ）」は概念情報識別情報「ｐ」に対応する概念情報への参照情報をもつ概念情報識別情報「ｑ」に対応する概念情報の存在確率であり、「Ｃ（ｑ）」は概念情報識別情報「ｑ」に対応する概念情報から他の概念情報への参照情報の数である。また、「Ｇ」は概念情報をノード、参照情報をエッジとしたグラフであり、「ｎ」はグラフ「Ｇ」のノードである概念情報の総数である。「ε」は予め設定される定数である。）
【００６４】
また、請求項１３に記載のプログラムは、請求項１１または１２に記載のプログラムにおいて、上記影響度算出ステップは、数式２に示すように、上記存在確率と上記無効後存在確率との差に基づいて、上記影響度情報を算出することを特徴とする。
【数１４】
2005275794.tif 000016
（ここで、「ｄＲ（ｐ）(s)」は評価概念情報識別情報「ｓ」に対応する評価概念情報が注目概念情報識別情報「ｐ」対応する注目概念情報に及ぼす影響度情報である。「Ｒ（ｐ）(s)」は、概念情報識別情報「ｓ」に対応する概念情報を参照情報と共に無効にしたときの、概念情報識別情報「ｐ」に対応する概念情報の存在確率である無効後存在確率であり、「Ｒ（ｐ）」は、概念情報識別情報「ｐ」に対応する概念情報の存在確率である。）
【００６５】
これは影響度算出ステップの一例を一層具体的に示すものである。このプログラムによれば、影響度算出ステップは、数式２に示すように、存在確率と無効後存在確率との差に基づいて、影響度情報を算出するので、簡単な四則演算により影響度情報を算出することができる。
【数１５】
2005275794.tif 000017
（ここで、「ｄＲ（ｐ）(s)」は評価概念情報識別情報「ｓ」に対応する評価概念情報が注目概念情報識別情報「ｐ」対応する注目概念情報に及ぼす影響度情報である。「Ｒ（ｐ）(s)」は、概念情報識別情報「ｓ」に対応する概念情報を参照情報と共に無効にしたときの、概念情報識別情報「ｐ」に対応する概念情報の存在確率である無効後存在確率であり、「Ｒ（ｐ）」は、概念情報識別情報「ｐ」に対応する概念情報の存在確率である。）
【００６６】
また、請求項１４に記載のプログラムは、請求項１１〜１３のいずれか一つに記載のプログラムにおいて、上記特定概念情報設定ステップにより設定された上記特定概念情報が複数の場合、上記影響度算出ステップにより算出された上記影響度情報に基づいて、上記特定概念情報の上記影響度情報からの第１主成分得点を算出する主成分得点算出ステップをさらに含むことを特徴とする。
【００６７】
このプログラムによれば、設定された特定概念情報が複数の場合、算出された影響度情報に基づいて、特定概念情報の影響度情報からの第１主成分得点を算出するので、例えば、複数の疾患を共に対象とした関連遺伝子を自動的に抽出することができ、新しい薬剤や治療法の開発におけるコスト削減を実現することができる。
【００６８】
また、請求項１５に記載のプログラムは、請求項１１〜１４のいずれか一つに記載のプログラムにおいて、上記影響度算出ステップにより算出された上記影響度情報および／または上記主成分得点算出ステップにより算出された上記第１主成分得点に基づいて、上記特定概念情報および／または上記評価概念情報を序列する序列ステップをさらに含むことを特徴とする。
【００６９】
このプログラムによれば、算出された影響度情報および／または算出された第１主成分得点に基づいて、特定概念情報および／または評価概念情報を、例えば昇順または降順に序列するので、例えば、影響度情報および／または第１主成分得点の絶対値の大きい遺伝子は疾患との深い関連性があるなどの知見を容易に得ることができる。
【００７０】
また、本発明は記録媒体に関するものであり、請求項１６に記載の記録媒体は、上記請求項１１〜１５のいずれか一つに記載されたプログラムを記録したことを特徴とする。
【００７１】
この記録媒体によれば、当該記録媒体に記録されたプログラムをコンピュータに読み取らせて実行することによって、請求項１１〜１５のいずれか一つに記載されたプログラムを、コンピュータを利用して実現することができ、これら各プログラムと同様の効果を得ることができる。
【発明の効果】
【００７２】
本発明によれば、概念（例えば、疾患や遺伝子など）に対応する情報であって他の概念に対応する情報を参照するための参照情報（例えば、Ｗｅｂページにおけるリンクなど）を少なくとも含む概念情報（例えば、Ｗｅｂページなど）を複数含む概念情報群を取得し、取得された概念情報群の中から、特定の概念情報である特定概念情報を設定し、取得された概念情報群について、各々の概念情報を単位に当該概念情報を相互に参照する参照情報があるとき、各々の概念情報を状態とみて各々の参照情報を通して概念情報を遷移するモデルで捉え、当該モデルをマルコフ過程と考えたときの特定概念情報の存在確率を算出し、取得された概念情報群に含まれる概念情報に基づいて、評価対象として予め設定された概念情報である評価概念情報を当該評価概念情報に含まれる参照情報と共に無効にしたときの特定概念情報の存在確率である無効後存在確率を算出し、算出された存在確率および算出された無効後存在確率に基づいて、評価概念情報が特定概念情報に及ぼす影響度に関する情報である影響度情報を算出するので、概念情報間の参照情報に基づく関係性の解析ができる、情報間関係性解析装置、情報間関係性解析方法、プログラムおよび記録媒体を提供することができる。
【００７３】
すなわち、本発明によれば、相互に参照される概念情報群において、ターゲットとする概念情報（評価概念情報）を無効にしたときの他への影響を算出するので、互いの関係性を定量化することができる、情報間関係性解析装置、情報間関係性解析方法、プログラムおよび記録媒体を提供することができる。また、本発明によれば、参照情報に基づいて互いの概念情報の関係性を定量化しているので、背景にある構造を反映しており、情報間の潜在的な関係性を発見することができる、情報間関係性解析装置、情報間関係性解析方法、プログラムおよび記録媒体を提供することができる。
【００７４】
また、本発明によれば、例えば疾患関連遺伝子の抽出において、遺伝子を個別にみず、全体を系とみているので、遺伝子間の相互作用が考慮された疾患関連遺伝子を局所的な範囲に限定せずに抽出することができる、情報間関係性解析装置、情報間関係性解析方法、プログラムおよび記録媒体を提供することができる。具体的には、本発明によれば、例えば、疾患と遺伝子群が医学生物学的結び付きによりリンク構造をもつネットワークを状態遷移モデルで捉え、或る状態（例えば、遺伝子）を無効にしたときの、疾患に相当する状態へ与える変動から、遺伝子の疾患に対する寄与を定量化することができる、情報間関係性解析装置、情報間関係性解析方法、プログラムおよび記録媒体を提供することができる。また、本発明によれば、例えば、潜在的な疾患関連遺伝子の発見の他、疾患の発症に強く関与する主効果遺伝子を抽出することができる、情報間関係性解析装置、情報間関係性解析方法、プログラムおよび記録媒体を提供することができる。
【００７５】
また、本発明によれば、例えば、トラフィックに応じて地点間を結ぶ工事の際の最適な迂回路の設計や、人から人への情報の伝達、企業間の資金の流れなどのコネクションに応じたコミュニティ生成に関する予測が可能となる、情報間関係性解析装置、情報間関係性解析方法、プログラムおよび記録媒体を提供することができる。
【００７６】
また、本発明によれば、既存の状態遷移モデルを用いて存在確率および無効後存在確率を算出することができる、情報間関係性解析装置、情報間関係性解析方法、プログラムおよび記録媒体を提供することができる。
【００７７】
また、本発明によれば、簡単な四則演算により影響度情報を算出することができる、情報間関係性解析装置、情報間関係性解析方法、プログラムおよび記録媒体を提供することができる。
【００７８】
また、本発明によれば、例えば、複数の疾患を共に対象とした関連遺伝子を自動的に抽出することができ、新しい薬剤や治療法の開発におけるコスト削減を実現することができる、情報間関係性解析装置、情報間関係性解析方法、プログラムおよび記録媒体を提供することができる。
【００７９】
さらに、本発明によれば、例えば、影響度情報および／または第１主成分得点の絶対値の大きい遺伝子は疾患との深い関連性があるなどの知見を容易に得ることができる、情報間関係性解析装置、情報間関係性解析方法、プログラムおよび記録媒体を提供することができる。
【発明を実施するための最良の形態】
【００８０】
以下に、本発明にかかる情報間関係性解析装置、情報間関係性解析方法、プログラムおよび記録媒体の実施の形態を図面に基づいて詳細に説明する。なお、この実施の形態によりこの発明が限定されるものではない。
【００８１】
そして、後述する実施例では、知識発見の顕著なニーズがあるライフサイエンス分野を対象にし、疾患関連遺伝子の探索へ応用した実験結果を述べ、本発明の有効性を示す。
【００８２】
［本発明の概要］
以下、本発明の概要について説明し、その後、本発明の構成および処理等について詳細に説明する。
図１は、本発明の基本原理を示す原理構成図である。
【００８３】
本発明は、概略的に、以下の基本的特徴を有する。すなわち、本発明は、まず、概念（例えば、疾患や遺伝子など）に対応する情報であって他の概念に対応する情報を参照するための参照情報（例えば、Ｗｅｂページへのリンクなど）を少なくとも含む概念情報（例えば、Ｗｅｂページなど）を複数含む概念情報群を取得し、概念情報群ファイルの所定の記憶領域に格納する（ステップＳ−１）。具体的には、例えば、ある疾患または遺伝子に対応し、他の疾患または遺伝子に対応するページに向けてリンクを張るＷｅｂページを複数含むページ群を取得する。つまり、相互に参照しあう構造を有するページ群を取得する。
【００８４】
ついで、ステップＳ−１にて取得された概念情報群の中から、特定の概念情報である特定概念情報を設定し、設定された特定概念情報を特定概念情報関係情報ファイルの所定の記憶領域に格納する（ステップＳ−２）。具体的には、例えば、ユーザが注目する特定の疾患に対応するＷｅｂページを設定する。
【００８５】
ついで、ステップＳ−１にて取得された概念情報群について、各々の概念情報を単位に当該概念情報を相互に参照する参照情報があるとき、各々の概念情報を状態とみて各々の参照情報を通して概念情報を遷移するモデルで捉え、当該モデルをマルコフ過程と考えたときの特定概念情報の存在確率を算出し、特定概念情報関係情報ファイルの所定の記憶領域に格納する（ステップＳ−３）。
【００８６】
ついで、ステップＳ−１にて取得された概念情報群に含まれる概念情報に基づいて、評価対象として予め設定された概念情報である評価概念情報を当該評価概念情報に含まれる参照情報と共に無効にしたときの特定概念情報の存在確率である無効後存在確率を算出し、特定概念情報関係情報ファイルの所定の記憶領域に格納する（ステップＳ−４）。具体的には、例えば、ある遺伝子に対応するＷｅｂページをリンクと共に無効にした際の、先に設定した特定の疾患に対応するＷｅｂページの存在確率を算出する。
【００８７】
ここで、ステップＳ−３およびステップＳ−４において、以下の数式１に基づいて、それぞれ、存在確率および無効後存在確率を算出してもよい。
【数１６】
2005275794.tif 000018
（ここで、「Ｒ（ｐ）」は概念情報識別情報「ｐ」に対応する概念情報の存在確率であり、「Ｒ（ｑ）」は概念情報識別情報「ｐ」に対応する概念情報への参照情報をもつ概念情報識別情報「ｑ」に対応する概念情報の存在確率であり、「Ｃ（ｑ）」は概念情報識別情報「ｑ」に対応する概念情報から他の概念情報への参照情報の数である。また、「Ｇ」は概念情報をノード、参照情報をエッジとしたグラフであり、「ｎ」はグラフ「Ｇ」のノードである概念情報の総数である。「ε」は予め設定される定数である。）
【００８８】
ついで、ステップＳ−３にて算出された存在確率およびステップＳ−４にて算出された無効後存在確率に基づいて、評価概念情報がステップＳ−２にて設定された特定概念情報に及ぼす影響度に関する情報である影響度情報を算出し、特定概念情報関係情報ファイルの所定の記憶領域に格納する（ステップＳ−５）。具体的には、例えば、ある遺伝子に対応するＷｅｂページをリンクと共に無効にした（ノックアウトした）ことによる先に設定した疾患に対応するＷｅｂページへの影響の度合いを算出する。
【００８９】
ここで、ステップＳ−５において、以下の数式２に示すように、存在確率と無効後存在確率との差に基づいて、影響度情報を算出してもよい。
【数１７】
2005275794.tif 000019
（ここで、「ｄＲ（ｐ）(s)」は評価概念情報識別情報「ｓ」に対応する評価概念情報が注目概念情報識別情報「ｐ」対応する注目概念情報に及ぼす影響度情報である。「Ｒ（ｐ）(s)」は、概念情報識別情報「ｓ」に対応する概念情報を参照情報と共に無効にしたときの、概念情報識別情報「ｐ」に対応する概念情報の存在確率である無効後存在確率であり、「Ｒ（ｐ）」は、概念情報識別情報「ｐ」に対応する概念情報の存在確率である。）
【００９０】
また、ステップＳ−２にて設定された特定概念情報が複数の場合、ステップＳ−５にて算出された影響度情報に基づいて、特定概念情報の影響度情報からの第１主成分得点などの主成分得点を算出してもよい。
【００９１】
また、ステップＳ−５にて算出された影響度情報および／または算出された第１主成分得点などの主成分得点に基づいて、特定概念情報および／または評価概念情報を例えば昇順または降順に序列してもよい。
【００９２】
［システム構成］
次に、本システムの構成について説明する。
図２は、本発明が適用される本システムの構成の一例を示すブロック図であり、該構成のうち本発明に関係する部分のみを概念的に示している。
【００９３】
図２において、ネットワーク３００は、情報間関係性解析装置１００と外部システム２００とを相互に接続する機能を有し、例えば、インターネット等である。
【００９４】
図２において、外部システム２００は、ネットワーク３００を介して、情報間関係性解析装置１００と相互に接続され、利用者に対して文書情報、論文情報等に関する外部データベースや各種の外部プログラムを実行するウェブサイトを提供する機能を有する。
【００９５】
ここで、外部システム２００は、ＷＥＢサーバやＡＳＰサーバ等として構成してもよく、そのハードウェア構成は、一般に市販されるワークステーション、パーソナルコンピュータ等の情報処理装置およびその付属装置により構成してもよい。また、外部システム２００の各機能は、外部システム２００のハードウェア構成中のＣＰＵ、ディスク装置、メモリ装置、入力装置、出力装置、通信制御装置等およびそれらを制御するプログラム等により実現される。
【００９６】
図２において、情報間関係性解析装置１００は、概略的に、情報間関係性解析装置１００の全体を統括的に制御するＣＰＵ等の制御部１０２、通信回線等に接続されるルータ等の通信装置（図示せず）に接続される通信制御インターフェース部１０４、各種のデータベースやファイルなどを格納する記憶部１０６、および、入力装置１１２や出力装置１１４に接続される入出力制御インターフェース部１０８を備えて構成されており、これら各部は任意の通信路を介して通信可能に接続されている。さらに、この情報間関係性解析装置１００は、ルータ等の通信装置および専用線等の有線または無線の通信回線を介して、ネットワーク３００に通信可能に接続されている。
【００９７】
図２の記憶部１０６に格納される各種のデータベースやテーブルやファイル（概念情報群ファイル１０６ａ〜序列結果ファイル１０６ｄ）は、固定ディスク装置等のストレージ手段であり、各種処理に用いる各種のプログラムやテーブルやファイルやデータベースやウェブページ用ファイル等を格納する。
【００９８】
これら記憶部１０６の各構成要素のうち、概念情報群ファイル１０６ａは、後述する概念情報群取得部１０２ａにより取得された、概念に対応する情報であって他の概念に対応する情報を参照するための参照情報を少なくとも含む概念情報を複数含む概念情報群を格納する概念情報群格納手段である。この概念情報群ファイル１０６ａに格納される情報は、複数の概念情報で構成されている。
【００９９】
また、特定概念情報関係情報ファイル１０６ｂは、後述する特定概念情報設定部１０２ｂにより設定された特定概念情報において、後述する存在確率算出部１０２ｃにより算出された存在確率と、後述する無効後存在確率算出部１０２ｄにより算出された無効後存在確率と、後述する無効後存在確率算出部１０２ｄにより無効後存在確率を算出する際に無効にした評価概念情報と、後述する影響度算出部１０２ｅにより算出された影響度情報を格納する特定概念情報関係情報格納手段である。図３は、特定概念情報関係情報ファイル１０６ｂに格納される情報の一例を示す図である。この特定概念情報関係情報ファイル１０６ｂに格納される情報は、図３に示すように、特定概念情報を一意に識別するための特定概念情報識別情報と、特定概念情報の存在確率と、無効にした評価概念情報を一意に識別するための評価概念情報識別情報と、各評価概念情報に対応する無効後存在確率および影響度と、を相互に関連付けて構成されている。
【０１００】
また、主成分得点ファイル１０６ｃは、設定された特定概念情報が複数の場合、後述する主成分得点算出部１０２ｆにより、影響度情報に基づいて算出した、特定概念情報の影響度情報からの第１主成分得点などの主成分得点を格納する主成分得点格納手段である。この主成分得点ファイル１０６ｃに格納される情報は、特定概念情報の影響度情報からの第１主成分得点などの主成分得点で構成されている。
【０１０１】
また、序列結果ファイル１０６ｄは、後述する序列部１０２ｇにより、特定概念情報および／または評価概念情報を、影響度情報および／または第１主成分得点などの主成分得点に基づいて序列した結果である序列結果を格納する序列結果格納手段である。この序列結果ファイル１０６ｄに格納される情報は、特定概念情報および／または評価概念情報を、影響度情報および／または第１主成分得点などの主成分得点に基づいて序列した結果である序列結果で構成されている。
【０１０２】
また、図２において、通信制御インターフェース部１０４は、情報間関係性解析装置１００とネットワーク３００（またはルータ等の通信装置）との間における通信制御を行う。すなわち、通信制御インターフェース部１０４は、他の端末と通信回線を介してデータを通信する機能を有する。
【０１０３】
また、図２において、入出力制御インターフェース部１０８は、入力装置１１２や出力装置１１４の制御を行う。ここで、出力装置１１４としては、モニタ（家庭用テレビを含む）の他、スピーカ等を用いることができる（なお、以下においては出力装置１１４をモニタとして記載する場合がある）。また、入力装置１１２としては、キーボード、マウス、および、マイク等を用いることができる。また、モニタも、マウスと協働してポインティングデバイス機能を実現する。
【０１０４】
また、図２において、制御部１０２は、ＯＳ（Ｏｐｅｒａｔｉｎｇ  Ｓｙｓｔｅｍ）等の制御プログラム、各種の処理手順等を規定したプログラム、および所要データを格納するための内部メモリを有し、これらのプログラム等により、種々の処理を実行するための情報処理を行う。制御部１０２は、機能概念的に、概念情報群取得部１０２ａ、特定概念情報設定部１０２ｂ、存在確率算出部１０２ｃ、無効後存在確率算出部１０２ｄ、影響度算出部１０２ｅ、主成分得点算出部１０２ｆおよび序列部１０２ｇを含んで構成されている。
【０１０５】
このうち、概念情報群取得部１０２ａは、概念に対応する情報であって他の概念に対応する情報を参照するための参照情報を少なくとも含む概念情報を複数含む概念情報群を取得する概念情報群取得手段である。
【０１０６】
また、特定概念情報設定部１０２ｂは、概念情報群の中から、特定の概念情報である特定概念情報を設定する特定概念情報設定手段である。
【０１０７】
また、存在確率算出部１０２ｃは、概念情報群について、各々の概念情報を単位に当該概念情報を相互に参照する参照情報があるとき、各々の概念情報を状態とみて各々の参照情報を通して概念情報を遷移するモデルで捉え、当該モデルをマルコフ過程と考えたときの特定概念情報の存在確率を算出する存在確率算出手段である。
【０１０８】
また、無効後存在確率算出部１０２ｄは、概念情報群に含まれる概念情報に基づいて、評価対象として予め設定された概念情報である評価概念情報を当該評価概念情報に含まれる参照情報と共に無効にしたときの特定概念情報の存在確率である無効後存在確率を算出する無効後存在確率算出手段である。
【０１０９】
また、影響度算出部１０２ｅは、存在確率および無効後存在確率に基づいて、評価概念情報が特定概念情報に及ぼす影響度に関する情報である影響度情報を算出する影響度算出手段である。
【０１１０】
また、主成分得点算出部１０２ｆは、特定概念情報が複数の場合、影響度情報に基づいて、特定概念情報の影響度情報からの第１主成分得点などの主成分得点を算出する主成分得点算出手段である。
【０１１１】
また、序列部１０２ｇは、影響度情報および／または第１主成分得点などの主成分得点に基づいて、特定概念情報および／または評価概念情報を序列する序列手段である。
【０１１２】
なお、これら各部によって行われる処理の詳細については、後述する。
【０１１３】
［システムの処理］
次に、このように構成された本実施の形態における本システムの処理の一例について、以下に図４を参照して詳細に説明する。図４は、本実施形態における本システムのメイン処理の一例を示すフローチャートである。
【０１１４】
まず、情報間関係性解析装置１００は、概念情報群取得部１０２ａの処理により、概念（例えば、疾患や遺伝子など）に対応する情報であって他の概念に対応する情報を参照するための参照情報（例えば、Ｗｅｂページへのリンクなど）を少なくとも含む概念情報（例えば、Ｗｅｂページなど）を複数含む概念情報群を取得し、概念情報群ファイル１０６ａの所定の記憶領域に格納する（ステップＳＡ−１）。具体的には、例えば、ある疾患または遺伝子に対応し、他の疾患または遺伝子に対応するページに向けてリンクを張るＷｅｂページを複数含むページ群を取得する。つまり、相互に参照しあう構造を有するページ群を取得する。
【０１１５】
ついで、情報間関係性解析装置１００は、特定概念情報設定部１０２ｂの処理により、ステップＳＡ−１にて取得された概念情報群の中から、特定の概念情報である特定概念情報を設定し（具体的には、例えば、ユーザが注目する特定の疾患に対応するＷｅｂページを設定し）、設定された特定概念情報を特定概念情報関係情報ファイル１０６ｂの所定の記憶領域に格納する（ステップＳＡ−２）。具体的には、設定された特定概念情報を一意に識別する特定概念情報識別情報を特定概念情報関係情報ファイル１０６ｂの所定の記憶領域に格納する。
【０１１６】
ついで、情報間関係性解析装置１００は、存在確率算出部１０２ｃの処理により、ステップＳＡ−１にて取得された概念情報群について、各々の概念情報を単位に当該概念情報を相互に参照する参照情報があるとき、各々の概念情報を状態とみて各々の参照情報を通して概念情報を遷移するモデルで捉え、当該モデルをマルコフ過程と考えたときの特定概念情報の存在確率を算出し、特定概念情報関係情報ファイル１０６ｂの所定の記憶領域に格納する（ステップＳＡ−３）。
【０１１７】
ついで、情報間関係性解析装置１００は、無効後存在確率算出部１０２ｄの処理により、ステップＳＡ−１にて取得された概念情報群に含まれる概念情報に基づいて、評価対象として予め設定された概念情報である評価概念情報を当該評価概念情報に含まれる参照情報と共に無効にしたときの特定概念情報の存在確率である無効後存在確率を算出し、特定概念情報関係情報ファイル１０６ｂの所定の記憶領域に格納する（ステップＳＡ−４）。具体的には、例えば、ある遺伝子に対応するＷｅｂページをリンクと共に無効にした際の、先に設定した特定の疾患に対応するＷｅｂページの存在確率を算出する。
【０１１８】
ここで、ステップＳＡ−３およびステップＳＡ−４において、以下の数式１に基づいて、それぞれ、存在確率および無効後存在確率を算出してもよい。
【数１８】
2005275794.tif 000020
（ここで、「Ｒ（ｐ）」は概念情報識別情報「ｐ」に対応する概念情報の存在確率であり、「Ｒ（ｑ）」は概念情報識別情報「ｐ」に対応する概念情報への参照情報をもつ概念情報識別情報「ｑ」に対応する概念情報の存在確率であり、「Ｃ（ｑ）」は概念情報識別情報「ｑ」に対応する概念情報から他の概念情報への参照情報の数である。また、「Ｇ」は概念情報をノード、参照情報をエッジとしたグラフであり、「ｎ」はグラフ「Ｇ」のノードである概念情報の総数である。「ε」は予め設定される定数である。）
【０１１９】
ついで、情報間関係性解析装置１００は、影響度算出部１０２ｅの処理により、ステップＳＡ−３にて算出された存在確率およびステップＳＡ−４にて算出された無効後存在確率に基づいて、評価概念情報がステップＳＡ−２にて設定された特定概念情報に及ぼす影響度に関する情報である影響度情報を算出し、特定概念情報関係情報ファイル１０６ｂの所定の記憶領域に格納する（ステップＳＡ−５）。具体的には、例えば、ある遺伝子に対応するＷｅｂページをリンクと共に無効にした（ノックアウトした）ことによる先に設定した疾患に対応するＷｅｂページへの影響の度合いを算出する。
【０１２０】
ここで、ステップＳＡ−５において、以下の数式２に示すように、存在確率と無効後存在確率との差に基づいて、影響度情報を算出してもよい。
【数１９】
2005275794.tif 000021
（ここで、「ｄＲ（ｐ）(s)」は評価概念情報識別情報「ｓ」に対応する評価概念情報が注目概念情報識別情報「ｐ」対応する注目概念情報に及ぼす影響度情報である。「Ｒ（ｐ）(s)」は、概念情報識別情報「ｓ」に対応する概念情報を参照情報と共に無効にしたときの、概念情報識別情報「ｐ」に対応する概念情報の存在確率である無効後存在確率であり、「Ｒ（ｐ）」は、概念情報識別情報「ｐ」に対応する概念情報の存在確率である。）
【０１２１】
また、情報間関係性解析装置１００は、ステップＳＡ−２にて設定された特定概念情報が複数の場合、主成分得点算出部１０２ｆの処理により、ステップＳＡ−５にて算出された影響度情報に基づいて、特定概念情報の影響度情報からの第１主成分得点などの主成分得点を算出し、主成分得点ファイル１０６ｃの所定の記憶領域に格納してもよい。
【０１２２】
また、情報間関係性解析装置１００は、序列部１０２ｇの処理により、ステップＳＡ−５にて算出された影響度情報および／または算出された第１主成分得点などの主成分得点に基づいて、特定概念情報および／または評価概念情報を例えば昇順または降順に序列し、序列結果を序列結果ファイル１０６ｄの所定の記憶領域に格納してもよい。
【０１２３】
以上、メイン処理が終了する。
【０１２４】
（他の実施の形態）
さて、これまで本発明の実施の形態について説明したが、本発明は、上述した実施の形態以外にも、上記特許請求の範囲に記載した技術的思想の範囲内において種々の異なる実施の形態にて実施されてよいものである。
【０１２５】
例えば、情報間関係性解析装置１００は、情報間関係性解析装置１００とは別筐体で構成されるクライアント端末からの要求に応じて処理を行い、その処理結果を当該クライアント端末に返却するように構成してもよい。
【０１２６】
また、実施形態において説明した各処理のうち、自動的に行なわれるものとして説明した処理の全部または一部を手動的に行うこともでき、あるいは、手動的に行なわれるものとして説明した処理の全部または一部を公知の方法で自動的に行うこともできる。
【０１２７】
この他、上記文書中や図面中で示した処理手順、制御手順、具体的名称、各種の登録データや検索条件等のパラメータを含む情報、画面例、データベース構成については、特記する場合を除いて任意に変更することができる。
【０１２８】
また、情報間関係性解析装置１００に関して、図示の各構成要素は機能概念的なものであり、必ずしも物理的に図示の如く構成されていることを要しない。
【０１２９】
例えば、情報間関係性解析装置１００の各部または各装置が備える処理機能、特に制御部１０２にて行なわれる各処理機能については、その全部または任意の一部を、ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）および当該ＣＰＵにて解釈実行されるプログラムにて実現することができ、あるいは、ワイヤードロジックによるハードウェアとして実現することも可能である。なお、プログラムは、後述する記録媒体に記録されており、必要に応じて情報間関係性解析装置１００に機械的に読み取られる。
【０１３０】
すなわち、ＲＯＭまたはＨＤなどの記憶部１０６などには、ＯＳ（Ｏｐｅｒａｔｉｎｇ  Ｓｙｓｔｅｍ）と協働してＣＰＵに命令を与え、各種処理を行うためのコンピュータプログラムが記録されている。このコンピュータプログラムは、ＲＡＭ等にロードされることによって実行され、ＣＰＵと協働して制御部１０２を構成する。また、このコンピュータプログラムは、情報間関係性解析装置１００に対して任意のネットワーク３００を介して接続されたアプリケーションプログラムサーバに記録されてもよく、必要に応じてその全部または一部をダウンロードすることも可能である。
【０１３１】
また、本発明にかかるプログラムを、コンピュータ読み取り可能な記録媒体に格納することもできる。ここで、この「記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ、ＥＰＲＯＭ、ＥＥＰＲＯＭ、ＣＤ−ＲＯＭ、ＭＯ、ＤＶＤ等の任意の「可搬用の物理媒体」や、各種コンピュータシステムに内蔵されるＲＯＭ、ＲＡＭ、ＨＤ等の任意の「固定用の物理媒体」、あるいは、ＬＡＮ、ＷＡＮ、インターネットに代表されるネットワークを介してプログラムを送信する場合の通信回線や搬送波のように、短期にプログラムを保持する「通信媒体」を含むものとする。
【０１３２】
また、「プログラム」とは、任意の言語や記述方法にて記述されたデータ処理方法であり、ソースコードやバイナリコード等の形式を問わない。なお、「プログラム」は必ずしも単一的に構成されるものに限られず、複数のモジュールやライブラリとして分散構成されるものや、ＯＳ（Ｏｐｅｒａｔｉｎｇ  Ｓｙｓｔｅｍ）に代表される別個のプログラムと協働してその機能を達成するものをも含む。なお、実施の形態に示した各装置において記録媒体を読み取るための具体的な構成、読み取り手順、あるいは、読み取り後のインストール手順等については、周知の構成や手順を用いることができる。
【０１３３】
記憶部１０６に格納される各種のファイル等（概念情報群ファイル１０６ａ〜序列結果ファイル１０６ｄ）は、ＲＡＭ、ＲＯＭ等のメモリ装置、ハードディスク等の固定ディスク装置、フレキシブルディスク、光ディスク等のストレージ手段であり、各種処理やウェブサイト提供に用いる各種のプログラムやテーブルやファイルやデータベースやウェブページ用ファイル等を格納する。
【０１３４】
また、情報間関係性解析装置１００は、既知のパーソナルコンピュータ、ワークステーション等の情報処理端末等の情報処理装置にプリンタやモニタやイメージスキャナ等の周辺装置を接続し、該情報処理装置に本発明の方法を実現させるソフトウェア（プログラム、データ等を含む）を実装することにより実現してもよい。
【０１３５】
さらに、情報間関係性解析装置１００の分散・統合の具体的形態は図示のものに限られず、その全部または一部を、各種の負荷等に応じた任意の単位で、機能的または物理的に分散・統合して構成することができる。例えば、各データベースを独立したデータベース装置として独立に構成してもよく、また、処理の一部をＣＧＩ（Ｃｏｍｍｏｎ  Ｇａｔｅｗａｙ  Ｉｎｔｅｒｆａｃｅ）を用いて実現してもよい。
【０１３６】
また、ネットワーク３００は、情報間関係性解析装置１００と外部システム２００とを相互に接続する機能を有し、例えば、インターネットや、イントラネットや、ＬＡＮ（有線／無線の双方を含む）や、ＶＡＮや、パソコン通信網や、公衆電話網（アナログ／デジタルの双方を含む）や、専用回線網（アナログ／デジタルの双方を含む）や、ＣＡＴＶ網や、ＩＭＴ２０００方式、ＧＳＭ方式またはＰＤＣ／ＰＤＣ−Ｐ方式等の携帯回線交換網／携帯パケット交換網や、無線呼出網や、Ｂｌｕｅｔｏｏｔｈ等の局所無線網や、ＰＨＳ網や、ＣＳ、ＢＳまたはＩＳＤＢ等の衛星通信網等のうちいずれかを含んでもよい。すなわち、本システムは、有線・無線を問わず任意のネットワークを介して、各種データを送受信することができる。
【実施例】
【０１３７】
次に、本発明の実施例について、図５から図１４等を参照して詳細に説明する。
【０１３８】
本実施例においては、ヒトの遺伝病に関するデータベースとしてよく知られるＯＭＩＭ（Ｏｎｌｉｎｅ  Ｍｅｎｄｅｌｉａｎ  Ｉｎｈｅｒｉｔａｎｃｅ  ｉｎ  Ｍａｎ）（例えば、「Ａ．Ｈａｍｏｓｈ，  ＡＦ．Ｓｃｏｔｔ，  Ｊ．Ａｍｂｅｒｇｅｒ，  Ｄ．Ｖａｌｌｅ，  ＶＡ．ＭｃＫｕｓｉｃｋ，  “Ｏｎｌｉｎｅ  Ｍｅｎｄｅｌｉａｎ  Ｉｎｈｅｒｉｔａｎｃｅ  ｉｎ  Ｍａｎ  （ＯＭＩＭ）”，  Ｈｕｍ．  Ｍｕｔａｔ．，  ｖｏｌ．１５，  Ｎｏ．１，  ｐｐ５７−６１，  ２０００」参照。）を用い、高血圧（例えば、「三木  哲郎，  “高血圧遺伝子解析プロジェクト”，  遺伝子医学，  Ｖｏｌ．５，  Ｎｏ．４，  ｐｐ．１２６−１３０，  ２００１」参照。）をはじめ、複数の遺伝子が複雑に関与する生活習慣病（例えば、「堀内  正嗣，  福田  恵一，  森下  竜一，  “生活習慣と遺伝子疾患”，  メディカルレビュー社，  ２００２」参照。）を疾患の例にした実験結果を通して本発明の有効性を示す。
【０１３９】
本実施例のねらいは、情報の相互参照構造から、系でみたそれら情報間の関係性を捉え、本発明により、１万規模の遺伝子間の相互作用を通して、疾患関連遺伝子の探索を実現する点を示すことにある。
【０１４０】
（１．実験に用いるデータ）
ヒトの遺伝病に関するデータベースとしてよく知られるＯＭＩＭ（Ｏｎｌｉｎｅ  Ｍｅｎｄｅｌｉａｎ  Ｉｎｈｅｒｉｔａｎｃｅ  ｉｎ  Ｍａｎ）（例えば、「Ａ．Ｈａｍｏｓｈ，  ＡＦ．Ｓｃｏｔｔ，  Ｊ．Ａｍｂｅｒｇｅｒ，  Ｄ．Ｖａｌｌｅ，  ＶＡ．ＭｃＫｕｓｉｃｋ，  “Ｏｎｌｉｎｅ  Ｍｅｎｄｅｌｉａｎ  Ｉｎｈｅｒｉｔａｎｃｅ  ｉｎ  Ｍａｎ  （ＯＭＩＭ）”，  Ｈｕｍ．  Ｍｕｔａｔ．，  ｖｏｌ．１５，  Ｎｏ．１，  ｐｐ５７−６１，  ２０００」参照。）を本実験では用いる。
【０１４１】
ＯＭＩＭは疾患や遺伝子（上述した本実施形態における概念に対応）について記述されたページ（上述した本実施形態における概念情報に対応）から構成され、各々のページはその疾患や遺伝子を単位に番号（上述した本実施形態における概念情報識別情報に対応）が付与されている。各ページの記述に、多型をはじめとする医学生物学的知見に基づく結び付きから、疾患や遺伝子のページがリンク（上述した本実施形態における参照情報に対応）により参照されており、疾患や遺伝子間の関係はそれらのページ中の番号により相互参照情報として取得した。
【０１４２】
疾患と遺伝子の計１３７２２ページ（２００１年１２月時点）の間で、全ての組み合せの約０．０２３％の４３９５６のリンクが張られている。なお、本実験では、重複する同じリンクは２度以上数えないこととする。
【０１４３】
（２．予備実験）
ここでは、疾患や遺伝子にそれぞれ対応するページ群における参照構造から得られる存在確率「Ｒ」の値（上述した本実施形態における数式１参照。）について述べる。
【０１４４】
図７は、遺伝子のみを対象に存在確率Ｒの値を示す図である。なお、遺伝子は、遺伝子名の得られる１０００１件を対象としている。また、遺伝子名は、各ページのタイトル部分から得た。図７では、存在確率Ｒの値を縦軸に、横軸にはその値の降順に遺伝子を並べている。
【０１４５】
図７に示すように、相対的に大きな値をもつ数百の遺伝子がある一方で、全体の約７割から８割は同等の小さな値となる構造となっていることがわかる。ここで、ネットワーク研究（例えば、「アルバート＝ラズロ・バラバシ，  “新ネットワーク思考−世界のしくみを読み解く”，  ＮＨＫ出版，  ２００２」参照。）によると、インターネットにおけるＷｅｂページのアクセス人気度やＷｅｂページ間を結ぶ構造などは、ベキ法則に従うことが知られている。ベキ法則では上限がなく、下限に近づくほど増える分布となり、上には上が下には下がいる関係がどこにも同じように現れる特性から、スケールフリーと呼ばれる。例えば、大腸菌の代謝ネットワークなどに、スケールフリー・トポロジーが見出され、生物学的ネットワークの考察の１つの捉え方として興味深い。
【０１４６】
図８は、図７で示した存在確率Ｒの値の上位３０の遺伝子を示す図である。図８において、被リンク数の多い遺伝子が挙がっていることから、これら遺伝子は、疾患との関係や生物学的な意義などの理由から研究対象として注目されるハブ的役割の遺伝子とみることができる。ここで、存在確率Ｒの値は、当該遺伝子のもつ「票」の数に相当し、ネットワーク全体に対する関わりの深さを意味する量となっている。例えば、図８において、第１位の「ＴＰ５３」という遺伝子は、ヒト腫瘍で変異のみられる癌抑制遺伝子（例えば、「医学大辞典，  南山堂，  １９９８」参照。）として、極めて著名な遺伝子である。
【０１４７】
（３．疾患関連遺伝子の探索）
ここでは、本発明による疾患関連遺伝子の抽出において、注目する疾患の存在確率Ｒの値を予めもって基準とし、遺伝子をノックアウトすることによる変動から、関係の序列を得る。具体的には、高血圧を疾患の例に、疾患関連遺伝子抽出の実験結果について述べる。
【０１４８】
ここで、高血圧は、発症に少なくとも１０数個は存在する（例えば、「勝谷  友宏，  槍垣  實男，  荻原  俊男，  “高血圧症の病因遺伝子”，  最新医学，  Ｖｏｌ．５５，  Ｎｏ．１１，  ２０００」参照。）とされる複数の遺伝子が複雑に関与する代表的な多遺伝子性疾患である。そして、複数の遺伝子のうち、発症に強く関与する１〜数個の遺伝子を指す主効果遺伝子の見極めが重要とされている（例えば、「佐藤  憲幸，  槍垣  實男，  荻原  俊男，  “高血圧”，  臨床検査，  Ｖｏｌ．４２，  Ｎｏ．１２，  ｐｐ．１５４３−１５４８，  １９９８」参照。）。
【０１４９】
高血圧は、ＯＭＩＭの番号が１４５５００の疾患のページで記述される（例えば、「三木  哲郎，  “高血圧遺伝子解析プロジェクト”，  遺伝子医学，  Ｖｏｌ．５，  Ｎｏ．４，  ｐｐ．１２６−１３０，  ２００１」参照。）。なお、当該ページは、２６の外に向けたリンク、１３の被リンクがあり、存在確率Ｒ（１４５５００）は、０．０００５３５が得られている。
【０１５０】
関連遺伝子の抽出は、遺伝子のページをノックアウト、即ち、リンク構造から除外したことによる、この１４５５００のページに関する変動ｄＲ（上述した本実施形態における数式３を参照。）を序列化することで実現される。
【０１５１】
本実験では、処理量の観点から、便宜的に、アルツハイマー病、癌、糖尿病、高血圧および喘息の５つの疾患に関連が知られる、それぞれ、５６、４２８、９４、３９および２４の遺伝子からの計６２０遺伝子をノックアウトの対象とした。
【０１５２】
ここで、この５つの疾患は、ミレニアム・ゲノム・プロジェクトで克服を目指すターゲットの５大疾患に相当する（例えば、「三木  哲郎，  “高血圧遺伝子解析プロジェクト”，  遺伝子医学，  Ｖｏｌ．５，  Ｎｏ．４，  ｐｐ．１２６−１３０，  ２００１」参照。）。また、遺伝子は、遺伝子情報ポータルとしてよく知られる「ＧｅｎｅＣａｒｄｓ（例えば、「Ｍ．Ｒｅｂｈａｎ  ｅｔ  ａｌ．，  “ＧｅｎｅＣａｒｄｓ：  ａ  ｎｏｖｅｌ  ｆｕｎｃｔｉｏｎａｌ  ｇｅｎｏｍｉｃｓ  ｃｏｍｐｅｎｄｉｕｍ  ｗｉｔｈ  ａｕｔｏｍａｔｅｄ  ｄａｔａ  ｍｉｎｉｎｇ  ａｎｄ  ｑｕｅｒｙ  ｒｅｆｏｒｍｕｌａｔｉｏｎ  ｓｕｐｐｏｒｔ”，  Ｂｉｏｉｎｆｏｒｍａｔｉｃｓ，  １４，  ｐｐ６５６−６６４，  １９９８」参照。）」において、“Ａｌｚｈｅｉｍｅｒ”、“ｃａｎｃｅｒ”、“ｄｉａｂｅｔｅｓ”、“ｈｙｐｅｒｔｅｎｓｉｏｎ”および“ａｓｔｈｍａ”の検索後で得られた（２００１年１２月時点）中から、ＯＭＩＭの番号のある遺伝子を関連遺伝子として得た。なお、２１の遺伝子については、５つのうちの複数の疾患で重複して関連する。
【０１５３】
図９は、６２０の各遺伝子のノックアウトによる高血圧のページに関するｄＲ（１４５５００）の値の分布を示す図である。図９において、ｄＲの値を縦軸に、横軸はその値の昇順に遺伝子を並べている。
【０１５４】
図９に示すように、両サイドに偏って値をもっていることから、関連が顕著に示される遺伝子が全体の一部に限られていることがわかる。
【０１５５】
また、任意の２つのページをリンクによりつなぐ経路（但し、向きは考えない。）の中で、最短で経由するページの数を「隣接長」と定義するとき、６２０の遺伝子と高血圧のページとの隣接長Ｌ（１４５５００）を、図９と同じ横軸で、図１０に示す。図１０は、高血圧に関するｄＲ（１４５５００）の値の序列に対する隣接長Ｌ（１４５５００）の分布を示す図である。ここで、高血圧のページにおける隣接長の最大値は５であり、図１０において、隣接長が５を超える遺伝子は、リンクでつながりを持たないことを意味している。
【０１５６】
図１０に示すように、図９と対応して、中央部に上に凸の形状が示され、遺伝子疾患へのリンク上の隔たりと関連の深さに相関が確認される。
【０１５７】
図１１および図１２は、それぞれ、ｄＲ（１４５５００）の値の負の側および正の側からの上位２０による高血圧の関連遺伝子を示す図である。なお、図１１および図１２には、各遺伝子の高血圧のページとの隣接長Ｌ、ＯＭＩＭ本文のページの中の“ｈｙｐｅｒｔｅｎｓｉｏｎ”の語の頻度、並びに、ＧｅｎｅＣａｒｄｓ（例えば、「Ｍ．Ｒｅｂｈａｎ  ｅｔ  ａｌ．，  “ＧｅｎｅＣａｒｄｓ：  ａ  ｎｏｖｅｌ  ｆｕｎｃｔｉｏｎａｌ  ｇｅｎｏｍｉｃｓ  ｃｏｍｐｅｎｄｉｕｍ  ｗｉｔｈ  ａｕｔｏｍａｔｅｄ  ｄａｔａ  ｍｉｎｉｎｇ  ａｎｄ  ｑｕｅｒｙ  ｒｅｆｏｒｍｕｌａｔｉｏｎ  ｓｕｐｐｏｒｔ”，  Ｂｉｏｉｎｆｏｒｍａｔｉｃｓ，  １４，  ｐｐ６５６−６６４，  １９９８」参照。）および高血圧候補遺伝子リスト（例えば、「ＭＫ．Ｈａｌｕｓｈｋａ  ｅｔ  ａｌ．，  “ＧＩＳＴ：  Ａ  ｗｅｂ  ｔｏｏｌ  ｆｏｒ  ｃｏｌｌｅｃｔｉｎｇ  ｇｅｎｅ  ｉｎｆｏｒｍａｔｉｏｎ”，  Ｐｈｙｓｉｏｌ．  Ｇｅｎｏｍｉｃｓ，  ｖｏｌ．１，  Ｎｏ．２，  ｐｐ７５−８１，  １９９９」参照。）で高血圧に関連するとされる遺伝子（高血圧候補遺伝子リストにおいて、実験対象とした６２０遺伝子では４１遺伝子が含まれる。）が、それぞれ、“Ｈ”および“ｈ”の印により、あわせて示されている。ここで、図１１および図１２において、“ｉ”および“ｏ”は、それぞれ、当該遺伝子のページから高血圧のページへ、および、高血圧のページから当該遺伝子のページへ、隣接長０の直接のリンクのあることを表している。
【０１５８】
ここで、上述したサーチエンジン（例えば、「山名  早人，  近藤  秀和，  “解説  サーチエンジン  Ｇｏｏｇｌｅ”，  情報処理，  Ｖｏｌ．４２，  Ｎｏ．８，  ｐｐ．７７５−７８０，  ２００１」参照。）では、検索語の含まれたページを対象に、Ｒの値の降順により、検索結果を得る方法が採られている。この方法によれば、本発明で抽出された図１１や図１２において語の頻度が０の遺伝子は含まれず、範囲が限定された結果となる。
【０１５９】
図１３は、６２０遺伝子において、高血圧のページ（ｈｙｐｅｒｔｅｎｓｉｏｎ）からの隣接長が０または１となる３４遺伝子を対象に、高血圧のページとのリンク構造を矢印の向きで示す図である。なお、図１３において、網掛けは、ＧｅｎｅＣａｒｄｓでの関連遺伝子（図１１および図１２の“Ｈ”の印に対応）を表している。ここで、図１３は、高血圧と３４遺伝子のページに関する部分的な構造であり、解析の基となる全遺伝子の構造とは厳密には様相が異なる。すなわち、図１１および図１２の隣接長は全ページから算出され、６２０遺伝子で算出した図１３での隣接長０または１と一致しない。なお、多数のページからなるリンク構造のビジュアライゼーションは、それ自体が新たな別の研究テーマといえる（例えば、「ＥＲ．Ｇａｎｓｎｅｒ，  ＳＣ．Ｎｏｒｔｈ，  “Ａｎ  ｏｐｅｎ  ｇｒａｐｈ  ｖｉｓｕａｌｉｚａｔｉｏｎ  ｓｙｓｔｅｍ  ａｎｄ  ｉｔｓ  ａｐｐｌｉｃａｔｉｏｎｓ  ｔｏ  ｓｏｆｔｗａｒｅ  ｅｎｇｉｎｎｅｒｉｎｇ”，  Ｓｏｆｔｗ．  Ｐｒａｃｔ．  Ｅｘｐｅｒ．，  Ｓ１，  ｐｐ１−５，  １９９９」参照。）。
【０１６０】
図１３に示すように、高血圧のページと直接リンクのあるページは、２８件ある中で対象の６２０遺伝子の範囲では１５件あり、ＧｅｎｅＣａｒｄｓによる関連遺伝子に含まれないＡＣＥ、ＴＮＦ、ＣＡＳＴ遺伝子も含めて図１１および図１２に全て挙げられている。
【０１６１】
図１１および図１２をみると、その直接リンクのある“ｉ”や“ｏ”の印の遺伝子の他に、直接リンクのない、即ち、高血圧の記述されるページをみる限りでは関連を知り得ない遺伝子が含まれていることがわかる。
【０１６２】
ここで、ＶＨＬ遺伝子は、リンクのある遺伝子に挟まれる形で、図１１の５番目の上位に挙げられているが、“ｈｙｐｅｒｔｅｎｓｉｏｎ”の語の頻度が４となっている通り、ＯＭＩＭ本文で高血圧に関連するとの記述があり、文献「安部  道子，  名倉  潤，  三木  哲郎，  “高血圧感受性遺伝子を探して”，  医学のあゆみ，  Ｖｏｌ．１９７，  Ｎｏ．１３，  ｐｐ．１１１７−１１２１，  ２００１」でも、単一遺伝子による高血圧原因遺伝子として関連が知られている。
【０１６３】
また、図１１の１２番目のＡＶＰＲ２遺伝子、２０番目のＩＡＰＰ遺伝子、図１２の５番目のＩＮＳ遺伝子、１６番目のＴＧＦＢ１遺伝子、２０番目のＡＤＲＢ２遺伝子は、図１３で描かれた遺伝子群に含まれないが、“ｈ”の印を通して妥当に挙がっていることが確認できる。また、他には、図１１の１９番目のＲＥＴ遺伝子、図１２の１４番目のＷＴ１遺伝子では、共にＯＭＩＭ本文中に“ｈｙｐｅｒｔｅｎｓｉｏｎ”の語がなく、また“Ｈ”や“ｈ”の印もなくてあまり知られていないが、前者は文献「安部  道子，  名倉  潤，  三木  哲郎，  “高血圧感受性遺伝子を探して”，  医学のあゆみ，  Ｖｏｌ．１９７，  Ｎｏ．１３，  ｐｐ．１１１７−１１２１，  ２００１」で単一遺伝子による高血圧原因遺伝子として、後者は文献「Ｒ．Ｄａｖｉｅｓ，  Ａ．Ｍｏｏｒｅ，  Ａ．Ｓｃｈｅｄｌ，  Ｅ．Ｂｒａｔｔ，  Ｋ．Ｍｉｙａｈａｗａ，  Ｍ．Ｌａｄｏｍｅｒｙ，  Ｃ．Ｍｉｌｅｓ，  Ａ．Ｍｅｎｋｅ，  Ｖ．ｖａｎ  Ｈｅｙｎｉｎｇｅｎ，  Ｎ．Ｈａｓｔｉｅ，  “Ｍｕｌｔｉｐｌｅ  ｒｏｌｅｓ  ｆｏｒ  ｔｈｅ  Ｗｉｌｍｓ'  ｔｕｍｏｒ  ｓｕｐｐｒｅｓｓｏｒ，  ＷＴ１”，  Ｃａｎｃｅｒ  Ｒｅｓ．，  ５９，  ｐｐ．１７４７−１７５０，  １９９９」で高血圧との関連の報告があり、妥当といえる。
【０１６４】
これら遺伝子は、リンク構造に基づく関係性の解析により発見されたとみることができ、本発明により有用な疾患関連遺伝子が得られることがわかる。また、本発明による疾患関連遺伝子の序列結果は、生物学的知見と対応するものであると考えられる。
【０１６５】
ここで、高血圧に密接に関わるレニン−アンジオテンシン系と呼ばれる血圧調節機構があり、レニン−アンジオテンシン系に関連する遺伝子が高血圧関連遺伝子として最も重要（例えば、「檜垣  實男，  勝谷  友宏，  荻原  俊男，  “遺伝疫学：高血圧”，  循環器科，  ４８，  ｐｐ．１５８−１６２，  ２０００」参照。）で、ＡＣＥ、ＡＧＴ、ＡＧＴＲ１遺伝子の３つがよく取り上げられる（例えば、「Ｂ．Ｂａｕｄｉｎ，  “Ａｎｇｉｏｔｅｎｓｉｎ  ＩＩ  ｒｅｃｅｐｔｏｒ  ｐｏｌｙｍｏｒｐｈｉｓｍ  ｉｎ  ｈｙｐｅｒｔｅｎｓｉｏｎ．  ｈａｒｍｏｃｏｇｅｎｏｍｉｃ  ｃｏｎｓｉｄｅｒａｔｉｏｎｓ．”，  Ｐｈａｒｍａｃｏｇｅｎｏｍｉｃｓ，  ３，  ｐｐ．６５−７３，  ２００２」参照。）とされている。
【０１６６】
ここで、レニン−アンジオテンシン系の関連遺伝子について簡単に説明する。レニン−アンジオテンシン系の関連遺伝子としては、ＡＧＴ、ＲＥＮ、ＡＣＥ、ＡＧＴＲ１遺伝子が知られている。ホルモンの前駆体であるアンジオテンシノーゲン（ＡＧＴ）は、レニン（ＲＥＮ）、アンジオテンシン変換酵素（ＡＣＥ）の作用によって活性化されてアンジオテンシンレセプター（ＡＧＴＲ１）に結合して血管壁を収縮させ、血圧が上昇するとされている。なお、ＲＥＮ遺伝子は、対象の６２０遺伝子に含まれていない。
【０１６７】
本実験結果は、それと符合して、それら３つが図１１と図１２の上位（図１１の２番目、３番目、図１２の１番目）に自動抽出されることは、特筆すべき点である。
【０１６８】
ここで、アルドステロン分泌によるナトリウム代謝を介したフィードバックにより血圧を調節する系は、医学生物学的知見としてよく知られている（例えば、「医学大辞典，  南山堂，  １９９８」参照。）。ＡＣＥ、ＡＧＴ遺伝子が昇圧に作用するアンジオテンシンに関係する一方で、ＳＣＮＮ１Ｂ、ＳＣＮＮ１Ｇ遺伝子はナトリウム、ＡＧＴＲ１、ＣＹＰ１１Ｂ１、ＣＹＰ１１Ｂ２遺伝子はアルドステロンに関係することを対比すると、図１１と図１２に医学生物学的知見との整合が考えられる。
【０１６９】
（４．複数の疾患における関連遺伝子の探索）
ここでは、複数の疾患に共通の関連遺伝子を抽出する実験結果について述べる。本発明を用いることにより、疾患に応じた複数の変動ｄＲ（上述した本実施形態における数式３参照。）から算出される主成分得点で対象疾患への総合的寄与を表し、その序列で共通に関連する遺伝子発見を図る。
【０１７０】
本実験では、生活習慣病として文献「佐々木  博巳，  “バイオ研究イラストマップ〜研究テーマ別の重要因子と要点がすぐわかる〜”，  羊土社，  ２００１」により互いに関わりの知られる高血圧、糖尿病、高脂血症および肥満の４つの疾患を対象とした。
【０１７１】
ここで、この４つの疾患に共通な基盤としてインスリン抵抗性が知られており、心筋梗塞や脳卒中などの動脈硬化性疾患へつながる医学的知見が得られている（例えば、「堀内  正嗣，  福田  恵一，  森下  竜一，  “生活習慣と遺伝子疾患”，  メディカルレビュー社，  ２００２」、「佐々木  博巳，  “バイオ研究イラストマップ〜研究テーマ別の重要因子と要点がすぐわかる〜”，  羊土社，  ２００１」参照。）。ここで、インスリンとは、筋肉の運動時などに必要なブドウ糖をエネルギーに変えるホルモンである。このインスリンが十分あるのに利用されず血糖値が高くなる状態、即ち、インスリンの作用不足が、インスリン抵抗性がある、と呼ばれる。インスリン抵抗性は、肥満者に多く、インスリン抵抗性の状態になると脂質合成が促進されたり、血圧も上がるなどの悪循環から、臨床的に危険とされている（例えば、「医学大辞典，  南山堂，  １９９８」参照。）。
【０１７２】
糖尿病、高脂血症および肥満の疾患に対応するページは、ＯＭＩＭの番号がそれぞれ、１２５８５３、１４４２５０および６０１６６５である。なお、各番号の存在確率Ｒ（１２５８５３）、Ｒ（１４４２５０）およびＲ（６０１６６５）は、それぞれ、０．００１００３、０．０００１９４および０．０００２１１が得られている。また、本実験で対象とする遺伝子は前節と同じ６２０遺伝子である。
【０１７３】
ここでは、対象の４つの疾患が相乗的に進展して引き起こされる心筋梗塞に着目し、その関連遺伝子（例えば、「Ｙ．Ｙａｍａｄａ  ｅｔ  ａｌ．，  “Ｐｒｅｄｉｃｔｉｏｎ  ｏｆ  ｔｈｅ  ｒｉｓｋ  ｏｆ  ｍｙｏｃａｒｄｉａｌ  ｉｎｆａｒｃｔｉｏｎ  ｆｒｏｍ  ｐｏｌｙｍｏｒｐｈｉｓｍｓ  ｉｎ  ｃａｎｄｉｄａｔｅ  ｇｅｎｅｓ”，  Ｎ．  Ｅｎｇｌ．  Ｊ．  Ｍｅｄ．，  ｖｏｌ．３４７，  Ｎｏ．２４，  ｐｐ１９１６−１９２３，  ２００２」参照。）のうち実験対象に該当する１６遺伝子の分布を通して評価する。
【０１７４】
第１主成分得点の昇順による序列において、得点の小さい遺伝子の方から、１〜２００、２０１〜４２０、４２１〜６２０の３つの範囲にわけるとき、心筋梗塞関連遺伝子は、順に、９、２および５件が含まれる分布が得られ、両サイドへの偏りから医学的知見との整合性が認められる。
【０１７５】
ここで、本実験と同様に、第１主成分得点に応じた６２０遺伝子の序列を３つにわけるとき、骨格筋に関連する１２遺伝子は、６、１および５件が含まれる分布となる。骨格筋は、糖全体の６０％以上を消費するためにインスリン作用で最も重要な組織であることから、インスリン抵抗性の成因に骨格筋が注目されており（例えば、「船橋  徹，  松澤  佑次，  “肥満症の分子機構”，  最新医学，  Ｖｏｌ．５５，  Ｎｏ．１１，  ｐｐ．６５−７４，  ２０００」参照。）、ここでも医学的知見との整合性がみられる。この骨格筋のような関連臓器は、蛋白質配列データベースＳＷＩＳＳ−ＰＲＯＴ（例えば、「Ａ．Ｂａｉｒｏｃｈ，  Ｒ．Ａｐｗｅｉｌｅｒ，  “Ｔｈｅ  ＳＷＩＳＳ−ＰＲＯＴ  ｐｒｏｔｅｉｎ  ｓｅｑｕｅｎｃｅ  ｄａｔａｂａｓｅ  ａｎｄ  ｉｔｓ  ｓｕｐｐｌｅｍｅｎｔ  ＴｒＥＭＢＬ  ｉｎ  ２０００”，  Ｎｕｃｌｅｉｃ．  Ａｃｉｄｓ  Ｒｅｓ．，  ｖｏｌ．２８，  Ｎｏ．１，  ｐｐ４５−４８，  ２０００」参照。）より、遺伝子毎に得られ（２００２年１１月時点）、具体的には各遺伝子の“Ｒｅｆｅｒｅｎｃｅｓ”のパートの引用文献に付与される“Ｃｏｍｍｅｎｔｓ”の項目にて、“ＴＩＳＳＵＥ＝”に続いて記されている臓器を関連臓器としている。
【０１７６】
図１４には、横軸を４つの変動ｄＲ（上述した本実施形態における数式３参照。）から算出される第１主成分得点、縦軸をＲの値（上述した本実施形態における数式１参照。）にして６２０遺伝子の分布を示している。図１４は、４つの疾患（高血圧、糖尿病、高脂血症および肥満）のｄＲからの第１主成分得点とＲの値との関係を示す図である。なお、深い関連が示される横軸の絶対値の大きい遺伝子には遺伝子名を記した。
【０１７７】
図１４に示すように、ここでも上述した図９と同様、深い関連の示される遺伝子は、比較的限定される傾向にあることがわかる。図１４において、遺伝子名が記されている遺伝子について、医学生物学的知見を照らし合わせると次の通りである。
【０１７８】
インスリン遺伝子（ＩＮＳ）およびインスリン受容体遺伝子（ＩＮＳＲ）は、両者の結合によりインスリン作用が発揮されることから、インスリン抵抗性の医学的知見との整合性が認められる。図１４では、ＴＮＦ遺伝子に関連が顕著に示され、続いてＴＧＦＢ１遺伝子、ＡＰＯＥ遺伝子に深い関連が示されている。これらの遺伝子は、いずれも心筋梗塞関連遺伝子（例えば、「Ｙ．Ｙａｍａｄａ  ｅｔ  ａｌ．，“Ｐｒｅｄｉｃｔｉｏｎ  ｏｆ  ｔｈｅ  ｒｉｓｋ  ｏｆ  ｍｙｏｃａｒｄｉａｌ  ｉｎｆａｒｃｔｉｏｎ  ｆｒｏｍ  ｐｏｌｙｍｏｒｐｈｉｓｍｓ  ｉｎ  ｃａｎｄｉｄａｔｅ  ｇｅｎｅｓ”，  Ｎ．  Ｅｎｇｌ．  Ｊ．  Ｍｅｄ．，  ｖｏｌ．３４７，  Ｎｏ．２４，  ｐｐ１９１６−１９２３，  ２００２」参照。）に含まれ、妥当と判断できる。ここで、ＴＮＦ遺伝子、ＴＧＦＢ１遺伝子は細胞間情報伝達の役割を担うサイトカインである。近年、動脈硬化へのサイトカインの関わりが注目されている。
【０１７９】
エネルギーバランスに関与する遺伝子として知られるＵＣＰ１遺伝子及び糖代謝に関与する遺伝子として知られるＨＮＦ４Ａ遺伝子は、それぞれ、文献「Ｊ．Ｋｌｅｉｎ，  Ｍ．Ｆａｓｓｈａｕｅｒ，  Ｍ．Ｂｅｎｉｔｏ，  ＣＲ．Ｋａｈｎ，  “Ｉｎｓｕｌｉｎ  ａｎｄａ  ｔｈｅ  ｂｅｔａ３−ａｄｒｅｎｏｃｅｐｔｏｒ  ｄｉｆｆｅｎｔｉａｌｌｙ  ｒｅｇｕｌａｔｅ  ｕｎｃｏｕｐｌｉｎｇ  ｐｒｏｔｅｉｎ−１  ｅｘｐｒｅｓｓｉｏｎ”，  Ｍｏｌ．Ｅｎｄｏｃｒｉｎｏｌ．，  ｖｏｌ．１４，  Ｎｏ．６，  ｐｐ７６４−７７３，  ２０００」および文献「武田  純，  “遺伝子異常と糖尿病”，  日医大誌，  Ｖｏｌ．６６，  Ｎｏ．５，  ｐｐ．４１−４３，  １９９９」でインスリン抵抗性との関連が報告され、妥当性が確認できる。
【０１８０】
また、ＧＹＳ１遺伝子にも、インスリン抵抗性との関連がよく知られ（例えば、「Ｈ．Ｓｈｉｍｏｍｕｒａ，  Ｔ．Ｓａｎｋｅ，  Ｋ．Ｕｅｄａ，  Ｔ．Ｈａｎａｂｕｓａ，  Ｓ．Ｓａｋａｇａｓｈｉｒａ，  Ｋ．Ｎａｎｊｏ，  “Ａ  ｍｉｓｓｅｎｓｅ  ｍｕｔａｔｉｏｎ  ｏｆ  ｔｈｅ  ｍｕｓｃｌｅ  ｇｌｙｃｏｇｅｎ  ｓｙｎｔｈａｓｅ  ｇｅｎｅ  （Ｍ４１６Ｖ）  ｉｓ  ａｓｓｏｃｉａｔｅｄ  ｗｉｔｈ  ｉｎｓｕｌｉｎ  ｒｅｓｉｓｔａｎｃｅ  ｉｎ  ｔｈｅ  Ｊａｐａｎｅｓｅ  ｐｏｐｕｌａｔｉｏｎ”，  Ｄｉａｂｅｔｏｌｏｇｉａ．，  ｖｏｌ．４０，  Ｎｏ．８，  ｐｐ９４７−９５２，  １９９７」参照。）、いずれも有効な抽出結果とみることができる。ここで、Ｒの値は、上述したように、系全体に与える影響の量と相関があり、副作用回避の観点からＲの値の小さい順を指標に創薬ターゲットを選定することが考えられる。
【０１８１】
本実験の４つの疾患のような生活習慣病を共に対象にして、遺伝子間の相互作用の考慮された疾患関連遺伝子を抽出した検討例は、これまでみられず（例えば、「鈴木  洋一，  “疾患関連遺伝子の総合的解析と相互作用”，  メディカル・サイエンス・ダイジェスト，  Ｖｏｌ．２９，  Ｎｏ．７，  ｐｐ．２１−２４，  ２００３」参照。）、疾患関連遺伝子の抽出結果は、医学生物学的知見と対応するものであると考えられる。しかも、ノックアウトスタディによる生物実験で、必ずしも症状が見出されない場合のあることが指摘されている（例えば、「“産総研シリーズ  ポストゲノム”，  丸善株式会社，  ３．５  タンパク質ネットワーク解析，  ｐｐ．７８−８９，  ２００２」参照。）。
【０１８２】
本発明は、実験の設計と検証による知見蓄積のサイクルにおいて、生物学的実験と相補的な位置付けにあり、上述したような有用な疾患関連遺伝子がリンク構造を基に計算機による自動処理を通して得られることが特長といえる。
【０１８３】
本発明は、情報間の参照構造に基づく関係性の解析を提案する。本実施例では、本発明の疾患関連遺伝子探索への応用について述べ、高血圧をはじめ生活習慣病を疾患の例に、有用な疾患関連遺伝子が得られることをヒトの遺伝病に関するデータベースＯＭＩＭを用いた実験結果を通して示した。
【０１８４】
生命現象は、遺伝子の発現により生成されるタンパク質が相互に作用して成り立っており、この相互作用を通して情報が伝達する現象は遺伝子カスケードと呼ばれる（例えば、「村松  正実，  “解読されたゲノム情報をどう活かすか”，  東京化学同人，  １２．  遺伝子カスケードと新しい創薬プラン，  ｐｐ．１０３−１１１，  ２００１」参照。）。
【０１８５】
疾患は、この遺伝子カスケードが病因により変調をきたし、この変調を乗り切って恒常性を保つために各遺伝子が防御的に働く状態といえ、健常状態に比べて多くの遺伝子の発現に変化が生じる。ここで、例えば、高血圧の発症において、ＡＧＴ遺伝子とＡＣＥ遺伝子が相乗的に関連すること（例えば、「ＣＴ．Ｔｓａｉ，  Ｄ．Ｆａｌｌｉｎ，  ＦＴ．Ｃｈｉａｎｇ，  ＪＪ．Ｈｗａｎｇ，  ＬＰ．Ｌａｉ，  ＫＬ．Ｈｓｕ，  ＣＤ．Ｔｓｅｎｇ，  ＣＳ．Ｌｉａｕ，  ＹＺ．Ｔｓｅｎｇ，  “Ａｎｇｉｏｔｅｎｓｉｎｏｇｅｎ  ｇｅｎｅ  ｈａｐｌｏｔｙｐｅ  ａｎｄ  ｈｙｐｅｒｔｅｎｓｉｏｎ：  ｉｎｔｅｒａｃｔｉｏｎ  ｗｉｔｈ  ＡＣＥ  ｇｅｎｅ  Ｉ  ａｌｌｅｌｅ”，  Ｈｙｐｅｒｔｅｎｓｉｏｎ，  ｖｏｌ．４１，  Ｎｏ．１，  ｐｐ９−１５，  ２００３」参照。）、ＡＣＥ遺伝子とＡＧＴＲ１遺伝子が共働すること（例えば、「高橋  規郎，  村上  秀子，  児玉  和紀，  笠置  文善，  山田  美智子，  西下  聡英，  稲上  正，  “１型アンギオテンシンＩＩ型受容体遺伝子の５'−領域の多型と高血圧症との相関”，  放影研報告書，  １−９９，  １９９９」参照。）をはじめとした相互作用（例えば、「ＳＭ．Ｗｉｌｌｉａｍｓ，  ＪＨ．Ａｄｄｙ，  ＪＡ３ｒｄ．Ｐｈｉｌｌｉｐｓ，  Ｍ．Ｄａｉ，  Ｊ．Ｋｐｏｄｏｎｕ，  Ｊ．Ａｆｆｕｌ，  Ｈ．Ｊａｃｋｓｏｎ，  Ｋ．Ｊｏｓｅｐｈ，  Ｆ．Ｅａｓｏｎ，  ＭＭ．Ｍｕｒｒａｙ，  Ｐ．Ｅｐｐｅｒｓｏｎ，  Ａ．Ａｄｕｏｎｕｍ，  ＬＪ．Ｗｏｎｇ，  ＰＡ．Ｊｏｓｅ，  ＲＡ．Ｆｅｌｄｅｒ，  “Ｃｏｍｂｉｎａｔｉｏｎｓ  ｏｆ  ｖａｒｉａｔｉｏｎｓ  ｉｎ  ｍｕｌｔｉｐｌｅ  ｇｅｎｅｓ  ａｒｅ  ａｓｓｏｃｉａｔｅｄ  ｗｉｔｈ  ｈｙｐｅｒｔｅｎｓｉｏｎ”，  Ｈｙｐｅｒｔｅｎｓｉｏｎ．，  ｖｏｌ．３６，  Ｎｏ．１，  ｐｐ．２−６，  ２０００」参照。）が明らかにされており、上述したレニン−アンジオテンシン系を代表例とするカスケード理解が発症メカニズム解明と治療や創薬支援への応用に向けて重要となる。
【０１８６】
遺伝子間の相互作用の考慮された本発明による解析は、全体を系としてみるシステム生物学の命題に対する１つの実現形といえる。ここで、これまで、例えば、高血圧に関する医学生物学的検討では、環境要因や人種などの集団の違いにより、幾つかの実験結果が必ずしも一致しないことから、決定的な関連遺伝子が特定できない状況下にある（例えば、「安部  道子，  名倉  潤，  三木  哲郎，  “高血圧感受性遺伝子を探して”，  医学のあゆみ，  Ｖｏｌ．１９７，  Ｎｏ．１３，  ｐｐ．１１１７−１１２１，  ２００１」、「Ｂ．Ｂａｕｄｉｎ，  “Ａｎｇｉｏｔｅｎｓｉｎ  ＩＩ  ｒｅｃｅｐｔｏｒ  ｐｏｌｙｍｏｒｐｈｉｓｍ  ｉｎ  ｈｙｐｅｒｔｅｎｓｉｏｎ．  Ｐｈａｒｍｏｃｏｇｅｎｏｍｉｃ  ｃｏｎｓｉｄｅｒａｔｉｏｎｓ．”，  Ｐｈａｒｍａｃｏｇｅｎｏｍｉｃｓ，  ３，  ｐｐ．６５−７３，  ２００２」参照。）とされ、生活習慣病における遺伝要因と環境要因の分離が期待される。
【０１８７】
本発明は、対象を疾患に依存せず解析可能であるが、本実施例で取り上げた生活習慣病への適用が効果的と考えられる。
【０１８８】
また、遺伝子カスケードのもたらすダイナミクスの多様性を考察し、食事、運動やストレスなどの環境要因を引き金に疾患が引き起こされる機序理解の解析が今後期待される。
【０１８９】
また、本発明は、全体像のふかんを図る遺伝子機能の解析（例えば、「松永  務，  “大規模知識構造化のための知識表現モデルと遺伝子機能解析への応用”，  信学論，  Ｖｏｌ．Ｊ８６−Ｄ−ＩＩ，  Ｎｏ．８，  ｐｐ．１１９６−１２０４，  ２００３」参照。）と統合させつつ、対象遺伝子規模を数万規模の全遺伝子へ拡大した実験が可能である。
【０１９０】
本発明を基にした遺伝子カスケード推定に関する検討も可能である。
【０１９１】
本実施例では、リンク情報を無効にするだけであったが、本発明は、例えば、トラフィックに応じて地点間を結ぶ工事の際の最適な迂回路の設計や、人から人への情報の伝達、企業間の資金の流れや倒産の影響などのコネクションに応じたコミュニティ生成に関する予測などへの展開も可能である。つまり、本発明は、上述した本実施例における疾患関連遺伝子探索などのバイオ、ライフサイエンス、製薬分野に限られず、幅広く他の分野においても適用可能である。
【産業上の利用可能性】
【０１９２】
以上のように、本発明にかかる情報間関係性解析装置、情報間関係性解析方法、プログラムおよび記録媒体は、相互に参照される情報（例えば、リンクを張るＷｅｂページ）群においてターゲットとする情報を無効にしたときの他への影響を算出し、互いの関係性を定量化することができる。
【０１９３】
また、本発明にかかる情報間関係性解析装置、情報間関係性解析方法、プログラムおよび記録媒体は、疾患関連遺伝子探索や創薬ターゲット候補選定などを行うバイオインフォマティクス分野において極めて有用である。
【０１９４】
さらに、本発明にかかる情報間関係性解析装置、情報間関係性解析方法、プログラムおよび記録媒体は、産業上の多くの分野、特に、医薬品や医療、ライフサイエンス等の分野で広く実施することができ、極めて有用である。
【図面の簡単な説明】
【０１９５】
【図１】本発明の基本原理を示す原理構成図である。
【図２】本発明が適用される本システムの構成の一例を示すブロック図である。
【図３】特定概念情報関係情報ファイル１０６ｂに格納される情報の一例を示す図である。
【図４】本実施形態における本システムのメイン処理の一例を示すフローチャートである。
【図５】７つのページとそれらの相互参照関係の例を示す図である。
【図６】図５におけるページのノックアウトの例を示す図である。
【図７】ＯＭＩＭにおける遺伝子のみを対象に存在確率Ｒの値を示す図である。
【図８】図７で示した存在確率Ｒの値の上位３０の遺伝子を示す図である。
【図９】６２０の各遺伝子のノックアウトによる高血圧のページに関するｄＲ（１４５５００）の値の分布を示す図である。
【図１０】高血圧に関するｄＲ（１４５５００）の値の序列に対する隣接長Ｌ（１４５５００）の分布を示す図である。
【図１１】ｄＲ（１４５５００）の値の負の側からの上位２０による高血圧の関連遺伝子を示す図である。
【図１２】ｄＲ（１４５５００）の値の正の側からの上位２０による高血圧の関連遺伝子を示す図である。
【図１３】高血圧と６２０遺伝子中で隣接長が０または１となる３４遺伝子によるリンク構造を示す図である。
【図１４】４つの疾患（高血圧、糖尿病、高脂血症および肥満）のｄＲからの第１主成分得点とＲの値との関係を示す図である。
【符号の説明】
【０１９６】
１００  情報間関係性解析装置
１０２  制御部
１０２ａ  概念情報群取得部
１０２ｂ  特定概念情報設定部
１０２ｃ  存在確率算出部
１０２ｄ  無効後存在確率算出部
１０２ｅ  影響度算出部
１０２ｆ  主成分得点算出部
１０２ｇ  序列部
１０４  通信制御インターフェース部
１０６  記憶部
１０６ａ  概念情報群ファイル
１０６ｂ  特定概念情報関係情報ファイル
１０６ｃ  主成分得点ファイル
１０６ｄ  序列結果ファイル
１０８  入出力制御インターフェース部
１１２  入力装置
１１４  出力装置
２００  外部システム
３００  ネットワーク
【図１】
2005275794.tif 000028
【図２】
2005275794.tif 000029
【図３】
2005275794.tif 000030
【図４】
2005275794.tif 000031
【図５】
2005275794.tif 000032
【図６】
2005275794.tif 000033
【図７】
2005275794.tif 000034
【図８】
2005275794.tif 000035
【図９】
2005275794.tif 000036
【図１０】
2005275794.tif 000037
【図１１】
2005275794.tif 000038
【図１２】
2005275794.tif 000039
【図１３】
2005275794.tif 000040
【図１４】
2005275794.tif 000041