作者azonazon (azon)
看板Statistics
標題[程式] R有沒有期望值的數值解?
時間Sun Nov 22 20:45:12 2009
------------------------------------------------------------------------
[軟體程式類別]: R
[程式問題]:
[軟體熟悉度]:
中(3個月到1年)
[問題敘述]: 我想要求 來自某分配的某個函數的期望值
例如: 欲求E(X) , X~ gamma(1,2)
在R中 我可以這樣做:
> X=rgamma(1000,1,2)
> mean(X)
但是不論生成隨機變數多寡 這個值是很不穩定的
但我想問的是
能不能給我一個穩定的近似值
有這方面的package 嗎?
[程式範例]:
-----------------------------------------------------------------------------
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.115.221.173
1F：→ bmka:how about increase the sample size? eg rgamma(10^8,1,2)? 11/22 20:57
2F：→ azonazon:能夠唯一嗎 不要每次都不一樣 有方法嗎? 11/22 21:58
3F：→ azonazon:bmka 感謝你回應我 11/22 21:58
4F：推 pran:可以設定seed，讓每次都一樣 11/22 22:22
5F：→ granzi:X ~ gamma(a,b), E(X) = ab, 為什麼要數值解? 11/23 00:29
6F：→ casella123:增加樣本數 11/23 00:47
7F：推 bugle:除了增加樣本數，那就是去推公式了，應該沒其他方法了吧 11/23 07:29
8F：→ azonazon:我只是舉例 實際上我碰到的問題是 一個很醜的函數 11/23 09:20
9F：→ granzi:你可以用HI裡面的arms來抽樣, 找empirical mean 11/23 22:39
10F：→ azonazon:granzi casella123 bugle pran 感謝回應我 11/24 23:09
11F：推 puwawa:期望值不存在就會發生這種事 11/25 20:05