不確実性の世界へようこそ：遷移強度の信頼区間
こんにちは！ALTAMカリキュラムのここまでたどり着いた皆さんなら、すでに遷移強度（「健康」から「病気」へといった、ある状態から別の状態へ人を押しやる「力」）の計算方法はご存知ですね。しかし、現実世界ではデータは完璧ではありません。通常、私たちは母集団全体ではなく、一部のサンプルを扱っているからです。
この章では、私たちが推定した値に対してどれくらい「確信」が持てるのか、その表現方法を学びます。単に「強度は0.05です」と言うのではなく、「真の強度は0.04から0.06の間にあると95%の確信をもって言えます」というように表現します。これはリスク管理において極めて重要です！
1. クイック復習：遷移強度のMLE（最尤推定量）
信頼区間を作る前に、まずは出発点が必要です。定数の遷移強度 \(\mu_{ij}\) （状態 \(i\) から状態 \(j\) へ移動する力）に対しては、最尤推定量（MLE）を用います。
公式：
\(\hat{\mu}_{ij} = \frac{d_{ij}}{V_i}\)
ここで：
- \(d_{ij}\) は、状態 \(i\) から状態 \(j\) への遷移が観測された総数。
- \(V_i\) は、すべての個人が状態 \(i\) に留まった合計時間（待ち時間や暴露時間とも呼ばれます）。
例え話：ホテルの回転ドアを想像してください。\(V_i\) は人々がロビーにいた合計時間であり、\(d_{ij}\) は実際にドアを通って通りに出た人数です。「強度」とは、ロビーに滞在した時間あたりの離脱率に過ぎないのです。
クイック復習ボックス：
常に覚えておきましょう：強度 = (イベント数) / (リスクにさらされた合計時間)。表記に惑わされないでくださいね！
2. 推定量の分散
信頼区間を作るには、推定量の「ばらつき」や分散を知る必要があります。データが少ない（\(V_i\) が小さい）場合、推定値は大きくばらつきます。逆に、膨大なデータがあれば、推定値は非常に安定します。
十分大きなサンプルの場合、強度の推定分散は以下のようになります：
\(\text{Var}(\hat{\mu}_{ij}) \approx \frac{\hat{\mu}_{ij}}{V_i} = \frac{d_{ij}}{(V_i)^2}\)
標準誤差（SE）は、分散の平方根をとるだけです：
\(SE(\hat{\mu}_{ij}) = \sqrt{\frac{d_{ij}}{(V_i)^2}} = \frac{\sqrt{d_{ij}}}{V_i}\)
豆知識：この分散の公式は、一定期間内の遷移回数のモデル化によく使われるポアソン分布の性質から導かれています。
3. 標準（線形）信頼区間
これは最もシンプルな区間です。（中心極限定理のおかげで）推定量が正規分布に従うと仮定します。
公式：
\(\hat{\mu}_{ij} \pm z_{\alpha/2} \cdot SE(\hat{\mu}_{ij})\)
95%信頼区間の場合、\(z_{\alpha/2}\) は 1.96 です。
問題点：遷移強度は必ず正の値でなければなりません。しかし、強度が小さく分散が大きい場合、この公式の「マイナス」計算部分が負の下限値を与えてしまうことがあります。死亡率などの強度は負になり得ないため、より良い方法が必要です！
4. 対数変換信頼区間
「負の値の問題」を避けるために、対数変換を用います。まず \(\ln(\hat{\mu}_{ij})\) の信頼区間を計算し、それを元に戻します。\(e^x\) は常に正であるため、最終的な区間も必ず0より大きくなります！
ステップ・バイ・ステップ：
1. 推定量の自然対数の標準誤差を計算します：
\(SE(\ln \hat{\mu}_{ij}) = \sqrt{\frac{1}{d_{ij}}}\)
（なんとシンプルなことでしょう！遷移数のみに依存します。）
2. 「誤差係数（Error Factor、\(EF\)）」を計算します：
\(EF = \exp \left( z_{\alpha/2} \cdot \frac{1}{\sqrt{d_{ij}}} \right)\)
3. 信頼区間は以下の通りです：
下限値： \(\frac{\hat{\mu}_{ij}}{EF}\)
上限値： \(\hat{\mu}_{ij} \cdot EF\)
暗記のコツ：対数変換版では、足し算・引き算をする代わりに、誤差係数を掛けたり割ったりします。これにより、結果が必ず正になることが保証されるのです！
例：25人の死亡（\(d_{ij} = 25\)）が500人年（\(V_i = 500\)）の間に観測された場合：
- \(\hat{\mu} = 25/500 = 0.05\)
- \(SE(\ln \hat{\mu}) = 1/\sqrt{25} = 0.2\)
- 95%信頼区間の場合：\(EF = \exp(1.96 \cdot 0.2) = 1.479\)
- 下限： \(0.05 / 1.479 = 0.0338\)
- 上限： \(0.05 \cdot 1.479 = 0.0739\)
5. まとめと重要ポイント
ALTAM試験で遷移強度の信頼区間を求められたら、以下の点を心に留めておいてください：
重要ポイント：
- MLE： まずは必ず \(\hat{\mu}_{ij} = d_{ij} / V_i\) を計算しましょう。
- 線形信頼区間： \(\hat{\mu} \pm z \cdot \frac{\sqrt{d}}{V}\) を使います。問題で対数変換の指定がない場合にのみ使用してください（実務では対数変換の方が好まれます）。
- 対数変換信頼区間： 境界値が正であることを保証するために使用します。重要な公式は \(SE(\ln \hat{\mu}) = 1/\sqrt{d}\) です。
- サンプルサイズ： 観測される遷移数（\(d\)）が多いほど、信頼区間は狭くなります。
避けるべきよくある間違い：
\(V_i\)（状態 \(i\) にいた時間）と研究対象の全人数を混同しないように。その特定の強度については、出発状態にいた時間のみが重要です！
結論： 信頼区間は、真の遷移強度として妥当な範囲を教えてくれます。対数変換は、これらの率が負になり得ないという現実を尊重してくれる、私たちの最高の味方です。
最初は数式が難しく見えるかもしれませんが、大丈夫です。要は「ベストな推測値」を見つけ、その推測がどれくらい外れる可能性があるかを把握するだけですよ！