/*BPS Problem 4.8*/
data archaeopteryx;
input femur humerus;
lines;
38 41
56 63
59 70
64 72
74 84
;
run;
proc plot data = archaeopteryx;
plot humerus*femur;
run;
proc corr data = archaeopteryx;
var femur humerus;
run;
proc reg data = archaeopteryx;
model humerus=femur/p;
output out=error residual=resid
run;
proc univariate data = error normal;
var resid;
run;
proc plot data = error;
plot resid * femur;
run;
quit;