一道概率题求助

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

一道经典概率题：一场足球比赛有22名运动员，有人生日相同的概率是多少。
这个如果用古典概率算法，算出来是48%。
我手头刚好有一份600人左右的身份证号码，已经从中提取了生日信息。
想用SAS来抽样验证以上概率
已经用pickit=ceil(rannui(0)*600)在样本中随机抽取观测，每个数据集抽取22个
这样可以通过生成很多个数据集
判断有几个数据集出现生日一样的，除以总数据集，应该接近48%的概率
不过对于每一个数据集，如何判断有生日一样？
还请各位大侠不吝赐教，谢谢。

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：身份证号码 Pick 各位大侠 身份证号 数据集 身份证号码 运动员 经典 如何 生日

48% 怎么算出来的啊？

Schizor 发表于 2012-2-24 16:33
48% 怎么算出来的啊？

1-365*364*...*344/365^22

1. data birth_dt;
   
2.         set id;
   
3.         if length(id)=18 then birth_dt=substr(id,11,4);
   
4.         if length(id)=15 then birth_dt=substr(id,9,4);
   
5. run;
   
6. 
   
7. 
   
8. %let total_num=10000;
   
9. %let count=0;
   
10. %macro Prob;
    
11.         %do i=1 %to &total_num;
    
12.                 proc surveyselect data =birth_dt method = srs n =22 out=Sample noprint;run ;
    
13. 
    
14.                 proc sql noprint;
    
15.                         create table tmp as
    
16.                         select count(*) as num from Sample
    
17.                         group by birth_dt
    
18.                         order by num desc
    
19.                         ;
    
20.                 select num into:max_num from tmp;
    
21.                 quit;
    
22.                 %if &max_num>1 %then %let count=%eval(&count+1);
    
23.         %end;
    
24. %mend;
    
25. %Prob;
    
26. %put %eval(&count/&total_num);

复制代码

看到高手了～

...
proc surveyselect data =a method = srs n =22 out=a2 reps =10000 noprint seed =11;
run ;
ods listing;
proc sql;
   select sum(sameB)*1e-4 label =' the probability is ' from (
      select (count(distinct monthD) <22) as sameB from (
         select replicate, put(month(birthday), z2.)||put(day(birthday), z2.) as monthD from a2)
         group by replicate
)
;
quit;

Ricky 认为程序的效率可以归结为尽量减少DO循环的次数。在我眼里使用%DO往往更糟糕。不过我并没有测试楼上程序的效率。京剧

jingju11 发表于 2012-2-25 06:30
...
proc surveyselect data =a method = srs n =22 out=a2 reps =10000 noprint seed =11;
run ;

恩，很赞！学习之！