"生物学的試料の画像の特徴量から、その画像が示す薬効を判定するための解析を行う画像解析方法において、生物学的試料の画像の特徴量と、前記画像に対応する薬効の程度との組み合わせの入力を受け付ける受付ステップと、前記受付ステップにより入力された組み合わせに基づいて、複数の前記特徴量を用いた演算により前記薬効を算出するためのアルゴリズムを取得するアルゴリズム取得ステップと、生物学的試料の画像の特徴量を取得する特徴量取得ステップと、前記アルゴリズム取得ステップにより算出された前記アルゴリズムを用いて、前記特徴量取得ステップにより特徴量が取得された画像が示す薬効を算出する薬効算出ステップと、をコンピュータが実行することを特徴とする画像解析方法。前記アルゴリズムはニューラルネットワークまたはサポートベクターマシンにより定義されることを特徴とする請求項1に記載の画像解析方法。前記特徴量を抽出する特徴量抽出ステップをコンピュータが実行し、前記アルゴリズムは前記特徴量抽出ステップにおける画像解析のパラメタの指定を含むことを特徴とする請求項1または2に記載の画像解析方法。生物学的試料の画像の特徴量から、その画像が示す薬効を判定するための解析を行う画像解析装置において、生物学的試料の画像の特徴量と、前記画像に対応する薬効の程度との組み合わせの入力を受け付ける受付手段と、前記受付手段により入力された組み合わせに基づいて、複数の前記特徴量を用いた演算により前記薬効を算出するためのアルゴリズムを取得するアルゴリズム取得手段と、生物学的試料の画像の特徴量を取得する特徴量取得手段と、前記アルゴリズム取得手段により算出された前記アルゴリズムを用いて、前記特徴量取得手段により特徴量が取得された画像が示す薬効を算出する薬効算出手段と、を備えることを特徴とする画像解析装置。図5の例のように、細胞の面積のような自明な特徴量で薬効が定量化できる場合は薬効を容易に判定できるが、一般には、薬効の有無は、画像の明るさに関するパラメタや画像の形状に関するパラメタの組み合わせなど、異なる特徴量の複雑な組み合わせになると考えられ、直感的に数式で表すことも困難である。画像の明るさとしては、細胞領域内の画素値の合計、平均、標準偏差、あるいはそれらと空間2次元画像のXY方向、あるいは3次元画像のXYZ方向の量を組み合わせた各種特徴量統計量が考えられる。形状に関するパラメタは、形状のモーメントやオイラー数等の幾何学的な特徴量や、Haralickのテクスチャの特徴量が考えられる。このように、創薬スクリーニングの画像解析では、薬効を表すパラメタ特徴量を抽出することが必要である。特徴量は自明でなく、また1種類とは限らない。一般には、複数の特徴量を最適に組み合わせることで、1種類の特徴量のみに依存するよりも高い信号対雑音比で薬効を定量できると考えられる。しかしこれまで、どのような特徴量をどう組み合わせるかは、研究者の経験や勘という主観に委ねられてきたのが実情であり、パラメタの最適化を自動的に系統立てて行う手法は存在しない。また、HCSやHCAは、薬学分野の研究者にとって馴染みが少ない画像処理の専門知識が必要であり、敷居が高い技術という印象が持たれている。本発明の目的は、画像処理の専門知識を要することなく、適切に薬効を判定することができる画像解析方法等を提供することにある。本発明の画像解析方法は、生物学的試料の画像の特徴量から、その画像が示す薬効を判定するための解析を行う画像解析方法において、生物学的試料の画像の特徴量と、前記画像に対応する薬効の程度との組み合わせの入力を受け付ける受付ステップと、前記受付ステップにより入力された組み合わせに基づいて、複数の前記特徴量を用いた演算により前記薬効を算出するためのアルゴリズムを取得、最適化するアルゴリズム取得ステップと、生物学的試料の画像の特徴量を取得する特徴量取得ステップと、前記アルゴリズム取得ステップにより算出された前記アルゴリズムを用いて、前記特徴量取得ステップにより特徴量が取得された画像が示す薬効を算出する薬効算出ステップと、をコンピュータが実行することを特徴とする。この画像解析方法によれば、画像の特徴量と、画像に対応する薬効の程度との組み合わせに基づいて、複数の特徴量を用いた演算により薬効を算出するためのアルゴリズムを取得するので、画像処理の専門知識を要することなく、適切に薬効を判定することができる。前記アルゴリズムはニューラルネットワークや、サポートベクターマシンSVMにより定義されてもよい。前記特徴量を抽出する特徴量抽出ステップをコンピュータが実行し、前記アルゴリズムは前記特徴量抽出ステップにおける画像解析のパラメタの指定を含んでもよい。本発明の画像解析装置は、生物学的試料の画像の特徴量から、その画像が示す薬効を判定するための解析を行う画像解析装置において、生物学的試料の画像の特徴量と、前記画像に対応する薬効の程度との組み合わせの入力を受け付ける受付手段と、前記受付手段により入力された組み合わせに基づいて、複数の前記特徴量を用いた演算により前記薬効を算出するためのアルゴリズムを取得、最適化するアルゴリズム取得手段と、生物学的試料の画像の特徴量を取得する特徴量取得手段と、前記アルゴリズム取得手段により算出された前記アルゴリズムを用いて、前記特徴量取得手段により特徴量が取得された画像が示す薬効を算出する薬効算出手段と、を備えることを特徴とする。この画像解析装置によれば、画像の特徴量と、画像に対応する薬効の程度との組み合わせに基づいて、複数の特徴量を用いた演算により薬効を算出するためのアルゴリズムを取得するので、画像処理の専門知識を要することなく、適切に薬効を判定することができる。本発明の画像解析方法によれば、画像の特徴量と、画像に対応する薬効の程度との組み合わせに基づいて、複数の特徴量を用いた演算により薬効を算出するためのアルゴリズムを取得するので、画像処理の専門知識を要することなく、適切に薬効を判定することができる。本発明の画像解析装置によれば、画像の特徴量と、画像に対応する薬効の程度との組み合わせに基づいて、複数の特徴量を用いた演算により薬効を算出するためのアルゴリズムを取得するので、画像処理の専門知識を要することなく、適切に薬効を判定することができる。以下、本発明による画像解析方法の一実施形態について説明する。図1は、本発明による画像解析方法を実行するための画像解析装置の構成例を示すブロック図である。図1に示すように、画像解析装置は、顕微鏡やカメラにより得られた細胞の画像データを取得する画像データ取得部1と、画像解析を行う機能ソフトウェアが実装され、画像データ取得部1から取得した画像データから、個々の細胞について多数の特徴量を演算して演算結果を記録する特徴量抽出部2と、ニューラルネットワークを構築する組み合わせ最適化部3と、ニューラルネットワークを用いて反応曲線を算出する反応曲線出力部4と、を備える。図2は、本発明による画像解析方法で利用するニューラルネットワークの例を示している。図2の例では、層の数は3層である。しかし、本発明による画像解析方法に適用されるニューラルネットワークにおける出力数nの値および層の数はこれに限定されない。次に、画像解析装置の動作について説明する。画像データ取得部1には、種々の薬剤について、複数多数の異なる濃度における画像データが取り込まれ、特徴量抽出部2は、それぞれの画像についての特徴量を演算、記録する。特徴量抽出部2において演算、記録される特徴量は、例えば、画像の明るさに関するパラメタや画像の形状に関するパラメタなどである。画像の明るさとしては、細胞領域内の画素値の合計、平均、標準偏差、あるいはそれらと空間2次元画像のXY方向、あるいは3次元画像のXYZ方向の量を組み合わせた各種特徴量統計量が考えられる。形状に関するパラメタは、形状のモーメントやオイラー数等の幾何学的な特徴量や、テクスチャに関する特徴量が考えられる。図3は、ニューラルネットワークを構築する際における組み合わせ最適化部3の動作を示すフローチャートである。図3のステップS1では、組み合わせ最適化部3は、特定の薬剤についてのネガティブコントロールおよびポジティブコントロールについての画像の特徴量を特徴量抽出部2から取得する。ここで、ネガティブコントロールは最低濃度または薬剤を投与しない場合であり、ポジティブコントロールは最大濃度の場合である。次に、ステップS2では、組み合わせ最適化部3においてニューラルネットワークを構成する。ネガティブコントロールおよびポジティブコントロールのデータは各々、1,0,0,1というベクトルで表し、出力が2ノードであるニューラルネットワークを構成する。各ノードの応答関数は、例えば、下記の1式で示すロジスティック関数等を用いることができるが、これと異なる関数でもよい。fx=a/1exp-x-T1式ただし、aは規格化定数、Tは閾値である。入力ノードは特徴量抽出部2から得られる各細胞の画像の特徴量であり、入力ノード数はこの特徴量の個数である。入力する特徴量は一般に異なる次元を持っており、また、数値の範囲も単位によって異なる。このため、ニューラルネットワークに入力する前段階で規格化を施す必要がある。規格化の方法としては、例えば、その特徴量の最大値と最小値を用いて値を0-1の範囲に規格化し、あるいは平均と標準偏差を用いて規格化してもよい。ステップS2では、複数多数の異なる濃度のうち、とくにネガティブコントロールとポジティブコントロールの2種類の特徴量の組に対して、それぞれ結果が1,0,0,1となるよう、ニューラルネットワークの構成、学習を行う。学習の方法としては、当該分野で公知の誤差逆伝播法Backpropagationを用い、学習係数の調整等を行う。入力する教師信号の数は、ネガティブコントロール、ポジティブコントロール各々のデータにおいて検出した細胞の数に該当するが、必要に応じて一部のデータのみを使用してもよい。ニューラルネットワークは特徴量を組み合わせて薬効を算出するアルゴリズムに対応しており、ニューラルネットワークの学習により、どの特徴量をどのような重みで関連付けるとネガティブコントロールとポジティブコントロールを明確に識別可能になるかが判明する。この学習過程が、熟練者が画像処理の知識に基づいて特徴量を組み合わせて薬効の定量化を試みる作業に対応している。ニューラルネットワークの中間層のノード数については何らかの経験則に従って決定してもよく、あるいは複数の異なる中間層のノード数のニューラルネットワークに対して繰り返し学習を試み、出力誤差が充分に小さくなるようなノード数を選択してもよい。図4は、反応曲線出力部4の動作を示すフローチャートである。反応曲線出力部4は、学習済みのニューラルネットワークを用いて、ネガティブコントロールおよびポジティブコントロール以外の濃度を含むすべての濃度での薬効を判定する処理を実行する。図4のステップS11では、判定の対象となる特定薬剤の特定濃度についての特徴量を特徴量抽出部2から取得する。次に、ステップS12では、特徴量抽出部2から取得された特徴量を学習済みのニューラルネットワークに入力し、出力ベクトルを取得する。次に、ステップS13では、ステップS12で得られた出力ベクトルをスカラー量に変換する。出力ベクトルは2次元のベクトル量であるため、反応曲線として表示しにくい。このため、ステップS13では、出力ベクトルx、yに対してy-x1/2またはy2-x21/212など、1,0が0に、0,1が1になるような演算を行ってベクトル量をスカラー量に変換する。このスカラー量を用いて反応曲線を描き、創薬スクリーニング分野で標準的に用いられているZ ́-factor最大薬効と最小薬効の差が誤差も考慮してどの程度有意かを示す評価量やEC50薬効が最大の半分になる薬剤濃度、SN薬効の最大と最小の比等の評価係数を出力することができる。しかし、スカラー量への変換を行わず、ベクトル量のまま評価を行ってもよい。なお、組み合わせ最適化部3によるニューラルネットワークの学習においては、反応曲線出力部4による算出結果をフィードバックして学習に反映させる処理を繰り返しており、この場合、反応曲線出力部4は学習中のニューラルネットワークを用いた演算を実行し、演算結果を組み合わせ最適化部3に向けて出力する。以上のように、本発明による画像解析方法を用いることにより、薬効が最小と考えられるネガティブコントロールと、薬効が最大と考えられるポジティブコントロールのデータを指定するだけで、自動的に薬効の判定を行うための特徴量の組み合わせが決定される。したがって、画像処理に関する専門知識を要することなく、適切な特徴量の組み合わせによる薬効の判定が可能となる。これにより、例えば、薬学の研究者がネガティブとポジティブの画像が定性的に異なることは認識できるが、定量化の手法が分からないような場合、本発明による画像解析方法を用いることにより、自動的に最適な解析方法を得ることができる。また、本発明による画像解析方法を用いて得られたニューラルネットワークの重み係数を分析することで、薬効による画像の変化に対して、経験等に頼った場合よりも深い知見を得ることができる。上記実施形態では、特徴量抽出部2における画像解析の内容については言及していないが、画像解析での結果得られる特徴量は一般的に画像解析を行う際のパラメタに依存する。このパラメタの単純な例としては、2値化を行うための閾値が挙げられる。滑らかな断面形状を持った細胞であれば、閾値を低くすれば細胞の面積が大きく見え、閾値を高くすれば細胞の面積が小さく見えることは自明である。このため特徴量は画像解析のパラメタの組み合わせに依存する。そこで、ニューラルネットワークの学習結果に基づいて結果をフィードバックし、薬効あるいは薬効を表すZ ́-factorやSNが最大になるように、またはEC50が最小になるように、画像解析のパラメタの最適化を行ってもよい。この場合、特徴量抽出部2で用いられる画像解析のパラメタを変更しながら、組み合わせ最適化部3における演算を繰り返すことにより、最適なパラメタを求めることができる。画像解析パラメタの最適化は、例えば、NelderMeadeのDownhillsimplexなどの手法を用いて行ってもよい。上記実施形態では、細胞の顕微鏡画像を対象とする場合について説明したが、本発明は、細胞以外の生物学的試料細胞の塊であるコロニー、プランクトン、藻類など、および顕微鏡以外の光学系の画像に対しても適用できる。本発明の適用範囲は上記実施形態に限定されることはない。本発明は、生物学的試料の画像の特徴量から薬効を判定するための解析を行う画像解析方法等に対し、広く適用することができる。"