【R】2次元の正規分布に従うデータ

2次元の正規分布に従うデータの生成方法です。例として、統計検定2級　2015年11月　問３　[8]　の散布図を作ってみます。

以下のデータを作ります。

　国語　：　平均56.0　分散236.6
　英語　：　平均59.1　分散170.1
　国語と英語の共分散　：　133.1

データを生成するために、MASSパッケージのmvrnorm関数を用います。mvrnorm関数は正規分布に従う2つの変数の平均値と分散共分散行列を入力して使います。

library(MASS)	　　  				                     # MASSパッケージの読込
Mu <- c(56, 59.1)					                     # 平均値ベクトル
Si <- matrix(c(236.6, 133.1, 133.1, 170.1), ncol=2)    # 分散共分散行列
dat <- mvrnorm(300, Mu, Si)			                     # データの生成
plot(x=dat[,1], y=dat[,2], xlab="国語", ylab="英語", xlim=c(0, 100), ylim=c(0, 100))              # プロット

こちらのサイトを参考にさせていただきました。