options ls=72;
data aortic1;
input bsa ava group $;
avabsa = ava / bsa;
lavabsa = log(1 + avabsa);
levene = 0;
if group = 'Healthy'  then levene = abs(lavabsa - 1.303);
if group = 'Stenotic' then levene = abs(lavabsa - 0.579);
label bsa = Body Surface Area
      ava = Aortic Valve Area
      group = Health Status
      avabsa = AVA to BSA Ratio
      levene = Levene on log(1 + AVA/BSA Ratio)
      lavabsa = log(1 + AVA/BSA Ratio);
cards;
  0.900   0.000 Healthy
  0.090   0.047 Healthy
  0.080   0.100 Healthy
  0.090   0.234 Healthy
  0.200   0.170 Healthy
  0.150   0.226 Healthy
  0.180   0.170 Healthy
  0.170   0.275 Healthy
  0.190   0.236 Healthy
  0.170   0.296 Healthy
  0.240   0.446 Healthy
  0.280   0.476 Healthy
  0.190   0.481 Healthy
  0.360   0.451 Healthy
  0.330   0.388 Healthy
  0.570   0.418 Healthy
  0.590   0.519 Healthy
  0.600   0.519 Healthy
  0.340   0.714 Healthy
  0.550   0.605 Healthy
  0.700   0.878 Healthy
  0.680   1.056 Healthy
  0.950   1.100 Healthy
  1.060   1.228 Healthy
  0.600   1.288 Healthy
  1.000   0.819 Healthy
  0.680   1.579 Healthy
  1.150   0.895 Healthy
  0.700   1.172 Healthy
  0.750   1.726 Healthy
  0.960   1.260 Healthy
  0.740   0.729 Healthy
  0.740   0.729 Healthy
  0.770   1.205 Healthy
  0.900   1.430 Healthy
  1.620   1.986 Healthy
  1.450   2.078 Healthy
  1.730   1.794 Healthy
  2.100   2.632 Healthy
  1.450   3.106 Healthy
  1.800   2.615 Healthy
  1.450   2.321 Healthy
  1.540   2.377 Healthy
  1.700   3.221 Healthy
  0.700   0.708 Healthy
  1.170   1.329 Healthy
  0.120   0.201 Healthy
  1.200   1.438 Healthy
  1.270   1.179 Healthy
  0.570   0.749 Healthy
  1.810   2.740 Healthy
  0.770   0.748 Healthy
  1.350   1.292 Healthy
  0.190   0.375 Healthy
  0.490   0.369 Healthy
  1.150   1.608 Healthy
  0.950   2.166 Healthy
  1.510   1.337 Healthy
  0.850   1.245 Healthy
  0.720   0.656 Healthy
  0.730   0.961 Healthy
  1.540   1.497 Healthy
  0.190   0.400 Healthy
  1.280   2.361 Healthy
  0.600   0.776 Healthy
  0.750   1.340 Healthy
  0.290   0.547 Healthy
  0.650   0.506 Healthy
  0.720   0.670 Healthy
  0.550   0.761 Healthy
  0.640   0.640 Stenotic
  1.140   1.730 Stenotic
  1.610   0.770 Stenotic
  1.050   0.670 Stenotic
  0.710   0.400 Stenotic
  0.530   0.400 Stenotic
  0.970   0.710 Stenotic
  1.050   0.850 Stenotic
  2.030   2.900 Stenotic
  1.670   0.760 Stenotic
  1.170   0.860 Stenotic
  2.050   2.260 Stenotic
  0.460   0.290 Stenotic
  0.200   0.070 Stenotic
  1.080   0.700 Stenotic
  1.260   0.730 Stenotic
  1.390   0.940 Stenotic
  0.970   1.310 Stenotic
  0.870   0.400 Stenotic
  1.540   0.530 Stenotic
  1.270   0.680 Stenotic
  1.910   0.740 Stenotic
  2.070   1.110 Stenotic
  1.600   1.120 Stenotic
  1.680   1.310 Stenotic
  2.240   1.470 Stenotic
  0.920   0.660 Stenotic
  0.600   0.340 Stenotic
  1.600   2.550 Stenotic
  1.470   0.970 Stenotic
  1.610   1.200 Stenotic
  1.550   0.660 Stenotic
  1.020   0.260 Stenotic
  1.980   1.390 Stenotic
  1.430   0.440 Stenotic
  0.900   0.580 Stenotic
  0.280   0.150 Stenotic
  0.770   0.270 Stenotic
  1.370   0.840 Stenotic
  1.600   1.010 Stenotic
  1.270   0.490 Stenotic
  2.000   0.810 Stenotic
  1.910   0.930 Stenotic
  1.390   0.550 Stenotic
  1.800   1.320 Stenotic
  0.390   0.080 Stenotic
  1.110   0.620 Stenotic
  0.820   0.520 Stenotic
  1.740   0.710 Stenotic
  0.260   0.070 Stenotic
  0.330   0.120 Stenotic
  2.320   0.580 Stenotic
  0.870   0.650 Stenotic
  0.250   0.060 Stenotic
  0.410   3.000 Stenotic
  1.020   0.340 Stenotic

proc sort; by group;

proc ttest cochran;
class group;
var avabsa lavabsa;
title Stenotic Data with the Outlier: Group=0 (Healthy) =1 (Stenotic);

proc npar1way wilcoxon;
class group;
var avabsa lavabsa;
title Stenotic Data with the Outlier: Group=0 (Healthy) =1 (Stenotic);


data aortic2;
set aortic1;
if avabsa < 5.0;

proc ttest cochran;
class group;
var avabsa lavabsa;
title Stenotic Data without the Outlier: Group=0 (Healthy) =1 (Stenotic);

proc npar1way wilcoxon;
class group;
var avabsa lavabsa;
title Stenotic Data without the Outlier: Group=0 (Healthy) =1 (Stenotic);

proc ttest;
class group;
var levene;
title Stenotic Data without the Outlier: Levene Test;

proc npar1way wilcoxon;
class group;
var levene;
title Stenotic Data without the Outlier: Levene Test;

run;