作者kria5304 (XenoMegaREENovaSaga)
看板Statistics
標題[問題] 連續型機率密度函數抽樣算均數
時間Wed Feb 6 16:27:02 2019
一個觀念問題 感覺有點基本但太久沒碰統計有點生疏(遮臉...
題目要從λ=2指數分配抽樣10000筆 然後算其均數 變異數
答案應該是趨近於(λ,λ^2)吧
原本我是直接取 x = (1,100,step=0.01) 丟進機率密度函數轉成機率f(x)
然後 sum(x*f(x))
結果就炸了 f(x)自己累加就超過1了...
自己想想 上面的作法是離散吧 連續感覺不能這樣
可是就一直卡在這跳不出去XDD
可以請大神開導一下嗎 謝謝!!
--
恩?你說去去武器走之類的很弱 想學強大一點的咒語?
霹靂卡霹靂拉拉波波莉納貝貝魯多
這個你覺得怎樣?
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 118.168.175.173
※ 文章網址: https://webptt.com/m.aspx?n=bbs/Statistics/M.1549441628.A.87E.html
1F：推 tnlomtt: 機率積分轉換？ 02/06 18:10
2F：→ kria5304: 是說*J嗎 好像有可能誒 回去再算算看 謝謝 02/06 18:57
3F：推 tnlomtt: 用PIT從uniform生成指數亂數 02/06 21:36
4F：→ Pieteacher: 用 -log(U) 為 Exp(1) where U ~ U(0,1) 在調一下 sca 02/06 22:08
5F：→ Pieteacher: le parameter 02/06 22:08
6F：→ recorriendo: 紹面兩樓說的是蒙地卡羅法 你做的是直接積分但是你忘 02/07 06:04
7F：→ recorriendo: 了成dx 可以想想如果你的step更精細 得到更多組f(x) 02/07 06:07
8F：→ recorriendo: 全部加起來不就更大? 所以顯然是缺一個normalization 02/07 06:08
9F：推 goshfju: 用excel就可以啦 先生成U(0,1)亂數:u1,...,u10000 02/07 12:28
10F：推 goshfju: u=F(x)=1-e^(-2x) , x = -ln(1-u)/2,可算出x1,...,x10000 02/07 12:30
11F：→ goshfju: 計算xbar,Sx^2 然後我看不懂你在幹嗎 你的命題不是要算 02/07 12:31
12F：→ goshfju: 樣本平均數跟樣本變異數? 02/07 12:31
13F：→ goshfju: 如果你要用pdf去算,根本不用模擬阿 02/07 12:32
14F：→ goshfju: E(X)=∫xf(x)dx=1/λ, E(X^2)=∫x^2f(x)dx 02/07 12:33
15F：→ goshfju: Var(X)=E(X^2)-(E(X))^2=1/λ^2 02/07 12:33