options ps=54 ls=80; title 'RPD - Example 12.1 - Box-tidwell Transformation'; data one; input family x y; if y=0 then y=0.01; if x=0 then x=1; cards; 1 0 .36 1 0 .47 1 0 .30 1 72 3.54 1 72 4.35 1 72 4.88 1 144 6.13 1 144 6.49 1 144 7.04 2 0 .32 2 0 .63 2 0 .51 2 72 4.95 2 72 4.45 2 72 1.50 2 144 6.46 2 144 4.35 2 144 2.18 3 0 .00 3 0 .43 3 0 .72 3 72 4.26 3 72 3.89 3 72 6.54 3 144 5.93 3 144 6.29 3 144 9.62 4 0 .54 4 0 .70 4 0 .49 4 72 3.69 4 72 2.81 4 72 4.08 4 144 5.68 4 144 4.68 4 144 5.79 5 0 .44 5 0 .42 5 0 .39 5 72 3.01 5 72 4.08 5 72 4.54 5 144 6.06 5 144 6.05 5 144 6.97 6 0 .55 6 0 .57 6 0 .45 6 72 2.32 6 72 3.57 6 72 3.59 6 144 4.32 6 144 6.11 6 144 6.49 7 0 .20 7 0 .51 7 0 .27 7 72 3.16 7 72 3.17 7 72 3.75 7 144 4.79 7 144 5.74 7 144 5.95 8 0 .31 8 0 .44 8 0 .84 8 72 2.80 8 72 2.96 8 72 2.04 8 144 10.58 8 144 4.44 8 144 1.70 ; run; * proc print; proc glm; class family; model y = family x / noint solution; run; data two; set one; z=0.03756*x*log(x); run; proc glm; class family; model y = family x z / noint solution; run; data three; set two; x1= x**.69517; run; proc glm; class family; model y = family x1 / noint solution; run; data four; set three; z1=0.17512*x1*log(x1); run; proc glm; class family; model y = family x1 z1 / noint solution; run; data five; set four; x2= x1**1.02917; run; proc glm; class family; model y = family x2 / noint solution; run; data six; set five; z2=0.15800*x2*log(x2); run; proc glm; class family; model y = family x2 z2 / noint solution; run; quit;