(57)【特許請求の範囲】
【請求項１】
生物学的データの複数の分子署名のセットの分析に基づき、複数のフィーチャを有する生物学的データの新しい署名を構成する、コンピュータにより実行される方法であって、
前記コンピュータのプロセッサが、Ｎ回、複数の実行の各実行において、生物学的データにパターン発見アルゴリズムを実行して、前記生物学的データの複数の署名のセットを求める段階と、
前記プロセッサが、サイズ（Ｆ×Ｆ）の第１の行列であって、前記第１の行列の各行列要素が実行ｒｉで生物学的データの前記分子署名における共起性に基づき、対応する各フィーチャペアを定量化する第１の行列を構成する段階と、
前記プロセッサが、サイズ（Ｆ×Ｎ）の第２の行列であって、前記第２の行列の各要素が対応する第１の行列の対応する行の和である第２の行列を構成する段階と、
前記プロセッサが、サイズ（Ｆ）のベクトルであって、前記ベクトルの各要素が前記第２の行列のすべての列の対応する行の平均であるベクトルを構成する段階と、
ベクトルの各値を降順にソートして、降順のベクトルの各要素により表される各フィーチャをランキングする段階と、
前記プロセッサが、行列のフィーチャペアの相関の計算に基づき、ランキングと相関に応じて分子署名の少なくとも２つのフィーチャを組み合わせて、ランク付けしたフィーチャに基づき、新しい分子署名のセットを構成する段階とを有する、方法。
【請求項２】
前記分子署名は、ヌクレオチドシーケンス、遺伝子変異、メチル化状態、または遺伝子発現よりなるグループから選択される、請求項１に記載の方法。
【請求項３】
生物学的データの新しい分子署名を構成する装置であって、
プロセッサを有し、前記プロセッサは、
生物学的データに対してＮ回パターン発見アルゴリズムを実行して、前記生物学的データの複数の署名のセットを求め、
複数の実行の各実行において、
サイズ（Ｆ×Ｆ）の第１の行列であって、前記第１の行列の各行列要素が実行ｒｉで生物学的データの分子署名における共起性に基づき、対応する各フィーチャペアを定量化する第１の行列を構成し、
サイズ（Ｆ×Ｎ）の第２の行列であって、前記第２の行列の各要素が対応する第１の行列の対応する行の和である第２の行列を構成し、
サイズ（Ｆ）のベクトルであって、前記ベクトルの各要素が前記第２の行列のすべての列の対応する行の平均であるベクトルを構成し、、
ベクトルの各値を降順にソートして、降順のベクトルの各要素により表される各フィーチャをランキングし、
行列のフィーチャペアの相関の計算に基づいて、ランキングと相関に応じて分子署名の少なくとも２つのフィーチャを組み合わせて、ランク付けしたフィーチャに基づき、新しい分子署名のセットを構成する、
装置。
【請求項４】
プロセッサに、
生物学的データに対してＮ回パターン発見アルゴリズムを実行して、前記生物学的データの複数の分子署名のセットを求める段階と、
複数の実行の各実行において、
サイズ（Ｆ×Ｆ）の第１の行列であって、前記第１の行列の各行列要素が実行ｒｉで生物学的データの分子署名における共起性に基づき、対応する各フィーチャペアを定量化する第１の行列を構成する段階と、
サイズ（Ｆ×Ｎ）の第２の行列であって、前記第２の行列の各要素が対応する第１の行列の対応する行の和である第２の行列を構成する段階と、
サイズ（Ｆ）のベクトルであって、前記ベクトルの各要素が前記第２の行列のすべての列の対応する行の平均であるベクトルを構成する段階と、
ベクトルの各値を降順にソートして、降順のベクトルの各要素により表される各フィーチャをランキングする段階と、
行列のフィーチャペアの相関の計算に基づき、ランニングと相関に応じて分子署名の少なくとも２つのフィーチャを組み合わせて、ランク付けしたフィーチャに基づき、新しい分子署名のセットを構成する段階とを実行させる、コンピュータプログラム。
【請求項５】
前記プロセッサに、
コンピュータ処理機能を有する装置により実行されたとき、請求項１または２に記載の方法のすべての段階を実行する段階をさらに実行させる、請求項４に記載のコンピュータプログラム。
【請求項６】
臨床データを統計的に分析する、請求項１に記載の方法。
【請求項７】
複数の分子署名に基づき検査を設計する、請求項１に記載の方法。
【請求項８】
複数の分子署名に基づき検査結果を解釈する、請求項１に記載の方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、概してバイオインフォマティックスの分野に関する。より具体的に、本発明は複数の生物学的データ署名のランキングに基づき、新しい生物学的データ署名を構成する方法に関する。また、本発明は、複数の分子署名（molecular  signature）を比較する装置、新しい署名を構成する装置、コンピュータプログラム製品、及びかかる方法と装置の使用に関する。
【背景技術】
【０００２】
生物学的サンプルの高スループット分子プロファイリングにより、サンプルをカテゴリに階層化するのに用いられる署名が得られることが知られている。かかるグループは、状態の診断（正常か、病気か）、疾病経過（悪性か、非悪性か）、治療選択（薬剤応答）などである。署名の例は分子署名である。
【０００３】
分子署名などの署名は、一般的に、臨床研究により発見され実証されるが、臨床における同じ疑問に答えるには、かかる署名のコンテンツには大きなばらつきがあることが分かってきた。これは主に、分析に利用できるデータの性質によるものである。高スループット測定により何万ないし何十万というフィーチャ（features）が得られるが、測定点（サンプル）の数は数百である。また、これは、フィーチャ（features）の数とサンプルの数の間に何桁も差がある次元の呪いとも呼ばれる。さらに、同じ臨床研究や同じ発見方法のコンテキストであれば、パラメータの選択にばらつきがあり、または異なるランダム化をした場合には、まったく同じプロシージャを繰り返しても同様だがそれでもなお異なる署名が得られることが多い。
【０００４】
本技術分野の発展は速く、患者をカテゴリに階層化する多くの署名が使われている。標準が無いため、このデータは一貫性のないバラバラなものであることが多い。また、生物学的システムの複雑性とその研究方法により、出所が異なるゲノムのアイデンティティのセットの比較は困難になる。
【０００５】
そのため、生物学的データの分析方法を改善すれば有利であり、柔軟性、コストエフェクティブネス、速さ、及び／または分析精度が高い方法が有利である。
【発明の概要】
【発明が解決しようとする課題】
【０００６】
従って、本発明は、好ましくは、単独でまたは組み合わされて、本技術分野における上記の１つ以上の欠点を緩和もしくは解消するものである。この目的は、添付した独立請求項による生物学的データの複数の署名をグループ化する方法、装置、コンピュータプログラム製品、及び使用を提供することにより、達成される。
【０００７】
本発明によるソリューションは、概して、同じ予測値の複数の署名の測定を可能とすることである。
【課題を解決するための手段】
【０００８】
本発明の一態様において、生物学的データの複数の署名のセットの分析に基づき、複数のフィーチャ
【０００９】
［外１］
0005714001.tif 000002
を有する生物学的データの新しい署名を構成する方法を提供する。該方法は、生物学的データに対してＮ回フィーチャサブセット選択アルゴリズムを実行して、生物学的データの複数の署名のセットを求める段階を有する。該方法は、複数の実行
【００１０】
［外２］
0005714001.tif 000003
のうちの各実行ｒｉについて、サイズ（Ｆ×Ｆ）の第１の行列（Ｍｉ）を構成する段階を有する。第１の行列（Ｍｉ）の各行列要素（ｉ，ｊ）は、実行ｒｉにおける分析する署名の共起性に基づき、対応する各フィーチャペア（ＦｉＦｊ）を定量化する。該方法は、サイズ（Ｆ×Ｎ）の第２の行列（Ｔ）を構成する段階も有する。第２の行列（Ｔ）の各要素（Ｔ（ｉ，ｊ））には、対応する第１の行列（Ｍｉ）の対応する行（Ｆｉ）の和が入っている。さらに、該方法は、サイズ（Ｆ）のベクトル（Ｃ）であって、ベクトル（Ｃ）の各要素（ＣＦｉ）が第２の行列（Ｔ）のすべての列（１ないしＮ）の対応する行（Ｔｉ）の平均であるベクトルを構成する段階を有する。該方法はさらに、ベクトル（Ｃ）の各値を降順にソートして、降順のベクトル（Ｃ）の各要素により表される各フィーチャ（Ｆｉ）をランキングする段階を有する。該方法は、ランキングが同じ署名の少なくとも２つのフィーチャを組み合わせて、ランク付けしたフィーチャに基づき、新しい署名のセットを構成する段階も含む。
【００１１】
本発明の第２の態様によると、互いに動作可能に接続されているとき、本発明の第１の態様によるステップを実行するように構成されたユニットを有する装置を提供する。
【００１２】
本発明の第３の態様によると、コンピュータにより処理されるコンピュータプログラムを含むコンピュータプログラム製品を提供する。このコンピュータプログラムは、本発明の第１の態様による方法を実行するコードセグメントを有する。
【００１３】
本発明の第４の態様によると、臨床データを統計的に分析する、第１の態様による方法または第２の態様による装置の使用を提供する。
【００１４】
本発明の第５の態様によると、複数の分子署名に基づき検査を設計する、第１の態様による方法または第２の態様による装置の使用を提供する。
【００１５】
本発明の第６の態様によると、複数の分子署名に基づき検査結果を解釈する、第１の態様による方法または第２の態様による装置の使用を提供する。
【００１６】
本発明の実施形態は従属項に規定した。
【００１７】
本発明は、署名発見プロセスからのより包括的な出力を得られるという点で、従来技術に対する優位性を有している。単一の署名出力の替わりに、一組の測定値の複数の臨床的見方を可能にする、より広い記述トレンドをキャプチャする。この記述は様々な方法で利用できる。臨床研究では、複数の署名により、調べている臨床的疑問の基礎にある生命現象の評価をよりよくすることができる。これを用いて、例えば、複数の薬剤ターゲットを探し、疾病メカニズムを評価することができる。さらにまた、診断検査の設計において、本発明により、予測値が同じ複数の署名が共存できる。これにより、例えば、ゲノムフィーチャ（genomic  features）の最適なセットの選択を可能とし、検査設計における問題を回避できる。臨床医は患者データを異なった見方で見て、ゲノムフィーチャの多くの組み合わせを利用し、生理的サンプルの生物学的状態と患者の臨床的状態を最もよく記述する最適なフィーチャを選択できる。
【００１８】
さらなる優位性は、本発明により患者データを異なる角度から見て、多くのゲノムフィーチャの組み合わせを利用し、遺伝子情報から結論を引き出せる最適なものを選択できることである。
【図面の簡単な説明】
【００１９】
本発明による上記その他の態様、特徴、利点は、添付した図面を参照して、本発明の実施形態に関する以下の説明により明らかになり、説明される。
【図１】一実施形態による方法を示すフローチャートである。
【図２】一実施形態によるフィーチャサブセット生成を示す図である。
【図３】一実施形態による装置を示すフローチャートである。
【図４】一実施形態によるコンピュータプログラム製品を示す図である。
【発明を実施するための形態】
【００２０】
本技術分野の当業者が本発明を実施できるように、添付した図面を参照しつつ本発明の実施形態を詳細に説明する。しかし、本発明は、多数の異なる形式で実施でき、ここに開示した実施形態だけに限定されると解釈してはならない。むしろ、これらの実施形態は、この開示が完全であり、当業者に本発明の範囲を完全に伝えるように提供するものである。実施形態は本発明を限定するものではなく、本発明は添付した特許請求の範囲によってのみ限定される。さらに、添付した図面に示した実施形態の詳細な説明で用いる用語は、本発明を限定するものではない。
【００２１】
図１の一実施形態において、複数の生物学的データの署名のセットの分析に基づき、複数のフィーチャ
【００２２】
［外３］
0005714001.tif 000004
を有する生物学的データの新しい署名を構成する方法１０を示す。該方法は、ステップ１１０において、生物学的データに対してＮ回フィーチャサブセット選択アルゴリズムを実行して、生物学的データの複数の署名のセットを求める段階を有する。複数の実行
【００２３】
［外４］
0005714001.tif 000005
のうちの各実行ｒｉについて、方法１０は、サイズ（Ｆ×Ｆ）の第１の行列
【００２４】
［外５］
0005714001.tif 000006
を構成する段階１２０ａを有する。第１の行列（Ｍｉ）の各行列要素（ｉ，ｊ）は、実行ｒｉにおける分析する署名の共起性に基づき、対応する各フィーチャペア（ＦｉＦｊ）を定量化する。方法１０は、サイズ（Ｆ×Ｎ）の第２の行列
【００２５】
［外６］
0005714001.tif 000007
を構成する段階１２０ｂも有する。第２の行列（Ｔ）の各要素（Ｔ（ｉ，ｊ））には、対応する第１の行列（Ｍｉ）の対応する行（Ｆｉ）の和が入っている。方法１０は、さらにステップ１３０において、サイズ（Ｆ）のベクトル
【００２６】
［外７］
0005714001.tif 000008
を構成する段階を有する。このベクトルの各要素（ＣＦｉ）には、第２の行列（Ｔ）のすべての列（１ないしＮ）の対応する行（Ｔｉ）の平均が入っている。この方法は、ベクトル（Ｃ）の各値を降順にソートして、降順のベクトル（Ｃ）の各要素により表される各フィーチャ（Ｆｉ）をランキングする段階１４０を有する。さらに、方法１０は、同じランキングの署名の少なくとも２つのフィーチャを組み合わせて、ランク付けしたフィーチャに基づき、新しい署名のセットを構成する段階１５０を有する。
【００２７】
署名を、測定した分子フィーチャ（molecular  features）のランキング等の特徴付けと、フィーチャの共起の発生などのトレンドにより説明する。これは、どのフィーチャをどんな順序で組み合わせるか判断する役に立つ。また、分子署名のファミリを構成（インスタンス化）できる。これらは、診断分析を設計し、ゲノムフィーチャのセットを変えて分子測定を行うのに用いられる。各実行において平均フィーチャサブセットサイズや（サンプルの分類などにおける）平均性能などの追加的特性を保存し、考慮してもよい。
【００２８】
一実施形態では、かかる署名は、ヌクレオチドシーケンス、ジェネティックバリエーション、メチル化状態、またはジェネティックエクスプレッションなどの分子署名でもよい。しかし、どんな署名を利用してもよい。
【００２９】
一実施形態では、前記フィーチャサブセット選択は、フィルタリング法、前方フィーチャ選択法、ジェネティックアルゴリズムなどを含み、これらは本技術分野で周知である。しかし、どんなフィーチャサブセット選択を利用してもよい。
【００３０】
【表１】
0005714001.tif 000009
【００３１】
プログラミング言語「Ｒ」に基づく、より具体的な実施形態では、表１の入力を用いて、以下のステップを実行する。
【表２】
0005714001.tif 000010
【００３２】
上記の擬似コードによると、Ｔに基づくフィーチャの「近傍」に基づき、最も相関が高いフィーチャと最も逆相関が高いフィーチャを同様の方法で求める。これを行う１つの方法は、Ｎ個の最も相関（または逆相関）が高いフィーチャを探して、Ｃを用いて、そのＮ個のフィーチャから最もランクの高いものを取ることである。この方法を用いると、相関（逆相関）しているが場合の数が少ない、関連が低いフィーチャを避けることができる。
【００３３】
図２を参照して、さらにフィーチャサブセット生成を説明する。
【００３４】
大文字の表記はサブセットに加えられるフィーチャを示す。その各フィーチャについて、逆相関する遺伝子を示す小文字のノードをマークした。これらのうち、１つのノードに大文字と小文字の表記がある場合のみに興味がある。点線で囲んだ領域は最も逆相関が高いと考えられるフィーチャである。相関が高く足し合わせるフィーチャは、図示を簡明にするため示していない。
【００３５】
レベル１では、サブセット＝｛｛Ａ｝｝であり、ＢとＣが加えられ（features.added.per.level=3）、Ｄ、Ｅ、Ｆは次のレベルで加えるようマークされる。４つの（逆）相関フィーチャの近傍を考える。レベル１の後、Ｓ＝｛｛Ａ，Ｂ，Ｃ｝｝である。今のところ逆相関コンフリクトがないからである。
【００３６】
レベル２において、Ｄをサブセットに加え、次のレベルでＧを加えることを考える。サブセット＝｛｛Ａ，Ｂ，Ｃ，Ｄ｝｝。次に、Ｅをサブセットに加え、次のレベルでＨを加えることを考える。サブセット＝｛｛Ａ，Ｂ，Ｃ，Ｄ，Ｅ｝｝。Ｆをサブセットに加え、次のレベルでＩを加えることを考える。しかし、ＥとＦは一緒には使えないので、サブセット＝｛｛Ａ，Ｂ，Ｃ，Ｄ，Ｅ｝，｛Ａ，Ｂ，Ｃ，Ｄ，Ｆ｝｝。
【００３７】
レベル３では、Ｇをサブセットに加え、サブセット={{A,B,C,D,E,G},  {A,B,C,D,F,G}}。Ｈをサブセットに加え、サブセット={{A,B,C,D,E,G,H},  {A,B,C,D,F,G,H}}。Ｉをサブセットに加えるが、ＩとＨは一緒には使えないので、サブセット={{A,B,C,D,E,G,H},  {A,B,C,D,F,G,H},  {A,B,C,D,E,G,I},  {A,B,C,D,F,G,I}}。
【００３８】
上記の例は、フィーチャ間の相関や逆相関（anti-correlations）など、ペアごとのトレンド値をどう使って、新しい署名を求めるかを示している。
【００３９】
一実施形態では、Ｎが４の場合、フィーチャサブセット選択アルゴリズムは、本技術分野で知られた任意のフィーチャサブセットアルゴリズムでもよいが、生物学的データの複数の署名のセットを求めるため、４回実行される（１１０）。表２に示すように、これらの署名に基づき、第１の行列（Ｍｉ）を構成する（１２０ａ）。
【００４０】
【表３】
0005714001.tif 000011
【００４１】
次に、表３に示すように、第２の行列（Ｔ）を構成する（１２０ｂ）。
【００４２】
【表４】
0005714001.tif 000012
【００４３】
第２行列（Ｔ）に基づき、第２行列（Ｔ）のすべての列（１ないしＮ）の対応する行（Ｔｉ）の平均を割り当てることにより、Ｃベクトルを構成する（１３０）。表４に示すように、Ｃベクトルも降順にソートする（１４０）。
【００４４】
【表５】
0005714001.tif 000013
【００４５】
次に、ランク付けしたフィーチャに基づき、新しい署名のセットを構成する（１５０）。この実施形態では、Ａがトップである。そこで、何をするか決めるため、ランキングが同様のフィーチャ、すなわち相関フィーチャが必要であり、関連のないフィーチャ、すなわち逆相関フィーチャも必要である。これはＴ行列に基づきフィーチャペアの相関を計算することにより求められる。結果を表５に示す。
【００４６】
【表６】
0005714001.tif 000014
【００４７】
フィーチャ｛Ａ｝から始めて、次のトップランクのフィーチャを加える。｛Ａ，Ｂ｝となる。次のトップランクのフィーチャ（Ｃ）は、ＢとＣとの逆相関が高いので、加えない。その替わり２つのセットを生成する｛｛Ａ，Ｂ｝，｛Ａ，Ｃ｝｝。以下、ＤとＥは逆相関が高いので、一緒にしておくことはできない。このように、この実施形態では、構成される新しい署名のセットは（１５０）、｛｛Ａ，Ｂ，Ｄ｝，｛Ａ，Ｃ，Ｅ｝｝である。
【００４８】
本技術分野で周知の、フィーチャサブセットアルゴリズムの各実行の追加的特性を、保存して考慮してもよい。一実施形態では、平均フィーチャサブセットサイズ、サンプルの分類の平均性能などを考慮してもよい。例えば、フィーチャサブセットにおいて、フィーチャの平均数などの統計を利用してもよい。一サブセットは２０個のフィーチャを有し、他の一サブセットは３５個のフィーチャを有してもよい。さらにまた、結果の予測力に関する各フィーチャサブセットの性能は異なってもよい。例えば、あるサブセットは９８％の予測精度を有し、他のサブセットは９５％の精度を有する。このように、フィーチャの特徴Ｃを計算しつつ、各サブセットの貢献をその性能に関して重み付けしてもよい。それゆえ、フィーチャのペアが署名において共起するが、予測性能が低いとき、そのペアごとの共起スコアを、署名の平均予測精度に反比例する量だけペナルティをかけて、これを調整することができる。例えば、共起ごとに１を加えるのではなく、性能がよくないサブセットには０．５を足してもよい。
【００４９】
臨床における一実施形態において、臨床医は患者データを異なった見方で見て、ゲノムフィーチャの多くの組み合わせを利用し、生理的サンプルの生物学的状態と患者の臨床的状態を最もよく記述する最適なフィーチャを選択する。例えば、分析により患者サンプルの１００個のゲノムフィーチャを測定するとき、複数の測定のセットで「ワンショット」を取るのではなく、診断テストは単一の研究に基づく複数の評価よりなる。そのため、単一のテストではなく、テストの結果を集合的に決定する署名の代表（committee）を用いてもよい。例えば、臨床研究に基づく場合、５つのゲノムフィーチャＡ、Ｂ、Ｃ、ＤおよびＥに基づくトレンドは３つのコンビネーションＡ、Ｂ、Ｃ；Ａ、Ｃ、Ｄ；Ａ、Ｃ、Ｅによりキャプチャできる。ノイズが大きいフィーチャＢの測定結果によりＡ、Ｂ、Ｃからの出力は歪んでしまうだろう。しかし、残る２つの新しい署名によりサンプルのシグナルを正しくキャッチできる。
【００５０】
図３に示した本発明の他の一実施形態では、生物学的データの新しい署名を構成する装置３０を示す。この装置３０は、生物学的データに対してＮ回フィーチャサブセット選択アルゴリズムを実行して、生物学的データの複数の署名のセットを求めるように構成された第１のユニット３１０を有する。
【００５１】
装置３０は、複数の実行
【００５２】
［外８］
0005714001.tif 000015
のうちの各実行ｒｉについて、サイズ（Ｆ×Ｆ）の第１の行列
【００５３】
［外９］
0005714001.tif 000016
を構成する（３２０ａ）ように構成された第２のユニット３２０を有する。第１の行列（Ｍｉ）の各行列要素（ｉ，ｊ）は、実行ｒｉにおいて分析する署名の共起性に基づき、対応する各フィーチャペア（ＦｉＦｊ）を定量化する。第２のユニット３２０は、サイズ（Ｆ×Ｎ）の第２の行列
【００５４】
［外１０］
0005714001.tif 000017
を構成する（３２０ｂ）ようにも構成されている。第２の行列（Ｔ）の各要素（Ｔ（ｉ，ｊ））には、対応する第１の行列（Ｍｉ）の対応する行（Ｆｉ）の和が入っている。
【００５５】
装置３０はさらに、サイズ（Ｆ）のベクトル
【００５６】
［外１１］
0005714001.tif 000018
を構成するように構成された第３のユニット３３０を有する。このベクトルの各要素（ＣＦｉ）には、第２の行列（Ｔ）のすべての列（１ないしＮ）の対応する行（Ｔｉ）の平均が入っている。
【００５７】
装置３０は、ベクトル（Ｃ）の各値を降順にソートして、降順のベクトル（Ｃ）の各要素により表される各フィーチャ（Ｆｉ）をランキングする第４のユニット３４０を有する。
【００５８】
装置３０はさらに、同じランキングの署名の少なくとも２つのフィーチャを組み合わせて、ランク付けしたフィーチャに基づき、新しい署名のセットを構成するように構成された第５のユニット３５０を有する。
【００５９】
ユニット３１０、３２０、３３０、３４０、３５０は互いに接続され動作している。ユニット３１０、３２０、３３０、３４０、３５０は、接続された、別々の物理的実体として実施してもよい。しかし、ユニット３１０、３２０、３３０、３４０、３５０は、単一の物理的実体として実施してもよい。ユニット３１０、３２０、３３０、３４０、３５０のどんな組み合わせも、異なる別々の物理的実体として実施してもよいし、単一の物理的実体として実施してもよい。この実体はさらにどのような設定で結合して、物理的実体間の接続を構成してもよい。
【００６０】
ユニット３１０、３２０、３３０、３４０、３５０は、メモリを有するプロセッサなどのハードウェアなどの、関連するタスクの実行に使える任意のユニットであり得る。
【００６１】
一実施形態では、装置３０またはシステム１００は、コンピュータ断層撮影（ＣＴ）システム、磁気共鳴画像化（ＭＲＩ）システム、または超音波画像化（ＵＳ）システムなどの医療用ワークステーションや医療用システムに含まれてもよい。
【００６２】
図４に示した一実施形態では、コンピュータプログラム製品は、コンピュータにより処理するコンピュータプログラム４０を有する。コンピュータプログラム４０は、Ｎ回、複数の実行
【００６３】
［外１２］
0005714001.tif 000019
の各実行において、生物学的データにフィーチャサブセット選択アルゴリズムを実行して、生物学的データの複数の署名のセットを取得する第１のコードセグメント４１０を有する。コンピュータプログラム４０は、複数の実行
【００６４】
［外１３］
0005714001.tif 000020
のうちの各実行（ｒｉ）について、サイズ（Ｆ×Ｆ）の第１の行列
【００６５】
［外１４］
0005714001.tif 000021
を構成する（４２０ａ）第２のコードセグメント４２０を有する。第１の行列（Ｍｉ）の各行列要素（ｉ，ｊ）は、実行ｒｉにおいて分析する署名の共起性に基づき、対応する各フィーチャペア（ＦｉＦｊ）を定量化する。第２のコードセグメント４２０は、サイズ（Ｆ×Ｎ）の第２の行列
【００６６】
［外１５］
0005714001.tif 000022
を構成する（４２０ｂ）。第２の行列（Ｔ）の各要素（Ｔ（ｉ，ｊ））には、対応する第１の行列（Ｍｉ）の対応する行（Ｆｉ）の和が入っている。
【００６７】
コンピュータプログラム４０はさらに、サイズ（Ｆ）のベクトル
【００６８】
［外１６］
0005714001.tif 000023
を構成するように構成された第３のコードセグメント４３０を有する。このベクトルの各要素（ＣＦｉ）には、第２の行列（Ｔ）のすべての列（１ないしＮ）の対応する行（Ｔｉ）の平均が入っている。
【００６９】
コンピュータプログラム４０は、ベクトル（Ｃ）の各値を降順にソートして、降順のベクトル（Ｃ）の各要素により表される各フィーチャ（Ｆｉ）をランキングする第４のコードセグメント４４０を有する。コンピュータプログラム４０はさらに、同じランキングの署名の少なくとも２つのフィーチャを組み合わせて、ランク付けしたフィーチャに基づき、新しい署名のセットを構成する第５のコードセグメント４５０を有する。
【００７０】
一実施形態では、コンピュータプログラム製品は、コンピュータにより処理されるコンピュータプログラム４０を有する。外コンピュータプログラム４０は、コンピュータ処理機能を有する装置により実行されると、方法１０で画定された方法ステップの全てを実行するように構成されている。
【００７１】
一実施形態では、方法１０または装置３０は、医師が患者を診断または治療する際に支援する、生物学的データの複数の分子署名などの情報を提供する。
【００７２】
このように、一実施形態では、臨床データを統計的に分析する方法１０または装置３０の使用を提供する。
【００７３】
例えば、１００人の患者のデータを含むデータセットにおいて、一部は癌の再発などの結果を生じ、一部は異なる結果を生じる。上記のように、複数の署名を発見し、特徴付ける。例えば、各患者の再発状態を予測する能力に基づき、署名を評価する。診断検査によりオリジナルな研究に基づく関連フィーチャをすべて測定する。新しい患者に対して、臨床医は、各々がその患者の再発を予測する複数の署名を生成し、その署名の特性に基づきそれぞれを評価する。例えば、臨床医は、特定の治療に関する生物学的経路に関連する遺伝子についての情報を含む署名のみを用いることができる。
【００７４】
一実施形態では、複数の分子署名に基づき検査を設計する、方法１０と装置３０の使用を提供する。
【００７５】
上記に基づき、生成したフィーチャの一部（smaller  subset）を含む検査（assays）を設計してもよい。この場合、フィーチャは、例えば、疾病に関連することが分かっている遺伝子や、その測定装置でよい特性が得られるフィーチャ（遺伝子）などの特性に基づいて選択する。
【００７６】
一実施形態では、複数の分子署名に基づき検査結果を解釈する、方法１０と装置３０の使用を提供する。
【００７７】
上記の検査結果により、臨床医は、その特徴を用いて、第１の検査結果からの出力を解釈する第２の検査を設計する。検査結果におけるフィーチャの信号の存否に基づき、臨床医は第１の検査では明示的に測られていないフィーチャについての追加情報を推論できる。例えば、分子署名においてフィーチャＡとＢが広く共起することが分かっており、第１の検査ではフィーチャＡの測定結果のみが得られている場合、臨床医は、フィーチャＢのありうる測定結果を推論し、患者に関する新しい臨床的に関連のある情報を得ることができる。
【００７８】
本発明は、ハードウェア、ソフトウェア、ファームウェアまたはこれらの組み合わせを含むいかなる好適な形式で実施することもできる。しかし、本発明は、１つ以上のデータプロセッサ及び／またはデジタル信号プロセッサ上で実行されるコンピュータソフトウェアとして実施することが好ましい。本発明の実施形態の構成要素は、いかなる好適な方法で物理的、機能的、論理的に実施してもよい。機能は単一のユニット、複数のユニット、または他の機能ユニットの一部として実施することもできる。このように、本発明は、単一ユニットで実施することもできるし、異なる複数のユニットやプロセッサに物理的かつ機能的に分散して実施することもできる。
【００７９】
好ましい実施形態を参照して本発明を説明したが、ここに記載した具体的な形式に限定することを意図したものではない。むしろ、本発明は添付した請求項によってのみ限定され、上記の実施形態以外の実施形態も添付した請求項の範囲に等しく入りうる。
【００８０】
請求項では、「有する」という用語は他の要素やステップの存在を排除するものではない。さらに、個別的に列挙されていても、複数の手段、要素、方法ステップは、例えば単一のユニットまたはプロセッサにより実施してもよい。また、個々の機能（feature）は異なる請求項に含まれていても、これらを有利に組み合わせることが可能であり、異なる請求項に含まれていても、機能を組み合わせられないとか、組み合わせても有利ではないということを示唆するものでもない。また、単数扱いをしても複数の場合を排除するものではない。よって、「１つの」、「第１の」、「第２の」等は複数の場合を排除するものではない。請求項中の参照符号は、明りょうにするために設けており、請求項の範囲を限定するものと解してはならない。
【図１】
0005714001.tif 000024
【図２】
0005714001.tif 000025
【図３】
0005714001.tif 000026
【図４】
0005714001.tif 000027