読む編集履歴表示

EMアルゴリズム

ウィキペディアから、無料の百科事典

応用数学 > 統計学 > EMアルゴリズム
データサイエンス > 機械学習 > EMアルゴリズム

機械学習およびデータマイニング

問題分類クラスタリング回帰異常検知相関ルール（英語版）強化学習構造化予測（英語版）特徴量設計（英語版）表現学習（英語版）オンライン学習（英語版）半教師あり学習（英語版）教師なし学習ランキング学習（英語版）文法獲得（英語版）
教師あり学習（分類 • 回帰）決定木（英語版）アンサンブル（バギング、ブースティング、ランダムフォレスト） k-NN 線形回帰単純ベイズニューラルネットワークロジスティック回帰パーセプトロン関連ベクトルマシン (RVM)（英語版）サポートベクトルマシン (SVM)
クラスタリング BIRCH（英語版）階層的（英語版） k平均法期待値最大化法 (EM) DBSCAN OPTICS（英語版）平均値シフト（英語版）
次元削減因子分析 CCA ICA LDA（英語版） NMF（英語版） PCA t-SNE
構造化予測（英語版）グラフィカルモデルベイジアンネットワーク CRF HMM
異常検知 k-NN 局所外れ値因子法
ニューラルネットワークオートエンコーダディープラーニング DeepDream 多層パーセプトロン RNN LSTM GRU 制約ボルツマンマシン（英語版） SOM CNN
強化学習 TD学習 Q学習 SARSA
理論偏りと分散のトレードオフ計算論的学習理論（英語版）経験損失最小化（英語版）オッカム学習（英語版） PAC学習統計的学習（英語版） VC理論（英語版）
学会・論文誌等 NIPS（英語版） ICML（英語版） ML（英語版） JMLR（英語版） ArXiv:cs.LG
全般統計学および機械学習の評価指標
Category:機械学習 Category:データマイニング
表話編歴

EMアルゴリズム（英: expectation–maximization algorithm）とは、統計学において、確率モデルのパラメータを最尤推定する手法の一つであり、観測不可能な潜在変数に確率モデルが依存する場合に用いられる。EM法、期待値最大化法（きたいちさいだいかほう）^[1]^[2]とも呼ばれる。その一般性の高さから、機械学習、音声認識、因子分析など、広汎な応用がある^[1]。

EMアルゴリズムは反復法の一種であり、期待値（英: expectation, E）ステップと最大化（英: maximization, M）ステップを交互に繰り返すことで計算が進行する。Eステップでは、現在推定されている潜在変数の分布に基づいて、モデルの尤度の期待値を計算する。Mステップでは、E ステップで求まった尤度の期待値を最大化するようなパラメータを求める。M ステップで求まったパラメータは、次の E ステップで使われる潜在変数の分布を決定するために用いられる。

概要

セッティング・目標

今、2値 $x$ 、 $z$ を取る確率分布があり、その確率分布の確率密度関数 $p(x,z|\theta )$ が未知の母数 $\theta \in \mathbb {R} ^{m}$ によりパラメトライズされているとする。ここで $\mathbb {R}$ は実数全体の集合を表す。

そして $p(x,z|\theta )$ に従って標本 $(x_{1},z_{1}),\ldots ,(x_{n},z_{n})$ を独立に抽出したものの、何らかの事情で $Z=(z_{1},\ldots ,z_{n})$ の値は観測できず、 $X=(x_{1},\ldots ,x_{n})$ だけが観測できたとする。実応用上は例えば、 $\theta =(\theta _{1},\theta _{2})$ という形をしており、まず観測不能な $z_{i}\sim p_{1}(z|\theta _{1})$ が選ばれた後、 $z_{i}$ に依存して観測可能な $x_{i}\sim p_{2}(x|\theta _{2},z_{i})$ が選ばれる、といったケースにEMアルゴリズムが使われる事が多いが、必ずしもこのケースにあてはまらなくてもよい。

簡単の為、記号を混用して $X$ 、 $Z$ の同時確率分布の確率密度関数も $p(X,Z|\theta )$ と書く。以下では $Z$ が離散変数の場合について説明するが、 $Z$ が連続変数の場合も総和を積分に置き換える以外は同様である^[3]。

このような状況において母数 $θ$ を最尤推定する事が我々の目標である。しかし $Z$ を知らない場合の $X=(x_{1},\ldots ,x_{n})$ に関する対数尤度

\ell (\theta |X):=\log p(X|\theta )=\log \sum _{Z}p(X,Z|\theta )

