回归分析实验程序-SAS版..docVIP

proc import out= xt49 /*使用import过程导入数据并输出到数据集xt4.9*/ datafile=E:\xt49.xls dbms=excel2000 replace; getnames=yes; /*首行为变量名*/ run; proc plot data=xt49;/*对xt49绘图*/ plot y*x=*;/*以x为横坐标,y为纵坐标,以*为各点,画散点图*/ run; proc corr pearson data=xt49;/*对xt49运行相关分析过程*/ var y x;/*计算y和x的Pearson相关系数*/ run; proc reg data=xt49;/*对xt4.9运行回归分析过程*/ model y=x;/*建立以y为因变量,以x为自变量的线性回归方程*/ model y=x/p r dw;/*建立以y为因变量,以x为自变量的线性回归方程，p是要求输出拟合值，r是要求输出残差值, dw是要求输出DW检验统计量的值*/ model y=x1-x4/vif；/*建立以y为因变量,以x1-x4为自变量的线性回归方程，vif是要求输出各自变量的VIF值*/ output out=res p=yhat r=residual;/*输出拟合值和残差值至数据集res，以便绘制残差图*/ run; -------------------------以下是绘制残差图的程序， data res_new;/*创建新数据集res_new*/ set res;/*先把res数据集复制过来*/ lag1residual=lag1(residual);/*lagn（n自定）函数可把一变量的各观测值移后n位；residual即,lag1residual即*/ t=_n_;/*_n_是data步内读取观测值的计数器变量，从1开始，每读取一观测值自加1，因此变量t的观测值即为期数1,2,...,n*/ run; proc plot data=res_new;/*绘制残差图*/ plot residual*lag1residual=*;/*以residual即残差值为纵坐标，以residual2即拟合值为横坐标*/ plot residual*t=*;/*以residual即残差值为纵坐标，以t即拟合值为横坐标*/ run; - ------------------以下是进行一阶差分后建立回归模型以及其自相关检验的程序 data et49_new;/*创建新数据集ch4_new*/ set xt49;/*先把ch4数据集复制过来*/ difx=x-lag1(x);/*lagn（n自定）函数可把一变量的各观测值移后n位；对x各观测值作一阶差分*/ dify=y-lag1(y);/*lagn（n自定）函数可把一变量的各观测值移后n位；对y各观测值作一阶差分*/ run; proc reg;/*对ex4.9_new运行回归分析过程*/ model dify=difx/p r dw;/*建立以y为因变量,以difx为自变量的线性回归方程，p是要求输出拟合值，r是要求输出残差值,dw是要求输出DW检验统计量的值*/ output out=res p=yhat r=residual;/*输出拟合值和残差值至数据集res，以便绘制残差图*/ run; data res_new;/*创建新数据集res_new*/ set res;/*先把res数据集复制过来*/ lag1residual=lag1(residual);/*lagn（n自定）函数可把一变量的各观测值移后n位；residual即 ,lag1residual即 */ t=_n_;/*_n_是data步内读取观测值的计数器变量，从1开始，每读取一观测值自加1，因此变量t的观测值即为期数1,2,...,n*/ run; proc plot data=res_new;/*绘制残差图*/ plot residual*lag1residual=*;/*以residual即残差值为纵坐标，以residual2即拟合值为横坐标*/ plot residual*t=*;/*以residual即残差值为纵坐标，以t即拟合值为横坐标*/ run; _ _ _ _ _ _ _ _ _ 异常值的的识别 假定有一数据集ch，因变量为y，自变量为x1-x2(或x)。 1. 求各观测体的学生化残差值的程序： proc reg;/*对ch运行回归分析过程*/ model y=x1-x2;/*建立以y为因变量,以x1-x2为自变量的线性回归方程*/ output out=abnormal student=sre;/*输出各观测体的学生化残差值至数据集abno