NSM-005

前進/後進 Kolmogorov 方程式 vs 逆問題

結論を先に:「後進方程式 (backward Kolmogorov) 」と「逆問題 (inverse problem) 」はまったく違うもの。混同しやすいが、前者は「方程式の書き方の違い」、後者は「データから未知パラメータを当てる枠組み」。

1. 何が問題か — 用語の混乱

イオンチャネルの状態数を実験データから推定する話を聞くとき、こんな問いが浮かぶかもしれない:

「マルコフ過程やQ行列を前進方程式で解析して未来の挙動を予測するのが普通だけど、
われわれがやりたいのは逆 (データから状態数を当てる) だから、後進方程式を使うべきでは?」

これは自然な発想だが、用語の選び方が少しズレている。整理しよう。

2. 前進方程式と後進方程式 — どちらも「順方向」の計算

まず、両者の式を並べる:

	前進方程式 (forward)	後進方程式 (backward)
式	$\dfrac{dP(t)}{dt} = P(t)\, Q$	$\dfrac{dP(t)}{dt} = Q\, P(t)$
動かすもの	終時刻 $t$ (未来の時刻)	初期状態 $i$ (どこから始めるか)
固定するもの	初期状態 $i$	終状態 $j$
得たいもの	「未来の状態確率」	「初期条件の関数として見た確率」
典型的用途	マクロ電流波形のシミュレーション	滞在時間・初到達時間 (first passage time) の計算

重要: 「前進」「後進」というのは時間の向きではない。両方とも時刻 $t$ について微分する順方向の微分方程式。「何を変数として見るか」が違うだけ。
解はどちらも同じ $P(t) = \exp(Qt)$。
（理由: $e^{Qt}$ は行列指数関数の定義 $\sum_n (Qt)^n/n!$ から $Q$ と可換 $\;Q e^{Qt} = e^{Qt} Q\;$ なので、$\dot P = PQ$ と $\dot P = QP$ の両方を同時に満たす。）

図解: 何を動かして何を固定するか

3. 逆問題 (inverse problem) とは別の話

われわれが本当にやりたいのは:

観測データ (open/closed の時系列) ⟹ Q 行列・状態数を推定

これは「方程式の向き」の話ではなく、パラメータ推定 (parameter estimation) の話。
方程式は前進でも後進でも、解 $P(t) = \exp(Qt)$ は同じ。逆問題ではこの解をモデルとして使い、「観測データの確率 (尤度) を最大化する Q を探す」。

図解: 順問題と逆問題

覚え方:

順問題: Q が既知 → 確率/データを予測 (前進・後進方程式どちらも使う)
逆問題: データが既知 → Q を推定 (尤度の中で順問題の答えを計算する道具として使う)

4. イオンチャネルの例で具体化

2状態モデル $\mathrm{C} \rightleftharpoons \mathrm{O}$ で考える ($\alpha$: 閉から開、$\beta$: 開から閉) 。

$Q = \begin{pmatrix} -\alpha & \alpha \\ \beta & -\beta \end{pmatrix}$

順問題のシナリオ

$\alpha, \beta$ が分かっている (例: 過去の研究で測定済み)
前進方程式で「t秒後に開状態にいる確率」を計算 → マクロ電流波形を予測
後進方程式で「閉状態から始めたとき開状態への平均初到達時間」を計算

逆問題のシナリオ (今回のテーマ)

実験データから滞在時間ヒストグラムを作る
$\alpha, \beta$ は未知
「データの尤度 $L(\alpha, \beta \mid \text{data})$ を最大化する $\alpha, \beta$ は?」を解く
さらに難しい問い: 「そもそも状態数は本当に2つか? 3つかも?」 → モデル選択 (BIC など)

5. なぜこの混乱が起きやすいか

「後進 (backward) 」という言葉から「データ → モデル」という逆向きの矢印を連想しがちだから。
しかし backward Kolmogorov の「backward」は初期状態を変数として見るという意味で、時間や因果の向きとは無関係。

歴史的注: Kolmogorov 自身が1931年の論文で「前進」「後進」と命名した。英語の forward / backward は数学者には自然な区別だが、応用分野の人には混乱を招くことがある。

6. プレゼンでの取り扱い指針

聴衆 (数学者) に対しては、「前進・後進方程式は両方とも順方向に確率を計算する道具」と明示し、後進方程式は滞在時間分布の導出のために必要だと位置付ける
「われわれのゴールは逆問題 (inverse problem)」と冒頭で宣言し、最尤推定の文脈で順問題の方程式が尤度計算の中で使われると説明する
「後進方程式 = 逆問題」と誤って混同しない

7. まとめ

用語	何の話か	使うとき
前進方程式	方程式の書き方 (初期状態固定)	未来の状態確率を計算
後進方程式	方程式の書き方 (終状態固定)	滞在時間・初到達時間を計算
逆問題	パラメータ推定の枠組み	データから Q を当てる

参考文献

Colquhoun, D. & Hawkes, A. G. (1995). The principles of the stochastic interpretation of ion channel mechanisms. In: Sakmann, B. & Neher, E. (eds.) Single-Channel Recording, 2nd edition, pp. 397–482. Plenum Press, New York. → プロジェクト全体の主要典拠：Kolmogorov 方程式と最尤推定枠組みの統合解説
Colquhoun, D. & Hawkes, A. G. (1982). On the stochastic properties of bursts of single ion channel openings and of clusters of bursts. Philosophical Transactions of the Royal Society of London B, 300, 1–59. → 尤度計算と Kolmogorov 後進方程式の応用：パラメータ推定（逆問題）の定式化
Colquhoun, D., Hawkes, A. G. & Srodzinski, K. (1996). Joint distributions of apparent open and shut times of single-ion channels and maximum likelihood fitting of mechanisms. Philosophical Transactions of the Royal Society of London A, 354, 2555–2590. → dead time 補正を含む尤度計算（HJC theory）：逆問題の実装における最重要論文