data t1;
   do i = 1 to 30;
      g='M';
      x=rannor(12356);
      dummy=1;dummyx=x;
      y=1+2*x+rannor(12356);
      output;
      g='F';
      x=rannor(12356);
      dummy=0;dummyx=0;
      y=1.5+1.5*x+rannor(12356);
      output;
   end;
run;

proc reg data=t1;
  model y= dummy dummyx x;
  output out=p p=py;
  s: test dummy=0;
  t: test dummyx=0;
  u: test dummy=0, dummyx=0;
  run;
  quit;

data t2;
  set p;
  if g='M' then do;
     yy=y;
     group='YMO';
     output;
     yy=py;
     group='YMP';
     output;
  end;
  else do;
     yy=y;
     group='YFO';
     output;
     yy=py;
     group='YFP';
     output;
  end;

SYMBOL1 i=none c=green v='$';
SYMBOL2 i=j c=green v=star;
SYMBOL3 i=none c=blue  v='&';
SYMBOL4 i=j c=blue  v='^';;
proc gplot data=t2;
  plot yy*x=group;
  run;
  quit;