第7讲多元线性回归3.pptxVIP

2.4 多元回归计算程序[例2.1] 某精神病学医生想知道精神病患者经过６个月治疗后疾病恢复的情况Y是否能通过精神错乱的程度X1、猜疑的程度X2两项指标来较为准确地预测。资料如下,试作分析。西南科技大学生命科学与工程学院周海廷制作DATA d4p15; INPUT y x1 x2; CARDS; 28 3.36 6.924 3.23 6.514 2.58 6.221 2.81 6.022 2.80 6.410 2.74 8.428 2.90 5.6 8 2.63 6.923 3.15 6.516 2.60 6.313 2.70 6.922 3.08 6.320 3.04 6.821 3.56 8.813 2.74 7.118 2.78 7.2;西南科技大学生命科学与工程学院周海廷制作PROC REG; MODEL y=x1 x2;RUN;以上是程序１的参数估计结果，不难看出截距项可以去掉。西南科技大学生命科学与工程学院周海廷制作计算各点上因变量的预测值DATA abc2;SET D4P15;PROC REG; MODEL y=x1 x2 / NOINT P R;RUN;将截距项去掉进行残差分析西南科技大学生命科学与工程学院周海廷制作西南科技大学生命科学与工程学院周海廷制作 这是程序２的方差分析和参数估计结果，方程与各参数的检验结果都有显著性意义，所求得的二元线性回归方程为:线性回归方程的估计标准误为：西南科技大学生命科学与工程学院周海廷制作西南科技大学生命科学与工程学院周海廷制作这是对程序２中的二元回归模型作残差分析的结果，从第④、⑤两列发现第８个观测点所对应的学生化残差的绝对值大于２(因STUDENT=-2.170)，故认为该点可能是异常点，需认真检查核对原始数据。西南科技大学生命科学与工程学院周海廷制作DATA abc3;SET D4P15;IF _N_=8 THEN DELETE;PROC REG; MODEL y=x1 x2 / NOINT P R;RUN;西南科技大学生命科学与工程学院周海廷制作西南科技大学生命科学与工程学院周海廷制作 这是程序３的方差分析和参数估计结果，方程与各参数的检验结果都有显著性意义，所求得的二元线性回归方程为∶线性回归方程的估计标准误为：西南科技大学生命科学与工程学院周海廷制作西南科技大学生命科学与工程学院周海廷制作 比较第８个观测点去掉前后预测平和Press的值从122.8819降为79.9550；对整个方程检验的F值从474.661上升为660.326；线性回归方程的估计标准误由2.53714下降为2.14515，表明该点对因变量预测值的影响是比较大的,值得注意。西南科技大学生命科学与工程学院周海廷制作[例2.2]　有人在某地抽样调查了29例儿童的血红蛋白与４种微量元素的含量， 资料如下，试问∶可否用４种微量元素(单位都是μmol/L)钙(X1)、镁(X2)、铁(X3)、 铜(X4)来较好地预测血红蛋白(Y，g/L)的含量？SAS程序：DATA d4p16a; INPUT y x1-x4; CARDS;135.0 13.70 12.68 80.32 0.16………………………………………….70.0 11.80 11.73 52.75 0.13 ;西南科技大学生命科学与工程学院周海廷制作PROC reg; MODEL y=x1-x4 / COLLIN COLLINOINT;RUN;共线性诊断选择项，对截距未进行校正共线性诊断选择项，对截距进行校正先用程序１作试探性分析，并用了共线性诊断的技术。西南科技大学生命科学与工程学院周海廷制作输出结果及其解释西南科技大学生命科学与工程学院周海廷制作程序１的参数估计结果∶截距项、X2、X4都无显著性意义，但不应过早将X2、X4从模型中去掉。最好等截距项从模型中去掉之后，重新拟合，视最后的结果再作决定。西南科技大学生命科学与工程学院周海廷制作这是共线性诊断的第１部分，即未对截距项校正的回归诊断结果∶从最后一行的条件数25.5585＞10可知,自变量之间有较强的共线性;从该行方差分量的数值可看出∶自变量之间的共线性主要表现在X2、X3两变量上。西南科技大学生命科学与工程学院周海廷制作这是共线性诊断的第２部分，即对截距项校正之后的回归诊断结果∶因本例的截距项无显著性意义，故用第１部分诊断结果就可以了。西南科技大学生命科学与工程学院周海廷制作用逐步回归法筛选自变量data d4p16b;set d4p16a;PROC REG; MODEL y=x1-x4 / NOINT SELECTION=STEPWISE SLE=0.30 SLS=0.10 STB;RUN;变量进入方程