options linesize=80 pagesize=65 nodate;

title1 'Factorial Design - No Replication';
title2 'Survival Times for Minnows';
title3 'Sokal & Rohlf - Biometry - 2nd ed. - pg. 376';

data minnows;
  do temps = 'A', 'B', 'C';
    do cnconc = 1 to 5;
      do o2conc = 'X', 'Y', 'Z';
        input time @@;
        output;
      end;
    end;
  end;
cards;
201 246 271 150 164 170 131 138 149 130 136 127  97 102  99
124 158 207 104 111 117  86  99  81  89  91  87  60  74  72
 79 129 142  63  54  93  50  51  62  51  52  51  32  46  52
;;;;
proc print;
proc glm;
  classes temps cnconc o2conc;
  model time = temps cnconc o2conc temps*cnconc temps*o2conc cnconc*o2conc;
  means temps cnconc o2conc;

proc chart;
  vbar temps / group=o2conc sumvar=time type=mean;
  vbar temps / group=cnconc sumvar=time type=mean;

proc sort data=minnows;
  by temps cnconc;
proc means data=minnows noprint;
  by temps cnconc;
  output out=tcaves mean=tcm;
  var time;
proc plot;
  plot tcm*cnconc=temps;
  plot tcm*temps=cnconc;

proc sort data=minnows;
  by temps o2conc;
proc means data=minnows noprint;
  by temps o2conc;
  output out=toaves mean=tom;
  var time;
proc plot;
  plot tom*o2conc=temps;
  plot tom*temps=o2conc;

proc sort data=minnows;
  by cnconc o2conc;
proc means data=minnows noprint;
  by cnconc o2conc;
  output out=coaves mean=com;
  var time;
proc plot;
  plot com*cnconc=o2conc;
  plot com*o2conc=cnconc;