UA MATH563 概率论的数学基础 中心极限定理5 Renyi定理
UA MATH563 概率論的數學基礎 中心極限定理5 Renyi定理
這是獨立隨機變量及其性質的一個應用,假設X1,?,Xn~iidFX_1,\cdots,X_n \sim_{iid} FX1?,?,Xn?~iid?F,其中FFF連續可微。
引理 P(Xi≠Xj)=1,?i≠jP(X_i \ne X_j) = 1,\forall i \ne jP(Xi??=Xj?)=1,?i?=j
證明 我們來說明X,YX,YX,Y獨立同分布于FFF,則P(X=Y)=0P(X=Y)=0P(X=Y)=0,這個結果非常類似在平面中,直線的Lebesgue測度為0,雖然結論非常直觀,但是我們要處理的對應也是可能取到無窮的,這在分析上會帶來一些困難,但我們可以用truncation技巧,任選一個正數M>0M>0M>0,
{X=Y}={X=Y,∣X∣,∣Y∣≤M}∩{X=Y,∣X∣,∣Y∣>M}\{X=Y\}=\{X=Y,|X| ,|Y| \le M\} \cap \{X=Y,|X|,|Y|>M\}{X=Y}={X=Y,∣X∣,∣Y∣≤M}∩{X=Y,∣X∣,∣Y∣>M}
先考慮第二項,根據獨立性
P({X=Y,∣X∣,∣Y∣>M})≤P(∣X∣,∣Y∣>M)=P(∣X∣>M)P(∣Y∣>M)=[1?F(M)+F(?M)]2P(\{X=Y,|X|,|Y|>M\}) \le P(|X|,|Y|>M) \\ = P(|X|>M)P(|Y|>M)=[1-F(M)+F(-M)]^2P({X=Y,∣X∣,∣Y∣>M})≤P(∣X∣,∣Y∣>M)=P(∣X∣>M)P(∣Y∣>M)=[1?F(M)+F(?M)]2
根據累積分布函數的性質,當MMM足夠大時,即?δ>0\exists \delta >0?δ>0,?M>δ\forall M > \delta?M>δ,
F(?M)<?,1?F(M)<?,??>0F(-M)<\epsilon,1-F(M)<\epsilon,\forall \epsilon>0F(?M)<?,1?F(M)<?,??>0
下面考慮第一項,對于∣X∣,∣Y∣≤M|X| ,|Y| \le M∣X∣,∣Y∣≤M的區域,我們可以把它分解成(2Mn)2(2Mn)^2(2Mn)2個小區域,并且只有對角線區域,也就是(a/n,(a+1)/n]×(a/n,(a+1)/n](a/n,(a+1)/n] \times (a/n,(a+1)/n](a/n,(a+1)/n]×(a/n,(a+1)/n],a=?Mn,?,Mna=-Mn,\cdots,Mna=?Mn,?,Mn這樣的區域才有概率,計算
P(X,Y∈(a/n,(a+1)/n]×(a/n,(a+1)/n])=(F(a+1n)?F(an))2P(X,Y \in (a/n,(a+1)/n] \times (a/n,(a+1)/n]) \\ = (F(\frac{a+1}{n})-F(\frac{a}{n}))^2P(X,Y∈(a/n,(a+1)/n]×(a/n,(a+1)/n])=(F(na+1?)?F(na?))2
于是
P(X=Y,∣X∣,∣Y∣≤M)≤P(?a(a/n,(a+1)/n]×(a/n,(a+1)/n])=∑a=?MnMn(F(a+1n)?F(an))2≤sup?a∣F(a+1n)?F(an)∣∑a=?MnMn[F(a+1n)?F(an))]=sup?a∣F(a+1n)?F(an)∣[F(M)?F(?M)]≤sup?a∣F(a+1n)?F(an)∣P(X=Y,|X|,|Y| \le M) \\ \le P(\bigcup_a(a/n,(a+1)/n] \times (a/n,(a+1)/n]) \\ = \sum_{a=-Mn}^{Mn}(F(\frac{a+1}{n})-F(\frac{a}{n}))^2 \\ \le \sup_a|F(\frac{a+1}{n})-F(\frac{a}{n})|\sum_{a=-Mn}^{Mn}[F(\frac{a+1}{n})-F(\frac{a}{n}))] \\ = \sup_a|F(\frac{a+1}{n})-F(\frac{a}{n})|[F(M)-F(-M)] \\ \le \sup_a|F(\frac{a+1}{n})-F(\frac{a}{n})|P(X=Y,∣X∣,∣Y∣≤M)≤P(a??(a/n,(a+1)/n]×(a/n,(a+1)/n])=a=?Mn∑Mn?(F(na+1?)?F(na?))2≤asup?∣F(na+1?)?F(na?)∣a=?Mn∑Mn?[F(na+1?)?F(na?))]=asup?∣F(na+1?)?F(na?)∣[F(M)?F(?M)]≤asup?∣F(na+1?)?F(na?)∣
因為FFF是連續的,事實上連續的累積分布函數在閉區間上會是一致連續的,所以當nnn足夠大時,
sup?a∣F(a+1n)?F(an)∣<?\sup_a|F(\frac{a+1}{n})-F(\frac{a}{n})|<\epsilonasup?∣F(na+1?)?F(na?)∣<?
綜上,P(X=Y)=0P(X=Y)=0P(X=Y)=0。
定義 稱XiX_iXi?是一個record如果Xi>max?(X1,?,Xi?1)X_i>\max(X_1,\cdots,X_{i-1})Xi?>max(X1?,?,Xi?1?);XjX_jXj?的rank為Rj=∑k=1j1Xi≥XjR_j = \sum_{k=1}^j1_{X_i \ge X_j}Rj?=∑k=1j?1Xi?≥Xj??,顯然Rj=1?XjisarecordR_j=1 \Leftrightarrow X_j\ is\ a\ recordRj?=1?Xj??is?a?record。
評注 我們可以這樣來理解,假設X1,?,XnX_1,\cdots,X_nX1?,?,Xn?表示某氣象站在某地一天內進行的nnn次溫度測量,但實際上該氣象站只需要披露最低溫與最高溫,Xi>max?(X1,?,Xi?1)X_i>\max(X_1,\cdots,X_{i-1})Xi?>max(X1?,?,Xi?1?)表示XiX_iXi?是到第iii次測量為止的最高溫,于是它形成一個新記錄,RjR_jRj?表示的是到第jjj次測量為止,第jjj次測量值的排名,顯然排名為1等價于它是一個新記錄。
Renyi定理 R1,?,RnR_1,\cdots,R_nR1?,?,Rn?獨立且服從同樣的分布:
P(Rj=k)=1j,k=1,?,jP(R_j=k)=\frac{1}{j},k=1,\cdots,jP(Rj?=k)=j1?,k=1,?,j
評述 定義Aj={Rj=1}A_j = \{R_j=1\}Aj?={Rj?=1}(XjX_jXj?是一個新記錄),Renyi定理說明1A1,?,1An1_{A_1},\cdots,1_{A_n}1A1??,?,1An??獨立,且1Aj~Ber(1/j)1_{A_j} \sim Ber(1/j)1Aj??~Ber(1/j),簡單來說也就是第一、二、。。。,n個測量成為新記錄的概率分別是1,1/2,?,1/n1,1/2,\cdots,1/n1,1/2,?,1/n。
證明
對于X1,?,XnX_1,\cdots,X_nX1?,?,Xn?,用X1(w),?,Xn(w)X_1(w),\cdots,X_n(w)X1?(w),?,Xn?(w)表示一組測量數據(也就是這個隨機變量序列的一個realization),則X1(w),?,Xn(w)X_1(w),\cdots,X_n(w)X1?(w),?,Xn?(w)可能有n!n!n!種次序(其實就是n!n!n!種排列),因為X1,?,XnX_1,\cdots,X_nX1?,?,Xn?獨立同分布,我們可以說明每一種排列的可能性相同,下面是一個簡單的思路:
計算某一種順序的概率,比如X1<X2<?<XnX_1 <X_2 < \cdots < X_nX1?<X2?<?<Xn?,
P(X1<X2<?<Xn)=∫X1<X2<?<Xndμ1×μ2×?×μnP(X_1 <X_2 < \cdots < X_n) = \int_{X_1 <X_2 < \cdots < X_n}d\mu_1 \times \mu_2 \times \cdots \times \mu_nP(X1?<X2?<?<Xn?)=∫X1?<X2?<?<Xn??dμ1?×μ2?×?×μn?
如果交換X1,X2X_1,X_2X1?,X2?的次序,
P(X2<X1<?<Xn)=∫X2<X1<?<Xndμ2×μ1×?×μnP(X_2 <X_1 < \cdots < X_n) = \int_{X_2 <X_1 < \cdots < X_n}d\mu_2 \times \mu_1 \times \cdots \times \mu_nP(X2?<X1?<?<Xn?)=∫X2?<X1?<?<Xn??dμ2?×μ1?×?×μn?
因為X1,X2X_1,X_2X1?,X2?獨立同分布,μ1=μ2\mu_1=\mu_2μ1?=μ2?,所以
dμ1×μ2×?×μn=dμ2×μ1×?×μnd\mu_1 \times \mu_2 \times \cdots \times \mu_n=d\mu_2 \times \mu_1 \times \cdots \times \mu_ndμ1?×μ2?×?×μn?=dμ2?×μ1?×?×μn?
上面兩個積分區域對應的都是Rn\mathbb{R}^nRn被平均分成n!n!n!份后的一份,根據對稱性,
P(X1<X2<?<Xn)=P(X2<X1<?<Xn)=1n!P(X_1 <X_2 < \cdots < X_n) = P(X_2 <X_1 < \cdots < X_n)=\frac{1}{n!}P(X1?<X2?<?<Xn?)=P(X2?<X1?<?<Xn?)=n!1?
對于每一種排列,我們都可以做這個操作,得到相同的概率。
經過簡單觀察,我們可以發現{R1,R2,?,Rn}\{R_1,R_2,\cdots,R_n\}{R1?,R2?,?,Rn?}的值與X1,?,XnX_1,\cdots,X_nX1?,?,Xn?的每一種排列之間是一一對應。比如R1=1,R2=2,R3=3R_1=1,R_2=2,R_3=3R1?=1,R2?=2,R3?=3對應的排列順序就是X1>X2>X3X_1>X_2>X_3X1?>X2?>X3?,再比如R1=1,R2=1,R3=1R_1=1,R_2=1,R_3=1R1?=1,R2?=1,R3?=1對應X1<X2<X3X_1<X_2<X_3X1?<X2?<X3?。所以
P(R1=r1,?,Rn=rn)=P(某種排列順序)=1n!P(R_1=r_1,\cdots,R_n=r_n)=P(某種排列順序)=\frac{1}{n!}P(R1?=r1?,?,Rn?=rn?)=P(某種排列順序)=n!1?
這是R1,?,RnR_1,\cdots,R_nR1?,?,Rn?的聯合分布,下面我們計算每一個RjR_jRj?的邊緣分布。
P(Rj=rj)=∑r1,?,rj?1,rj+1,?,rnP(R1=r1,?,Rn=rn)=1×2×?(j?1)?(j+1)×?×n/n!=1jP(R_j=r_j)=\sum_{r_1,\cdots,r_{j-1},r_{j+1},\cdots,r_n} P(R_1=r_1,\cdots,R_n=r_n) \\ = 1 \times 2 \times \cdots (j-1) \cdots (j+1) \times \cdots \times n/n! = \frac{1}{j}P(Rj?=rj?)=r1?,?,rj?1?,rj+1?,?,rn?∑?P(R1?=r1?,?,Rn?=rn?)=1×2×?(j?1)?(j+1)×?×n/n!=j1?
因此,我們不難驗證
∏j=1nP(Rj=rj)=P(R1=r1,?,Rn=rn)=1n!\prod_{j=1}^n P(R_j = r_j) = P(R_1=r_1,\cdots,R_n = r_n) = \frac{1}{n!}j=1∏n?P(Rj?=rj?)=P(R1?=r1?,?,Rn?=rn?)=n!1?
所以R1,?,RnR_1,\cdots,R_nR1?,?,Rn?獨立。
總結
以上是生活随笔為你收集整理的UA MATH563 概率论的数学基础 中心极限定理5 Renyi定理的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: UA MATH567 高维统计I 概率不
- 下一篇: [概统]本科二年级 概率论与数理统计 第