OPTIONS NODATE NONUMBER PS=60 LS=64;

************************************;
*** USING THE %INCLUDE STATEMENT ***;
************************************;

FILENAME NPPLOT 'NORMPLOT.SAS';
FILENAME PRPLOT 'PREDPLOT.SAS';
FILENAME RNPLOT 'RUNPLOT.SAS';

DATA in; INPUT order voltage exhaust pressure @@; 
CARDS;
1  127  60 51  5  127  60 44   9  127  60 46
2  127 150  8  4  127 150 15  11  127 150 13
3  220  60 48  8  220  60 43  10  220  60 38
6  220 150 28  7  220 150 25  12  220 150 31
;
PROC GLM DATA=in NOPRINT;
     CLASS voltage exhaust;
     MODEL pressure = voltage|exhaust / SS3;
     OUTPUT OUT=diag P=pred R=resid;

%INCLUDE PRPLOT;
%INCLUDE RNPLOT;
%INCLUDE NPPLOT;
RUN;