作者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