[論文レビュー] Analysis of the Conradi-Kahle Algorithm for Detecting Binomiality on Biological Models

Alexandru Iosif, Hamid Rahkooy|arXiv (Cornell University)|Dec 14, 2019

Computational Drug Discovery Methods参考文献 12被引用数 4

ひとこと要約

この論文は、生物学的モデルから導出された多項式イデアルにおける二項性の検出に向けたConradi-Kahleアルゴリズムを分析している。本研究では、Macaulay2およびMapleにそのアルゴリズムを実装し、BioModelsリポジトリの20のバイオモデルに対してグレブナー基底計算と性能を比較した。主な発見は、Conradi-Kahleアルゴリズムは高速（大多数のケースで1秒未塔）であるが、非同次イデアルに対して完全性に欠けるため、正しく二項性を特定できるのは3つのモデルにとどまり、一方グレブナー基底は遅いものの、600秒以内に終了する場合には非二項性を検出できる。

ABSTRACT

We analyze the Conradi-Kahle Algorithm for detecting binomiality. We present experiments using two implementations of the algorithm in Macaulay2 and Maple on biological models and assess the performance of the algorithm on these models. We compare the two implementations with each other and with Gröbner bases computations up to their performance on these biological models.

研究の動機と目的

• 生物学的モデルからの多項式イデアルにおける二項性検出のためのConradi-Kahleアルゴリズムの効率性および有効性を分析すること。
• 実際の生物学的モデルに対して、Macaulay2およびMapleに実装されたConradi-Kahleアルゴリズムの性能をグレブナー基底計算と比較すること。
• 生物学的反応ネットワークで一般的な非同次イデアルに適用した場合のConradi-Kahleアルゴリズムの限界を調査すること。
• グレブナー基底が合理的な時間内に終了しない場合でも、アルゴリズムが生化学的ネットワークの定常状態イデアルにおける二項性を信頼性を持って同定できるかどうかを評価すること。

提案手法

• 20のバイオモデル（BioModelsリポジトリより）を対象とし、反応速度パラメータを有理関数体上の非同次多項式イデアルとして扱う。
• アルゴリズムは、係数行列の簡約行階段形を繰り返し計算し、3つ以上の非ゼロ成分を含む行があるかどうかをチェックすることで二項性をテストする。
• ベンチマークとしてグレブナー基底計算を用い、EisenbudとSturmfelsの補題1.2を応用して、グレブナー基底の還元による二項性のテストを実施する。
• 両システム（Macaulay2およびMaple）に実装し、実行時間を記録し、両システム間での結果を比較する。
• 元の非同次イデアルの同次化版にもアルゴリズムを適用し、同次イデアルに限って完全性が保証されることを認識する。
• CPU時間と二項性同定の成功確率を測定することで性能を評価し、グレブナー基底には600秒のタイムアウトを設定する。

実験結果

リサーチクエスチョン

• RQ1 生物学的モデルからの非同次多項式イデアルにおける二項性検出において、Conradi-Kahleアルゴリズムの速度と正確性はいかがであるか？
• RQ2 実際の生物学的モデルにおいて、Conradi-Kahleアルゴリズムはグレブナー基底を上回るか、下回るか、その程度は？
• RQ3 効率性に欠かず、Conradi-Kahleアルゴリズムが大多数のバイオモデルで二項性を同定できないのはなぜか？
• RQ4 Conradi-Kahleアルゴリズムは非二項イデアルを検出できるのか、それとも二項性の確認のみに限定されるのか？
• RQ5 Macaulay2とMapleにおけるConradi-Kahleアルゴリズムの計算性能は、生物学的モデルにおいてどのように比較されるか？

主な発見

• テストしたすべてのバイオモデルに対して、Conradi-Kahleアルゴリズムは1秒未塔で実行され、Macaulay2では0.01〜0.7秒、Mapleでは0.02〜6秒の実行時間であった。
• アルゴリズムは3つのバイオモデル（9, 262, 264）でのみ二項性を正しく同定し、Bin (C-K) 列に「Yes」と出力したが、残りのすべてのモデルでは「No」と出力した。
• グレブナー基底計算は20のバイオモデルのうち10つで600秒以内に終了し、最も長かったのはMapleでのバイオモデル153で500秒であった。
• バイオモデル46, 187, 243, 335では、グレブナー基底が非二項性を正しく同定したが、Conradi-Kahleアルゴリズムは非同次イデアルにおける不完全性のため、同定に失敗した。
• 両手法が成功したケースでは実行時間がほぼ同等であった：バイオモデル264ではC-K (M2) が0.7秒、GB (M2) が2秒、C-K (Maple) が0.03秒、GB (Maple) が0.04秒であった。
• Conradi-Kahleアルゴリズムおよびグレブナー基底計算において、Macaulay2とMapleの実装間に顕著な性能差は認められなかった。