モデルが勝てないときに追加する観測

Daiki Kumakura · 2026-10-25 · Blog · 1.0

摂動予測とVirtual Cellの最新研究を読み、情報の不足、モデルの未活用、評価の低感度を切り分ける。観測時点と信号の希釈の数式例から、次に測るものを考える。

文献確認日：2026年10月9日。以下では査読済み論文、プレプリント、開発元の報告を区別する。図と数値表は本記事の数式例であり、実測の摂動データを再解析した結果ではない。

新しいモデルが、平均値を返すだけの予測に勝てない。そこで層を増やすか、事前学習を大きくするか。それより先に、確かめたいことがある。予測したい違いがデータに含まれているか。その違いを評価が拾えているか。

この二つは別の問題だ。情報のない観測から推定できないことはある。一方、情報があってもモデルが使えていない、あるいは評価が改善を見落としていることもある。「勝てなかった」という結果だけでは区別できない。

Eric Kernfeldは2025年5月のブログ記事で、遺伝子制御ネットワークの推論ではモデルの工夫に先立ち、データの情報量と必要な追加観測を考えるべきだと論じた。本稿もこの問いから始める。ただし、これは研究方針についての意見であり、すべての摂動予測に対する不可能性の証明ではない。

細胞数と区別できる条件数

単一細胞データでは細胞数が大きくなりやすい。しかし、同じ条件から細胞を増やすことと、別の時点・摂動・背景を測ることは同じではない。

同じ条件の細胞を増やせば、その条件の平均や分布は詳しくなる。だが、一つの時点しかなければ、その前に何が起きたかは直接にはわからない。同じ細胞株しかなければ、別の細胞株への移行は検証できない。多数の細胞があっても、独立した実験バッチやドナーの反復が少なければ、その不確実性は細胞数だけでは解消しない。

例えば、一つの培養から得た1万細胞を、そのまま1万回の独立した生物学的反復とは数えない。細胞内・条件内のばらつきと、培養・ドナー・実験日を跨ぐばらつきを分け、後者を一般化したいなら独立した反復を用意する。

これは「細胞を増やしても無意味」という話ではない。追加する観測が、どの不確実性を減らすのかを分けて考えるという話だ。

観測条件による情報量の違い

減衰する応答を考える。

\[ y(t)=A\exp(-kt)+\varepsilon,\qquad \varepsilon\sim N(0,\sigma^2). \]

ここでは各測定の誤差は独立、分散は既知で一定と仮定する。振幅\(A\)が既知なら、減衰速度\(k\)についてのFisher情報は

\[ I_k=\frac{A^2}{\sigma^2}\sum_r t_r^2\exp(-2kt_r) \]

になる。時刻0を何回測っても、\(k\)の情報は0だ。一観測あたりの情報は\(t=1/k\)で最大になるが、これは振幅既知・独立な一定分散誤差・このモデルの条件下での話である。未知パラメータが増えれば、最適な時点も変わる。

振幅も未知なら、一つの正の時点\(t_0\)を繰り返し測ってわかるのは\(A\exp(-kt_0)\)だけだ。同じ値を与える振幅と速度の組は無数にある。反復数を増やせばその時点の平均は精密になるが、この二つは分離できない。

図1．観測時点と識別可能性。 正規化した単位で計算した説明用の例。上段は\(A=k=1,\sigma=0.1\)。下段の二曲線は時刻1の平均が一致する。実際の細胞応答に指数減衰を当てはめたものではない。

同じ10測定でも設計を変えると、振幅と速度についての情報行列の階数が変わる。以下はパラメータ\((\log A,\log k)\)についての局所計算である。

出力：results/designs.csv。振幅を固定した条件付き情報と、振幅を未知として扱う情報を混同しないことが重要だ。二時点ならこの例では局所的に分離できるが、実データで十分に精密に推定できるかは別途確認する必要がある。

実際のモデルでは、パラメータを一つずつ固定し、他を再推定するprofile likelihoodが役立つ。Raueらの論文は、構造上の識別不能と、有限データによる実用上の識別困難を扱っている。平坦な尤度だけで一般的な不可能性を断定するのではなく、モデルの対称性、観測設計、誤差モデル、探索範囲を併せて調べたい。

摂動予測の論文が比較するもの

近年の議論を「深層学習は勝てる／勝てない」の二択にすると、肝心な違いを見失う。比較するベースライン、評価する出力、未知条件の定義が違うからだ。

Millerらの論文は特に重要だ。同じ摂動の細胞を二分したtechnical duplicateと、弱い応答を平均値へ寄せるinterpolated duplicateを用い、陽性・陰性対照の差を指標がどれだけ拾うかを調べた。前者は有限細胞から作る推定なので、弱い信号では平均ベースラインよりばらつくこともある。真の条件付き平均を知る仮想的な予測器とは違い、独立した生物学的実験の再現性を保証するものでもない。

線形・加算ベースラインとの比較も含み、Wessels23の組合せ予測では加算を上回るモデルが報告された。一方、Norman19では加算が強かった。学習範囲だけでなく、応答の強さやデータ構成、モデル世代も異なるため、差の原因を一つに決められない。「どの設定でも勝てない」も「常に勝つ」も適切でなく、この研究は未知の細胞背景への一般化を検証したものでもない。

ここでのcalibrationは評価指標の感度の点検であり、予測確率の較正と同義ではない。また、評価用の実測応答を使う陽性対照は、未知条件に投入できる予測器ではない。試験データに基づく対照を採点のために使うことと、その情報を学習やモデル選択へ戻すことを分ける必要がある。

