options ls=72;
data armspan1;
input height armspan;
differ = height - armspan;
label height = Height
      armspan = Armspan
      differ  = Height - Armspan;
cards;
   73.0    75.0 
   69.0    69.0 
   74.0    72.0 
   71.0    70.0 
   70.0    70.0 
   71.5    72.5 
   72.0    75.5 
   71.5    71.0 
   74.0    76.0 
   68.5    71.0 
   75.0    76.0 
   73.0    71.0 
   66.0    65.0 
   74.0    74.0 
   73.0    70.0 
   70.0    69.0 
   72.0    75.5 
   72.0    74.0 
   70.0    71.0 
   69.0    66.0 
   74.5    74.5 
   73.5    73.5 
   73.0    73.0 
   72.5    72.5 
   72.0    72.0 
   67.5    67.0 
   70.0    70.5 
   72.0    72.5 
   70.0    69.5 
   75.0    75.5 
   72.0    72.0 
   73.5    73.5 
   76.5    76.0 
   72.5    72.5 
   72.0    72.0 
   72.5    72.5 
   73.5    73.5 
   75.0    75.0 
   72.0    72.0 
   76.5    76.5 
   73.0    75.0 
   69.0    69.5 
   69.5    73.0 
   69.0    66.5 
   71.5    72.5 
   70.5    73.0 
   78.0    80.0 
   73.0    77.0 
   71.0    70.0 
   68.0    64.0 




proc reg data = armspan1;
model differ = height;
title 'Armspan Data';
output out = conmd cookd = cookd residual=resid;
data conmd2; set conmd;
aresid = abs(resid);
label cookd = Cooks Distance
      aresid = Absolute Residual: SVL on Tail Length
      resid = Residual of SVL on Tail Length;
proc plot data = conmd2;
plot cookd * height ;
title 'Armspan Data, Males, Cook Distance';
proc reg data = conmd2;
model aresid = height;
title 'Armspan Data, Males, Absolute Residuals regressed on height';


proc means;

symbol1 value=dot height=2;
symbol1 w=5 i=rl;
symbol2 value=dot height=2;
proc gplot data = conmd2;
plot aresid * height aresid * height / overlay;
title 'Armspan Data, Males, Absolute Residuals';


run;