怎样用sas 蒙特卡罗方法 simulate gamma distribution

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

我有一组数据，通过proc histogram 得出他符合gamma分布：
                              Threshold   Theta           0
                                 Scale       Sigma    4.529809
                                 Shape       Alpha    62.27631
                                 Mean                 282.0998
                                 Std Dev              35.74714
我想用蒙特卡罗的方法模拟输出：
1. 通过给定的gamma distribution 参数（如上），确定一个gamma分布，求出gamma分布的cdf，并通过这个函数求反函数。
2.产生随机数0－1之间，通过gamma函数的反函数， 得到模拟的值。

想请教有没有高人做过类似问题，sas有没有写好的gamma distribution 反函数的function。
我在gamma distribution求反函数的地方就卡壳了，目前找资料中。希望各位指点。

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：distribution simulate 蒙特卡罗方法 gamma Istr 蒙特卡罗 反函数 Sigma

有点深奥，没看明白~~呵呵~

gamma distribution的反函数很难求出来，看来我原来的想法存在问题，希望有做过simulation的人指点一下，如何在已知distribution 的条件下 进行 蒙特卡罗simualtion. 貌似要用到x = RAND('GAMMA',a) 函数。

我自己上网搜到了，希望对以后的人有帮助

data gammasimulation;
seed=1234;
n=100;alpha=62.27631;beta=4.529809;
do i=1 to n;
x=beta*rangam(seed,alpha);
output;
end;
keep x;
run;

limuqi 发表于 2011-9-19 10:30
我自己上网搜到了，希望对以后的人有帮助

data gammasimulation;

楼主，想请教一下，在SAS中怎么确定变量的分布，比如是伽马分布还是对数正态分布还是其他什么分布？谢谢你。