data frogs;
   label x="body length(cm)";
   label y="mouth width(cm)";
   array sname{2} $20 ("LeptodactylusOcellatus","BufaMariunus");
   do i=1 to 2;
      input x y @;
      species=sname{i};
      output;
   end;
cards;
3.8 1.0 3.8 1.6
4.0 1.2 4.3 1.7
4.9 1.7 6.2 2.3
7.1 2.0 6.3 2.5
8.1 2.7 7.8 3.2
8.5 2.5 8.5 3.0
8.9 2.4 9.0 3.5
9.1 2.9 10.0 3.8
9.8 3.2 . . 
;
run;
proc sort data=frogs;
   by species;
run;

symbol value=dot i=rl;

proc gplot;
   by species;
   plot y*x;
run;
proc gplot;
   plot y*x=species;
run;
proc glm data=frogs;
   class species;
   *model y=species|x  / solution;
   model y=species x  / solution;
   output out=predz p=p r=r;
run;
proc gplot data=predz; 
   plot r*p;
   plot p*y;
   plot p*x y*x/overlay;
run;