资源描述：
《回归分析实验程序-sas版》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、procimportout=xt49/*使用import过程导入数据并输出到数据集xt4.9*/datafile="E:xt49.xls"dbms=excel2000replace;getnames=yes;/*首行为变量名*/run;procplotdata=xt49;/*对xt49绘图*/ploty*x='*';/*以x为横坐标,y为纵坐标,以*为各点,画散点图*/run;proccorrpearsondata=xt49;/*对xt49运行相关分析过程*/varyx;/*计算y和x的Pearson相关系数*/run;proc
2、regdata=xt49;/*对xt4.9运行回归分析过程*/modely=x;/*建立以y为因变量,以x为自变量的线性回归方程*/modely=x/prdw;/*建立以y为因变量,以x为自变量的线性回归方程，p是要求输出拟合值，r是要求输出残差值,dw是要求输出DW检验统计量的值*/modely=x1-x4/vif；/*建立以y为因变量,以x1-x4为自变量的线性回归方程，vif是要求输出各自变量的VIF值*/outputout=resp=yhatr=residual;/*输出拟合值和残差值至数据集res，以便绘制残差图*/ru
3、n;-------------------------以下是绘制残差图的程序，datares_new;/*创建新数据集res_new*/setres;/*先把res数据集复制过来*/lag1residual=lag1(residual);/*lagn（n自定）函数可把一变量的各观测值移后n位；residual即,lag1residual即*/t=_n_;/*_n_是data步内读取观测值的计数器变量，从1开始，每读取一观测值自加1，因此变量t的观测值即为期数1,2,...,n*/run;procplotdata=res_new;/
4、*绘制残差图*/plotresidual*lag1residual='*';/*以residual即残差值为纵坐标，以residual2即拟合值为横坐标*/plotresidual*t='*';/*以residual即残差值为纵坐标，以t即拟合值为横坐标*/run;-------------------以下是进行一阶差分后建立回归模型以及其自相关检验的程序dataet49_new;/*创建新数据集ch4_new*/setxt49;/*先把ch4数据集复制过来*/difx=x-lag1(x);/*lagn（n自定）函数可把一变量的各
5、观测值移后n位；对x各观测值作一阶差分*/dify=y-lag1(y);/*lagn（n自定）函数可把一变量的各观测值移后n位；对y各观测值作一阶差分*/run;procreg;/*对ex4.9_new运行回归分析过程*/modeldify=difx/prdw;/*建立以y为因变量,以difx为自变量的线性回归方程，p是要求输出拟合值，r是要求输出残差值,dw是要求输出DW检验统计量的值*/outputout=resp=yhatr=residual;/*输出拟合值和残差值至数据集res，以便绘制残差图*/run;datares_n
6、ew;/*创建新数据集res_new*/setres;/*先把res数据集复制过来*/lag1residual=lag1(residual);/*lagn（n自定）函数可把一变量的各观测值移后n位；residual即,lag1residual即*/t=_n_;/*_n_是data步内读取观测值的计数器变量，从1开始，每读取一观测值自加1，因此变量t的观测值即为期数1,2,...,n*/run;procplotdata=res_new;/*绘制残差图*/plotresidual*lag1residual='*';/*以residua
7、l即残差值为纵坐标，以residual2即拟合值为横坐标*/plotresidual*t='*';/*以residual即残差值为纵坐标，以t即拟合值为横坐标*/run;_________异常值的的识别假定有一数据集ch，因变量为y，自变量为x1-x2(或x)。1.求各观测体的学生化残差值的程序：procreg;/*对ch运行回归分析过程*/modely=x1-x2;/*建立以y为因变量,以x1-x2为自变量的线性回归方程*/outputout=abnormalstudent=sre;/*输出各观测体的学生化残差值至数据集abno
8、rmal*/run;procprintdata=abnormal;/*打印出abnormal数据集中的数据，查看各观测体的学生化残差值*/run;运行程序后查看数据集abnormal，StudentizedResidual列即为各观测体的学生化残差
