(57)【特許請求の範囲】
【請求項１】
分解系係数を含み分解速度を示す項と、生成系係数を含み生成速度を示す項とからなる遺伝子の速度式を前提として、遺伝子間の相互作用を表す遺伝子ネットワークを求める遺伝子発現相互作用推定装置であって、以下の要素を有することを特徴とする遺伝子発現相互作用推定装置
（１）非破壊株及び破壊株のタイム・コース値を入力装置により入力する非破壊株・破壊株タイム・コースデータ入力処理部
（２）非破壊株・破壊株タイム・コースデータ入力処理部が入力した非破壊株及び破壊株のタイム・コース値に対してタイム・コース値の対数値を演算装置により求め、求めたタイム・コース値の対数値を基に、タイム・コース値の対数値のスプライン補間係数を演算装置により算出し、
スプライン補間係数を用いて非破壊株及び破壊株のタイム・コース値を演算装置によりスプライン補間し、スプライン補間されたタイム・コース値を補間後タイム・コース値として演算装置により算出し、
補間後タイム・コース値に対して補間後タイム・コース値の対数値を演算装置により求め、求めた補間後タイム・コース値の対数値の微係数を演算装置により算出する前処理部
（３）前処理部が算出した補間後タイム・コース値のうち、破壊株の補間後タイム・コース値と、
前処理部が算出した補間後タイム・コース値の対数値の微係数のうち、破壊株の補間後タイム・コース値の対数値の微係数と
を使用して、最小二乗法により、分解系係数の最適値を演算装置により算出する分解系係数推定部
（４）分解系係数推定部が算出した分解系係数の最適値と、
前処理部が算出した補間後タイム・コース値のうち、非破壊株の補間後タイム・コース値と、
前処理部が算出した補間後タイム・コース値の対数値の微係数のうち、非破壊株の補間後タイム・コース値の対数値の微係数と
を使用して、最小二乗法により、生成系係数の最適値を演算装置により算出する生成系係数推定部。
【請求項２】
生成系係数推定部は、非破壊株の補間後タイム・コース値の対数値の微係数から、生成速度プロファイルを演算装置により算出し、算出した生成速度プロファイルを利用して、最小二乗法により生成系係数の最適値を演算装置により算出することを特徴とする請求項１記載の遺伝子発現相互作用推定装置。
【請求項３】
遺伝子発現相互作用推定装置は、更に、算出した生成系係数及び分解系係数に基づいて、遺伝子発現についての遺伝子間の相互作用を表すパスウェイを表示装置に描画するパスウェイ描画部を有することを特徴とする請求項１記載の遺伝子発現相互作用推定装置。
【請求項４】
分解系係数を含み分解速度を示す項と、生成系係数を含み生成速度を示す項とからなる遺伝子の速度式を前提として、遺伝子間の相互作用を表す遺伝子ネットワークを求める遺伝子発現相互作用推定装置であって演算装置と入力装置とを備える遺伝子発現相互作用推定装置による遺伝子発現相互作用推定方法において、
（１）非破壊株・破壊株タイム・コースデータ入力処理部により、
非破壊株及び破壊株のタイム・コース値を入力装置により入力する工程と、
（２）前処理部により、
非破壊株・破壊株タイム・コースデータ入力処理部が入力した非破壊株及び破壊株のタイム・コース値に対してタイム・コース値の対数値を演算装置により求め、求めたタイム・コース値の対数値を基に、タイム・コース値の対数値のスプライン補間係数を演算装置により算出し、
スプライン補間係数を用いて非破壊株及び破壊株のタイム・コース値を演算装置によりスプライン補間し、スプライン補間されたタイム・コース値を補間後タイム・コース値として演算装置により算出し、
補間後タイム・コース値に対して補間後タイム・コース値の対数値を演算装置により求め、求めた補間後タイム・コース値の対数値の微係数を演算装置により算出する工程と、
（３）分解系係数推定部により、
前処理部が算出した補間後タイム・コース値のうち、破壊株の補間後タイム・コース値と、
前処理部が算出した補間後タイム・コース値の対数値の微係数のうち、破壊株の補間後タイム・コース値の対数値の微係数と
を使用して、最小二乗法により、分解系係数の最適値を演算装置により算出する工程と、
（４）生成系係数推定部により、
分解系係数推定部が算出した分解系係数の最適値と、
前処理部が算出した補間後タイム・コース値のうち、非破壊株の補間後タイム・コース値と、
前処理部が算出した補間後タイム・コース値の対数値の微係数のうち、非破壊株の補間後タイム・コース値の対数値の微係数と
を使用して、最小二乗法により、生成系係数の最適値を演算装置により算出する工程と
を有する遺伝子発現相互作用推定方法。
【請求項５】
分解系係数を含み分解速度を示す項と、生成系係数を含み生成速度を示す項とからなる遺伝子の速度式を前提として、遺伝子間の相互作用を表す遺伝子ネットワークを求める遺伝子発現相互作用推定装置となるコンピュータであって演算装置と入力装置とを備えるコンピュータに、以下の処理を実行させるためのプログラムを記録したコンピュータ読み取り可能な記録媒体
（１）非破壊株・破壊株タイム・コースデータ入力処理部により、
非破壊株及び破壊株のタイム・コース値を入力装置により入力する処理
（２）前処理部により、
非破壊株・破壊株タイム・コースデータ入力処理部が入力した非破壊株及び破壊株のタイム・コース値に対してタイム・コース値の対数値を演算装置により求め、求めたタイム・コース値の対数値を基に、タイム・コース値の対数値のスプライン補間係数を演算装置により算出し、
スプライン補間係数を用いて非破壊株及び破壊株のタイム・コース値を演算装置によりスプライン補間し、スプライン補間されたタイム・コース値を補間後タイム・コース値として演算装置により算出し、
補間後タイム・コース値に対して補間後タイム・コース値の対数値を演算装置により求め、求めた補間後タイム・コース値の対数値の微係数を演算装置により算出する処理
（３）分解系係数推定部により、
前処理部が算出した補間後タイム・コース値のうち、破壊株の補間後タイム・コース値と、
前処理部が算出した補間後タイム・コース値の対数値の微係数のうち、破壊株の補間後タイム・コース値の対数値の微係数と
を使用して、最小二乗法により、分解系係数の最適値を演算装置により算出する処理
（４）生成系係数推定部により、
分解系係数推定部が算出した分解系係数の最適値と、
前処理部が算出した補間後タイム・コース値のうち、非破壊株の補間後タイム・コース値と、
前処理部が算出した補間後タイム・コース値の対数値の微係数のうち、非破壊株の補間後タイム・コース値の対数値の微係数と
を使用して、最小二乗法により、生成系係数の最適値を演算装置により算出する処理。
【請求項６】
分解系係数を含み分解速度を示す項と、生成系係数を含み生成速度を示す項とからなる遺伝子の速度式を前提として、遺伝子間の相互作用を表す遺伝子ネットワークを求める遺伝子発現相互作用推定装置となるコンピュータであって演算装置と入力装置とを備えるコンピュータに、以下の手順を実行させるためのプログラム
（１）非破壊株・破壊株タイム・コースデータ入力処理部により、
非破壊株及び破壊株のタイム・コース値を入力装置により入力する手順
（２）前処理部により、
非破壊株・破壊株タイム・コースデータ入力処理部が入力した非破壊株及び破壊株のタイム・コース値に対してタイム・コース値の対数値を演算装置により求め、求めたタイム・コース値の対数値を基に、タイム・コース値の対数値のスプライン補間係数を演算装置により算出し、
スプライン補間係数を用いて非破壊株及び破壊株のタイム・コース値を演算装置によりスプライン補間し、スプライン補間されたタイム・コース値を補間後タイム・コース値として演算装置により算出し、
補間後タイム・コース値に対して補間後タイム・コース値の対数値を演算装置により求め、求めた補間後タイム・コース値の対数値の微係数を演算装置により算出する手順
（３）分解系係数推定部により、
前処理部が算出した補間後タイム・コース値のうち、破壊株の補間後タイム・コース値と、
前処理部が算出した補間後タイム・コース値の対数値の微係数のうち、破壊株の補間後タイム・コース値の対数値の微係数と
を使用して、最小二乗法により、分解系係数の最適値を演算装置により算出する手順
（４）生成系係数推定部により、
分解系係数推定部が算出した分解系係数の最適値と、
前処理部が算出した補間後タイム・コース値のうち、非破壊株の補間後タイム・コース値と、
前処理部が算出した補間後タイム・コース値の対数値の微係数のうち、非破壊株の補間後タイム・コース値の対数値の微係数と
を使用して、最小二乗法により、生成系係数の最適値を演算装置により算出する手順。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、遺伝子間の相互作用を表す遺伝子ネットワークを求める遺伝子発現相互作用推定装置及び方法に係り、多くの遺伝子を、より早く、より小さな計算機資源で推定できる技術に関する。
【背景技術】
【０００２】
本発明は、遺伝子間の相互作用を表す遺伝子ネットワークを求めるためのものである。このネットワークの解析に用いられる遺伝子発現データを取得する方法に関して、多くの技術が開発され、近年そのデータの精度や取得数が増加している。
【０００３】
同時に、試験結果を解析する技術も研究されてきた。しかし、例えばＧＡ（遺伝的アルゴリズム）を使用した幾つかの試みでは（非特許文献１及び２）、一度に推定できる遺伝子数は５程度で、多くの計算機資源を必要とするように、遺伝子ネットワークの推定には多くの困難が伴う。
【０００４】
そのため、ＧＡやＳＡ（シミュレーイテットアニーリング）のような発見的な最適化手法の場合よりも、多くの遺伝子を、より早く、より小さな計算機資源で推定できる方法が求められている。
【非特許文献１】トミナガ・Ｄ、オカモト・Ｍ、マキ・Ｙ、ワタナベ・Ｓ、エグチ・Ｙ，「ノンリニアー・ニューメリカル・オプティマイゼイション・テクニック・ベースド・オン・ア・ジェネチィック・ネットワークス：トゥワーズ・ザ・インファレンス・オブ・ジェネティック・ネットワーク（Ｎｏｎｌｉｎｅａｒ  ｎｕｍｅｒｉｃａｌ  ｏｐｔｉｍｉｚａｔｉｏｎ  ｔｅｃｈｎｉｑｕｅ  ｂａｓｅｄ  ｏｎ  ａ  ｇｅｎｅｔｉｃ  ａｌｇｏｒｉｔｈｍ  ｆｏｒ  ｉｎｖｅｒｓｅ  ｐｒｏｂｌｅｍｓ：  ｔｏｗａｒｄｓ  ｔｈｅ  ｉｎｆｅｒｅｎｃｅ  ｏｆ  ｇｅｎｅｔｉｃ  ｎｅｔｗｏｒｋｓ）」，（ドイツ），プロック・ジャーマン・コンファレンス・オン・バイオインフォマティックス，１９９９年，ｐ．１２７−１４０
【非特許文献２】ヴォイト・Ｅ・Ｏ、アルメイダ・Ｊ，「ダイナミック・プロファイリング・アンド・キャノニカル・モデリング：パワフル・パートナーズ・イン・メタボリック・パスウェイ・アイデンティフィケイション（Ｐｏｗｅｒｆｕｌ  Ｐａｒｔｎｅｒ  ｉｎ  Ｍｅｔａｂｏｌｉｃ  Ｐａｔｈｗａｙ  Ｉｄｅｎｔｉｆｉｃａｔｉｏｎ）」，（オランダ），クラワー・アカデミック・パブリッシング（Ｋｌｕｗｅｒ  Ａｃａｄｅｍｉｃ  Ｐｕｂｌｉｓｈｉｎｇ），２００３年
【非特許文献３】松原謙一、榊佳之監修「ゲノム情報生物学」中山書店，２０００年１０月３０日，ｐ．１６５−１８７
【発明の開示】
【発明が解決しようとする課題】
【０００５】
本発明は、上記した従来技術の欠点を除くためになされたものであって、その目的とするところは、遺伝子発現相互作用推定手法として、Ｓ−Ｓｙｓｔｅｍモデルを適用することにより、対象遺伝子数及び計算量の問題を改善することである。これにより、遺伝子ネットワークの研究の進展に役立つ。
【課題を解決するための手段】
【０００６】
本発明に係る遺伝子発現相互作用推定装置は、
分解系係数を含み分解速度を示す項と、生成系係数を含み生成速度を示す項とからなる遺伝子の速度式を前提として、遺伝子間の相互作用を表す遺伝子ネットワークを求める遺伝子発現相互作用推定装置であって、以下の要素を有することを特徴とする
（１）野生株及び破壊株のタイム・コースデータを入力する野生株・破壊株タイム・コースデータ入力処理部
（２）野生株及び破壊株のタイム・コースデータについて、タイム・コース値の対数値を基に、スプライン補間係数を算出し、算出したスプライン補間係数から、タイム・コース値及び対数値データの微係数を算出する前処理部
（３）破壊株の前記タイム・コース値を対数空間に変換した値、及び破壊株の前記対数値データの微係数を使用して、最小二乗法により、分解系係数の最適値を算出する分解系係数推定部
（４）算出した分解系係数の最適値を前提とし、非破壊株データの前記タイム・コース値を対数空間に変換した値、及び非破壊株データの前記対数値データの微係数を使用して、最小二乗法により生成系係数の最適値を算出する生成系係数推定部。
【０００７】
生成系係数推定部は、非破壊株データの前記対数値データの微係数から生成速度プロファイルを算出し、算出した生成速度プロファイルを利用して、最小二乗法により生成系係数の最適値を算出することを特徴とする。
【０００８】
遺伝子の数を、ｎとし、ｉ番目遺伝子の発現量を、Ｘｉとし、ｉ番目遺伝子の生成系係数を、αｉ及びｇｉｊとし、ｉ番目遺伝子の分解系係数を、βｉ及びｈｉｊとした場合に、前記遺伝子の速度式は、以下の式で示されることを特徴とする
【０００９】
【数２】
0004611656.tif 000002
遺伝子発現相互作用推定装置は、更に、算出した生成系係数及び分解系係数に基づいて、遺伝子発現についての遺伝子間の相互作用を表すパスウェイを描画するパスウェイ描画部を有することを特徴とする。
【００１０】
本発明に係る遺伝子発現相互作用推定方法は、
分解系係数を含み分解速度を示す項と、生成系係数を含み生成速度を示す項とからなる遺伝子の速度式を前提として、遺伝子間の相互作用を表す遺伝子ネットワークを求める遺伝子発現相互作用推定装置による方法であって、以下の要素を有することを特徴とする
（１）野生株及び破壊株のタイム・コースデータを入力する工程
（２）野生株及び破壊株のタイム・コースデータについて、タイム・コース値の対数値を基に、スプライン補間係数を算出し、算出したスプライン補間係数から、タイム・コース値及び対数値データの微係数を算出する工程
（３）破壊株の前記タイム・コース値を対数空間に変換した値、及び破壊株の前記対数値データの微係数を使用して、最小二乗法により、分解系係数の最適値を算出する工程
（４）算出した分解系係数の最適値を前提とし、非破壊株データの前記タイム・コース値を対数空間に変換した値、及び非破壊株データの前記対数値データの微係数を使用して、最小二乗法により生成系係数の最適値を算出する工程。
【００１１】
本発明に係るプログラムを記録したコンピュータ読み取り可能な記録媒体は、
分解系係数を含み分解速度を示す項と、生成系係数を含み生成速度を示す項とからなる遺伝子の速度式を前提として、遺伝子間の相互作用を表す遺伝子ネットワークを求める遺伝子発現相互作用推定装置となるコンピュータに、以下の処理を実行させるためのプログラムを記録したコンピュータ読み取り可能な記録媒体であることを特徴とする
（１）野生株及び破壊株のタイム・コースデータを入力する処理
（２）野生株及び破壊株のタイム・コースデータについて、タイム・コース値の対数値を基に、スプライン補間係数を算出し、算出したスプライン補間係数から、タイム・コース値及び対数値データの微係数を算出する処理
（３）破壊株の前記タイム・コース値を対数空間に変換した値、及び破壊株の前記対数値データの微係数を使用して、最小二乗法により、分解系係数の最適値を算出する処理
（４）算出した分解系係数の最適値を前提とし、非破壊株データの前記タイム・コース値を対数空間に変換した値、及び非破壊株データの前記対数値データの微係数を使用して、最小二乗法により生成系係数の最適値を算出する処理。
【００１２】
本発明に係るプログラムは、
分解系係数を含み分解速度を示す項と、生成系係数を含み生成速度を示す項とからなる遺伝子の速度式を前提として、遺伝子間の相互作用を表す遺伝子ネットワークを求める遺伝子発現相互作用推定装置となるコンピュータに、以下の手順を実行させるためのプログラムであることを特徴とする
（１）野生株及び破壊株のタイム・コースデータを入力する手順
（２）野生株及び破壊株のタイム・コースデータについて、タイム・コース値の対数値を基に、スプライン補間係数を算出し、算出したスプライン補間係数から、タイム・コース値及び対数値データの微係数を算出する手順
（３）破壊株の前記タイム・コース値を対数空間に変換した値、及び破壊株の前記対数値データの微係数を使用して、最小二乗法により、分解系係数の最適値を算出する手順
（４）算出した分解系係数の最適値を前提とし、非破壊株データの前記タイム・コース値を対数空間に変換した値、及び非破壊株データの前記対数値データの微係数を使用して、最小二乗法により生成系係数の最適値を算出する手順。
【発明の効果】
【００１３】
本発明においては、ＧＡやＳＡのような発見的な最適化手法の場合よりも多くの遺伝子を、より早く、より小さな計算機資源で推定する。
【発明を実施するための最良の形態】
【００１４】
実施の形態１．
以下本発明を図面に示す実施例に基づいて説明する。遺伝子発現相互作用を推定するためには、遺伝子毎に分解系及び生成系を解析する必要がある。
【００１５】
本発明では、分解系及び生成系を同時に動作させることができるＳ−Ｓｙｓｔｅｍモデルによる推定を行う。本手法におけるＳ−Ｓｙｓｔｅｍモデルの微分方程式を次式のように表す。
【００１６】
【数３】
0004611656.tif 000003
ｉ＝１，２，…，ｎ
但し、
Ｘｉ：ｉ番目遺伝子の発現量
αｉ，ｇｉｊ：ｉ番目遺伝子の生成系係数
βｉ，ｈｉｊ：ｉ番目遺伝子の分解系係数
処理の概要について説明する。本装置及び方法では、破壊株のタイムコースデータを利用し、その対数空間に変換した値及び補間データから得た微係数を使用して最小二乗法により分解系係数βｉ，ｈｉｊの最適値を算出する。
【００１７】
また、その分解系を前提とし、非破壊株データのタイムコースデータを利用し、その対数空間に変換した値及び補間データから生成速度プロファイルを算出し、それを利用して最小二乗法により生成系係数αｉ，ｇｉｊの最適値を算出する。
【００１８】
装置のモジュール構成を示す。図１は、遺伝子発現相互作用推定装置のモジュール構成を示す図である。遺伝子発現相互作用推定装置は、入力装置１、野生株・破壊株タイムコース・データ入力処理部２、前処理部３、分解系係数推定部４、生成系係数推定部５、解析データ・手法記憶部６、パスウェイ描画部７、及び表示装置８を有している。また、前処理部３は、スプライン補間処理部３１、タイムコースデータ算出部３２、及び対数値データの微係数算出部３３を備えている。
【００１９】
次に処理について説明する。図２は、全体処理のフローを示す図である。まず、入力装置１を介して、野生株のタイムコースデータ及び破壊株のタイムコースデータを入力し（Ｓ２０１）、前処理部３により前処理（Ｓ２０２）を行う。本処理では、野生株及び破壊株のタイムコースデータに対する補間等の処理を行う。
【００２０】
図３は、前処理のフローを示す図である。野生株及び破壊株のタイム・コースデータを読込み、スプライン補間処理部３１により、タイム・コース値の対数値を基にスプライン補間係数を算出する（Ｓ３０１）。そして、スプライン補間係数から、タイムコースデータ算出部３２によりタイム・コース値を算出し（Ｓ３０２）、更に、対数値データの微係数算出部３３により対数値データの微係数を求める（Ｓ３０３）。スプライン補間係数、タイム・コース値、及び対数値データの微係数は、それぞれ記憶領域に記憶される。
タイム・コース値：
【００２１】
【数４】
0004611656.tif 000004
対数値データの微係数：
【００２２】
【数５】
0004611656.tif 000005
次に、分解系係数推定部４により分解系係数の推定処理（Ｓ２０３）を行う。本処理では、式１の右辺の第二項に示す分解系係数の推定を行う。図４は、分解系係数の推定処理のフローを示す図である。Ｓ４０２〜Ｓ４０４の処理を、全遺伝子に対して処理し、全遺伝子の分解系係数を推定する。
【００２３】
その為に、すべての遺伝子について以下の処理を繰り返す（Ｓ４０１）。遺伝子ｉの破壊株（株番号ｋ：株番号は、カッコつきの上付き添え字で表示する。）のデータに対して、Ｓ３０２で算出したタイム・コース値と、Ｓ３０３で算出した対数値データの微係数から、ベクトルｂと行列Ａを求める（Ｓ４０２）。ベクトルｂと行列Ａは、記憶領域に記憶される。
【００２４】
ベクトルｂは、以下に示す成分からなる。
【００２５】
【数６】
0004611656.tif 000006
行列Ａは、以下に示す成分からなる。
【００２６】
【数７】
0004611656.tif 000007
そして、算出したＡ行列とｂベクトルを使用して、最小二乗法によりベクトルｃを推定する（Ｓ４０３）。このとき、下記の関係にあることを前提としている。
ｂベクトル＝Ａ行列×ｃベクトル
ｃベクトルは、以下の成分から構成される。
【００２７】
【数８】
0004611656.tif 000008
但し、ｈ´ｉｊ＝ｈｉｊ−δｉｊである（δｉｊは、クロネッカーのデルタである。）。従って、各成分から分解系係数βｉ，ｈｉｊを算出できる（Ｓ４０４）。ベクトルｃ、分解系係数βｉ，ｈｉｊは、記憶領域に記憶される。
【００２８】
そして、すべての遺伝子について処理した時点で終了する（Ｓ４０５）。
【００２９】
続いて、生成系係数推定部５により生成系係数の推定処理（Ｓ２０４）を行う。本処理では、本処理では、式１の右辺の第一項に示す生成系係数の推定を行う。図５は、生成系係数の推定処理のフローを示す図である。全遺伝子について、非破壊株（株番号ｋ＿１〜ｋ＿ｎ）のデータに対してＳ５０２〜Ｓ５０４の処理を行い、全遺伝子の生成系係数を推定する。
【００３０】
その為に、全遺伝子について以下の処理を繰り返す（Ｓ５０１）。Ｓ３０２で算出したタイム・コース値と、Ｓ３０３で算出した対数値データの微係数から、ベクトルｂと行列Ａを求める（Ｓ５０２）。ベクトルｂと行列Ａは、記憶領域に記憶される。
【００３１】
ベクトルｂは、以下に示す成分からなる。
【００３２】
【数９】
0004611656.tif 000009
行列Ａは、以下に示す成分からなる。
【００３３】
【数１０】
0004611656.tif 000010
但し
【００３４】
【数１１】
0004611656.tif 000011
算出したＡ行列とｂベクトルを使用して、最小二乗法によりベクトルｃを推定する（Ｓ５０３）。このとき、下記の関係にあることを前提としている。
ｂベクトル＝Ａ行列×ｃベクトル
ｃベクトルは、以下の成分から構成される。
【００３５】
【数１２】
0004611656.tif 000012
但し、ｇ´ｉｊ＝ｇｉｊ−δｉｊである。従って、各成分から生成系係数αｉ，ｇｉｊを算出できる（Ｓ４０４）。ベクトルｃ、生成系係数αｉ，ｇｉｊは、記憶領域に記憶される。
【００３６】
そして、すべての遺伝子について処理した時点で終了する（Ｓ５０５）。
【００３７】
最後に、パスウェイ描画部７によりパスウェイ描画処理（Ｓ２０５）を行う。本処理では、算出した生成系係数及び分解系係数に基づいて、遺伝子発現についての遺伝子間の相互作用を表すパスウェイを描画し、表示装置に表示させる。
【００３８】
遺伝子発現相互作用推定装置は、コンピュータであり、各要素はプログラムにより処理を実行することができる。また、プログラムを記憶媒体に記憶させ、記憶媒体からコンピュータに読み取られるようにすることができる。
【００３９】
図６は、遺伝子発現相互作用推定装置のハードウェア構成例を示す図である。バスに、演算装置６０１、データ記憶装置６０２、メモリ６０３が接続されている。データ記憶装置６０２は、例えばＲＯＭ（Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）やハードディスクである。メモリ６０３は、通常ＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）である。前述した記憶領域は、データ記憶装置６０２の記憶回路や記憶媒体に形成されるものである。
【００４０】
プログラムは、通常データ記憶装置６０２に記憶されており、メモリ６０３にロードされた状態で、順次演算装置６０１に読み込まれ処理を行う。
【００４１】
また、入力装置６０４と出力装置６０５にも接続している。入力装置６０４は、通信媒体や記憶媒体からデータを入力し、出力装置６０５は、データの表示、印字、または送信等の処理を行う。
【図面の簡単な説明】
【００４２】
【図１】遺伝子発現相互作用推定装置のモジュール構成を示す図である。
【図２】全体処理のフローを示す図である。
【図３】前処理のフローを示す図である。
【図４】分解系係数の推定処理のフローを示す図である。
【図５】生成系係数の推定処理のフローを示す図である。
【図６】遺伝子発現相互作用推定装置のハードウェア構成例を示す図である。
【符号の説明】
【００４３】
１  入力装置、２  野生株・破壊株タイムコース・データ入力処理部、３  前処理部、４  分解系係数推定部、５  生成系係数推定部、６  解析データ・手法記憶部、７  パスウェイ描画部、８  表示装置、３１  スプライン補間処理部、３２  タイムコースデータ算出部、３３  対数値データの微係数算出部。
【図１】
0004611656.tif 000013
【図２】
0004611656.tif 000014
【図３】
0004611656.tif 000015
【図４】
0004611656.tif 000016
【図５】
0004611656.tif 000017
【図６】
0004611656.tif 000018