を最大値を直接計算するのは一般には簡単ではない。

EMアルゴリズムは、反復法により、数列 ${\hat {\theta }}^{(t)}$ で対数尤度 $\ell ({\hat {\theta }}^{(t)}|X)$ が単調非減少であるものを作るアルゴリズムである。最尤推定量を ${\hat {\theta }}_{\mathrm {MLE} }$ とすると、

\ell ({\hat {\theta }}_{\mathrm {MLE} }|X)\geq \ell ({\hat {\theta }}^{(t)}|X)

である事から、 $\ell ({\hat {\theta }}_{\mathrm {MLE} }|X)$ が有限であれば $\ell ({\hat {\theta }}^{(t)}|X)$ の単調性より $\ell ({\hat {\theta }}^{(t)}|X)$ は必ず収束する。

アルゴリズム

EMアルゴリズムでは、以下の手順により数列 ${\hat {\theta }}^{(0)},{\hat {\theta }}^{(1)},\ldots$ を作る^[3]。

初期値 ${\hat {\theta }}^{(0)}$ を（何らかの方法で）選ぶ。
$t=0,1,\ldots$ $t=0,1,\ldots$ に対して以下を実行する
- E ステップ: $p(Z|X,{\hat {\theta }}^{(t)})$ を求める。
- M ステップ: ${\hat {\theta }}^{(t+1)}={\underset {\theta }{\operatorname {arg\,max} }}\ Q(\theta |{\hat {\theta }}^{(t)})\,$ を求める。

ここでQは対数尤度関数 $\log p(X,Z|\theta )$ の $Z$ に関する条件付き期待値

Q(\theta |\theta ^{(t)}):=\operatorname {E} _{Z|X,{\hat {\theta }}^{(t)}}{\big [}\log p(X,Z|\theta ){\big ]}=\sum _{Z}p(Z|X,{\hat {\theta }}^{(t)})\log p(X,Z|\theta )\,

である。実応用上は、 ${\hat {\theta }}^{(t)}$ の値が十分小さくなったと判定する何らかの条件を事前に定めておき、その条件を満たしたら上述のループを終了する。ループを終了する条件は、パラメータ値や対数尤度関数を使って定められる^[3]。

留意点

EステップとMステップの切れ目は書籍により異なるので注意が必要である。本項では次節の議論と整合性をとる為に文献^[3]の切れ目に従ったが、文献^[4]では $Q(\theta |{\hat {\theta }}^{(t)})$ を計算する所までがEステップであり、 $Q(\theta |{\hat {\theta }}^{(t)})$ の $\operatorname {arg\,max}$ を取るところだけがMステップである。

ステップの名称「E」と「M」はそれぞれExpectation（期待値）、Maximization（最大化）の略であり^[4]、文献^[4]のようにEステップで $Q(\theta |{\hat {\theta }}^{(t)})$ を求める為に期待値を計算し、Mステップで $Q(\theta |{\hat {\theta }}^{(t)})$ の $\operatorname {arg\,max}$ を取るところに名称の由来がある。

動作原理

EMアルゴリズムで我々が求めたいのは、 $X=(x_{1},\ldots ,x_{n})$ を観測した際における対数尤度

\ell (\theta |X):=\log p(X|\theta )

を最大化する母数 $\theta$ であった。EMアルゴリズムの動作原理を説明する為、以下のような汎関数を考える：

{\mathcal {L}}(q,\theta ):=\sum _{Z}q(Z)\log {p(X,Z|\theta ) \over q(Z)}

　　...(Eq.1)

ここで $q(Z)$ は任意の確率密度関数である。 $p_{X,\theta }(Z):=p(Z|X,\theta )$ とすると、 $p(Z|X,\theta )p(X|\theta )=p(X,Z|\theta )$ より、カルバック・ライブラー情報量

\mathrm {KL} (q||p_{X,\theta })=-\sum _{Z}q(Z)\log {p(Z|X,\theta ) \over q(Z)}

を使って

{\mathcal {L}}(q,\theta )=\ell (\theta |X)-\mathrm {KL} (q||p_{X,\theta })

　...(Eq.2)

と書ける事が分かる。カルバック・ライブラー情報量が常に非負である事（ギブスの不等式）から、

\ell (\theta |X)\geq {\mathcal {L}}(q,\theta )

