
X_1,…,X_nを正規母集団から取った大きさnの標本とします。
簡単のため、母集団の平均は0、分散は1と仮定します。
このとき標本平均X=(X_1+…+X_n)/nと
標本(不偏)分散s=((X_1-X)^2+…+(X_n-X)^2)/(n-1)
を考えます。
Xは平均0、分散1/nの正規分布に、
(n-1)sは自由度n-1のχ^2分布に従うと思いますが、
このXとsの独立性の証明はどうやったらよいのでしょうか?
結合分布の計算にX_i^2が混じるので大変に面倒です。
非芯χ^2分布の特性関数の計算などを使うのでしょうか。
方針は立つものの、あまりに煩雑な計算になりそうで尻込みしています。
簡便な計算法をご存知であれば教えていただきたく思います。
よろしくお願いします。
No.2ベストアンサー
- 回答日時:
申し訳ありませんが、行列を使って書きます。
あと、ちょっと難しめです。u, v を正規分布に従う平均0の確率変数(スカラー)とします。この2つが独立である条件は、E[uv]=0 であることはご存じだと思います。
今、x=(x1,x2,...,xn)' を標準正規分布に従う確率変数ベクトルとします。したがって、E[x]=0, E[xx']=I です。また、i=(1,..,1)' とします。すると
平均 X = 1/n i'x
分散 s = 1/(n-1) (x-iX)'(x-iX)
と書けることが分かります。ここで
x-X = x - 1/n ii'x
= (I - 1/n ii')x
であり、
(I-1/n ii')'(I-1/n ii') = (I-1/n ii')
であることから
(n-1)s = x'(I-1/n ii')'(I-1/n ii')x
と書くことが出来ます。ここで u=i'x, v=(I-1/n ii')x とおくと、u, v 共に正規分布に従い、
E[uv'] = E[i'xx' (I-1/n ii')']
= i'E[xx'](I-1/n ii')'
= i'(I-1/n ii')' = 0
となるため、u, v は独立であることが分かります。
したがって、独立な確率変数から計算された統計量はやはり独立なので、X と s は独立になります。
u,vが多次元正規だから少し心配でしたが、ガウス系だから共分散行列0で独立が出てくるのですね。♯1様の回答を眺めて、何らかの線形変換で独立な正規確率ベクトルuとvがあって、標本平均はuとある定数ベクトルの内積、標本分散はvのノルムという感じになっているんだなあと思いましたが、まさにこの方法でそれが現れているのですね。大変参考になりました。ありがとうございます。
No.1
- 回答日時:
x=(x1,x2,…,xn)' xi ~ N(μ,σ^2) i.i.d, i=1,2,…,n
X=Σxi/n
s^2=Σ(xi-X)^2/(n-1)
とします。
第1行が(n^(-1/2),…,n^(-1/2))であるような直交変換Gを用いて y=Gxとします。
ただしy=(y1,y2,…,yn)'
すると以下が成り立ちます。
(1) y1=√n・X
(2) Σyi^2=Σxi^2=Σ(xi-X)^2-nX^2, Σyi^2[i=2 to n]=(n-1)s^2
(3) (y1-√n・μ)^2+Σyi^2[i=2 to n]=Σ(xi-μ)^2
Gは直交変換なのでヤコビアンは1ですから、xiの同時密度関数
C・exp(-Σ(xi-μ)^2/2σ^2)
はyiの同時密度関数
C・exp{(-(y1-μ√n)^2-Σyi^2[i=2 to n])/2σ^2}
=C・exp{-(y1-μ√n)^2/2σ^2}・Πexp(-yi^2/2σ^2)[i=2 to n]
に変換されます((3)より)。したがってy1,…,ynは独立であり、
√n・X=y1 ~ N(√n・μ,σ^2)
(n-1)s^2/σ^2=Σyi^2[i=2 to n]/σ^2 ~ χ^2(n-1)
となります((1)(2)より)。またXはy1のみ、s^2はy2,…,ynのみの関数ですから互いに独立になります。
(証明おわり)
どうもありがとうございました。実は今日、稲垣先生の数理統計学を借りてきて読んでみてエレガントだなあ、と感心してたのです。こんなうまい方法があるんですね。独立性の証明はいつも結合分布を求めようとばかりしていたので目から鱗でした。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 統計学 母集団分布を平均 μ, 分散 σ2 の正規分布と想定し, 母集団から無作為抽出した標本のデータ(標本 4 2023/01/30 20:25
- 統計学 不偏分散について 3 2022/03/29 15:57
- 統計学 統計学 最大値の標準偏差 15 2023/02/02 18:36
- 大学・短大 大学 統計学 1 2022/09/14 11:27
- 統計学 ある集団の平均身長を調べたい. 集団全体を母集団とし, 母集団から無作為抽出した15人の平均身長が1 1 2023/02/03 15:26
- 統計学 t値の計算方法 1 2022/11/29 18:37
- 統計学 統計学が分かりません!詳しい解説と回答を教えてくださる方お願いいします! 5 2022/08/23 03:10
- 統計学 統計学の問題です。教えてください(_ _) 数万人の有権者がいる選挙区で, 無作為に400人の標本を 2 2023/02/03 15:27
- 統計学 確率統計の問題です。 3 2022/04/07 04:39
- 統計学 統計学の質問【帰無仮説】 昨年度のA大学卒業者の平均初任給(月額・万円)について調べるために、昨年度 1 2023/05/25 23:30
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
平均偏差MDと標準偏差SDの理想...
-
任意の面内にある点の座標から...
-
最小二乗法の分母について
-
分散共分散行列は半正定値?
-
アップ率の求め方について
-
小数点以下
-
1週間当たりの労働時間の計算方法
-
パーセントの平均の計算式。 42...
-
パーセントの合計と平均について
-
x*y=log(e^x+e^y)と定義すると...
-
関数電卓を使った3乗根の計算方...
-
相加平均と相乗平均
-
数学(算数)教えてください 平...
-
算数の問題です。
-
52-8を工夫して計算してください
-
中3でFカップは大きいですか? ...
-
対数平均について
-
【座標の変換】日本測地系から...
-
エクセルで過去5年間のデータ...
-
Excelで平方2乗平均を計算するには
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
おすすめ情報