options nocenter ls=72;

*** generate some data for the example ***;
*** real data would go here ***;

data a;
  do i=1 to 1000;
    a=normal(0);
    x1=a+normal(0);
    x2=a+normal(0);
    y=a+normal(0);
    output;
   end;
run;

data b;set a;
  if x1>0 then x1=1;else x1=0;
  if x2>0 then x2=1;else x2=0;
  if y>0  then y=1;else y=0;
run;

*** perform the regression analyses, writing out parameter estimates
    and standard errors.;

proc logistic outest=c covout;model y=x1 x2;
run;

proc print;run;

data d;set c(keep =   intercep x1 x2);
  file 'temp.dat';put intercep x1 x2;
run;

data e;infile 'temp.dat';input
    intercep x1 x2
    seint    x  x
    x     sex1  x
    x     x   sex2;

  keep intercep x1 x2 seint sex1 sex2;
run;

data f;set e;
  file 'result1.dat';put
    intercep seint/
    x1 sex1/
    x2 sex2;
run;