をとして証明する


8

母集団の平均の信頼区間を含む統計的問題は、次の重み関数の観点から組み立てることができます。

w(α,n)tn1,α/2nfor 0<α<1 and n>1.

たとえば、無限の超母集団の平均の標準的な古典的なレベルの信頼区間は、次のように書くことができます。1α

CI(1α)=[x¯n±w(α,n)sn].

次の関数を使用して、限界および\ lim _ {\ alpha \ uparrow 1} w(\ alpha、n)= 0を確立することは簡単ですT分布。信頼区間のコンテキストでは、これは、信頼レベルを下げると区間が1つのポイントに縮小し、信頼レベルを上げると実際の線全体に拡大することを示しています。保持する必要があるもう1つの直感的なプロパティは、データを取得するにつれて間隔が1つのポイントに縮小することです。つまり、次のことを意味します。limα0w(α,n)=limα1w(α,n)=0

limnw(α,n)=0.

質問:後者の重み関数の特性を証明してください。


詳細情報:T分布の臨界点に慣れていない数学の読者にとって、値は、陰的な方程式で定義された関数です。 ntn1,α/2n

α2=1(n1)πΓ(n2)Γ(n12)tn1,α/2(1+r2n1)n/2dr.

2
これは、生徒が全体的な変動で正常に収束しないために、些細なことではありませんか?分子は定数行きます。zα/2

1
@guy:「些細な」を構成するものは、あなたが尋ねている人のレベルに依存します。私がこの質問をする理由は、これが私が教えている数人の大学院生に提起した問題であり、このサイトの人々がどのような種類の証拠を思いついたかを見たいからです(これを学生によって与えられた種類の答え)。卒業生は通常、証明にどのようにアプローチするかについて大まかな考えを持っていますが、詳細に苦労しています。ほとんどはを表示してそこから移動しようとしますが、積分の境界を確立することで解決しようとします。tn1,α/2zα/2<
ベン-モニカを復活させる

(私はまた、結果を提示する良い方法として、最もエレガントな証明を期待しています。)
ベン-モニカを復活させる

2
私は誰かがこれをどれほど優雅に得ることができるかに興味があります。個人的な意見ですが、密度の使用は、不要な追加の構造をもたらすため、エレガントではないと主張します。私の回答の証明には2つの重要な教訓があります。まず、命題は確率論の有用なビットであり、DFがどのように動作するかについていくつかのアイデアを与えることができます。第2に、証明手法(数値をカバーするという点で言い換えることができます)は、多くのことを証明するのに非常に役立ちます。

回答:


4

チェビシェフの不等式による証明

これは、チェビシェフの不等式 を使用した証明です。Pr(|T|kσ)1k2

我々は埋める場合とセットその後、制限がありますσtν=νν21/k2=α=Pr(|T|tν,α/2)

Pr(|T|νν21α)Pr(|T|tν,α/2)

したがって、は上記の境界になりますtν,α/2

tν,α/2νν21α

明らかな下限を追加し、によって分割しν+1

0tn1,α/2ν+1νν+1(ν2)1α

これはをから0に ntn1,α/2/nn


私は通常「いい答え」のコメントを避けますが、これは本当にエレガントです!(明らかに+1!)
usεr1185218

1
@usεr11852t統計の分布がn空間の角度の分布に関連しているという私の幾何学的証明が好きかもしれません。
Sextus Empiricus

4

これを行う簡単な方法があると私は確信していますが、結果は以下からすぐに得られます:

命題:レッツ連続分布関数となりそのような分布関数のシーケンス(分布すなわち、)弱いです。次に、を均一にします。F n F nF F nx F x xFFnFnFFn(x)F(x)x

証明:連続性と単調性を使用して、任意の自然数に対して、(および)となるを選択できます。弱い収束とが連続であることにより、。いずれかのために、間隔見つけるを含有するとメモをその。したがって、と、X 0xは1... xはmは F XのJ= J / M X 0 = - X M = F F NX JF XのJY [ XのJ - 1Xをj ] y | F ny Fmx0,x1,,xmF(xj)=j/mx0=xm=FFn(xj)F(xj)y[xj1,xj]y¯ LIMのnのSUPY| FnyFy| 1|Fn(y)F(y)|supj|Fn(xj)F(xj)|+|F(xj)F(xj1)|1m msupy| FnyFy| 0lim¯nsupy|Fn(y)F(y)|1mm任意でした。supy|Fn(y)F(y)|0

次に、が標準正規分布に分布で収束するのは、Slutskyの定理のよく知られたアプリケーションです。前の結果は、、つまり。通常の分位関数を両側に適用すると、ます。 F NT N - 1 α- F T N - 1 α0 F T N - 1 αα T N - 1 αZ αtn1Fn(tn1,α)F(tn1,α)0F(tn1,α)αtn1,αzα

