多次元正規乱数に従うサンプルを生成したく,以下の資料を参考にし,プログラムを作成しました。生成したサンプルが本当に要求した多次元正規乱数にしたがっているか確認したいのですが,以下の手順で正しいといえるでしょうか?
例:
平均μ=[1; 2]
分散協分散行列
Σ=[1,0;
0,1];
に従う多次元正規乱数y=[y1,y2]を全部でn=100サンプル生成した.
このとき,n=100個のy1の平均が1,n=100個のy2の平均が2となり,
n=100個のy1の分散が1,n=100個のy2の分散が1,y1とy2の共分散が0となれば正しく生成できているといってよいでしょうか?
この方法で確認したのですが,生成したサンプルが元の平均,共分散から0.1~0.2のずれがあり正しいのか不安です。
参考にした資料:
「多次元正規乱数の作成法」
http://ci.nii.ac.jp/naid/110001186428/ja
「共分散の計算の仕方」
http://www.ge.fukui-nct.ac.jp/~nagamizu/f-3-s.pdf
No.2
- 回答日時:
追加です。
ANo.1の例は2次元正規分布でデータ数が100個ですが、p次元でデータ数がn個の場合は、 (Yi-μ)'Σ^(-1)(Yi-μ)と自由度pのカイ二乗分布の(i-0.5)/n×100%点をプロットすることになります。
ご質問の例の
平均μ=[1; 2]
分散共分散行列
Σ=[1,0;
0,1];
という場合なら、
(y1 - 1)^2 + (y2 - 2)^2
を計算し、昇順に並び替えてプロットすることになります。
No.1ベストアンサー
- 回答日時:
Q-Qプロットの一つのカイ二乗確率プロットで確認するのはどうでしょうか?
100組のデータをベクトルYi(i=1~100)と表し、X'を行列Xの転置行列、X^(-1)をXの逆行列とすると、Q-Qプロットの作成方法は、
1. (Yi-μ)'Σ^(-1)(Yi-μ)を計算する。
2. (Yi-μ)'Σ^(-1)(Yi-μ)を昇順に並び替える。
3. 小さいほうからi番目のデータと自由度が2のカイ二乗分布の(i-0.5)/100×100%点を組にして、散布図を作成する。
4. プロットされた点が極端に直線から外れていれば、サンプルが本当に要求した多次元正規乱数にしたがっていないと結論付けられます。
例として、
μ=(1, 2)
Σ =
[1 0.5]
[0.5 1]
である二次元正規分布に従う乱数を100個生成し、散布図とQ-Qプロットを作成してみました。
特に大きく外れた点がないので、うまく乱数が生成できていると考えられます。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 統計学 Excelによるサンプルの拡大について 6 2023/08/22 16:03
- 統計学 統計学の問題です。 数学 51 49 23 77 78 56 44 37 7 29 80 61 36 1 2023/02/03 15:24
- 数学 数学の答えと解き方を教えてください。 問:ある(人数の非常に多い)集団から無作為に6名を選んで身長を 4 2022/12/14 10:06
- 統計学 加重最小二乗法=①「変数を自然対数変換」=②「誤差項の分散の逆数を重み付け」? 8 2022/11/26 11:15
- 数学 以下の数学の問題を教えてください。 確率変数Xは標準正規分布N(0、1)に確率変数Yは平均3のポアソ 3 2022/12/02 19:13
- 統計学 母集団分布を平均 μ, 分散 σ2 の正規分布と想定し, 母集団から無作為抽出した標本のデータ(標本 4 2023/01/30 20:25
- 大学・短大 大学 統計学 1 2022/09/14 11:27
- 統計学 確率統計の問題です。 3 2022/04/07 04:39
- 統計学 統計学の問題です よろしくお願いします 回帰直線 次のデータから集計表を作成し,以下の問いに答えよ。 2 2023/01/31 23:36
- 統計学 統計学の問題です よろしくお願いします 回帰直線 次のデータから集計表を作成し,以下の問いに答えよ。 1 2023/01/31 18:55
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
正規確率紙
-
QQプロット
-
標準偏差
-
統計学における有効数字につい...
-
数学Bの「統計的な推測」の【復...
-
ヒトの可聴周波数は20Hz~20kHz...
-
偏差平方和なぜ2乗
-
顔面偏差値100点満点中何点ぐら...
-
相対標準偏差についてですが…
-
ランダム出力の数値をコントロ...
-
標準偏差バーをグラフに入れた...
-
推定(統計)です。有効数字の...
-
高校偏差値67ってどのくらい字...
-
統計学でいうRSD%とは何ですか。
-
中央値の求め方ってどのように...
-
吹田市に出来た大和大学は、経...
-
統計学の問題についてです。
-
3者見積の異常値
-
数bです。 分散の求め方二つは...
-
正規分布に従わないと標準偏差...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
おすすめ情報