OPTIONS PS=66 LS=72 NODATE NONUMBER; ***********************************************; *** LACK OF FIT ANALYSIS FOR TREADWEAR DATA ***; ***********************************************; DATA IN; INPUT WEAR TEMP MILES @@; CARDS; 2.25847 53.2 388 2.19915 53.2 388 2.19068 53.2 388 1.99153 53.2 388 2.27837 53.3 438 2.22698 53.3 438 2.03854 53.3 438 2.05139 53.3 438 2.68891 66.0 58 2.63003 66.0 58 2.59078 66.0 58 2.58587 66.0 58 2.35556 70.3 7 2.41333 70.3 7 2.34222 70.3 7 2.42667 70.3 7 2.15361 76.9 28 2.20181 76.9 28 2.22892 76.9 28 2.16867 76.9 28 2.09884 78.4 25 2.25000 78.4 25 2.04360 78.4 25 2.08721 78.4 25 2.07979 88.1 275 2.23404 88.1 275 2.07713 88.1 275 2.21011 88.1 275 2.01934 89.6 324 2.26796 89.6 324 2.10221 89.6 324 2.01105 89.6 324 ; DATA IN; SET IN; TMPMLS = TEMP*MILES; TEMP2 = TEMP**2; MILES2 = MILES**2; *** CODE THE VARIABLE LEVELS ***; T = (TEMP-71.4)/18.2; M = (MILES-222.5)/215.5; TM = T*M; T2 = T**2; M2 = M**2; /* PROC RSREG DATA=IN OUT=CONF; MODEL WEAR = T TM M / LACKFIT COVAR=2 PREDICT U95 L95 U95M L95M COOKD; TITLE 'ANOVA WITH T2 REMOVED USING PROC RSREG'; DATA CONF2; SET CONF; IF _TYPE_='D'; PROC PRINT DATA=CONF2; DATA CONF; SET CONF; IF _TYPE_='D' THEN DELETE; PROC SORT DATA=CONF; BY _TYPE_ T M; PROC MEANS DATA=CONF NOPRINT; BY _TYPE_ T M; VAR WEAR; OUTPUT OUT=WSET MEAN=WEAR; PROC PRINT DATA = WSET; BY _TYPE_; */ DATA IN; SET IN END=EOF; OUTPUT; IF EOF THEN DO; DO T= -1 TO 1 BY .05; DO M= -1 TO 1 BY .05; WEAR =.; OUTPUT; END; END; END; PROC RSREG DATA=IN OUT=CONT NOPRINT; MODEL WEAR = T TM M / PREDICT L95M U95M COVAR=2; PROC SORT DATA=CONT; BY _TYPE_; PROC PLOT DATA=CONT; BY _TYPE_; PLOT T*M=WEAR / CONTOUR=8; TITLE 'CONTOUR PLOTS BY CODED TEMPERATURE AND MILES'; RUN;