问题请教

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

大家好，我的问题
    根据指数分布参数theta的先验分布GAMMA（a,b）随机生成一组theta值，如何从这一组theta值中任取一个生成一组指数分布的数据？
请大家帮忙看看，用R软件怎么做？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：Theta gamma 指数分布 我的问题 先验分布 如何 软件

先生成一组（例如10个）Gamma随机数：x=rgamma(10,a,b)
随机抽取一个：y=sample(x)
再生成指数随机数（例如10个）：rexp(10,y)

谢谢楼上，学习了

hugebear 发表于 2011-10-12 21:48
先生成一组（例如10个）Gamma随机数：x=rgamma(10,a,b)
随机抽取一个：y=sample(x)
再生成指数随机数（例 ...

明白了，谢谢啊。再请教：上面那过程需要重复1000次，怎么做呢

高手们，请帮帮忙吧。我需要将那过程重复1000次

levyangel 发表于 2011-10-13 09:00
高手们，请帮帮忙吧。我需要将那过程重复1000次

那就循环呗，不知道速度怎么样。。。

a.star=rep(1,1000);a.hat=rep(1,1000)
for (i in 1:1000)
{
a=rgamma(100,3.5,1.1)

a.star[i]=sample(a)



x=rexp(10,a.star)

T=sum(x)+10*x[10]

a.hat[i]=13.5/(T+1.1)
}
ER=sum((a.star-a.hat)^2)/1000
ER
我这样循环有问题，帮我看看，而且我要求的ER值应该比较小，麻烦了

levyangel 发表于 2011-10-13 09:59
a.star=rep(1,1000);a.hat=rep(1,1000)
for (i in 1:1000)
{

对不起，我有个地方没说清楚，应该是sample(a,1)，1用来指明要随机抽取的size是1，默认的话就会给出一次完整不放回抽样的序列。

哦，问题解决了，太感谢你了，真是大好人