options ps=54 ls=76; FILENAME website HTTP "http://www.stat.ufl.edu/~winner/data/stature_hand_foot.dat"; data one; infile website; input id gender stature hand foot; female = gender - 1; femhand = female*hand; femfoot = female*foot; handfoot = hand*foot; hand2 = hand*hand; foot2 = foot*foot; run; ods PDF file="proj1_part3.pdf"; ods graphics on; proc reg; model stature = female hand foot femhand femfoot handfoot hand2 foot2 / selection = cp aic bic best=12; run; proc reg; model stature = ?????; output out=regout p=yhat r=e; run; proc univariate normal; var e; run; proc model; parms ???? stature = ??????; fit stature / breusch = (1 ????); run; ods graphics off; ods PDF close; quit;