であるので、 ${\mathcal {L}}(q,\theta )$ は $\ell (\theta |X)$ の下限になっている。EMアルゴリズムはこの下限 ${\mathcal {L}}(q,\theta )$ を逐次的に改善していくことで、 $\ell (\theta |X)$ を可能な限り最大化するアルゴリズムである。すなわち、EステップとMステップは以下のように書き換えられる事を示す事ができる^[3]：

E ステップ: ${\hat {q}}^{(t)}={\underset {q}{\operatorname {arg\,max} }}{\mathcal {L}}(q,{\hat {\theta }}^{(t)})$ を求める。
M ステップ: ${\hat {\theta }}^{(t+1)}={\underset {\theta }{\operatorname {arg\,max} }}{\mathcal {L}}({\hat {q}}^{(t)},\theta )$ を求める。

この事実から対数尤度 $\ell ({\hat {\theta }}^{(t)}|X)$ の単調非減少性が明らかに従う。（但し反復法の常として、初期値しだいでは尤度の最大点ではない極大点に到達してそこで停止する可能性がある。）

証明

本節ではEステップ、Mステップが上述のように書き換えられることを示す。本節の証明は文献^[3]を参考にした。

Eステップの証明

カルバック・ライブラー情報量 $\mathrm {KL} (q||p_{X,\theta })$ が最小値0になるのは $q=p_{\theta ,X}$ の場合だけであった事から、(Eq.2)より ${\mathcal {L}}(q,\theta )$ は

q(Z)=p(Z|X,\theta )

が満たされる場合に最大値を取る。すなわちEMアルゴリズムにおけるEステップは、 $\theta ={\hat {\theta }}^{(t)}$ を固定したままの状態で、 ${\mathcal {L}}(q,\theta )$ を最大化する $q$ である

{\hat {q}}^{(t)}:=p_{X,{\hat {\theta }}^{(t)}}={\underset {q}{\operatorname {arg\,max} }}{\mathcal {L}}(q,{\hat {\theta }}^{(t)})

を求めるステップである。

Mステップの証明

${\mathcal {L}}(q,\theta )$ の定義式(Eq.1)に ${\hat {q}}^{(t)}=p_{X,{\hat {\theta }}^{(t)}}$ を代入すると、

{\mathcal {L}}({\hat {q}}^{(t)},\theta )=\sum _{Z}p(Z|X,\theta ^{(t)})\log {p(X,Z|\theta ) \over p(Z|X,\theta ^{(t)})}=Q(\theta |\theta ^{(t)})-H_{X,\theta ^{(t)}}(Z)

が成立し（ここで $H_{X,\theta ^{(t)}}(Z)=\textstyle \sum _{Z}p(Z|X,\theta ^{(t)})\log p(Z|X,\theta ^{(t)})$ は条件付きエントロピー）、上式右辺第二項は $θ$ に依存しないので、

{\hat {\theta }}^{(t+1)}={\underset {\theta }{\operatorname {arg\,max} }}Q(\theta |{\hat {\theta }}^{(t)})={\underset {\theta }{\operatorname {arg\,max} }}{\mathcal {L}}(p_{X,{\hat {\theta }}^{(t)}},\theta )

が成立する。

一般化

EMアルゴリズムは観測データの対数尤度を、E ステップとM ステップの繰り返しにより最大化するアルゴリズムであるので、正確にはlog-EMアルゴリズムというべきものである。log関数にはα-logとよばれる一般化された対数があるので、それを用いるとlog-EMを特例として含むアルゴリズムを作り上げることができる。ただし、この場合は尤度ではなくてα-log尤度比とαダイバージェンスを用いて基本等式を導くことになる。このようにして得られたものがα-EMアルゴリズム ^[5] であり、log-EMアルゴリズムをサブクラスとして含んでいる。α-EMアルゴリズムは適切なαを選ぶことにより、log-EMアルゴリズムよりも高速になる。また、log-EMが隠れマルコフモデル推定アルゴリズム（Baum-Welchアルゴリズム）を含んでいるように、α-EMアルゴリズムから高速なα-HMMアルゴリズムを得ることができる。 ^[6]

歴史

EMアルゴリズムは、アーサー・デンプスター（英語版）、ナン・レアード（英語版）、ドナルド・ルービンによる1977年の論文^[7]で導入され、その名が付けられた。彼らは、EMアルゴリズムがほかの複数の著者によって「特殊な文脈でなんども提案されてきた」("proposed many times in special circumstances") ことを述べた上で、EMアルゴリズムの一般化を行い、その背後にある理論を追求した。

