options nodate nonumber ps=55 ls=80;

title 'Kuehl -- Box_Cox Transformation';

data crab;
do site=1 to 6;
do rep=1 to 25;
input count @@;
tcount=count+(1/6);
tcm10=((tcount**(-1))-1)/(-1);
tcm09=((tcount**(-0.9))-1)/(-0.9);
tcm08=((tcount**(-0.8))-1)/(-0.8);
tcm07=((tcount**(-0.7))-1)/(-0.7);
tcm06=((tcount**(-0.6))-1)/(-0.6);
tcm05=((tcount**(-0.5))-1)/(-0.5);
tcm04=((tcount**(-0.4))-1)/(-0.4);
tcm03=((tcount**(-0.3))-1)/(-0.3);
tcm02=((tcount**(-0.2))-1)/(-0.2);
tcm01=((tcount**(-0.1))-1)/(-0.1);
tcp10=((tcount**(1))-1)/(1);
tcp09=((tcount**(0.9))-1)/(0.9);
tcp08=((tcount**(0.8))-1)/(0.8);
tcp07=((tcount**(0.7))-1)/(0.7);
tcp06=((tcount**(0.6))-1)/(0.6);
tcp05=((tcount**(0.5))-1)/(0.5);
tcp04=((tcount**(0.4))-1)/(0.4);
tcp03=((tcount**(0.3))-1)/(0.3);
tcp02=((tcount**(0.2))-1)/(0.2);
tcp01=((tcount**(0.1))-1)/(0.1);
tcmp0=log(tcount);
output;
end; end;
datalines;
0 0 22 3 17 0 0 7 11 11 73 33 0 65 13 44 20 27 48 104 233 81 22 9 2
415 466 6 14 12 0 3 1 16 55 142 10 2 145 6 4 5 124 24 204 0 0 56 0 8
0 0 4 13 5 1 1 4 4 36 407 0 0 18 4 14 0 24 52 314 245 107 5 6 2
0 0 0 4 2 2 5 4 2 1 0 12 1 30 0 3 28 2 21 8 82 12 10 2 0
0 1 1 2 2 1 2 29 2 2 0 13 0 19 1 3 26 30 5 4 94 1 9 3 0
0 0 0 2 3 0 0 4 0 5 4 22 0 64 4 4 43 3 16 19 95 6 22 0 0
;
run;


proc glm noprint;
class site;
model tcm10-tcm01 tcmp0 tcp01-tcp10=site;
output out=boxcox r=em10-em01 emp0 ep01-ep10;
run;

proc means uss;
var em10-em01 emp0 ep01-ep10;
title2 'This gives the Error SS for each transformation';
run;


quit;