[程式] R有沒有期望值的數值解?

看板Statistics作者 (azon)時間16年前 (2009/11/22 20:45), 編輯推噓3(3推 0噓 9→)
留言12則, 8人參與, 7年前最新討論串1/1
------------------------------------------------------------------------ [軟體程式類別]: 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
→
11/22 20:57, , 1F
how about increase the sample size? eg rgamma(10^8,1,2)?
11/22 20:57, 1F
→
11/22 21:58, , 2F
能夠唯一嗎 不要每次都不一樣 有方法嗎?
11/22 21:58, 2F
→
11/22 21:58, , 3F
bmka 感謝你回應我
11/22 21:58, 3F
推
11/22 22:22, , 4F
可以設定seed，讓每次都一樣
11/22 22:22, 4F
→
11/23 00:29, , 5F
X ~ gamma(a,b), E(X) = ab, 為什麼要數值解?
11/23 00:29, 5F
→
11/23 00:47, , 6F
增加樣本數
11/23 00:47, 6F
推
11/23 07:29, , 7F
除了增加樣本數，那就是去推公式了，應該沒其他方法了吧
11/23 07:29, 7F
→
11/23 09:20, , 8F
我只是舉例 實際上我碰到的問題是 一個很醜的函數
11/23 09:20, 8F
→
11/23 22:39, , 9F
你可以用HI裡面的arms來抽樣, 找empirical mean
11/23 22:39, 9F
→
11/24 23:09, , 10F
granzi casella123 bugle pran 感謝回應我
11/24 23:09, 10F
推
11/25 20:05, , 11F
期望值不存在就會發生這種事
11/25 20:05, 11F
→
01/02 15:00, 7年前 , 12F
可以設定seed，讓每 http://yofuk.com
01/02 15:00, 12F
‣ 更多 azonazon 的文章
文章代碼(AID): #1B2J9QZ9 (Statistics)