options ls=80;
data consump;
input region energy;
cards;
1 15
1 10
1 13
1 14
1 13
2 17
2 12
2 18
2 13
2 15
2 12
3 11
3 7
3 9
3 13
4 10
4 12
4 8
4 7
4 9
;

proc univariate data=consump;
by region;

proc glm data=consump;
class region;
model energy=region;
output out=next r=resid p=yhat;

proc print data=next;
  
proc rank normal=blom;
    var resid;
    ranks nscore;
    
proc plot;
    plot resid*nscore;
    plot resid*yhat;
run;