options ls=72;
data ester1;
input esterase count;
logcount = log(count);
logester = log(esterase);
label esterase = Esterase
      logester = Log(Esterase)
      count = Count
      logcount = Log of Count;
cards;
 14.4 255.0
 20.8 459.0
 13.1 215.0
 8.8  153.0
 29.0 595.0
 13.7 249.0
 21.8 416.0
 10.8 167.0
 18.1 342.0
 23.7 466.0
 18.8 262.0
 41.2 597.0
 12.3 256.0
 15.8 271.0
 27.4 646.0
  6.4  84.0
  8.6 105.0
  8.1  96.0
 12.8 219.0
  6.7  86.0
 40.9 239.0
  8.2 130.0
 30.8 438.0
 21.8 296.0
 46.6 599.0
 16.4 256.0
 21.4 270.0
 52.4 697.0
 16.0 148.0
 14.6 239.0
  9.0  79.0
 25.2 418.0
 26.9 472.0
 20.8 270.0
 15.2 193.0
 22.1 317.0
  6.5  85.0
 13.8 154.0
 14.8 266.0
  8.0 104.0
  7.8 127.0
 20.5 354.0
 35.2 416.0
 25.5 435.0
 10.3 159.0
 15.9 250.0
 24.2 412.0
 39.1 1042.0
 41.7 1006.0
 14.6 317.0
 19.0 486.0
 24.4 566.0
 23.2 376.0
 11.2 97.0
 23.8 529.0
 17.5 136.0
 17.0 137.0
 16.9 126.0
 12.6 115.0
 11.6 170.0
 8.0 117.0
 15.2 389.0
 8.6 79.0
 13.1 144.0
 13.8 226.0
 17.7 276.0
 11.8 131.0
 8.4 124.0
 10.9 100.0
 16.0 103.0
 29.0 474.0
 25.2 531.0
 10.6 93.0
 16.1 315.0
 15.0 301.0
 14.1 288.0
 13.3 268.0
 33.6 635.0
 45.0 778.0
 13.9 97.0
 20.5 393.0
 23.5 381.0
 20.9 343.0
 52.0 789.0
 44.5 718.0
 19.2 336.0
 10.5 191.0
 9.2 100.0
 12.8 201.0
 14.6 216.0
 27.2 208.0
 20.8 260.0
 15.2 278.0
 17.1 340.0
 11.1 116.0
 13.3 139.0
 9.5 203.0
 21.2 474.0
 12.1 233.0
 14.0 329.0
 38.6 717.0
 38.2 695.0
 23.0 409.0
 24.6 369.0
 27.7 412.0
 30.3 527.0
 40.5 695.0
 52.1 921.0


proc plot;
plot count * esterase;
title 'Count versus Esterase';

proc plot;
plot logcount * logester;
title 'Log of Count versus Log Esterase';

proc reg data = ester1;
model count = esterase;
title 'Count versus Esterase';
output out = temp01 residual = resid;
data temp02; set temp01;
aresid = abs(resid);
label  resid = Residual from Fit for Count
      aresid = Absolute Residual for Count;
proc plot data = temp02;
plot aresid * esterase;
title 'Absolute residuals for Count';
proc reg data = temp02;
model aresid = esterase;
title 'Absolute residuals for versus esterase';


proc reg data = ester1;
model logcount = logester;
title 'Log of Count versus log Esterase';
output out = temp03 residual = resid;
data temp04; set temp03;
aresid = abs(resid);
label  resid = Residual from Fit for log Count
      aresid = Absolute Residual for log count;
proc plot data = temp04;
plot aresid * logester;
title 'Absolute residuals for Log Count';
proc reg data = temp04;
model aresid = logester;
title 'Absolute residuals for Log Count versus log esterase';