DATA example7_1;
INPUT X Y;
CARDS;
6.8  546
7.8  553
8.7  562
8.7  563
8.9  570
9.5  575
10.1 581
10.2 605
10.3 607
10.4 621
11.1 624
12.4 626
13.3 632
13.1 640
13.2 656
; 
PROC REG;              /*调用线性回归分析程序*/
MODEL y=x;             /*拟合自变量对因变量的关系*/
OUTPUT out=res r=r;  /*将残差输出到数据集RES*/
RUN;
PROC UNIVARIATE normal data=res;  
/*调用univariate程序，normal选项用于残差的正态性检验*/
VAR r;
RUN; 
PROC REG;       /*调用线性回归分析程序*/
MODEL y=x;      /*拟合自变量对因变量的关系*/
plot student.*p./nomodel nostat;     
/*绘制学生化残差散点图，student.为学生化残差，p.为因变量预测值*/
RUN;
PROC REG;
MODEL y=x;                  /*拟合自变量对因变量的关系*/
OUTPUT out=res r=r;       /*将残差输出到数据集RES，残差命名为R*/
RUN;
DATA res1;
SET res;                    /*建立数据集RES1，将数据集RES的数据导入*/
res=abs(r);                 /*将残差R取绝对值，并赋值给变量RES*/
PROC CORR spearman ;      
/*调用相关分析程序，选项spearman表示进行spearman等级相关分析*/
VAR x res;
RUN;
