科学の箱

科学・IT・登山の話題

機械学習

Pythonではじめる機械学習 – GMM

投稿日:

k-Meansの課題

  • クラスタ間で微妙な位置関係にあるデータについて特定クラスタに分類された際の不確実の度合いを知ることができない。
  • 例えば51%の確率でクラスタAであり、49%の確率でクラスタBかもしれないがそれはモデルからは判明しない。

クラスタリングの精度を知るための一つの方法として、クラスタの中心から円を描き、その円に入らないデータはクラスタから除外する

  • この方法ではクラスタの形が円であることが前提である。
  • 楕円等になった場合には、正しくクラスタリングできない

PCAにより軸を変える方法が考えられるが、しかしPCAによりデータが円状になることは保証されていない。

k-meansの問題は2つである。

  • クラスタの形状は柔軟性がない
  • 確率によるデータのクラスタへの割り当てができない

この問題を解決する方法としては2つある。

  1. クラスタの形状に楕円を用いる
  2. データについてすべてのクラスタ中心からの距離を測定する

ガウス分布と混合ガウス分布

ガウス分布とは正規分布のことである。正規分布は統計で最も基本となる確率分布となり、中心となる平均と広がりである分散となる2つのパラメータをもつ曲線である。

混合ガウス分布とは平均と分散が異なる複数のガウス分布を組み合わせた分布となる。データの散らばり具合、その確率を混合ガウス分布で表現する。複数の正規分布を用いるので正規分布より柔軟な表現が可能になる。

混合ガウスモデル

混合ガウスモデルの学習では、「データはすでに存在する確率分布=混合ガウス分布により生成したとみなす」を前提にする。その前提に立ち、データから混合ガウス分布を求める。

ガウスモデルは平均と分散をパラメータとする。分散を変えることで楕円形のクラスタを生成することが可能になる。

 

混合ガウスモデルの求め方は以下のようになる。

  1. 各クラスタごとに初期のガウス分布を分布を生成する(初期パラメータを平均と分散に設定する)
  2. 各データについてそれぞれのガウス分布の重みを計算する
  3. 得られた重みによりガウス分布を更新する。
  4. ガウス分布の変化が十分小さくなるまで繰り返す。

実装

covariance_typeによるクラスタ形状の指定

デフォルト値は”diag”。これは楕円であるが軸は軸に依存する。”spherial”ではクラスタの形状は球形となり、k-meansと似た結果を返す。”full”は最も柔軟な楕円を返す。ただし計算量は多くなる。

n_componentsの大きさ

n_componentsによりデータポイントに合わせた柔軟なクラスタを得ることができる。ただし大きすぎても意味がなくなる。AICおよびBICにより最適なn_componentsを取得する

メタ情報

inarticle



メタ情報

inarticle



-機械学習

執筆者:


comment

メールアドレスが公開されることはありません。 * が付いている欄は必須項目です

関連記事

no image

データサイエンス超入門

シンプソンのパラドックス レコメンドエンジン コンテンツベースフィルタリング 強調 アイテム ユーザー ビジネスにおけるデータ分析の手順 課題→ビジネスケース→仮説→分析→用途 データ分析の手順 デー …

no image

売り上げデータの分析

利益 = 売り上げ – コスト 売り上げ = 客数 x 客単価 コスト ≒ 人件費 + 廃棄コスト 客数 客単価 = Σ 品物i x 購入数 客数を増やす方法 来てもらう方法 安売りキャン …

no image

スピアマンの順位相関

スピアマンの順位相関を利用して県のデータを使った相関係数を出してみる。 持家普及率を自動車普及率のランキングは、土地がせまいほど低く、有効な土地が広いほど高くなるために順位相関があると考えられる。 持 …

no image

仕事で始める機械学習 – 3.学習結果を評価しよう – 指標

主な指標 モデル構築後に確認する主な指標は4つある。 正解率 適合率 再現率 F値 正解率 $$ 正解率 = \frac{TP + TN}{TP+FP+TN+FN} $$ 正解率は全データ数に対する正 …

no image

Excelによる単回帰分析

エクセルでの単回帰分析をについて手順をまとめる。エクセルの使い方にあわせて分析結果の解釈の仕方も書いておく。 データは勉強時間に対する点数とする。これは作ったデータであるが、線形回帰にフィットするよう …

2019年11月
« 10月   12月 »
 123
45678910
11121314151617
18192021222324
252627282930  

side bar top



アーカイブ

カテゴリー