sas生成确定概率的随机数

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

各位大侠好，现在有一个变量可以取1000个值，分别对应1000个概率值，储存在一个数据集里面，想随机生成一列数，符合上述的概率分布，sas怎么样可以做到啊？

谢谢哈！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：随机数 概率分布 各位大侠 数据集 怎么样

1. data a;
   
2. do a=1 to 4;
   
3. p=a/10;
   
4. output;
   
5. end;
   
6. run;
   
7. 
   
8. data aa;
   
9. set a;
   
10. pp+p;
    
11. run;
    
12. 
    
13. %macro ff(n);
    
14. %do i=1 %to &n.;
    
15. %let x=%sysfunc(ranuni(0));
    
16.   data aa;
    
17.   set aa;
    
18.   retain xx s 0;
    
19.   xx=&x.;
    
20.   if xx<pp then s=1;
    
21.   run;
    
22.   data aa;
    
23.   set aa;
    
24.   by s;
    
25.   if first.s and s=1 then do;
    
26.     put xx=  a=  / ;
    
27.   end;
    
28.   drop xx s;
    
29.   run;
    
30. %end;
    
31. %mend ff;
    
32. 
    
33. %ff(10)

复制代码

先来看个简单的情况：

dataset a 里面 a 是变量值，p 是概率值。在 log 里 put 出来的 a 就是你要的符合上述概率分布的随机变量值。

tabled distn， 对应一千个取值和概率，概率和为1.
proc sql; select prob into :plist separated by ',' from ...; quit;
data data2;
do i =1 to 1e5;
x =rand('table', &plist)
;...; run;
生成的x即为原数据的行号_N_.
根据行号进行索引。京剧

zhangzachary 发表于 2012-2-2 11:34
先来看个简单的情况：

谢谢，这种想法很不错！

jingju11 发表于 2012-2-2 13:23
应该是tabled distn， 对应一千个取值和概率，概率和为1.
proc sql; select prob into :plist separated b ...