推 tnlomtt: 機率積分轉換？ 02/06 18:10
→ kria5304: 是說*J嗎 好像有可能誒 回去再算算看 謝謝 02/06 18:57
推 tnlomtt: 用PIT從uniform生成指數亂數 02/06 21:36
→ Pieteacher: 用 -log(U) 為 Exp(1) where U ~ U(0,1) 在調一下 sca 02/06 22:08
→ Pieteacher: le parameter 02/06 22:08
→ recorriendo: 紹面兩樓說的是蒙地卡羅法 你做的是直接積分但是你忘 02/07 06:04
→ recorriendo: 了成dx 可以想想如果你的step更精細 得到更多組f(x) 02/07 06:07
→ recorriendo: 全部加起來不就更大? 所以顯然是缺一個normalization 02/07 06:08
推 goshfju: 用excel就可以啦 先生成U(0,1)亂數:u1,...,u10000 02/07 12:28
推 goshfju: u=F(x)=1-e^(-2x) , x = -ln(1-u)/2,可算出x1,...,x10000 02/07 12:30
→ goshfju: 計算xbar,Sx^2 然後我看不懂你在幹嗎 你的命題不是要算 02/07 12:31
→ goshfju: 樣本平均數跟樣本變異數? 02/07 12:31
→ goshfju: 如果你要用pdf去算,根本不用模擬阿 02/07 12:32
→ goshfju: E(X)=∫xf(x)dx=1/λ, E(X^2)=∫x^2f(x)dx 02/07 12:33
→ goshfju: Var(X)=E(X^2)-(E(X))^2=1/λ^2 02/07 12:33