[求助]一道关于随机数的题.

1645

收藏 2009-04-02

老师给了一道题,要用值序法生成密度为3/2*x^2,在(-1,1)内的随机数.

但是这个密度不是服从BETA分布的.

我是觉得不能做..

后来又想先值序法做密度为3*x^2的,然后再变换..

可这样感觉又怪怪的..

请高手告之解决方法..

还是说这题真的不能做..

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

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

全部回复

2009-4-2 11:12:00

这种方法我怎么没有见过???

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

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

2009-4-4 06:39:00

/*inverse transformation method*/
data test;
do i = 1 to 1000;
x = ranuni(-1);
output;
end;
run;

data test;
set test;
y = abs(2*x-1)**(1/3);
y = y*sign(2*x-1);
py = 3*y**2/2;
run;

proc gplot;
plot py*y;
run;

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

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

分享