Create the scatter plot. On the COMPARE statement, the GROUP= option groups the data by the SPECIES variable.

proc sgscatter data=sashelp.iris;
  title "Iris Data: Length and Width";
  compare x=(sepallength petallength)
          y=(sepalwidth petalwidth)
          / group=species;
run;