* file: plotMeans.wolves.sas
*
* Plots of least squares means for Wolf skull measurements
* NOTE: least squares means taken from the output of file
*       manova.wolves.sas;
*;
proc format;
  value facfmt
     1='Length'
     2='Width1'
     3='Width2';
run;
data sexes;
  input sex $ factor lsmean;
  format factor facfmt.;
  label lsmean='Least Squares Mean';
cards;
female   1   -.53
female   2   -1.03
female   3   -.43
male     1   .61
male     2   .49
male     3   .37
;
data area;
  input location $ factor lsmean;
  format factor facfmt.;
  label lsmean='Least Squares Mean';
cards;
Rockies    1   .74
Rockies    2  -.53
Rockies    3   .21
Arctic     1  -.66
Arctic     2  -.01
Arctic     3  -.27
;
data all;
  input wolf $ factor lsmean;
  format factor facfmt.;
  label lsmean='Least Squares Mean';
cards;
RM_female    1   .09
RM_female    2  -1.60
RM_female    3  -.57
RM_male    1    1.40
RM_male    2     .54
RM_male    3     .99
AR_female    1  -1.14
AR_female    2  -.46
AR_female    3  -.28
AR_male    1   -.18
AR_male    2   .44
AR_male    3   -.25
;
*
* plots
*;
title 'Wolf Skull Measurements';
title2 'Least Square Means: Sex';
proc gplot data=sexes;
   symbol1 c=black l=1 w=3 i=join v=* h=5;
   symbol2 c=black l=2 w=3 i=join v=> h=5;
   plot lsmean*factor=sex;
run;
title2 'Least Square Means: Location';
proc gplot data=area;
   symbol1 c=black l=1 w=3 i=join v=A h=3;
   symbol2 c=black l=2 w=3 i=join v=R h=3;
   plot lsmean*factor=location;
run;
title2 'Least Square Means: Location and Sex';
proc gplot data=all;
   symbol1 c=black l=1 w=3 i=join v=* h=5;
   symbol2 c=black l=2 w=3 i=join v=> h=5;
   symbol3 c=black l=3 w=3 i=join v=X h=3;
   symbol4 c=black l=4 w=3 i=join v=Y h=3;
   plot lsmean*factor=wolf;
run;
quit;