/* ch 5 untransformed */
/* ch 5 plot look for transform */
options ls=72;
data beads;
  input size finglnth fingcat time;
  LINES;
   1 64 1 100
   1 72 2 140
   1 76 3 93
   2 70 1 85
   2 71 2 45
   2 79 3 105
   3 70 1 40
   3 75 2 37
   3 77 3 35
;
data; set beads;
proc sort; by size;
proc means noprint mean var;
  var time; by size;
  output out=bdmnvr mean=avtime var=vartime;
proc print; var size avtime vartime;
data; set bdmnvr;
  ln_av=log(avtime); ln_var=log(vartime);
proc print;
proc plot;
  plot ln_var*ln_av/vpos=19 hpos=50;
run;quit;
/* ch 5 plot look for transform */
/* class should look at residual analysis for htime */