元データ分析の会社で働いていた人の四方山話スマートフォン用の表示で見る元データ分析の会社で働いていた人の四方山話この広告は、90日以上更新していないブログに表示しています。論文でm次元ガウス分布を例にして評価実験が行われていたので、それを試してみるために色々と調べていたら、やっぱりありました。自分でゴリゴリ書くのが若干めんどくさかったのでありがたいですね。多変量解析用パッケージのMASSのmvrnormを使います。第2引数に各次元の平均(ベクトル)第3引数に共分散行列(行列)を与えてやります。第3引数が若干混乱しそうですが、共分散行列は対角成分に分散が、それ以外は例えば1行2列目には1番目と2番目の次元の相関関係が入ります。各次元が独立している状況を考えたければ単位行列を与えればよいわけです。第2引数の長さと、共分散行列なので当然正方行列である第3引数の行と列は等しくないといけません。論文では、N(0,I)とあったので、第2引数は全て0のベクトル、第3引数は第2引数の長さと同じ行と列を持つ単位行列を用意します。Sigma)これで5次元のガウス分布に従うデータができます。-0.6248058こんな感じの結果になりました。これが正しいのか調べてみます。1.239451という感じで、全ての次元がだいたい平均0、標準偏差1くらいになっていますね。変数間の相関関係も気になるので相関係数を求めるcor()を使って調べてみます。0.178124336924986"ということで若干相関関係出とるやないかい!となるのですが、よく考えると今データは10個と統計量を求めるにはちょっと少ないですね。データ数を1000にしたときの各次元間の相関係数です。-0.0500139428051309"とほぼ0になっており、次元間が無相関であるという感じになってますね。Entriesはてなブログをはじめよう!rindai87さんは、はてなブログを使っています。ブログを報告するスターをつけました引用をストックしましたストック一覧を見る閉じる引用するにはまずログインしてくださいログイン閉じる引用をストックできませんでした。再度お試しください閉じる限定公開記事のため引用できません。