行列の推測統計をしたい。

Question

平均[168, 71], 分散共分散⾏列[[81,68],[68,250]]の2次元正規分布に従う乱数を1000個⽣成した上で
• X軸1列⽬,Y軸2列⽬データの散布図を描け．
• 1列, 2列⽬の平均を求めよ．
• 1列,2列⽬の相関係数を求めよ．
この上記の三つを示すプログラムを作りたいのですが、作り方がさっぱりわからないので、教えていただきたいです。

Answer 1

2次元正規分布はscipyに関数があるのでこれを使うと実装できます。

import scipy.stats as stats

sample = stats.multivariate_normal.rvs(mean=[168, 71],cov=[[81,68],[68,250]],size=1000)

上記で平均[168, 71], 分散共分散⾏列[[81,68],[68,250]]の2次元正規分布に従う乱数を1000個生成ができるので、あとはsampleから必要なデータを取り出して散布図や系列の平均を計算すればいいかと思います。