作者WANG3213 (WANG3213)
看板Statistics
标题Re: [问题] 数理统计的问题
时间Thu Apr 6 14:24:25 2006
※ 引述《WANG3213 (WANG3213)》之铭言:
: ※ 引述《taldy ()》之铭言:
: : U1,U2,...~iid U(0,1)
: : 设N为一个随机变数,其值为0,1,2...
: : n n+1
: : 定义N=n iff Π Ui >_ e^(-λ) > Π Ui
: : i=1 i=1
: : 求N的分布
把过程写完整好了........
n n+1
P(N=n) = P(Π Ui ≧ e^(-λ) > Π Ui) By 题目的定义
i=1 i=1
n n+1
= P(Σ -log(Ui) ≦ λ < Σ -log(Ui)) By part a
i=1 i=1
n n+1
= P(Σ Xi ≦ λ < Σ Xi) Xi iid follow exp(1) for i=1~n+1
i=1 i=1
n
= P( Gn ≦ λ < Gn + Xn+1) Gn = Σ Xi ~ gamma(n,1)
i=1
= P( Gn ≦ λ and Gn + Xn+1 > λ)
注意到 Gn 与 Xn+1 仍然独立, 接下来你可以用图解 Gn 与 Xn+1
的范围或条件机率的方式把答案求出来。
方法一:
原式 = P( Gn ≦ λ) - P( Gn + Xn+1 ≦ λ) 注意到 Gn + Xn+1 ~ gamma(n+1,1)
λ 1 n-1 -x λ 1 n -x
= ∫ ----- x e dx - ∫ ------- x e dx 使用分部积分
0 Γ(n) 0 Γ(n+1)
1 n -λ 1 n -λ
= ------- λ e = --- λ e
Γ(n+1) n!
方法二:
λ
原式 = ∫ P( Gn = y and Gn + Xn+1 > λ) dy
0
λ
= ∫ P( Gn + Xn+1 > λ | Gn = y ) P( Gn = y )dy
0
λ
= ∫ P( Xn+1 > λ- y | Gn = y ) P( Gn = y )dy
0
λ
= ∫ P( Xn+1 > λ - y) P( Gn = y )dy
0
λ -(λ-y) 1 n-1 -y
= ∫ e ----- y e dy
0 Γ(n)
λ 1 n-1 -λ
= ∫ ----- y e dy
0 Γ(n)
1 n -λ λ
= ------- y e |
n*Γ(n) 0
1 n -λ 1 n -λ
= ------- λ e = --- λ e 以上好像没啥好说的
Γ(n+1) n!
-λ
还有, N=0 时, P( G0 ≦ λ) = 1,要另外算,得 P(N=0) = e
所以 N ~ Poisson(λ)
--
.﹒‧∴
*˙﹒
★ 甜。心。小。北。鼻
* ◢◣◢◣ http://www.wretch.cc/album/u504053 .﹒‧
◥██◤
◥◤ 统。计。生。活。馆
* http://home.educities.edu.tw/rebecca0924/
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 219.71.62.141
1F:推 taldy:GAMMA参数写法颠倒了 然怪我会算错@@ 04/06 14:55
2F:推 taldy:感谢 你耐心回答我的问题@@ 04/06 15:18
※ 编辑: WANG3213 来自: 219.71.62.141 (04/06 15:22)