options ps=55 ls=76;
data a;
input y x1 x2 x3 x4;
datalines;
10 1.0 0.0 0.0 0.0 
20 1.1 1.2 3.5 0.0 
22 1.3 1.1 3.2 0.0 
29 1.5 2.0 3.3 1.6 
30 0.0 2.0 3.0 4.1 
33 0.0 2.1 3.2 4.1 
31 0.0 2.0 3.2 4.1 
40 0.0 0.0 3.6 4.5 
36 0.0 0.0 3.2 4.1 
43 0.0 0.0 3.6 4.0 
50 3.1 4.5 0.0 0.0 
50 3.0 4.7 0.0 0.0
proc reg;
model y = x1 x2 x3 x4;
output out=new1 p=pred r=resid;
proc univariate normal plot;
 var resid;
run;