2019-07-01から1ヶ月間の記事一覧

Rで見る主成分分析と特異値分解の関係

多次元主成分分析を行うにあたって、主成分分析と特異値分解の関係を復習してみた。 まず、データを作成する。 # x:1000個の正規分布N(0, 0.2)に従う乱数# y : xを3倍し、そこにN(0, 0.3)に従うノイズを加えたもの # x2, y2: x, yそれぞれについて平均を引…

jagsで潜在混合モデル

ベイス統計で実戦モデリング(北大路書房)の第6章 潜在混合モデルの練習で、次のコードを作成してみた。 8個のデータが4個づつ違う正規分布に従うとして、データがどちらの正規分布に従うか、またその時の正規分布の平均を推測する。ただし、標準偏差はど…