弱弱的问一下，w取值0或1，是随机的，还是有一定的数量要求，例如有40个取1,60个取0等。
还有你的N(8.0.25)指的是均值是8，方差是0.25吗，B(0.3)呢？

这里每个数据集w=0和w=1个数相等，各为50个。N(8.0.25)指的是均值是8，方差是0.25，B(0.3)指的是发生概率为0.3的贝努力分布。谢谢您的关注和回复，祝好！

还请大虾们多多帮助，谢谢了

%let ss=20;
%let seed1=12345;
%let seed2=23456;
%let seed3=34567;

%macro ds;

data ith;
do m=1 to 100;
   output;
end;
run;

%do k =1 %to 100;

data ds&k;
   do i =1 to 100;
      if i<=50 then w=1;
          else w=0;
      e=rannor(&seed1);
      x1=8+0.5*rannor(&seed2);
      x2=ranbin(&seed2,1,0.3);
      y=5*w+2*x1+3*x2+e;
      output;
   end;
run;

proc surveyselect data=ith
  method=srs
  rep=1
  sampsize=&ss
  seed=&seed3
  out=samds;
  id _all_;
run;

data samds;
  set samds;
  n=_n_;
run;

%do w =1 %to &ss;
proc sql noprint;
  select m into :m
  from samds
  where n=&w;
quit;

data ds&k;
  set ds&k;
  if i=&m then y=.;
  if missing(y) then r=1;
  else r=0;
run;

%end;


%let seed1=%eval(12345+&k);
%let seed2=%eval(23456+&k);
%let seed3=%eval(34567+&k);

%end;


%mend;
%ds;

谢谢大侠！程序很好用，谢谢您！如果每个数据集想增加一个变量，给出条件概率P（r=1|x1，x2，x3），怎样修改程序呀，谢谢大侠了！