data in;
	input TEMP SOAKING SOAKPCT DIFFTIME DIFFPCT PITCH;
datalines;
1650 0.58 1.10 0.25 0.90 0.013 
1650 0.66 1.10 0.33 0.90 0.016
1650 0.66 1.10 0.33 0.90 0.015
1650 0.66 1.10 0.33 0.95 0.016
1600 0.66 1.15 0.33 1.00 0.015
1600 0.66 1.15 0.33 1.00 0.016
1650 1.00 1.10 0.50 0.80 0.014
1650 1.17 1.10 0.58 0.80 0.021 
1650 1.17 1.10 0.58 0.80 0.018
1650 1.17 1.10 0.58 0.80 0.019
1650 1.17 1.10 0.58 0.90 0.021
1650 1.17 1.10 0.58 0.90 0.019
1650 1.17 1.15 0.58 0.90 0.021
1650 1.20 1.15 1.10 0.80 0.025
1650 2.00 1.15 1.00 0.80 0.025
1650 2.00 1.10 1.10 0.80 0.026
1650 2.20 1.10 1.10 0.80 0.024
1650 2.20 1.10 1.10 0.80 0.025
1650 2.20 1.15 1.10 0.80 0.024
1650 2.20 1.10 1.10 0.90 0.025
1650 2.20 1.10 1.10 0.90 0.027
1650 2.20 1.10 1.50 0.90 0.026
1650 3.00 1.15 1.50 0.80 0.029
1650 3.00 1.10 1.50 0.70 0.030
1650 3.00 1.10 1.50 0.75 0.028
1650 3.00 1.15 1.66 0.85 0.032
1650 3.33 1.10 1.50 0.80 0.033
1700 4.00 1.10 1.50 0.70 0.039
1650 4.00 1.10 1.50 0.70 0.040
1650 4.00 1.15 1.50 0.85 0.035
1700 12.5 1.00 1.50 0.70 0.056
1700 18.5 1.00 1.50 0.70 0.068
;
run;

** Create a scatterplot of the data;
* from http://www.ats.ucla.edu/stat/sas/faq/sas9_stat_plots.htm ; 
ods html style=journal; /*turn on html output*/
ods listing close; /*turn off list (regular) output window, optional*/
ods graphics on; /*turn on ods graphics*/
proc corr data = in plots;
  var PITCH TEMP SOAKING SOAKPCT DIFFTIME DIFFPCT;
run;
ods graphics off;
ods html close;
ods listing;

* Perform a regression with many diagnostic plots;
ods html; ods listing close; ods graphics on;
proc reg data = in; 
  model PITCH = TEMP SOAKING SOAKPCT DIFFTIME DIFFPCT; 
run;
ods graphics off; ods html close; ods listing;


/*
proc reg data=in;
 model PITCH=TEMP SOAKING SOAKPCT DIFFTIME DIFFPCT/alpha=.95;
 output out=M26out cookd=d h=h lcl=lpred ucl=upred lclm=lest uclm=uest 
press=press p=fit r=resid student=standres rstudent=studres stdr=se_resid
stdi=se_pred stdp=se_mean;
run;
*/

PROC GLM DATA= in;
   MODEL PITCH = TEMP SOAKING SOAKPCT DIFFTIME DIFFPCT / SS3 TOLERANCE;
TITLE 'LACK OF FIT ANALYSIS FOR HEAT TREATING (mm2.6) DATA';
run;
PROC SORT DATA=IN; BY TEMP MILES;
PROC MEANS DATA=IN CSS; BY TEMP MILES;
     VAR WEAR;
     OUTPUT OUT=CSSSET CSS=CSS;

PROC MEANS DATA=CSSSET SUM;
     VAR CSS;
RUN;