... procsurveyselect data =a method = srs n =22 out=a2 reps =10000 noprint seed =11; run ;
ods listing; procsql;
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;