本来のEMアルゴリズムでは、期待値の評価において潜在変数のとりうる値すべてを列挙することが必要なため、効率的に扱える分布が限られていた。しかしその後、マルコフ連鎖モンテカルロ法や変分ベイズ法（英語版）が考案されたことにより、より一般の分布でも現実的な時間での計算が可能になった^[1]^[8]。

脚注

[脚注の使い方]

1 2 3 計算統計I, p. 130.
↑ 計算統計I, p. 157.
1 2 3 4 5 6 #PRML pp.156, 164-171
1 2 3 #ESL pp.316-317.
↑ Matsuyama, Yasuo (2003). “The α-EM algorithm: Surrogate likelihood maximization using α-logarithmic information measures”. IEEE Transactions on Information Theory 49 (3): 692-706.
↑ Matsuyama, Yasuo (2011). “Hidden Markov model estimation based on alpha-EM algorithm: Discrete and continuous alpha-HMMs”. International Joint Conference on Neural Networks: 808-816.
↑ Dempster, A.P., Laird, N.M., Rubin, D.B., (1977). “Maximum Likelihood from Incomplete Data via the EM Algorithm”. Journal of the Royal Statistical Society. Series B (Methodological) 39 (1): 1–38. JSTOR2984875. MR 0501537.
↑ 計算統計I, p. 163.

参考文献

引用文献

Trevor Hastie; Robert Tibshirani; Jerome Friedman (2014/6/25). 統計的学習の基礎―データマイニング・推論・予測―. 共立出版. ISBN 978-4-320-12362-5
C.M. ビショップ (2012/2/29). パターン認識と機械学習下 (ベイズ理論による統計的予測). 丸善出版. ISBN 978-4621061244
汪金芳、手塚集、上田修功、田栗正章、樺島祥介、甘利俊一、竹村彰通、竹内啓、伊庭幸人『計算統計 I ―確率計算の新しい手法』 11巻〈統計科学のフロンティア〉、2003年。ISBN 4000068512。

その他の参考文献

Robert Hogg, Joseph McKean and Allen Craig. Introduction to Mathematical Statistics. pp. 359-364. Upper Saddle River, NJ: Pearson Prentice Hall, 2005.
The on-line textbook: Information Theory, Inference, and Learning Algorithms, by David J.C. MacKay includes simple examples of the E-M algorithm such as clustering using the soft K-means algorithm, and emphasizes the variational view of the E-M algorithm.
A Gentle Tutorial of the EM Algorithm and its Application to Parameter Estimation for Gaussian Mixture and Hidden Markov Models, by Jeff Bilmes includes a simplified derivation of the EM equations for Gaussian Mixtures and Gaussian Mixture Hidden Markov Models.
Variational Algorithms for Approximate Bayesian Inference, by M. J. Beal includes comparisons of EM to Variational Bayesian EM and derivations of several models including Variational Bayesian HMMs.
The Expectation Maximization Algorithm, by Frank Dellaert, gives an easier explanation of EM algorithm in terms of lowerbound maximization.
The Expectation Maximization Algorithm: A short tutorial, A self contained derivation of the EM Algorithm by Sean Borman.
The EM Algorithm, by Xiaojin Zhu.
Geoffrey J. McLachlan and Thriyambakam Krishnan: The EM Algorithm and Extensions, Wiley series in probability and statistics, John Wiley & Sons, Inc., ISBN 0-471-12358-7 (1997).
Geoffrey J. McLachlan and Thriyambakam Krishnan: The EM Algorithm and Extensions, 2nd Edition, Wiley & Sons Inc., ISBN 978-0-471-20170-0 (February 2008).　上記の改訂第2版。
Kenneth Lange: MM Optimization Algorithms, SIAM, ISBN 978-1-611974-39-3 (2016). ※ "EM" アルゴリズムの一般化として　"MM algorithm" を提唱している．

和書

小西貞則・越智義道・大森裕浩:「計算統計学の方法 ―ブートストラップ，EMアルゴリズム，MCMC―」、朝倉書店(シリーズ：予測と発見の科学、5)、ISBN 978-4-254-12785-0　(2008年3月25日)。
関原謙介：「ベイズ信号処理」、共立出版、ISBN 978-4-320-08574-9 (2015年4月11日)。
関原謙介：「ベイズ推論の基礎と信号処理への応用」
黒田正博:「EMアルゴリズム」、共立出版(シリーズ：統計学One Point、18巻)、ISBN 978-4-320-11269-8、(2020年7月31日)。

この項目は、統計学に関連した書きかけの項目です。この項目を加筆・訂正などしてくださる協力者を求めています（プロジェクト:数学／Portal:数学）。