したがって暗示いずれかに対する(特に、)。T N - 1 αtn1,αzαGNGN=tn1,αg(n)0g(n)g(n)=n


3

幾何学的証明

幾何学的ビュー

観測されたサンプルをn次元ユークリッド空間の点と見なし、平均の推定を観測のモデル線への投影と。x 1 = x 2 = = X N = ˉ Xx1,x2,...,xnx1=x2=...=xn=x¯

tスコアは、この空間の2つの距離の比率として表すことができます

  • 投影された点と母集団間の距離
    n(x¯μ)
  • この点と観測の間の距離
    i=1n(x^xi)2

これは、観測とそれが投影される線との間の角度の正接に関連しています。

tn1=n(x¯μ)i=1n(x^xi)2=1tanθ

幾何学的スケッチ

等価t分布と角度分布

この幾何学的ビューでは、tスコアが特定のよりも高い確率は、角度が特定の値よりも小さい確率と同等です。

Pr(|T|>tn1,α/2)=2Pr(θθν,α)=α

または

tn1,α/2n1=1tanθν,α

tスコアは、理論モデルの線と観測の角度に関連していると言えます。信頼区間の外側の点の場合(はから離れていて、角度は小さくなります)、角度はいくつかの制限下回ります。この制限は、より多くの観測で変化します。この角度の制限が大きいに対して90度になる(円錐の形状がより平坦になる、つまり、先のとがらなくなって長くなります)、これは信頼区間のサイズが小さくなり、近づくことを意味しますゼロ。μx¯θν,αθν,αn

角度vs t

n球のキャップの相対面積としての角度分布

独立した正規分布変数の同時確率分布の対称性により、すべての方向は確率が等しく、角度が特定の領域内にある確率は、n球のキャップの相対面積に等しくなります。

このnキャップの相対面積は、n錐台の面積を積分することによって求められます

2Pr(θθc)=211+tan(θc)21(1x2)n32B(12,n12)dx=11+tan(θc)21t0.5(1t)n32B(12,n12)dt=I11+tan(θc)2(12,n12)

ここで、は、上位の正規化された不完全ベータ関数です。Ix(,)

角度の限界

場合 90度に移行次に ゼロになります。θn,αntn1,α/2/n

または逆のステートメント:90度より小さい角度の場合、n球上のその角度の相対面積は、が無限大になるとゼロに減少します。n

直感的には、これは、次元が無限大に増加するにつれて、n球のすべての領域が赤道に集中することを意味します。n

定量的には、式を使用してこれを示すことができます

a1t0.5(1t)n32B(12,n12)dt<a1(1a)n32B(12,n12)dt=(1a)n12B(12,n12)=L(n)

との違いを考慮してください。L(n+2)L(n)

ある時点で、分母の減少は分子の減少に引き継がれますあり、関数はから無限大までゼロに減少します。

B(12,x+1)B(12,x)=xx+12
(1a)n+12(1a)n12=1a
L(n)n

1

我々は持っています

α2=tn1,α/2limn1(n1)πΓ(n2)Γ(n12)(1+r2n1)n/2dr=tn1,α/212πe12r2dr=1Φ(tn1,α/2)1[12+φ(tn1,α/2)(tn1,α/2+(tn1,α/2)33+(tn1,α/2)515+)]

これは、最大のがであるため、囲み括弧内の2番目の項は最大でであることを意味します。は正規分布のあることに注意してください。この近似もこれに基づいてます。12α1φ(x)

したがって、

0<α1+2φ(tn1,α/2)(tn1,α/2+(tn1,α/2)33+(tn1,α/2)515+)<1

@guerhuerhのサイトへようこそ。ここでアカウントの登録と通常のメンバーになることを検討してください(ヘルプセンターの「マイアカウント」セクションに情報があります)。ぜひお会いしましょう。必要に応じて、新規ユーザー向けの情報が記載されているツアーに参加できます。
-モニカの復活

2
これは興味深いアプローチですが、最初のステップを正当化する方法(積分の内部に制限操作を追加した場合)はわかりません。これは偽である暗黙の方程式を与えます。tn1,α/2
ベン-モニカを復活させる

@ベン:ませんか?limnα2=α2
PEV 2018

@PEV:はい。ただし、右側には、積分の内側ではなく、積分の外側に制限があるはずです。積分の制限は依存するため、現在のところ右側は関数であり、と等しくありません。、N α / 2nnα/2
ベン-モニカを復活

@PEV:下限をで置き換える部分と、残りを与える部分の2つに積分を分離することで、この方法を変更できると思います。前者はすでに表示されているものに収束し、後者はゼロに収束することを示す必要があります。私が言ったように、これは興味深いアプローチだと思います。これはメソッドとしてのメリットがあり、正しいものに変更できる可能性がありますが、現在は正しくありません。zα/2
ベン-モニカを復活させる
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.