options ls=75 nodate ps=1000;

data one;
   input GENDER $ bp DIET $;
   cards; 
                           female    104    LV 
                           female    112    LV 
                           female    110    LV 
                           female    102    LV 
                           female    122    LV 
                           male      124    LV 
                           male      120    LV 
                           male      112    LV 
                           male      108    LV 
                           male      112    LV 
                           female    114    NOR
                           female    120    NOR
                           female    116    NOR
                           female    102    NOR
                           female    124    NOR
                           male      134    NOR
                           male      112    NOR
                           male      128    NOR
                           male      134    NOR
                           male      128    NOR
                           female    114    SV 
                           female    106    SV 
                           female     88    SV 
                           female    108    SV 
                           female     96    SV 
                           male      122    SV 
                           male      104    SV 
                           male      114    SV 
                           male      128    SV 
                           male      122    SV 
;
run;

proc glm;
   class DIET GENDER;
   model bp=DIET|GENDER;
   means DIET|GENDER;
   lsmeans DIET/stderr;
   means DIET/scheffe cldiff;
   estimate "theta4" diet -1 2 -1/divisor=2;
   output out=two p=p r=r;
run;

data probh;  /* for Scheffe intervals */
   fvalue=sqrt((3-1)*finv(0.95,3-1,24));
run;

proc print;run;

proc gplot data=two;
   plot bp*diet=gender;
run;
symbol1 i=join line=1;
symbol2 i=join line=2;

proc gplot data=two;
   plot p*diet=gender;
run;

/* one more problem: estimate the simple gender effect separately
for each diet using the SLICE command in an LSMEANS statement */

PROC GLM;
   title "check out the output from the slice command";
   CLASS diet gendER;
   MODEL BP=diet|gender;
   LSMEANS diet*gender/SLICE=diet;
run;