实验五__回归分析sas过程

《实验五__回归分析sas过程》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、实验五回归分析SAS过程（2）实验目的：1.会对实际问题建立有效的多元回归模型，能对回归模型进行残差分析；2.掌握SAS输出结果用于判别回归方程优良性的不同统计量，能对回归模型进行运用，对实际问题进行预测或控制.实验要求:编写程序，结果分析．实验内容：1.误差的正态性检验有几种方法，何时认为误差项服从正态分布？答：1.学生化残差2.残差正态性的频率检验3.残差的正态QQ图检验判断若散点大致在一条直线上相关系数：认为来自正态分布，接受误差正态性检验．2.回归方程的选取的穷举法中，评价回归方程优良性的准则有哪些？

2、根据准则何时方程最优？答：1）修正的复相关系数准则或均方残差准则（或准则）2）准则3）预测平方和准则（准则）拟合所有可能的个回归方程，画出图:,在图中选取最接近参考直线的点所对应的回归方程为最优方程．3.简述逐步回归方法的思想和步骤．基本思想：逐个引入自变量建立回归方程，每次引入对Y影响最显著的自变量,并对方程中旧变量逐个进行检验,把变为不显著的变量逐个从方程中剔除掉,最终得到的方程中,既不漏掉对Y显著影响的变量,又不包含对Y影响不显著的变量.添加或删除某个自变量的准则是用残差平方和的相对减少或增加来衡量．步

3、骤：（1）修正的复相关系数准则、准则选择模型（2）预测平方和准则选择PRESSp最优回归方程（3）最优模型的拟合检验4.做2.62.8（选作）2.9注意：可以选课外综合题目。2.6程序：dataexamp2_6;inputx1x2y;cards;8.37010.38.66510.38.86310.210.57216.410.78118.810.88319.711.06615.611.07518.211.18022.611.27519.911.37924.211.47621.011.47621.411.7692

4、1.312.07519.112.97422.212.98533.813.38627.413.77125.713.86424.914.07834.514.28031.714.57436.316.07238.316.37742.617.38155.417.58255.717.98058.318.08051.518.08051.020.68777.0;run;procregdata=examp2_6;modely=x1-x2;outputout=ap=predictr=residh=hstudent=r;run;d

5、atab;seta;dropx1-x2;run;procprintdata=b;run;proccapabilitygraphicsnoprintdata=a;/*对数据集a调用capability过程，高分辨图，不打印输出*/qqplotr/normal;/*作student数据的正态QQ图*/run;goptionsreset=all;/*将图形的设置恢复为默认状态*/procgplotdata=a;/*对数据集a作出画高分辨的散点图或曲线图*/plotresid*predict;/*画纵坐标为残差、横坐

6、标为yi散点图*/symbolv=doti=none;/*散点表示符号圆点•，不画连线*/run;/*此处至Quit是计算学生化残差对应的标准正态分布的分位数*/procsortdata=a;byr;/*按r排序*/prociml;/*调用iml矩阵分析模块，计算数据*/usea;/*打开数据集a*/readallvar{r}intorr;/*读入集a中变量r(学生化残差)各观测值到矩阵rr中*/doi=1to31;/*此循环计算*/qi=probit((i-0.375)/54.25);q=q//qi;/*矩

7、阵qi上下连接而成，即得54*1阶矩q=（q(1),q(2),¼,q(30)）T*/end;rq=rr

8、

9、q;/*表示矩阵rq=（rrq）*/createcorrelvar{rq};/*创建数据集correl,变量为r、q*/appendfromrq;/*从矩阵rq读取数据（）*/quit;/*iml过程结束*/proccorrdata=correl;/*计算学生化残差与对应的标准正态分布的分位数的相关系数*/run;（a）学生化残差的正态QQ图（b）拟合值的残差图结果分析:1)由学生化残差的正态QQ图可知，

10、其点明显不在一条直线上;2)求得有序学生化残差与相应正态分布的分位数的相关系数=0.94091与1相差较大．因此，若拟合线性回归模型，则误差分布与正态分布有较大的偏离;3)Y拟合值的残差图也表明与不满足线性关系,且两个拟合值还为负数．由此知,直接假定体积与直径和树干高度之间的线性回归关系是不恰当的．（2）对因变量作Box-Cox变换第一步：确定变换参数鉴于（1）中的残差分析结果，我们对作Box-Co