評価が改善を隠す状況は、簡単な期待値でも説明できる。評価する\(G\)遺伝子のうち10遺伝子だけが1だけ変化し、各遺伝子の独立な観測誤差の分散を1とする。変化なしを予測するモデルの期待MSEは\(1+10/G\)、真の条件付き平均を知る仮想的な予測器は1になる。

図2．固定した信号の希釈。 両軸は対数。出力：results/metrics.csv。この例では順位は逆転せず、真の平均を知る予測器が常に良い。差が小さいことと統計的に検出不能なことも同じではない。十分な独立反復があれば小さな差を検出できる可能性がある。

これはMillerらの評価を再現した図ではなく、全体指標における信号の薄まりを示す例だ。実データで「変化した遺伝子だけ」を選ぶ場合も、選択方法と対照を明示しなければ、新たな評価の偏りが入る。

細胞数以外の情報を増やす方法

RNAだけでは区別しにくい関係が、別の観測層や背景情報で区別できることがある。LINGERはRNAとクロマチンの情報、外部の情報を利用して制御ネットワークを推論する。ベンチマーク上の改善は、推定された全ての辺が因果的だという保証ではないが、入力を増やすという方向を具体化している。

2026年には、観測の質と背景の広がりを扱う報告も出た。Wangらのプレプリントは、摂動が十分に効いているか、共通応答と特異的応答を区別できるかを点検する。そこでの分類は研究固有の定義に依存するため、「公開データの一定割合は使えない」という一般則にはしない。

Stateの査読論文は、細胞の集合と背景間の移行を扱う。10月5日公開のPIEは、既存の生物学的知識と摂動データから集団レベルの応答を予測するプレプリントだ。PIEの開発元による説明では、細胞背景と摂動が両方未知になると精度が落ちることも示されている。これは独立した追試ではなく、著者による報告として読む。

PIEのように入力する知識や出力の定義を変えた改善を、単純に「大きなモデルが勝った」とまとめるのは適切でない。比較するときは外部知識の出所、利用可能な時点、事前学習の対象、評価条件との重なりも確認したい。外部知識の利用そのものは不正ではなく、その情報を使える用途かどうかが問題になる。

また、Arcの2026年のデータ作成報告では、複数の細胞株で手順を揃え、ノックダウンと対照を点検している。背景を広げ、測定の一貫性を上げる試みだ。ただし、時点を揃えることと時系列を測ることは違う。こうしたデータから時間変化まで検証できるとは限らない。

追加実験を選ぶ目的

追加実験を選ぶ目的も一つではない。モデルの不確実性を減らすことと、関心のある状態へ細胞を動かすことは違う。

パラメータ\(\theta\)、既存データ\(D\)、実験\(e\)の結果\(Y_e\)について、期待情報利得を使うなら

\[ e^*=\arg\max_e I(\theta;Y_e\mid D,e) \]

と書ける。これは、結果を見る前に、事後分布が平均的にどれだけ変わるかを評価する考え方だ。ただし、事前分布、誤差モデル、実験候補、関心パラメータを決めて初めて意味を持つ。モデルが間違っていれば、そのモデルの中でだけ効率の良い実験を選ぶ可能性がある。

期待情報利得と効用に基づく設計の一般的な整理は、ChalonerとVerdinelliの総説を参照できる。候補モデルが違う応答を予測するだけでは十分でなく、その差を測定誤差の中で検出でき、実施できることも必要だ。

Zhangらは、望む応答に近づく介入の選択を扱う。SandorとAntalは、ネットワーク構造の不確実性を減らすBayesian active learningを扱い、シミュレーション等で検討した。後者を、前向きの実験で普遍的な有効性が確立した方法とは扱わない。

NAIADは組合せ摂動の応答予測と探索を扱う。大きな効果を持つ組合せを見つけたい場合と、ネットワークを区別したい場合では、同じ実験を選ぶとは限らない。情報量が多い実験が、実務上最も有用とも限らない。費用、実施可能性、意思決定の損失まで含めた効用が必要になる。

モデルを変える前の確認表

私なら、次のように切り分ける。これは各論文から導いた研究上の提案であり、この手順全体を検証した実験結果ではない。用途によっては、最後の「何の判断を改善したいか」を最初に固定する。

識別不能なパラメータでも、関心のある予測は安定することがある。逆に、既存条件でよく予測できても、新しい時点や背景では崩れることがある。したがって、情報量の監査はモデルの価値を否定するためではなく、どの用途なら使え、どの観測を追加すれば次へ進めるかを決めるために行う。

実験への提案としては、まず少数の標的・背景・時点を使ったパイロットで、摂動効率、条件内の精度、独立反復間の変動を確かめる。その結果から、細胞数を増やすか、独立反復を増やすか、時点や背景を広げるかを選ぶ。ここで配分の最適値を実測データから求めたわけではなく、必要な情報を見積もるための出発点である。

モデルが勝てないとき、すぐに「もっと良いモデル」や「もっと多い細胞」へ進まない。まず、比べたい違い、測っているもの、採点しているものを揃える。そのうえで、候補間の違いを測定誤差に対して見分けられ、実施可能な観測を選ぶ。候補同士が同意していても、共通の誤りを見つける独立検証は必要だ。

再実行

図と数値表はinformation_demo.pyから生成した。PythonとNumPy、Matplotlibが必要で、乱数や外部データは使わない。python information_demo.pyで生成し、python -m unittest test_demo.py -qで数式と情報行列の性質を検査する。文献の大規模モデルを再実装・検証した資料ではない。

図表と計算コードをダウンロード