「続?わかりやすいパターン認識」 第12章 ディリクレ過程混合モデルによるクラスタリング（前半 : 12.1 ）
5 likes5,301 views
2015.09.29 「続?わかりやすパターン認識」読書会 #13 の 発表資料。
More Related Content
What's hot (20)
Similar to 「続?わかりやすいパターン認識」 第12章 ディリクレ過程混合モデルによるクラスタリング（前半 : 12.1 ） (16)
「続?わかりやすいパターン認識」 第12章 ディリクレ過程混合モデルによるクラスタリング（前半 : 12.1 ）
- 2. 自己紹介 ? Twitter : @aich_08_ ? 所属組織 : 匿名知的集団ホクソエム（ピンク担当） ? 職種 : コンサル ? 趣味 : トレーニング，登山 ? 全然痩せない????orz 2 ラテアートではなく， ただのカプチーノだよ。
- 3. これまでのおはなし ? 第11章では、分割の確率モデルの具体例として， ホップの壺モデルと中華料理店過程（Chinese restaurant process : CRP)を取り上げた． ? この2モデルの分割ルールは， ディリクレ分布を無限次元に拡張することで導出できることを確認した． ? 第12章では，クラスタ数が未知の場合のクラスタリング方法として， ディリクレ過程混合モデルによるクラスタリングを紹介する． ? ディリクレ過程混合モデルでは，クラスタ数 ? はディリクレ過程により決定される． 3
- 4. これまでのおはなし（つづき） 表記のおさらい（p. 214） パターン集合 : ? = {?1, … , ? ?} パターン集合 ? に対応するクラスタ集合 : ? = {?1, … , ? ?} 各クラスタに対するパラメータ : ? = {?1, … , ? ?} 4
- 5. ディリクレ過程混合モデルによる生成モデル ? ディリクレ過程混合モデルによる生成モデルは以下のように 記される．ここで，? は集中度パラメータ， ?0(?)は基底分布を表す. 5
- 6. CRPに基づくディリクレ過程混合モデル ? 中華料理店過程（Chinese Restaurant Process : CRP）によるディリクレ過程 混合モデルの生成プロセスは以下のように記される． ? ここで、CPR(?)は集中度パラメータ ? のCPRを、 ?は?0 ? ハイパーパラ メーを表す． 6
- 7. 生成モデルに対する学習 ? いま、生成モデルからある観測パターン ? が得られたとする． ? 以下に示される事後分布を最大化することにより，クラスタ ? およびパラ メータ ? を決定する． ■ 所属クラスタとパラメータを決定する場合 7
- 8. 生成モデルに対する学習（つづき） ■ 所属クラスタのみを決定する場合 → ? = {?1, … , ? ?} の可能なすべての組合せを評価するのは，その組合せ が膨大ゆえ，直接求めるのは困難である（テキスト p.215， 問題 2） ? ギブスサンプリング（A.5 参照） 8
- 9. この問題に対するギブスサンプリングの考え方 ① ? 個のパターンのうち，任意の1コ ? ?を選択し，これを最後に入力 されたパターンとみなす． ② 当該パターンの所属クラスタ ? ?を未定とみなし，確率的に決定さ れるクラスタを，当該パターンの新たな所属クラスタとして更新す る． ③ ①，②の操作を ? 個すべてのパターンについて繰り返すことにより， 確率分布?(?1, … , ? ?|?1, … , ? ?, ?) に従う ?1, … , ? ? が得られる． 9
- 10. ディリクレ過程混合モデルにおける ギブスサンプリングの更新式の導出 10 ? 更新式の導出に先立ち，? から ? ? を除いた集合 ???および ? から ? ? を 除いた集合 ???を定義する.
- 11. ディリクレ過程混合モデルにおける ギブスサンプリングの更新式の導出（つづき） ? ギブスサンプリングにより ? ? を確率的に生成するためには，???， ???，? が既知の下で ? ? が与えられたときに，? ? = ??となる事後確 率を計算すれば良い． 11
- 12. ? 所属クラスタとパラメータを決定する方法，式（11.5）の解法を示す. ? まず，この場合の更新式（12.10）の計算方法を整理する． ? 右辺の第1項は，既存クラスタである場合 ? ? = ?? ，新規クラスタで ある場合 ? ? = ? ??? それぞれについて以下のように書き換えられる． クラスタリング法 1 : 所属クラスタとそのパラメータの決定 12
- 13. ? 式（12.10）の右辺の第2項は，式（11.42）から以下で書き換えられる． ? ここで，?? ′ は ? 番目のパターンを除く，すべてのパターンのうち，クラ スタ ?? に所属するパターンの数を表す． クラスタリング法 1 : 所属クラスタとそのパラメータの決定（つづき） 13
- 14. ? 式（12.11），（12.12），（12.13）を式（12.10）に代入することにより，次式を得る． ? パターン ? ? の所属クラスタを更新するためには，式（12.14）の上段を現時点 のクラスタ総数 ? に渡り計算し，下段の計算と合わせて正規化することにより， ? ? = ?? ，? ? = ? ??? となる事後確率を各々計算する．そして，その結果を用いて ? ? の値を決定する． クラスタリング法 1 : 所属クラスタとそのパラメータの決定（つづき） 14
- 15. ? 以上の操作により， ?1, … , ? ? の値が決まれば，事前分布?0(?) およ びクラスタ ?? に所属するパターン対する尤度 ?(??|{? ?; ? ? ∈ ?? }) を用いて，ベイズの定理により次式を得，この事後分布に基づいて ? 個のクラスタのパラメータ ??(? = 1, … , ?) を確率的に決定し，更新 する． クラスタリング法 1 : 所属クラスタとそのパラメータの決定（つづき） 15
- 16. ? 以上の操作（所属クラスタの更新，パラメータの更新）を繰り返し実 行し，式（11.4）の事後分布が最大となったときの ?, ? の値が最適解 である． クラスタリング法 1 : 所属クラスタとそのパラメータの決定（つづき） 16
- 17. クラスタリング法 1 のアルゴリズム 入力 : ? 個の観測パターン ?1, … , ? ? 出力 : 各観測パターンの所属クラスタを表す潜在変数の値 ?1, … , ? ?， 総クラスタ数 ?， クラスタのパラメータ ?1, … , ? ?． 17
- 18. クラスタリング法 1 のアルゴリズム（つづき） Step 1 : 初期設定 ? 潜在変数 ?1, … , ? ? を初期化する． このときのクラスタ数を ? とし，クラスタ ?? に所属する パターンの数を ?? とする（? = 1, … , ?）． ? クラスタパラメータ ?1, … , ? ? を初期化する． ? 事後確率の最大値 ???? ← 0 と初期化する． 18
- 19. クラスタリング法 1 のアルゴリズム（つづき） Step 2 : 所属クラスタの更新 以下を? = 1, … , ?に対して実行する． ? パターン ? ? の所属クラスタ ? ? （ = ?? と仮定する）を更新するため， ? ? を現在の所属クラ スタから除外し， ?? ← ?? ? 1 とする． ? これにより ?? = 0 となり，空きクラスタが発生した場合は，このからクラスタを除去すべく， ? ← ? ? 1 とするとともに， ? 以降のクラスタインデックスを全て一つ減らす．また，それに対 応するパラメータのインデックスも更新する． ? 式（12.14） を用いて ? ? の値を確率的に決定する．具体的には， ? ? = ?? ? = 1, … , ? に 対 して式（12.14）の右辺上段の式を計算し， ? ? = ? ??? に対して下段の式を計算する．これ ら(? + 1)個の値の比の確率で ? ? を決定する．また， ? ? = ? ???なら， ? ?+1 ← 1とするとと もに， ? ← ? + 1として総クラスタ数を更新する． 19
- 20. クラスタリング法 1 のアルゴリズム（つづき） Step 3 : 各クラスタのパラメータの更新 ? Step 2で得た ? = {?1, … , ? ?} の値に基いて，式（12.15）により，各ク ラスタのパラメータ ??を確率的に決定し，更新する． 20
- 21. クラスタリング法 1 のアルゴリズム（つづき） Step 4: 事後確率最大化 ? 現時点での ? = {?1, … , ? ?} ，および ?1, … , ? ? の値を用いて式 （11.4）により，事後確率 ? を計算する．具体的には，ベイズの定理 より次式（12.16）を計算すればよい． ?(?) は式（11.11）のイーウェン スの抽出公式で計算する． 21
- 22. クラスタリング法 1 のアルゴリズム（つづき） Step 4: 事後確率最大化（つづき） ? 以下の処理を行う． (1) ? > ???? なら，以下の更新を行い， ? ???? ← ? ? ? ← {?1, … , ? ?} (2) それ以外なら，これまでの????および ? をそのまま保持する． 22
- 23. クラスタリング法 1 のアルゴリズム（つづき） Step 5: 終了判定 ? 以上の処理を繰り返し， ????の更新されない状態が十分継続した 場合，ギブスサンプリングの処理を打ち切り， ?，?，および現時点の クラスタ総数 ? を出力して終了する．さもなければ，Step 2に戻る． 23
- 24. ? 所属クラスタのみを決定する方法，式（11.7）の解法を示す. ? まず?この場合の更新式（12.10）の計算方法を整理する． ? 式（12.8），（12.9），（12.10）から ? を削除することにより次式を得る． クラスタリング法 2 : 所属クラスタのみ決定 24
- 25. ? 式（12.17）右辺の第1項は，次式で計算できる． ? ここで，事後確率 ?(??|???) はベイズの定理より次式となる． クラスタリング法 2 : 所属クラスタのみ決定（つづき） 25
- 26. ? ここで，??? のうち ?? に所属しないパターンは ?? に依存しないので， ? 式（12.19）に代入すると，次式を得る． クラスタリング法 2 : 所属クラスタのみ決定（つづき） 26
- 27. ? ? ? = ? ??? のとき，式（12.21）は次式になる． ? 一方，更新式（12.17）の右辺の第2項 ?(? ?|???) は ，式（12.13）によ り求まる． クラスタリング法 2 : 所属クラスタのみ決定（つづき） 27
- 28. ? 以上を整理すると，次式を得る． ? パターン ? ? の所属クラスタを更新するためには，式（12.23）の上段を現時点 のクラスタ総数 ? に渡り計算し，下段の計算と合わせて正規化することにより， ? ? = ?? ，? ? = ? ??? となる事後確率を各々計算する．そして，その結果を用い て ? ? の値を決定する． クラスタリング法 2 : 所属クラスタのみ決定（つづき） 28
- 29. クラスタリング法 2 のアルゴリズム 入力 : ? 個の観測パターン ?1, … , ? ? 出力 : 各観測パターン所属クラスタを表す潜在変数の値 ?1, … , ? ?， 総クラスタ数 ? 29
- 30. クラスタリング法 2 のアルゴリズム（つづき） Step 1 : 初期設定 ? 潜在変数 ?1, … , ? ? を初期化する． このときのクラスタ数を ? とし，クラスタ ?? に所属する パターンの数を ?? とする（? = 1, … , ?）． ? 事後確率の最大値 ???? ← 0 と初期化する． 30
- 31. クラスタリング法 2 のアルゴリズム（つづき） Step 2 : 所属クラスタの更新 以下を? = 1, … , ?に対して実行する． ? パターン ? ? の所属クラスタ ? ? （ = ?? と仮定する）を更新するため， ? ? を現在の所属クラ スタから除外し， ?? ← ?? ? 1 とする． ? これにより ?? = 0 となり，空きクラスタが発生した場合は，このからクラスタを除去すべく， ? ← ? ? 1 とするとともに， ? 以降のクラスタインデックスを全て一つ減らす．また，それに対 応するパラメータのインデックスも更新する． ? 式（12.23） を用いて ? ? の値を確率的に決定する．具体的には， ? ? = ?? ? = 1, … , ? に 対 して式（12.23）の右辺上段の式を計算し， ? ? = ? ??? に対して下段の式を計算する．これ ら(? + 1)個の値の比の確率で ? ? を決定する．また， ? ? = ? ???なら， ? ?+1 ← 1とするとと もに， ? ← ? + 1として総クラスタ数を更新する． 31
- 32. クラスタリング法 2 のアルゴリズム（つづき） Step 3: 事後確率最大化 ? 現時点での ? = {?1, … , ? ?} の値を用いて式（11.6）により，事後確率 ? を計算する．具体的には，ベイズの定理より次の式（12.24）を計算 すればよい． ?(?) は式（11.11）のイーウェンスの抽出公式で計算す る． 32
- 33. クラスタリング法 2 のアルゴリズム（つづき） Step 3: 事後確率最大化（つづき） ? 以下の処理を行う． (1) ? > ???? なら，以下の更新を行い， ? ???? ← ? ? ? ← {?1, … , ? ?} (2) それ以外なら，これまでの????および ? をそのまま保持する． 33
- 34. クラスタリング法 2 のアルゴリズム（つづき） Step 4: 終了判定 ? 以上の処理を繰り返し， ????の更新されない状態が十分継続した 場合，ギブスサンプリングの処理を打ち切り， ? および現時点のクラ スタ総数 ? を出力して終了する．さもなければ，Step 2に戻る． 34