尤度 観測した結果 B を固定したうえで「もし原因が A だったとしたらこの B が出る確率」を表す量。条件付き確率 P(B|A) と同じ式だが、A を動かす関数として読む視点の違いから別名を持つ。ベイズの定理 P(A|B) = P(B|A) P(A) / P(B) では事後確率の分子に現れ、最大尤度推定では推定値を選ぶ基準として最大化される。 P ( B ∣ A ) P(B \mid A) P ( B ∣ A ) とは、観測した結果 B B B を固定したうえで「もし原因が A A A だったとしたらこの B B B が出る確率」を表す量です。
定義
尤度の式は条件付き確率 事象 B が起きたとわかったときに事象 A も起きている割合。A と B の同時確率を B の確率で割って定義され、全体を B の縦帯に絞り直してから A の比を測る操作として読める。 とまったく同じです。
ここで B B B は「観測した結果」、A A A は「考えられる原因や仮説」を指します。
条件付き確率 では A A A を固定して B B B の確率を計算しますが、尤度では B B B を固定して A A A を動かす関数 として読みます。この視点の違いが「尤度」という別名を生む理由です。
A A A を動かす関数として明示するときは次の記法も使います。
セミコロン ; の左が「動かすもの」、右が「固定するもの」です。
性質
条件付き確率 と同じ式、視点が違う : 条件付き確率 は A A A を固定して B B B の確率を見ます。尤度は B B B を固定して A A A を動かします。式は共通ですが、何を変数とするかが正反対です。
尤度は確率分布ではない : A A A について総和を取ると ∑ A P ( B ∣ A ) \sum_{A} P(B \mid A) ∑ A P ( B ∣ A ) は 1 になりません。1 を超えることも 1 未満になることもあります。個々の値 P ( B ∣ A ) P(B \mid A) P ( B ∣ A ) は確率ですが、A A A を動かす関数としての尤度は確率分布ではありません。
ベイズの定理 観測された結果から原因の確率を逆向きに計算する公式。事前確率と尤度から事後確率を導く、確率更新の基本道具。 の分子 : 事前確率 観測データを受け取る前に事象に対して割り当てておく確率。ベイズの定理 P(A|B) = P(B|A) P(A) / P(B) の中で事後確率を計算する起点になる。 P ( A ) P(A) P ( A ) と掛け合わせた積 P ( B ∣ A ) ⋅ P ( A ) P(B \mid A) \cdot P(A) P ( B ∣ A ) ⋅ P ( A ) が事後確率 観測したデータを受け取った後の、ある事象が成り立つ条件付き確率。事前確率 × 尤度 / 周辺確率の比でベイズの定理から計算され、観測前の確率を観測後の確率に更新する量を表す。 の分子を作ります。
尤度の比較が判断の基準になる : P ( B ∣ A 1 ) P(B \mid A_1) P ( B ∣ A 1 ) と P ( B ∣ A 2 ) P(B \mid A_2) P ( B ∣ A 2 ) の比 Λ = P ( B ∣ A 1 ) / P ( B ∣ A 2 ) \Lambda = P(B \mid A_1) / P(B \mid A_2) Λ = P ( B ∣ A 1 ) / P ( B ∣ A 2 ) は「観測 B B B が仮説 A 1 A_1 A 1 と A 2 A_2 A 2 のどちらをより強く支持するか」の指標になります。
視覚的に見る
医療検査の例で確認します。疾患の有病率は 1%(P ( 病気 ) = 0.01 P(\text{病気}) = 0.01 P ( 病気 ) = 0.01 )とします。検査の感度(病気のとき陽性が出る確率)は 0.90、偽陽性率(健康のとき陽性が出る確率)は 0.09 です。
縦軸を「陽性か陰性か」、横軸を「病気か健康か」で描くと、尤度は 各帯の中の縦比率 として現れます。
病気 1% 健康 99%
尤度: 各帯の中での陽性比率 横幅は病気1%、健康99%。それぞれの帯の中で、緑の陽性部分は病気90%、健康9%を占める。
横幅は事前確率 。各帯の中で、緑は陽性・灰色は陰性を表します。
病気のとき陽性となる比率 90% P(陽性 | 病気) = 0.90
健康のとき陽性となる比率 9% P(陽性 | 健康) = 0.09
帯の幅(横方向)が事前確率 です。病気帯は幅 0.01、健康帯は幅 0.99 です。
尤度はこの帯の幅とは無関係で、帯の中の縦方向の比率 です。病気帯の中で陽性が占める比率が 0.90、健康帯の中で陽性が占める比率が 0.09 です。この 2 つの値を「A A A を動かして比較する」操作が尤度を使うことの中身です。
A A A について和を取ると P ( 陽性 ∣ 病気 ) + P ( 陽性 ∣ 健康 ) = 0.90 + 0.09 = 0.99 P(\text{陽性} \mid \text{病気}) + P(\text{陽性} \mid \text{健康}) = 0.90 + 0.09 = 0.99 P ( 陽性 ∣ 病気 ) + P ( 陽性 ∣ 健康 ) = 0.90 + 0.09 = 0.99 で、1 になりません。尤度が確率分布ではないことが数字で確認できます。
実世界での使われ方
医療検査の感度と特異度 : 感度(se nsitivity)は「病気のときに陽性が出る確率」で、尤度 P ( 陽性 ∣ 病気 ) P(\text{陽性} \mid \text{病気}) P ( 陽性 ∣ 病気 ) そのものです。国立感染症研究所(感染症発生動向調査 )が公表する検査性能指標はこの尤度を基準に設計されています。特異度(specificity)から 1 − 特異度 = P ( 陽性 ∣ 健康 ) 1 - \text{特異度} = P(\text{陽性} \mid \text{健康}) 1 − 特異度 = P ( 陽性 ∣ 健康 ) が逆方向の尤度を与えます。
不良率の最大尤度推定 : 製品の不良率 p p p を推定するとき、n n n 個の検査結果 B = ( X 1 , … , X n ) B = (X_1, \ldots, X_n) B = ( X 1 , … , X n ) から「観測 B B B を最もよく説明する p p p 」を L ( p ; B ) L(p;\, B) L ( p ; B ) の最大化として求めます。総務省「統計的品質管理 」の基準で標準的に使われる手法です。
GARCH モデルのボラティリティ 推定 : 株式リターンのボラティリティ 金融でいうボラティリティは、資産価格の変動の大きさを示す尺度です。上昇・下落の方向や、投資リスク全体をそのまま表すものではありません。 を推定する GARCH モデルでは、観測リターン系列の尤度を最大化してパラメータを求めます。日本銀行「金融市場レポート 」のボラティリティ 推定に使われる標準手法です。
機械学習の損失関数 : ロジスティック回帰やニューラルネットの学習は、訓練データの対数尤度 ∑ i = 1 n log P ( y i ∣ x i ; θ ) \sum_{i=1}^{n} \log P(y_i \mid x_i;\, \theta) ∑ i = 1 n log P ( y i ∣ x i ; θ ) を最大化するパラメータ θ \theta θ を勾配法で探す操作です。「クロスエントロピー損失」は対数尤度の符号反転と等価で、損失の最小化は尤度の最大化と同じです(scikit-learn の実装解説: scikit-learn ロジスティック回帰 、理論的背景は Bishop Pattern Recognition and Machine Learning , Springer, 2006)。
深掘り
尤度比の議論は仮説検定を学んだ段階で戻ってくれば十分で、初読では飛ばせます。
尤度比と仮説の支持度
2 つの仮説 A 1 , A 2 A_1, A_2 A 1 , A 2 に対して、観測 B B B の 尤度比 を次のように定義します。
医療検査の例では Λ = 0.90 / 0.09 = 10 \Lambda = 0.90 / 0.09 = 10 Λ = 0.90/0.09 = 10 です。陽性という観測は病気仮説を健康仮説の 10 倍強く支持します。
ベイズの定理 を整理すると 事後オッズ = 事前オッズ × 尤度比 という形になります。
観測 B B B が事前オッズを何倍に更新するかを示す乗数が尤度比です。Neyman と Pearson は 1933 年に「サイズ固定の検定で最大の検出力を持つのは尤度比に基づく検定である」(Neyman-Pearson 補題)を証明しました。
対数尤度と独立 観測
独立 2 つの事象 A, B について、B が起きたという情報を得ても A の確率が変わらない関係。数学的には P(A ∩ B) = P(A) P(B) または同値に P(A|B) = P(A) で定義する。排反性とは正反対の概念。 な観測 B 1 , … , B n B_1, \ldots, B_n B 1 , … , B n に対して、結合尤度は積になります。
サンプルサイズが大きくなると積は数値的に 0 に近づいて扱いにくくなります。対数を取って和に直すと計算が安定します。
対数は単調増加関数なので、尤度の最大点と対数尤度の最大点は一致します。最大尤度推定では対数尤度の最大化が標準です(Fisher, 1922)。
Fisher 情報量と推定精度
A A A を連続パラメータとするとき、対数尤度関数の 2 階微分の期待値 確率変数が取りうる各値に、その確率を重みとして掛けて足し合わせた値。長期的に試行を繰り返したときの観測平均の収束先を表す。離散型では値と確率の積の総和、連続型では密度関数による積分で定義される。 の符号を反転した量を Fisher 情報量 I ( A ) I(A) I ( A ) と呼びます。
Cramér(1946)と Rao(1945)は独立 に、任意の不偏推定量の分散 分散は、データが平均の周りにどれほど散らばっているかを表す統計量です。各データと平均との差(偏差)を二乗して平均し、元データの単位の二乗で表します。 が 1 / I ( A ) 1/I(A) 1/ I ( A ) で下から押さえられることを証明しました(Cramér-Rao 下限)。Fisher 情報量が大きいほど観測から精度よくパラメータを推定できます。
直感的には、対数尤度のグラフが鋭く尖っているほど Fisher 情報量が大きく、最大点の位置を精度よく特定できます。平坦な対数尤度は情報量が少なく、推定精度が落ちます。
頻度主義とベイズでの使い方の違い
頻度主義と呼ばれる立場では、尤度 L ( A ; B ) L(A;\, B) L ( A ; B ) を A A A の関数として最大化します(最大尤度推定)。A A A は固定された真の値であり、確率を持ちません。
ベイズの立場では、尤度 L ( A ; B ) L(A;\, B) L ( A ; B ) を事前確率 P ( A ) P(A) P ( A ) と掛け合わせて事後確率 を求めます。
A A A にも確率分布を割り当てる点が頻度主義と異なります。Fisher(頻度主義の最大尤度推定の提唱者)と Jeffreys(ベイズ統計の体系化者)の論争が代表的な対立です。現代の機械学習では変分推論・経験ベイズなど両者の中間的な手法が広く使われており、実用上の対立は薄れています。
関連する用語
条件付き確率 (尤度 P ( B ∣ A ) P(B \mid A) P ( B ∣ A ) は条件付き確率 と同じ式。何を固定して何を動かすかの視点が異なります)
事前確率 (ベイズの定理 で尤度と掛け合わせて事後確率 の分子を作ります)
事後確率 (事前確率 × 尤度 / 周辺確率 同時分布から特定の変数を足し合わせて消したあとに残る、ある変数単独の確率。ベイズの定理の分母(evidence)として現れる。 で計算される観測後の確率)
ベイズの定理 (事前確率 ・尤度・事後確率 を結びつける公式)
周辺確率 (ベイズの定理 の分母 P ( B ) = ∑ A P ( B ∣ A ) ⋅ P ( A ) P(B) = \sum_{A} P(B \mid A) \cdot P(A) P ( B ) = ∑ A P ( B ∣ A ) ⋅ P ( A ) )
この用語を扱う本編記事
stats-07: 条件付き確率とベイズの面積 (医療検査の例で P ( 陽性 ∣ 病気 ) = 0.90 P(\text{陽性} \mid \text{病気}) = 0.90 P ( 陽性 ∣ 病気 ) = 0.90 を感度として導入し、事前確率 と尤度を動かしたときの事後確率 の変化を扱います)
よくある誤解
「尤度は確率だから A A A について足すと 1 になる」という誤解があります。個々の値 P ( B ∣ A ) P(B \mid A) P ( B ∣ A ) は確率ですが、A A A を動かす関数 L ( A ; B ) L(A;\, B) L ( A ; B ) は確率分布ではありません。∑ A P ( B ∣ A ) \sum_{A} P(B \mid A) ∑ A P ( B ∣ A ) は 1 になりません。医療検査の例では 0.90 + 0.09 = 0.99 0.90 + 0.09 = 0.99 0.90 + 0.09 = 0.99 です。B B B を固定して A A A を動かしているため、A A A の空間での総和が 1 になる保証はありません。