2019-01-01から1年間の記事一覧

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

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

jagsで潜在混合モデル

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

rTensorで多次元主成分分析をやってみた

RのパッケージのrTensorを使って多次元主成分分析を行ってみた。 多次元主成分分析のやり方は、 http://www.iri.pref.kumamoto.jp/sgk/2011/cdrom/thesis/04_hall/04_session/444.pdf からダウンロードできる日本語文献のpdfファイル ”多次元主成分分析によ…