options ls=80;

data a;
input time hardness;
timesq=time*time;
label timesq=time*time;
cards;
10 216
15 222
25 237
30 250
40 262
45 271
50 292
60 308
70 324
85 369
90 375
95 379
;
proc print data=a;
proc means data=a;
proc reg data=a;
 model hardness=time;
output out=b r=resid;

proc plot data=b;
 plot hardness*time resid*time;

proc capability data=b noprint;
  probplot resid;

proc univariate data=b normal;
 var hardness;

proc reg data=a;
 model hardness=time timesq;
options ls=80;

proc reg data=a;
 model hardness=time;
output out=b r=resid;

data c; set b;
residsq=resid*resid;

proc print data=c;

proc reg data=c; 
 model residsq=time;
title 'Breusch-Pagen test';

proc sort data=b; by resid;
proc print data=b;
title 'Sorted data by residual';


proc rank data=b normal=blom out=b;
 var resid;
 ranks nor_perc;    

data b; set b;
new=sqrt(41.16359)*nor_perc;    

proc print data=b;
title 'Sorted data by residual';    

proc corr data=b;
 var resid nor_perc new;
title 'Correlation for correlation test for normality';